6 |
CBOP |
CBOP |
7 |
C !ROUTINE: FIZHI_INIT_VEGSURFTILES |
C !ROUTINE: FIZHI_INIT_VEGSURFTILES |
8 |
C !INTERFACE: |
C !INTERFACE: |
9 |
subroutine fizhi_init_vegsurftiles(nymd,nhms,prec,myThid) |
subroutine fizhi_init_vegsurftiles(globalArr,xsize,ysize, |
10 |
|
& nymd,nhms,prec,myThid) |
11 |
|
|
12 |
C !DESCRIPTION: |
C !DESCRIPTION: |
13 |
C Read in grid space values of the land state |
C Read in grid space values of the land state |
14 |
C and then convert to vegetation tile space |
C and then convert to vegetation tile space |
15 |
|
|
16 |
C !USES: |
C !USES: |
25 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
26 |
#include "PARAMS.h" |
#include "PARAMS.h" |
27 |
#ifdef ALLOW_EXCH2 |
#ifdef ALLOW_EXCH2 |
28 |
|
#include "W2_EXCH2_SIZE.h" |
29 |
#include "W2_EXCH2_TOPOLOGY.h" |
#include "W2_EXCH2_TOPOLOGY.h" |
|
#include "W2_EXCH2_PARAMS.h" |
|
30 |
#endif /* ALLOW_EXCH2 */ |
#endif /* ALLOW_EXCH2 */ |
31 |
|
|
32 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
33 |
|
integer xsize, ysize |
34 |
|
Real*8 globalArr(xsize,ysize,8) |
35 |
CHARACTER*1 prec |
CHARACTER*1 prec |
36 |
INTEGER nhms,nymd |
INTEGER nhms,nymd |
37 |
INTEGER myThid |
INTEGER myThid |
47 |
integer bislot,bjslot,iunit |
integer bislot,bjslot,iunit |
48 |
integer recl |
integer recl |
49 |
integer bi,bj,fileprec |
integer bi,bj,fileprec |
|
integer xsize, ysize |
|
|
#if defined(ALLOW_EXCH2) |
|
|
PARAMETER ( xsize = exch2_domain_nxt * sNx ) |
|
|
PARAMETER ( ysize = exch2_domain_nyt * sNy ) |
|
|
#else |
|
|
PARAMETER ( xsize = Nx ) |
|
|
PARAMETER ( ysize = Ny ) |
|
|
#endif |
|
|
Real*8 globalarr(xsize,ysize,8) |
|
50 |
_RL tempgrid(sNx,sNy) |
_RL tempgrid(sNx,sNy) |
51 |
_RL temptile(nchp) |
_RL temptile(nchp) |
52 |
_RL fracland(sNx,sNy,Nsx,Nsy) |
_RL fracland(sNx,sNy,Nsx,Nsy) |
84 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
85 |
|
|
86 |
#if defined(ALLOW_EXCH2) |
#if defined(ALLOW_EXCH2) |
87 |
bislot = exch2_txglobalo(W2_myTileList(bi))-1 |
bislot = exch2_txglobalo(W2_myTileList(bi,bj))-1 |
88 |
bjslot = exch2_tyglobalo(W2_myTileList(bi))-1 |
bjslot = exch2_tyglobalo(W2_myTileList(bi,bj))-1 |
89 |
#else |
#else |
90 |
bislot = myXGlobalLo-1+(bi-1)*sNx |
bislot = myXGlobalLo-1+(bi-1)*sNx |
91 |
bjslot = myYGlobalLo-1+(bj-1)*sNy |
bjslot = myYGlobalLo-1+(bj-1)*sNy |