/[MITgcm]/MITgcm/pkg/fizhi/update_earth_exports.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/update_earth_exports.F

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

revision 1.21 by molod, Wed Jul 28 22:08:40 2004 UTC revision 1.22 by molod, Thu Aug 12 15:21:22 2004 UTC
# Line 84  C*              Get Leaf-Area-Index and Line 84  C*              Get Leaf-Area-Index and
84  C***********************************************************************  C***********************************************************************
85    
86         if( alarm('turb') .or. alarm('radsw') ) then         if( alarm('turb') .or. alarm('radsw') ) then
87         call getlgr (sec,month,day,chlt,ityp,nchpland,nchp,nSx,nSy,bi,bj,         call getlgr (sec,month,day,chlt,ityp,nchpland(bi,bj),
88       .                                                       alai,agrn )       .       nchp,nSx,nSy,bi,bj,alai,agrn )
89         endif         endif
90    
91  C **********************************************************************  C **********************************************************************
# Line 95  C ************************************** Line 95  C **************************************
95         if( alarm('radsw') ) then         if( alarm('radsw') ) then
96          call astro(nymd,nhms,lats,lons,im2*jm2,cosz,radius)          call astro(nymd,nhms,lats,lons,im2*jm2,cosz,radius)
97          call getalb(sec,month,day,cosz,snodep,fraci,fracl,im2,jm2,nchp,          call getalb(sec,month,day,cosz,snodep,fraci,fracl,im2,jm2,nchp,
98       .    nchptot,nchpland,nSx,nSy,bi,bj,igrd,ityp,chfr,chlt,alai,agrn,       .    nchptot(bi,bj),nchpland(bi,bj),nSx,nSy,bi,bj,igrd,ityp,
99         .    chfr,chlt,alai,agrn,
100       .    albvisdr,albvisdf,albnirdr,albnirdf )       .    albvisdr,albvisdf,albnirdr,albnirdf )
101         endif         endif
102    
# Line 104  C                      Compute Surface E Line 105  C                      Compute Surface E
105  C **********************************************************************  C **********************************************************************
106    
107         if( alarm('radlw') ) then         if( alarm('radlw') ) then
108          call grd2msc(fraci,im2,jm2,igrd,ficetile,nchp,nchptot)          call grd2msc(fraci,im2,jm2,igrd,ficetile,nchp,nchptot(bi,bj))
109          call getemiss(fracl,im2,jm2,nchp,nchptot,nSx,nSy,bi,bj,          call getemiss(fracl,im2,jm2,nchp,nchptot(bi,bj),nSx,nSy,bi,bj,
110       .   igrd,ityp,chfr,snodep,ficetile,emiss)       .   igrd,ityp,chfr,snodep,ficetile,emiss)
111         endif         endif
112    
# Line 119  C*************************************** Line 120  C***************************************
120          tmpij(i,j) = 0.          tmpij(i,j) = 0.
121         enddo         enddo
122         enddo         enddo
123         do i = 1,nchpland         do i = 1,nchpland(bi,bj)
124          tmpchp(i) = tcanopy(i,bi,bj)          tmpchp(i) = tcanopy(i,bi,bj)
125         enddo         enddo
126         call msc2grd(igrd(1,bi,bj),chfr(1,bi,bj),tmpchp,         call msc2grd(igrd(1,bi,bj),chfr(1,bi,bj),tmpchp,
127       .                           nchp,nchpland,fracl,tmpij,im2,jm2)       .                    nchp,nchpland(bi,bj),fracl,tmpij,im2,jm2)
128         do j = jm1,jm2         do j = jm1,jm2
129         do i = im1,im2         do i = im1,im2
130          tgz(i,j,bi,bj) = tmpij(i,j)          tgz(i,j,bi,bj) = tmpij(i,j)
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22