8 |
C !ROUTINE: FIZHI_MNC_INIT |
C !ROUTINE: FIZHI_MNC_INIT |
9 |
|
|
10 |
C !INTERFACE: |
C !INTERFACE: |
11 |
SUBROUTINE FIZHI_MNC_INIT( |
SUBROUTINE FIZHI_MNC_INIT( myThid ) |
|
I sNx,sNy, OLx,OLy, nSx,nSy, nPx,nPy, NrPhys, |
|
|
I myThid ) |
|
12 |
|
|
13 |
C !DESCRIPTION: |
C !DESCRIPTION: |
14 |
C Create the pre-defined grid types and variable types. |
C Create some pre-defined MNC grid types and variable types useful |
15 |
|
C for the FIZHI package. Borrows heavily from MNC_CW_INIT(). |
|
C The grid type is a character string that encodes the presence and |
|
|
C types associated with the four possible dimensions. The character |
|
|
C string follows the format |
|
|
C \begin{center} |
|
|
C \texttt{H0\_H1\_H2\_\_V\_\_T} |
|
|
C \end{center} |
|
|
C where the terms \textit{H0}, \textit{H1}, \textit{H2}, \textit{V}, |
|
|
C \textit{T} can be almost any combination of the following: |
|
|
C \begin{center} |
|
|
C \begin{tabular}[h]{|ccc|c|c|}\hline |
|
|
C \multicolumn{3}{|c|}{Horizontal} & Vertical & Time \\ |
|
|
C \textit{H0}: location & \textit{H1}: dimensions & \textit{H2}: halo |
|
|
C & \textit{V}: location & \textit{T}: level \\\hline |
|
|
C \texttt{-} & xy & Hn & \texttt{-} & \texttt{-} \\ |
|
|
C U & x & Hy & i & t \\ |
|
|
C V & y & & c & \\ |
|
|
C Cen & & & & \\ |
|
|
C Cor & & & & \\\hline |
|
|
C \end{tabular} |
|
|
C \end{center} |
|
16 |
|
|
17 |
|
C !USES: |
18 |
implicit none |
implicit none |
19 |
|
#include "SIZE.h" |
20 |
|
#include "fizhi_SIZE.h" |
21 |
|
#include "fizhi_land_SIZE.h" |
22 |
|
|
23 |
C !INPUT PARAMETERS: |
C !INPUT PARAMETERS: |
24 |
integer myThid |
integer myThid |
|
integer sNx,sNy, OLx,OLy, nSx,nSy, nPx,nPy |
|
|
integer NrPhys |
|
25 |
|
|
26 |
C !USES: |
C !USES: |
27 |
#ifdef ALLOW_MNC |
#ifdef ALLOW_MNC |
60 |
data time_dat / |
data time_dat / |
61 |
& '- ', 't ', ' ', ' ', ' ' / |
& '- ', 't ', ' ', ' ', ' ' / |
62 |
|
|
63 |
|
C Define a few "standard" or "convenience" types |
64 |
ncomb = 0 |
ncomb = 0 |
65 |
DO ihorz = 1,5 |
DO ihorz = 1,5 |
66 |
DO is = 1,3 |
DO is = 1,3 |
207 |
ENDDO |
ENDDO |
208 |
ENDDO |
ENDDO |
209 |
|
|
210 |
|
C Define a "fizhi_veg" type to be used in FIZHI_INIT_VEG() |
211 |
|
C |
212 |
|
C name 1234567890 |
213 |
|
ndim = 3 |
214 |
|
dn(1)(1:10) = 'X ' |
215 |
|
dim(1) = sNx + 2*OLx |
216 |
|
ib(1) = OLx + 1 |
217 |
|
ie(1) = OLx + sNx |
218 |
|
dn(2)(1:10) = 'Y ' |
219 |
|
dim(2) = sNy + 2*OLy |
220 |
|
ib(2) = OLy + 1 |
221 |
|
ie(2) = OLy + sNy |
222 |
|
dn(3)(1:10) = 'VegType ' |
223 |
|
dim(3) = maxtyp |
224 |
|
ib(3) = 1 |
225 |
|
ie(3) = maxtyp |
226 |
|
|
227 |
|
CALL MNC_CW_ADD_GNAME('fizhi_veg', ndim, dim, dn,ib,ie, myThid) |
228 |
|
CALL MNC_CW_ADD_VNAME('surftype', 'fizhi_veg', 4,5, myThid) |
229 |
|
CALL MNC_CW_ADD_VNAME('tilefrac', 'fizhi_veg', 4,5, myThid) |
230 |
|
|
231 |
|
|
232 |
#endif /* ALLOW_MNC */ |
#endif /* ALLOW_MNC */ |
233 |
|
|
234 |
RETURN |
RETURN |