44 |
integer bislot,bjslot,iunit |
integer bislot,bjslot,iunit |
45 |
integer recl |
integer recl |
46 |
integer bi,bj,fileprec |
integer bi,bj,fileprec |
47 |
Real*8 globalarr(Nx,Ny,8) |
integer xsize, ysize |
48 |
|
#if defined(ALLOW_EXCH2) |
49 |
|
PARAMETER ( xsize = exch2_domain_nxt * sNx ) |
50 |
|
PARAMETER ( ysize = exch2_domain_nyt * sNy ) |
51 |
|
#else |
52 |
|
PARAMETER ( xsize = Nx ) |
53 |
|
PARAMETER ( ysize = Ny ) |
54 |
|
#endif |
55 |
|
Real*8 globalarr(xsize,ysize,8) |
56 |
_RL tempgrid(sNx,sNy) |
_RL tempgrid(sNx,sNy) |
57 |
_RL temptile(nchp) |
_RL temptile(nchp) |
58 |
_RL fracland(sNx,sNy,Nsx,Nsy) |
_RL fracland(sNx,sNy,Nsx,Nsy) |
89 |
#endif /* ALLOW_EXCH2 */ |
#endif /* ALLOW_EXCH2 */ |
90 |
|
|
91 |
call get_landfrac(sNx,sNy,Nsx,Nsy,bi,bj,maxtyp, |
call get_landfrac(sNx,sNy,Nsx,Nsy,bi,bj,maxtyp, |
92 |
. surftype,tilefrac,fracland) |
. surftype,tilefrac,fracland(1,1,bi,bj)) |
93 |
|
|
94 |
do j = 1,sNy |
do j = 1,sNy |
95 |
do i = 1,sNx |
do i = 1,sNx |