27 |
#include "fizhi_ocean_coms.h" |
#include "fizhi_ocean_coms.h" |
28 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
29 |
|
|
30 |
integer myTime, myIter, myThid |
integer myIter, myThid |
31 |
|
_RL myTime |
32 |
|
|
33 |
logical alarm |
logical alarm |
34 |
external alarm |
external alarm |
85 |
C*********************************************************************** |
C*********************************************************************** |
86 |
|
|
87 |
if( alarm('turb') .or. alarm('radsw') ) then |
if( alarm('turb') .or. alarm('radsw') ) then |
88 |
call getlgr (sec,month,day,chlt,ityp,nchpland,nchp,nSx,nSy,bi,bj, |
call getlgr (sec,month,day,chlt,ityp,nchpland(bi,bj), |
89 |
. alai,agrn ) |
. nchp,nSx,nSy,bi,bj,alai,agrn ) |
90 |
endif |
endif |
91 |
|
|
92 |
C ********************************************************************** |
C ********************************************************************** |
96 |
if( alarm('radsw') ) then |
if( alarm('radsw') ) then |
97 |
call astro(nymd,nhms,lats,lons,im2*jm2,cosz,radius) |
call astro(nymd,nhms,lats,lons,im2*jm2,cosz,radius) |
98 |
call getalb(sec,month,day,cosz,snodep,fraci,fracl,im2,jm2,nchp, |
call getalb(sec,month,day,cosz,snodep,fraci,fracl,im2,jm2,nchp, |
99 |
. nchptot,nchpland,nSx,nSy,bi,bj,igrd,ityp,chfr,chlt,alai,agrn, |
. nchptot(bi,bj),nchpland(bi,bj),nSx,nSy,bi,bj,igrd,ityp, |
100 |
|
. chfr,chlt,alai,agrn, |
101 |
. albvisdr,albvisdf,albnirdr,albnirdf ) |
. albvisdr,albvisdf,albnirdr,albnirdf ) |
102 |
endif |
endif |
103 |
|
|
106 |
C ********************************************************************** |
C ********************************************************************** |
107 |
|
|
108 |
if( alarm('radlw') ) then |
if( alarm('radlw') ) then |
109 |
call grd2msc(fraci,im2,jm2,igrd,ficetile,nchp,nchptot) |
call grd2msc(fraci,im2,jm2,igrd,ficetile,nchp,nchptot(bi,bj)) |
110 |
call getemiss(fracl,im2,jm2,nchp,nchptot,nSx,nSy,bi,bj, |
call getemiss(fracl,im2,jm2,nchp,nchptot(bi,bj),nSx,nSy,bi,bj, |
111 |
. igrd,ityp,chfr,snodep,ficetile,emiss) |
. igrd,ityp,chfr,snodep,ficetile,emiss) |
112 |
endif |
endif |
113 |
|
|
121 |
tmpij(i,j) = 0. |
tmpij(i,j) = 0. |
122 |
enddo |
enddo |
123 |
enddo |
enddo |
124 |
do i = 1,nchpland |
do i = 1,nchpland(bi,bj) |
125 |
tmpchp(i) = tcanopy(i,bi,bj) |
tmpchp(i) = tcanopy(i,bi,bj) |
126 |
enddo |
enddo |
127 |
call msc2grd(igrd(1,bi,bj),chfr(1,bi,bj),tmpchp, |
call msc2grd(igrd(1,bi,bj),chfr(1,bi,bj),tmpchp, |
128 |
. nchp,nchpland,fracl,tmpij,im2,jm2) |
. nchp,nchpland(bi,bj),fracl,tmpij,im2,jm2) |
129 |
do j = jm1,jm2 |
do j = jm1,jm2 |
130 |
do i = im1,im2 |
do i = im1,im2 |
131 |
tgz(i,j,bi,bj) = tmpij(i,j) |
tgz(i,j,bi,bj) = tmpij(i,j) |