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

Diff of /MITgcm_contrib/jscott/igsm/src/daily_new.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:30 2006 UTC revision 1.4 by jscott, Tue Sep 1 22:03:56 2009 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "ctrparam.h"  #include "ctrparam.h"
5    
# Line 37  C**** Line 39  C****
39       &  sstmin(12,2)       &  sstmin(12,2)
40       &       ,miceo(JM0)       &       ,miceo(JM0)
41        common/qfl/QFLUX(JM0,0:13),ZOAV(JM0),QFLUXT(JM0)        common/qfl/QFLUX(JM0,0:13),ZOAV(JM0),QFLUXT(JM0)
42        common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTSURF(JM0)  !     common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTSURF(JM0)
43    #include "TSRF.COM"
44        common/fixcld/cldssm(JM0,LM0,0:13),cldmcm(JM0,LM0,0:13),        common/fixcld/cldssm(JM0,LM0,0:13),cldmcm(JM0,LM0,0:13),
45       &  CLDSST(JM0,LM0),       &  CLDSST(JM0,LM0),
46       &  CLDMCT(JM0,LM0)       &  CLDMCT(JM0,LM0)
47        common/surps/srps(JM0+3),nsrps        common/surps/srps(JM0+3),nsrps
48  #if ( defined OCEAN_3D || defined ML_2D)  #if ( defined OCEAN_3D || defined ML_2D)
49  #include "AGRID.COM"  #include "AGRID.h"
50    #endif
51    #if ( defined CLM )
52    #include "CLM.h"
53  #endif  #endif
54        LOGICAL HPRNT        LOGICAL HPRNT
55        common/conprn/HPRNT,JPR,LPR        common/conprn/HPRNT,JPR,LPR
# Line 146  c     MONTH=7 Line 152  c     MONTH=7
152   7241  continue   7241  continue
153        endif        endif
154  #if (defined OCEAN_3D || defined ML_2D)  #if (defined OCEAN_3D || defined ML_2D)
       do 725 j=1,JM  
        DTSURF(j)=TSURFD(j)-TSURFT(j)  
        TSURFD(j)=0.  
   725  continue  
155         if(JDATE.le.16)then         if(JDATE.le.16)then
156         do 723 j=1,JM         do 723 j=1,JM
157           TSURFT(j)=((16-JDATE)*TSURFC(j,MONTH-1)+           TSURFT(j)=((16-JDATE)*TSURFC(j,MONTH-1)+
158       *   (JDATE+15)*TSURFC(j,MONTH))/31.       *   (JDATE+15)*TSURFC(j,MONTH))/31.
159             TLANDT(j)=((16-JDATE)*TLANDC(j,MONTH-1)+
160         *   (JDATE+15)*TLANDC(j,MONTH))/31.
161   723  continue   723  continue
162        else        else
163         do 724 j=1,JM         do 724 j=1,JM
164          TSURFT(j)=((JDATE-16)*TSURFC(j,MONTH+1)+          TSURFT(j)=((JDATE-16)*TSURFC(j,MONTH+1)+
165       *            (31-JDATE+16)*TSURFC(j,MONTH))/31.       *            (31-JDATE+16)*TSURFC(j,MONTH))/31.
166            TLANDT(j)=((JDATE-16)*TLANDC(j,MONTH+1)+
167         *            (31-JDATE+16)*TLANDC(j,MONTH))/31.
168    724  continue    724  continue
169        endif        endif
170    !       print *,'From daily_new TSURFD'
171    !       print *,TSURFD
172    !       print *,'TSURFT'
173    !       print *,TSURFT
174    !       print *,'From daily_new TLANDD'
175    !       print *,TLANDD
176    !       print *,'TLANDT'
177    !       print *,TLANDT
178          do 725 j=1,JM
179           DT2MGL(j)=TSURFD(j)-TSURFT(j)
180           DT2MLD(j)=TLANDD(j)-TLANDT(j)
181           TSURFD(j)=0.
182           TLANDD(j)=0.
183      725  continue
184    #if ( defined CLM )
185            DT2MLAND=0.
186            if(PRTREND) then
187            AREAL=0.
188            do j=1,jm
189             DT2MLAND=DT2MLAND+DT2MLD(J)*DXYP(j)*FDATA(1,j,2)
190             AREAL=AREAL+DXYP(j)*FDATA(1,j,2)
191            end do  !j
192            DT2MLAND=DT2MLAND/AREAL
193    !       print *,'DT2MLD'
194    !       print *,DT2MLD
195           if(JDATE.eq.1)then
196            print *,'JDATE=',JDATE,' DT2MLAND=',DT2MLAND
197           endif
198           endif
199    !       print *,'AREAL=',AREAL
200    #endif
201    
202  #endif  #endif
203    
204        RETURN                                                            1108.5        RETURN                                                            1108.5
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22