-
Notifications
You must be signed in to change notification settings - Fork 0
/
parameter.inc.n96e
128 lines (128 loc) · 4.82 KB
/
parameter.inc.n96e
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
INTEGER NZ,NZM1,NZP1,NDIM,NX,NY
INTEGER*8 npts
PARAMETER(NZ = 99 , NZM1 = NZ-1 , NZP1 = NZ+1 )
#ifdef CFS
PARAMETER(NDIM = 1 , NX = 87, NY = 34, NPTS = NX * NY )
#else
#ifdef OASIS2
PARAMETER(NDIM = 1 , NX = 129, NY = 73, NPTS = NX * NY )
#else
#ifdef OASIS3
PARAMETER(NDIM = 1, NX = 192, NY = 115, NPTS = NX*NY)
c PARAMETER(NDIM = 1, NX = 81, NY = 52, NPTS = NX*NY)
#else
PARAMETER(NDIM = 1 , NX = 192, NY = 115, NPTS = NX * NY )
c PARAMETER(NDIM = 1, NX = 10, NY = 10, NPTS = NX*NY)
#endif /*OASIS3*/
#endif /*OASIS2*/
#endif /*CFS*/
INTEGER NVEL,NSCLR,NVP1,NSP1,NSB
PARAMETER(NVEL = 2 , NSCLR= 2 , NVP1 = NVEL+1,NSP1=NSCLR+1)
PARAMETER(NSB = 1 ) !NSCLR-2
INTEGER itermax
REAL hmixtolfrac
PARAMETER(hmixtolfrac = 0.1 , itermax = 200 )
c temporary grid
INTEGER NGRID,NZL,NZU,NZDIVmax,NZtmax,NZP1tmax,igridmax
PARAMETER(NGRID = NZ )
PARAMETER(NZL = 1 , NZU = 2 , NZDIVmax = 8 )
PARAMETER(NZtmax = NZ +(NZL+NZU)*(NZDIVmax-1), NZP1tmax=NZtmax+1)
PARAMETER(igridmax= 5)
c fluxes and forcing
INTEGER NSFLXS,NJDT,NSFLXSM1,NSFLXSP2,NDHARM
PARAMETER(NSFLXS = 9 , NJDT = 1 , NSFLXSM1 = NSFLXS-1 ,
+ NSFLXSP2 = NSFLXS+2 )
PARAMETER(NDHARM = 5 )
c ocean advection
INTEGER maxmodeadv
PARAMETER(maxmodeadv = 6 )
c richardson mixing
INTEGER MR,MRP1
PARAMETER(MR =100 , MRP1 = MR+1)
c rad/conv model
* INTEGER NPLEV,NPSAVE
* PARAMETER(NPLEV = 18 )
* PARAMETER(NPSAVE = 10 )
c output buffer
* INTEGER NDOUT,NBUFF
* PARAMETER(NDOUT = 10+NZP1)
* PARAMETER(NBUFF = NZP1*(NVp1+NSP1) + NZP1*(NVEL+NSCLR)
* + + NDOUT + 3*NZ + 5*NSFLXS )
c + + 2*NPLEV + NPSAVE ) ! to store rad/conv output
c
c Parameters for regional coupling
INTEGER NX_GLOBE,NY_GLOBE
#ifdef CFS
PARAMETER (NX_GLOBE=192, NY_GLOBE=94)
#else
#ifdef OASIS2
PARAMETER (NX_GLOBE=288, NY_GLOBE=217)
#else
#ifdef OASIS3
PARAMETER(NX_GLOBE=192, NY_GLOBE=144)
c PARAMETER (NX_GLOBE=182, NY_GLOBE=152)
#else
PARAMETER (NX_GLOBE=192, NY_GLOBE=144)
#endif /*OASIS3*/
#endif /*OASIS2*/
#endif /*CFS*/
INTEGER*8 NPTS_GLOBE
PARAMETER (NPTS_GLOBE= NX_GLOBE*NY_GLOBE)
c-----------------------------------------------------------------------
c Note in this version:
c -albedo for ocean is set to 0.06, which is used for QSW from fcomp
c or fread when rad/conv is not running:
c albocn=0.06 (in subroutine init cnsts)
c-----------------------------------------------------------------------
c
c main (permanent) grid
c
c NZ : number of layers
c NZP1 : number of grid points
c NDIM & NX & NY : dimension of the model(not used in this version)
c NVEL : number of velocity components, i.e. 2 for U and V
c NSCLR : number of scalars, i.e. T, S, and additional scalars
c NSB : number of biological scalars, used in "biocommon.inc"
c hmixtolfrac : convergence tolerance for hmix(new)-hmix(old)
c iteration in ocnstep: fraction of layer thickness hm(kmix)
c itermax : maximum number of hmix iterations (on main or temporary
c grids.
c
c temporary grid
c
c NGRID : number of grids = permanent grid + temporary grids,
c if only p-grid used: NGRID = 1, if t-grid used: NGRID = NZ
c NZL & NZU : refinement interval: it is defined on permanent-grid
c from (kmix+NZU) to (kmix-NZL)
c NZDIVmax: maximum number of fine(temporary) grid intervals
c per permanent grid interval. It is needed to dimension the
c necessary arrays; the number to be used is being read in
c as NZDIV.
c NZtmax: maximum number of layers on temporary grid,
c the number to be used in this run is read in as NZT.
c igridmax : maximum number of new grid refinements in ocntgrid
c
c fluxes and forcing
c
c NSFLXS: number of fluxes: sflux(NSFLXS,5,0:NJDT)
c NJDT : number of older flux values used to extrapolate new fluxes
c NDHARM: maximum number of harmonics used to specify forcing
c
c ocean advection
c maxmodeadv: maximum number of different modes for advection
c richardson mixing
c
c MR : dimension of "F of Ri" in ri_mix routine
c
c rad/conv model
c
c NPLEV : number of levels in atmospheric model.
c (Note: NPLEV=plev necessary in "rad.par")
c NPSAVE: number of extra atmospheric variables to be saved
c
c dimension of output array
c
c NDOUT : dimension of output array "dout", so that a number of
c scalar parameters, as well as one extra profile
c (on layer or grid) can be stored for diagnostic purposes.
c NBUFF : dimension of data array "buffer"