/[MITgcm]/MITgcm_contrib/jscott/igsm/src/surf_ocean.F
ViewVC logotype

Diff of /MITgcm_contrib/jscott/igsm/src/surf_ocean.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by jscott, Tue Aug 22 20:25:52 2006 UTC revision 1.4 by jscott, Wed Sep 2 15:43:37 2009 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "ctrparam.h"  #include "ctrparam.h"
5    
# Line 44  C**** Line 46  C****
46  #include "BD2G04.COM"  #include "BD2G04.COM"
47    
48  #if ( defined CLM )  #if ( defined CLM )
49  #include "CLM.COM"  #include "CLM.h"
50  #endif  #endif
51    
52        COMMON/SPEC2/KM,KINC,COEK,C3LAND(IO0,JM0),C3OICE(IO0,JM0)         5808.1          COMMON/SPEC2/KM,KINC,COEK,C3LAND(IO0,JM0),C3OICE(IO0,JM0)         5808.1  
# Line 63  C**** Line 65  C****
65        DIMENSION WMGMIN(JM0)        DIMENSION WMGMIN(JM0)
66        LOGICAL POLE,PRNT,HPRNT        LOGICAL POLE,PRNT,HPRNT
67        common/conprn/HPRNT                  common/conprn/HPRNT          
68        common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)  !     common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)
69    #include "TSRF.COM"
70        common/SURRAD/TRSURF(JM0,4),SRSURF(JM0,4)        common/SURRAD/TRSURF(JM0,4),SRSURF(JM0,4)
71  c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV                                    5818.    c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV                                    5818.  
72        COMMON/CWMG/WMGEA(JM0),NWMGEA(JM0),CHAVER(JM0),DTAV(JM0),DQAV(JM0)        COMMON/CWMG/WMGEA(JM0),NWMGEA(JM0),CHAVER(JM0),DTAV(JM0),DQAV(JM0)
# Line 72  c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV Line 75  c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV
75       & WSSL(JM0),T2ML(JM0),       & WSSL(JM0),T2ML(JM0),
76       & TSSL(JM0),QSSL(JM0),USSL(JM0),VSSL(JM0),TAUSL(JM0),BLJ(JM0,50)       & TSSL(JM0),QSSL(JM0),USSL(JM0),VSSL(JM0),TAUSL(JM0),BLJ(JM0,50)
77       & ,ELHTG(JM0),SHTG(JM0),TAUXG(JM0),TAUYG(JM0)       & ,ELHTG(JM0),SHTG(JM0),TAUXG(JM0),TAUYG(JM0)
78           real T2MZ(JM0)
79  C  C
80  #if ( defined OCEAN_3D || defined ML_2D)  #if ( defined OCEAN_3D || defined ML_2D)
81  #include "AGRID.h"  #include "AGRID.h"
# Line 186  C     SRCORX=1. Line 190  C     SRCORX=1.
190        IQ2=IM/2+1                                                        5882.          IQ2=IM/2+1                                                        5882.  
191        IQ3=3*IM/4+1                                                      5883.          IQ3=3*IM/4+1                                                      5883.  
192        DTSURF=NDYN*DT/NSURF                                              5884.          DTSURF=NDYN*DT/NSURF                                              5884.  
193          NCLMPERDAY=(24.*3600.)/DTSURF
194        print *,' DTSURF=',DTSURF        print *,' DTSURF=',DTSURF
195          print *,' NCLMPERDAY=',NCLMPERDAY
196           DTSRCE=DT*NDYN                                                 5885.             DTSRCE=DT*NDYN                                                 5885.  
197        SHA=RGAS/KAPA                                                     5886.          SHA=RGAS/KAPA                                                     5886.  
198        RVX=0.                                                            5887.          RVX=0.                                                            5887.  
# Line 780  c     print *,j,DU1(1,j),DV1(1,j) Line 786  c     print *,j,DU1(1,j),DV1(1,j)
786        BLDATA(I,J,6)=USS                                                 6593.          BLDATA(I,J,6)=USS                                                 6593.  
787        BLDATA(I,J,7)=VSS                                                 6594.          BLDATA(I,J,7)=VSS                                                 6594.  
788        BLDATA(I,J,8)=TAUS                                                6595.          BLDATA(I,J,8)=TAUS                                                6595.  
789         T2ML(J)=T2MS         T2MZ(J)=T2MS
790  c     print *,j,T(I,J,1),Q(I,J,1)  c     print *,j,T(I,J,1),Q(I,J,1)
791  c     print *,(TGRND(I,J,k),k=1,4)  c     print *,(TGRND(I,J,k),k=1,4)
792  c     print *,(EVAPOR(I,J,k),k=1,4)  c     print *,(EVAPOR(I,J,k),k=1,4)
# Line 1200  C**** MIX MOMENTUM THROUGHOUT THE BOUNDA Line 1206  C**** MIX MOMENTUM THROUGHOUT THE BOUNDA
1206        endif        endif
1207        enddo        enddo
1208   9000 CONTINUE                                                          6794.     9000 CONTINUE                                                          6794.  
1209  c     print *,' From surf_ocean T2ML'  !     print *,' From surf_ocean T2MZ'
1210  c     print *,T2ML  !     print *,T2MZ
1211        do 9001 J=1,JM        do 9001 J=1,JM
1212  c      TSURFD(J)=TSURFD(J)+(BLDATA(1,J,2)-273.16)/24.  c      TSURFD(J)=TSURFD(J)+(BLDATA(1,J,2)-273.16)/24.
1213         TSURFD(J)=TSURFD(J)+(T2ML(J)-273.16)/24.  !      TSURFD(J)=TSURFD(J)+(T2MZ(J)-273.16)/24.
1214           TSURFD(J)=TSURFD(J)+(T2MZ(J)-273.16)/NCLMPERDAY
1215         T2ML(J)=0.0         T2ML(J)=0.0
1216   9001 continue   9001 continue
1217        RETURN                                                            6795.          RETURN                                                            6795.  
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22