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

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

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

revision 1.1 by jscott, Fri Aug 11 19:35:32 2006 UTC revision 1.2 by jscott, Mon Apr 23 21:20:18 2007 UTC
# Line 44  C**** Line 44  C****
44  #include "BD2G04.COM"  #include "BD2G04.COM"
45    
46  #if ( defined CLM )  #if ( defined CLM )
47  #include "CLM.COM"  #include "CLM.h"
48  #endif  #endif
49    
50        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 62  C**** Line 62  C****
62        DIMENSION SINI(72),COSI(72)                                       5816.          DIMENSION SINI(72),COSI(72)                                       5816.  
63        LOGICAL POLE,PRNT,HPRNT        LOGICAL POLE,PRNT,HPRNT
64        common/conprn/HPRNT                  common/conprn/HPRNT          
65        common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)  !     common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)
66    #include "TSRF.COM"
67        common/SURRAD/TRSURF(JM0,4),SRSURF(JM0,4)        common/SURRAD/TRSURF(JM0,4),SRSURF(JM0,4)
68  c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV                                    5818.    c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV                                    5818.  
69        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 139  c     print *,'surface TAU=',TAU Line 140  c     print *,'surface TAU=',TAU
140        print *,' SURFACE FOR LAND ONLY'        print *,' SURFACE FOR LAND ONLY'
141        print *,' ZGS=30 m for LAND '        print *,' ZGS=30 m for LAND '
142        WMGMIN=8.        WMGMIN=8.
143        WMGMIN=4.        WMGMIN=5.
144        print *,'WMGMIN 4 LAND=',WMGMIN        print *,'WMGMIN 4 LAND=',WMGMIN
145          print *,'over land  WMG=max(WMG0,WMGMIN)'
146        IFIRST=0                                                          5877.          IFIRST=0                                                          5877.  
147        print *,' WMGE'        print *,' WMGE'
148          print 258,(WMGE(1,J),J=1,JM)          print 258,(WMGE(1,J),J=1,JM)
# Line 315  C**** ZERO OUT SURFACE DIAGNOSTICS WHICH Line 317  C**** ZERO OUT SURFACE DIAGNOSTICS WHICH
317        IM1=IM                                                            5969.          IM1=IM                                                            5969.  
318  #if ( defined CLM )  #if ( defined CLM )
319           if(NS.eq.1)then           if(NS.eq.1)then
320             tsl4clm(j)=0.0             i=1
321             qs4clm(j)=0.0             tsl4clm(i,j)=0.0
322             ps4clm(j)=0.0             qs4clm(i,j)=0.0
323             ws4clm(j)=0.0             ps4clm(i,j)=0.0
324             us4clm(j)=0.0             ws4clm(i,j)=0.0
325             vs4clm(j)=0.0             us4clm(i,j)=0.0
326               vs4clm(i,j)=0.0
327           endif           endif
328  #endif  #endif
329        DO 6000 I=1,IMAX                                                  5970.          DO 6000 I=1,IMAX                                                  5970.  
# Line 518  C**** LOOP OVER GROUND TIME STEPS Line 521  C**** LOOP OVER GROUND TIME STEPS
521        W1=SQRT(U1*U1+V1*V1)        W1=SQRT(U1*U1+V1*V1)
522        W1=SQRT(UG*UG+VG*VG)      ! 07/14/2006        W1=SQRT(UG*UG+VG*VG)      ! 07/14/2006
523        WS0=W1        WS0=W1
524    !     WMG=WMG0+WMGMIN
525    !     07/17/2006
526        WMG=WMG0+WMGMIN        WMG=WMG0+WMGMIN
527        WS=SQRT(W1**2+WMG)        WS=SQRT(W1**2+WMG)
528        if(J.eq.JPR)then        if(J.eq.JPR)then
# Line 582  c       print *,' ' Line 587  c       print *,' '
587    
588  #if ( defined CLM )  #if ( defined CLM )
589            if(ITYPE.EQ.4.or.ITYPE.EQ.3)then            if(ITYPE.EQ.4.or.ITYPE.EQ.3)then
590             tsl4clm(j)=tsl4clm(j)+TS*PTYPE/PLAND             tsl4clm(i,j)=tsl4clm(i,j)+TS*PTYPE/PLAND
591             qs4clm(j)=qs4clm(j)+QS*PTYPE/PLAND             qs4clm(i,j)=qs4clm(i,j)+QS*PTYPE/PLAND
592             ps4clm(j)=ps4clm(j)+PS*PTYPE/PLAND             ps4clm(i,j)=ps4clm(i,j)+PS*PTYPE/PLAND
593             ws4clm(j)=ws4clm(j)+WS*PTYPE/PLAND             ws4clm(i,j)=ws4clm(i,j)+WS*PTYPE/PLAND
594             us4clm(j)=us4clm(j)+US*PTYPE/PLAND             us4clm(i,j)=us4clm(i,j)+US*PTYPE/PLAND
595             vs4clm(j)=vs4clm(j)+VS*PTYPE/PLAND             vs4clm(i,j)=vs4clm(i,j)+VS*PTYPE/PLAND
596           endif           endif
597  #endif  #endif
598    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22