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

Diff of /MITgcm_contrib/jscott/igsm/src/surface.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.3 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 63  C**** Line 63  C****
63        DIMENSION WMGMINO(JM0)        DIMENSION WMGMINO(JM0)
64        LOGICAL POLE,PRNT,HPRNT        LOGICAL POLE,PRNT,HPRNT
65        common/conprn/HPRNT                  common/conprn/HPRNT          
66        common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)  !     common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)
67    #include "TSRF.COM"
68        common/SURRAD/TRSURF(JM0,4),SRSURF(JM0,4)        common/SURRAD/TRSURF(JM0,4),SRSURF(JM0,4)
69  c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV                                    5818.    c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV                                    5818.  
70        COMMON/CWMG/WMGEA(JM0),NWMGEA(JM0),CHAVER(JM0),DTAV(JM0),DQAV(JM0)        COMMON/CWMG/WMGEA(JM0),NWMGEA(JM0),CHAVER(JM0),DTAV(JM0),DQAV(JM0)
71       & ,Z0AV(JM0),WSAV(JM0),WS0AV(JM0),TAUAV(JM0)       & ,Z0AV(JM0),WSAV(JM0),WS0AV(JM0),TAUAV(JM0)
72  C  C
73  #if ( defined OCEAN_3D || defined ML_2D)  #if ( defined OCEAN_3D || defined ML_2D)
74  #include "AGRID.COM"  #include "AGRID.h"
75  #endif  #endif
76  c  c
77        DATA RVAP/461.5/                                                  5819.          DATA RVAP/461.5/                                                  5819.  
# Line 323  C**** ZERO OUT SURFACE DIAGNOSTICS WHICH Line 324  C**** ZERO OUT SURFACE DIAGNOSTICS WHICH
324        IM1=IM                                                            5969.          IM1=IM                                                            5969.  
325  #if ( defined CLM )  #if ( defined CLM )
326           if(NS.eq.1)then           if(NS.eq.1)then
327             tsl4clm(j)=0.0             i=1
328             qs4clm(j)=0.0             tsl4clm(i,j)=0.0
329             ps4clm(j)=0.0             qs4clm(i,j)=0.0
330             ws4clm(j)=0.0             ps4clm(i,j)=0.0
331             us4clm(j)=0.0             ws4clm(i,j)=0.0
332             vs4clm(j)=0.0             us4clm(i,j)=0.0
333               vs4clm(i,j)=0.0
334           endif           endif
335  #endif  #endif
336        DO 6000 I=1,IMAX                                                  5970.          DO 6000 I=1,IMAX                                                  5970.  
# Line 595  c     WS=SQRT(W1*W1+0.8*WMG) Line 597  c     WS=SQRT(W1*W1+0.8*WMG)
597          print *,'TG=',TG,' QG=',QG          print *,'TG=',TG,' QG=',QG
598          print *,'RVX=',RVX,' TG1=',TG1          print *,'RVX=',RVX,' TG1=',TG1
599        endif        endif
600    !      if(TAU.gt.3623.0.and.TAU.lt.3744.0) then
601    !        if(ITYPE.eq.4)then
602    !         if(J.eq.29)then
603    !           write (329) ,TAU,TOFDAY,TG1
604    !         endif
605    !         endif
606    !       endif
607    
608  #if ( defined CPL_OCEANCO2  || defined OCEAN_3D)  #if ( defined CPL_OCEANCO2  || defined OCEAN_3D)
609        if(ITYPE.eq.1)then        if(ITYPE.eq.1)then
# Line 677  c     TRHEAT=TRHR(I,J,1)-STBO*(TG*TG)*(T Line 686  c     TRHEAT=TRHR(I,J,1)-STBO*(TG*TG)*(T
686  #if ( defined CLM )  #if ( defined CLM )
687           if(NS.eq.1)then           if(NS.eq.1)then
688            if(ITYPE.EQ.4.or.ITYPE.EQ.3)then            if(ITYPE.EQ.4.or.ITYPE.EQ.3)then
689             tsl4clm(j)=tsl4clm(j)+TS*PTYPE/PLAND             tsl4clm(i,j)=tsl4clm(i,j)+TS*PTYPE/PLAND
690             qs4clm(j)=qs4clm(j)+QS*PTYPE/PLAND             qs4clm(i,j)=qs4clm(i,j)+QS*PTYPE/PLAND
691             ps4clm(j)=ps4clm(j)+PS*PTYPE/PLAND             ps4clm(i,j)=ps4clm(i,j)+PS*PTYPE/PLAND
692             ws4clm(j)=ws4clm(j)+WS*PTYPE/PLAND             ws4clm(i,j)=ws4clm(i,j)+WS*PTYPE/PLAND
693             us4clm(j)=us4clm(j)+US*PTYPE/PLAND             us4clm(i,j)=us4clm(i,j)+US*PTYPE/PLAND
694             vs4clm(j)=vs4clm(j)+VS*PTYPE/PLAND             vs4clm(i,j)=vs4clm(i,j)+VS*PTYPE/PLAND
695           endif           endif
696          endif          endif
697  #endif  #endif
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22