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

Diff of /MITgcm_contrib/jscott/igsm/src/sur4clm.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    
45  #include "BD2G04.COM"  #include "BD2G04.COM"
46    
47  #include "CLM.COM"  #include "CLM.h"
48    
49        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  
50       *  ,C3LICE(IO0,JM0),WMGE(IO0,JM0),TSSFC(1,JM0,4)                   5808.2         *  ,C3LICE(IO0,JM0),WMGE(IO0,JM0),TSSFC(1,JM0,4)                   5808.2  
# Line 60  C**** Line 60  C****
60        COMMON/RDATA/ROUGHL(IO0,JM0)                                      5815.          COMMON/RDATA/ROUGHL(IO0,JM0)                                      5815.  
61        LOGICAL POLE,PRNT,HPRNT        LOGICAL POLE,PRNT,HPRNT
62        common/conprn/HPRNT                  common/conprn/HPRNT          
63        common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)  !     common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)
64    #include "TSRF.COM"
65        common/SURRAD/TRSURF(JM0,4),SRSURF(JM0,4)        common/SURRAD/TRSURF(JM0,4),SRSURF(JM0,4)
66  c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV                                    5818.    c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV                                    5818.  
67        integer IQ1,IQ2,IQ3        integer IQ1,IQ2,IQ3
# Line 136  C     SRCORX=1. Line 137  C     SRCORX=1.
137        IQ1=IM/4+1                                                        5881.          IQ1=IM/4+1                                                        5881.  
138        IQ2=IM/2+1                                                        5882.          IQ2=IM/2+1                                                        5882.  
139        IQ3=3*IM/4+1                                                      5883.          IQ3=3*IM/4+1                                                      5883.  
140        DTSURF=NDYN*DT/NSURF                                              5884.    !     DTSURF=NDYN*DT/NSURF                                              5884.  
141        print *,' DTSURF=',DTSURF  !     print *,' From SRF4CLM DTSURF=',DTSURF
142           DTSRCE=DT*NDYN                                                 5885.    !        DTSRCE=DT*NDYN                                                 5885.  
143        SHA=RGAS/KAPA                                                     5886.          SHA=RGAS/KAPA                                                     5886.  
144        RVX=0.                                                            5887.          RVX=0.                                                            5887.  
145     50 CONTINUE     50 CONTINUE
# Line 202  c     print *,'After 100' Line 203  c     print *,'After 100'
203        BWS=0.        BWS=0.
204        BWMG=0.        BWMG=0.
205        IM1=IM                                                            5969.          IM1=IM                                                            5969.  
206         tsl4clm(j)=0.0         i=1
207         qs4clm(j)=0.0         tsl4clm(i,j)=0.0
208         ps4clm(j)=0.0         qs4clm(i,j)=0.0
209         ws4clm(j)=0.0         ps4clm(i,j)=0.0
210         us4clm(j)=0.0         ws4clm(i,j)=0.0
211         vs4clm(j)=0.0         us4clm(i,j)=0.0
212           vs4clm(i,j)=0.0
213        DO 6000 I=1,IMAX                                                  5970.          DO 6000 I=1,IMAX                                                  5970.  
214  C****                                                                   5971.    C****                                                                   5971.  
215  C**** DETERMINE SURFACE CONDITIONS                                      5972.    C**** DETERMINE SURFACE CONDITIONS                                      5972.  
# Line 330  c     QS=QS+X*(SHA/ELHX) Line 332  c     QS=QS+X*(SHA/ELHX)
332  c3500 CONTINUE  c3500 CONTINUE
333    
334            if(ITYPE.EQ.4.or.ITYPE.EQ.3)then            if(ITYPE.EQ.4.or.ITYPE.EQ.3)then
335             tsl4clm(j)=tsl4clm(j)+TS*PTYPE/PLAND             tsl4clm(i,j)=tsl4clm(i,j)+TS*PTYPE/PLAND
336             qs4clm(j)=qs4clm(j)+QS*PTYPE/PLAND             qs4clm(i,j)=qs4clm(i,j)+QS*PTYPE/PLAND
337             ps4clm(j)=ps4clm(j)+PS*PTYPE/PLAND             ps4clm(i,j)=ps4clm(i,j)+PS*PTYPE/PLAND
338             ws4clm(j)=ws4clm(j)+WS*PTYPE/PLAND             ws4clm(i,j)=ws4clm(i,j)+WS*PTYPE/PLAND
339  !          us4clm(j)=us4clm(j)+US*PTYPE/PLAND  !          us4clm(i,j)=us4clm(i,j)+US*PTYPE/PLAND
340  !          vs4clm(j)=vs4clm(j)+VS*PTYPE/PLAND  !          vs4clm(i,j)=vs4clm(i,j)+VS*PTYPE/PLAND
341             us4clm(j)=us4clm(j)+RW*US*PTYPE/PLAND             us4clm(i,j)=us4clm(i,j)+RW*US*PTYPE/PLAND
342             vs4clm(j)=vs4clm(j)+RW*VS*PTYPE/PLAND             vs4clm(i,j)=vs4clm(i,j)+RW*VS*PTYPE/PLAND
343           endif           endif
344    
345        TSSFC(I,J,ITYPE)=TS                                               6521.5          TSSFC(I,J,ITYPE)=TS                                               6521.5  
# Line 383  C**** QUANTITIES ACCUMULATED FOR SURFACE Line 385  C**** QUANTITIES ACCUMULATED FOR SURFACE
385           BLJ(J,28)=BWMG           BLJ(J,28)=BWMG
386           BLJ(J,23)=BTS             BLJ(J,23)=BTS  
387   7000 CONTINUE                                                          6677.     7000 CONTINUE                                                          6677.  
388    !     print *,' From SRF4CLM TAU=',TAU
389    !          i=1
390    !          j=1
391    !          print *,'TS(1),TS(2)'
392    !          print *,tsl4clm(1,1),tsl4clm(1,2)
393    !          print *,dsw4clm(i,j),dlw4clm(i,j)
394    !          print *,swinr4clm(i,j), swvis4clm(i,j)
395  C****                                                                   6678.    C****                                                                   6678.  
396  #endif  #endif
397        RETURN                                                            6795.          RETURN                                                            6795.  
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22