/[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.9 by molod, Thu Jun 10 21:50:33 2004 UTC revision 1.10 by molod, Wed Jun 16 19:19:49 2004 UTC
# Line 75  C*************************************** Line 75  C***************************************
75  C*              Get Leaf-Area-Index and Greenness Index                *  C*              Get Leaf-Area-Index and Greenness Index                *
76  C***********************************************************************  C***********************************************************************
77    
78        if( alarm('turb') .or. alarm('radsw') ) then         if( alarm('turb') .or. alarm('radsw') ) then
79        call getlgr (sec,month,day,chlt,ityp,nchpland,nSx,nSy,bi,bj,         call getlgr (sec,month,day,chlt,ityp,nchpland,nSx,nSy,bi,bj,
80       .                                                       alai,agrn )       .                                                       alai,agrn )
81        endif         endif
82    
83  C **********************************************************************  C **********************************************************************
84  C                      Compute Surface Albedo  C                      Compute Surface Albedo
85  C **********************************************************************  C **********************************************************************
86    
87        if( alarm('radsw') ) then         if( alarm('radsw') ) then
88         call astro(nymd,nhms,lats,lons,im2*jm2,cosz,ra)          call astro(nymd,nhms,lats,lons,im2*jm2,cosz,ra)
89         call getalb(sec,month,day,cosz,snodep,fraci,fracl,im2,jm2,nchp,          call getalb(sec,month,day,cosz,snodep,fraci,fracl,im2,jm2,nchp,
90       .             nchpland,nSx,nSy,bi,bj,igrd,ityp,chfr,chlt,alai,agrn,       .             nchpland,nSx,nSy,bi,bj,igrd,ityp,chfr,chlt,alai,agrn,
91       .             albvisdr,albvisdf,albnirdr,albnirdf )       .             albvisdr,albvisdf,albnirdr,albnirdf )
92        endif         endif
93    
94    
95  C **********************************************************************  C **********************************************************************
96  C                      Compute Surface Emissivity  C                      Compute Surface Emissivity
97  C **********************************************************************  C **********************************************************************
98    
99        if( alarm('radlw') ) then         if( alarm('radlw') ) then
100         call grd2msc(fraci,im2,jm2,igrd,ficetile,nchp,nchp)          call grd2msc(fraci,im2,jm2,igrd,ficetile,nchp,nchp)
101         call getemiss(fracl,im2,jm2,nchp,nSx,nSy,bi,bj,igrd,ityp,chfr,          call getemiss(fracl,im2,jm2,nchp,nSx,nSy,bi,bj,igrd,ityp,chfr,
102       .      snodep,ficetile,emiss)       .      snodep,ficetile,emiss)
103        endif         endif
104    
105    
106  C*********************************************************************  C*********************************************************************
107  C            Ground Temperature Over Ocean is from SST array,  C            Ground Temperature Over Ocean is from SST array,
108    C               Over land is from tcanopy
109  C*********************************************************************  C*********************************************************************
110    
111        do j = 1,jm         do j = jm1,jm2
112        do i = 1,im         do i = im1,im2
113        if(fracl(i,j).lt.0.3.and.sea_ice(i,j).eq.0.0)tgz(i,j) = sst(i,j)          tgz(i,j,bi,bj) = 0.
114        endif         enddo
115        enddo         enddo
116        enddo         call msc2grd(igrd(1,bi,bj),chfr(1,bi,bj),tcanopy(1,bi,bj),
117         .                           nchp,nchp,fracl,tgz(1,bi,bj),im2,jm2)
118           do j = jm1,jm2
119           do i = im1,im2
120            if(fracl(i,j).lt.0.3.and.sea_ice(i,j,bi,bj).eq.0.0)
121         .                                 tgz(i,j,bi,bj) = sst(i,j,bi,bj)
122           enddo
123           enddo
124    
125        enddo        enddo
126        enddo        enddo
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22