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

Diff of /MITgcm_contrib/jscott/igsm/src/gr_clm.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.3 by jscott, Mon Apr 23 21:20:18 2007 UTC
# Line 1  Line 1 
1    c source sokolov users 14712 Aug 25 14:32 /home/sokolov/IGSM2/SRC/gr_clm.F
2  #include "ctrparam.h"  #include "ctrparam.h"
3    
4  !       ==========================================================  !       ==========================================================
# Line 28  C**** Line 28  C****
28  #if ( defined OCEAN_3D || defined ML_2D )  #if ( defined OCEAN_3D || defined ML_2D )
29  #include "AGRID.h"  #include "AGRID.h"
30  #endif  #endif
31  #include "CLM.COM"  #include "CLM.h"
32        dimension DZSOI(nlevsoi)        dimension DZSOI(nlevsoi)
33    
34    
# Line 193  c     print *,'From gr_clm J=',j,' TAU=' Line 193  c     print *,'From gr_clm J=',j,' TAU='
193  c     print *,tsoiclm(i,j,1),tsoiclm(i,j,nlevsoi)  c     print *,tsoiclm(i,j,1),tsoiclm(i,j,nlevsoi)
194  c     print *,TG1,TG2,TG11  c     print *,TG1,TG2,TG11
195  c     print *,WTR1,WTR2,SNOW  c     print *,WTR1,WTR2,SNOW
196        albinr=0.7*aldirclm(j)+0.3*aldifclm(j)        albinr=0.7*aldirclm(i,j)+0.3*aldifclm(i,j)
197        albvis=0.7*asdirclm(j)+0.3*asdifclm(j)        albvis=0.7*asdirclm(i,j)+0.3*asdifclm(i,j)
198        swuclm=albinr*swinr4clm(j)+albvis*swvis4clm(j)        swuclm=albinr*swinr4clm(i,j)+albvis*swvis4clm(i,j)
199    !        if(J.eq.29)then
200    !          if(TAU.ge.3644.0.and.TAU.le.3743.0) then
201    !         print *,'Form gr_clm J=',j,' TAU=',TAU
202    !         print *,' albinr=',albinr,' albvis=',albvis
203    !         print *,' swuclm=',swuclm
204    !         endif
205    !        endif
206    
207        GDATA(I,J,2)=SNOW        GDATA(I,J,2)=SNOW
208        GDATA(I,J,4)=TG1                                                  7294.5          GDATA(I,J,4)=TG1                                                  7294.5  
# Line 241  C**** LONGITUDINALLY INTEGRATED QUANTITI Line 248  C**** LONGITUDINALLY INTEGRATED QUANTITI
248           BJ(J,19)=BJ(J,19)+BEVAP                                        7320.5             BJ(J,19)=BJ(J,19)+BEVAP                                        7320.5  
249  c        BJ(J,40)=BJ(J,40)+BERUN0                                       7323.5    c        BJ(J,40)=BJ(J,40)+BERUN0                                       7323.5  
250  c        BJ(J,41)=BJ(J,41)+BEDIFS                                       7324.    c        BJ(J,41)=BJ(J,41)+BEDIFS                                       7324.  
251           BJ(J,42)=BJ(J,42)+BF1DT                                        7325.    !        BJ(J,42)=BJ(J,42)+BF1DT                                        7325.  
252  c        BJ(J,43)=BJ(J,43)+BERUN2                                       7326.5    c        BJ(J,43)=BJ(J,43)+BERUN2                                       7326.5  
253           BJ(J,43)=BJ(J,43)+BSWUP           BJ(J,43)=BJ(J,43)+BSWUP
254  c        BJ(J,45)=BJ(J,45)+BDIFS                                        7327.5    c        BJ(J,45)=BJ(J,45)+BDIFS                                        7327.5  
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22