/[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.18 by molod, Mon Jul 26 20:23:16 2004 UTC revision 1.21 by molod, Wed Jul 28 22:08:40 2004 UTC
# Line 99  C ************************************** Line 99  C **************************************
99       .    albvisdr,albvisdf,albnirdr,albnirdf )       .    albvisdr,albvisdf,albnirdr,albnirdf )
100         endif         endif
101    
   
102  C **********************************************************************  C **********************************************************************
103  C                      Compute Surface Emissivity  C                      Compute Surface Emissivity
104  C **********************************************************************  C **********************************************************************
105    
106         if( alarm('radlw') ) then         if( alarm('radlw') ) then
107          call grd2msc(fraci,im2,jm2,igrd,ficetile,nchp,nchp)          call grd2msc(fraci,im2,jm2,igrd,ficetile,nchp,nchptot)
108          call getemiss(fracl,im2,jm2,nchp,nchptot,nSx,nSy,bi,bj,          call getemiss(fracl,im2,jm2,nchp,nchptot,nSx,nSy,bi,bj,
109       .   igrd,ityp,chfr,snodep,ficetile,emiss)       .   igrd,ityp,chfr,snodep,ficetile,emiss)
110         endif         endif
111    
   
112  C*********************************************************************  C*********************************************************************
113  C            Ground Temperature Over Ocean is from SST array,  C            Ground Temperature Over Ocean is from SST array,
114  C               Over land is from tcanopy  C               Over land is from tcanopy
# Line 121  C*************************************** Line 119  C***************************************
119          tmpij(i,j) = 0.          tmpij(i,j) = 0.
120         enddo         enddo
121         enddo         enddo
122         do i = 1,nchp         do i = 1,nchpland
123          tmpchp(i) = tcanopy(i,bi,bj)          tmpchp(i) = tcanopy(i,bi,bj)
124         enddo         enddo
125         call msc2grd(igrd(1,bi,bj),chfr(1,bi,bj),tmpchp,         call msc2grd(igrd(1,bi,bj),chfr(1,bi,bj),tmpchp,
126       .                           nchp,nchptot,fracl,tmpij,im2,jm2)       .                           nchp,nchpland,fracl,tmpij,im2,jm2)
127         do j = jm1,jm2         do j = jm1,jm2
128         do i = im1,im2         do i = im1,im2
129          tgz(i,j,bi,bj) = tmpij(i,j)          tgz(i,j,bi,bj) = tmpij(i,j)
# Line 892  C and now call sibalb Line 890  C and now call sibalb
890    
891  C finally some transformations back to grid space for albedos  C finally some transformations back to grid space for albedos
892    
       print *,' In getalb, chfr: '  
       print *,(chfr(i,1,1),i=1,nchptot)  
   
893        DO I=1,IM        DO I=1,IM
894        DO J=1,JM        DO J=1,JM
895         tmpij(i,j) = albvr(i,j,bi,bj)         tmpij(i,j) = albvr(i,j,bi,bj)
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22