Namelist parameters for NHM

Update $Date: 2010-12-06 09:00:00 $

NAMMSW: mode switch for the decision basic conditions of the model
(IC = initial condition, BC = boundary condition)

name of parameter kind of parameter value setting note
mswsys(1) lower BC 0 none flux on the land surface  
1st digit is 1 computation of bulk coefficient at the land:Sommeria (Jul. 7, 2003)
1st digit is 2 computation of bulk coefficient at the land:Businger
1st digit is 3 computation of bulk coefficient at the land:Kader and Yaglom
1st digit is 4 computation of bulk coefficient at the land:Louis
1st digit is 5 computation of bulk coefficient at the land:Beljaars and Holtslag (Feb. 26, 2004)
1st digit is 9 use MRI/JMA-SIB as LSM (Oct. 7, 2005)(Jan. 30, 2006)
10th digit is 0 computation of bulk coefficient at the sea surface:Kondo (Dec. 26, 2003)
10th digit is 1 computation of bulk coefficient at the sea surface:Sommeria
10th digit is 2 computation of bulk coefficient at the sea surface:Businger
10th digit is 3 computation of bulk coefficient at the sea surface:Kader and Yaglom
10th digit is 4 computation of bulk coefficient at the sea surface:Louis
10th digit is 5 computation of bulk coefficient at the sea surface:Beljaars and Holtslag (Feb. 26, 2004)
10th digit is 6 computation of bulk coefficient at the sea surface:Kondo/ Donelan/ Fairall (COAREversion3)  
100th digit is 0 not consider stomatal resistance and time, wetness variation of wetness (Jul. 18, 2003)
100th digit is 1 consider stomatal resistance
100th digit is 2 consider stomatal resistance and time variation of wetness (Mar. 30, 2004)
1000th digit is 0 computation of roughness at the sea surface:Kondo (Dec. 26, 2003)
1000th digit is 1 computation of roughness at the sea surface:Beljaars
1000th digit is 2 computation of roughness at the sea surface:simple Kondo (Jul. 2005)
1000th digit is 3 computation of roughness at the sea surface:Kondo/ Donelan/ Fairall (COAREversion3)  
10000th digit is 1 revise the speed used computation of surface flux by free-scale (Jul. 2005)
mswsys(2) outflow lateral BC for normal component 0 no radiation condition* (Sep. 18, 2003) available when all RATIOx in NAMVAL equal 1.0
1 use Orlanski type radiation condition  
mswsys(3) horizontal spacing 0 eigen function is read from stored file (not available)
1 variable eigen function is calculated by Jacobi method
2 uniform eigen function is calculated by trigonometrical function
mswsys(4) outflow lateral BC for tangential component 0 use flux adjusted values of the outer model (Oct. 7, 2003) These are able to use in the case of all RATIOx of NAMVAL is "1.0".
1 Orlanski type  
2 extrapolation in time and space a kind of the condition of radiation
3 inner closest point value (This option is "mswsys(4)=0" before Oct. 7, 2003)
mswsys(5) lateral BC for potential temperature, turbulence kinetic energy and variables for cloud physics 0 inner closest point value  
1 Orlanski type  
2 extrapolation in time and space a kind of the condition of radiation
mswsys(6) definition of density -2 fully compressible (with map factor)  
-1 fully compressible (without map factor)  
0 use values of the reference atmosphere non-elastic model, quasi-compressive model
1 Bousinesq approximation  
mswsys(7) Bousinesq approximation 0 flux-adjusted values of the outer model (2003/09/30)
1 given by time-integration  
2 not time-integrated, using previous values (2003/09/30)
3 not time-integrated, with consideration of wind direction (2003/09/30)
mswsys(8) Coriolis parameter
(11 recommended)		 0 without Coriolis parameter  
1 vertical component (2 w sin f) considered f=(0, 0, 2ωsinφ)
2 full evaluation  
+10 1th digit is 0 : consider treatment of curvature terms
1st digit is 1-9 : omit
( for computational efficiency)		 (2003/10/03)
mswsys(9) number of iteration in pressure eqn. solver 1 no iteration in the case of "the elastic model" and in the case of no topography or having friction
3 three time iteration  
mswsys(10) dimension 0 3D  
1 2D (x-z, may not work appropriately) not available
Please, use "Nhm/Ss/Ideal/squall2d" for 2-dimensional experiments
mswsys(11) upper BC 0 free slip, rigid wall no other choice
mswsys(12) start-up procedure 1 wind given weakly, then gradually increased until it=ITGROW (this option not available)  
2 pre-existing wind and mountain (for stand-alone ; not nest)  
3 same as 2, given from a file (this option not available)
4 given by nesting (w=0 at all levels)  
5 given by nesting (w converted from w)  
6 given by nesting (w computed using continuity equation)  
7 given by nesting (w computed using continuity equation, but zeroed at lateral boundary)  
mswsys(13) ground temperature and atmospheric radiation (12 is recommended) 0 no heat and moisture fluxes
Note that sensible and latent fluxes forced in SCM are considered. 		 
1 varying in accordance with sinusoidal curve  
2 predicted  
3 same as 2, but ground inclination considered  
4 same as 2, but orographic shadows considered  
5 predicted, ground inclination and orographic shadows considered  
6 predicted, contributions to atmosphere considered, cloud amounts diagnosed * method of RSM(Sugi and Tada 1988)
7 same as 6, but ground inclination considered  
8 same as 6, but cloud water and ice used  
9 same as 8, but cloud water and ice used  
10 same as 6, but radiation scheme of GSM0103 (2005/03/31)Cloud radiation: Kitagawa(2000), Clear sky radiation: Kitagawa(2000).
11 same as 10, considering ground inclination
12 same as 10, including more detailed clear sky radiation scheme. Radiation scheme of GSM0412 (2006/04/28)Cloud Radiation: Kitagawa(2000), Clear sky radiation: Yabu et al.(2005).
13 same as 11, considering ground inclination
mswsys(14) lateral boundary condition -2 free-slip rigid wall  
-1 open in x-direction, free-slip rigid wall for y-direction  
0 open in x- and y-direction *  
1 open in x-direction, periodic in y-direction  
2 periodic in x- and y-direction  
3 periodic in x-direction, open in y-direction  
mswsys(15) buoyancy term 0 linearized, computed from potential temperature perturbation anelastic
1 not linearized, computed from potential temperature perturbation (HE-VI / HI-VI)  
2 computed from density perturbation (HE-VI / HI-VI) for the HI-VI cases (usually )
mswsys(16) wind component (IC and BC) 0 multiply density and metric term (r G1/2) *  
1 use as it is, without modification for the self-nested cases
2 use as it is in IC, multiply "r G1/2" in BC  
3 multiply "r G1/2" in IC, use as it is in BC  
mswsys(17) splitting in HE-VI
(3 or 5 recommended, 7 is recommended for coarser resolution such as dx >= 10km, 0 recommended for very high resolution such as dx < 1km)		 0 no splitting for the HI-VI cases
1 splitting gravity-wave terms (forward scheme)  
2 same as 1 (leapfrog scheme)  
3 splitting gravity-wave and momentum terms  
5 splitting momentum and potential temperature terms  
7 3 for former half of leapfrog time step and 5 for latter half of leapfrog time step (2005/10/27)(2005/10/25)
mswsys(18) cloud microphysics 3 warm rain, qc, qr, Nc, and Nr predicted/td> [correspondence]
cloud microphysics mode : val_micro = only 1 or 3 ,scheme for control saturation : ipccnd = only 10 or 12 , velocity falling raindrop : irain = only 1 (2005/01/20)
2 dry  
1 warm rain, qc and qr predicted  
0 ice phase included, qc, qr, qi, qs, qg, and Ni predicted  
-1 same as 0, but qc, qr, qi, qs, qg, Ni and Ns predicted  
-2 same as 0, but qc, qr, qi, qs, qg, Ni, Ns and Ng predicted  
-3 -2 + Nr, Nc predicted [correspondence]
cloud microphysics mode : val_micro = only 1 or 3 , scheme for control saturation : ipccnd = only 10 or 12 , velocity falling raindrop : irain = only 1 (2005/01/20)
-4 ice phase included, qc, qr, qi, qs, qg, qh, Nc, Nr, Ni, Ns, Ng, and Nh predicted reservation for the kit arrangement
-10 same as 0, but Ni not predicted  
-30 same as 0 but with formation process of graupel simplified
-40 same as 0, but Ni not predicted but with formation process of graupel simplified
-50 same as 0, but qg and Ni not predicted not include graupel
-60 same as 0, but qg and Ni not predicted, and over-cooled phase ommited (mixed phase not considered)
-100 same as mswsys(18) = 0, but use MM5(ver2)  
-101 same as mswsys(18) = -1, but use MM5(ver2)  
-102 same as mswsys(18) = -2, but use MM5(ver2)  
-110 same as mswsys(18) = -10, but use MM5(ver2)  
mswsys(19) turbulence closure model (31001 recommended) 1st digit is 0 explicit solver of vertical turbulent diffusion  
1st digit is 1 implicit solver of vertical turbulent diffusion  
1th digit is 0 computations of mixing length:l_x = l_y = l_z, ds = (dx * dy * dz) ** (1 / 3) (2003/12/26)
not available for more than 20000
10th digit is 1 computations of mixing length:l_x = l_y >> l_z, ds = dz
10th digit is 2 computations of mixing length:l_x = l_y = l_z, ds = dz
10th digit is 3 computations of mixing length:l_x = l_y = 0, ds = dz
100th digit is 0 computations of vertical mixing length (l_z) :Blackadar ,1962 (2003/12/26)
not available for more than 20000
100th digit is 1 computations of vertical mixing length (l_z) :non local PBL scheme (Sun and Chang, 1986)
100th digit is 2 computations of vertical mixing length (l_z) :non local PBL scheme (Sun and Chang, 1986) + Blackadar
1000th digit is 0 computation of PBL height = 0 [m] (2003/08/28)
1000th digit is 1 computation of PBL height:calculate from vertical virtual PT profile
1000th digit is 2 computation of PBL height:Hong and Pan(1996)
1000th digit is 3 computation of PBL height:calculate from profile of coefficient of diffusion  
10000th digit is 0 prognosis treatment of turbulent energy : Deardorff(1973) (2003/10/03)
10000th digit is 1 diagnosis treatment of turbulent energy : Deardorff(1973)
10000th digit is 2 Improved Mellor-Yamada Level2.5  
10000th digit is 3 Improved Mellor-Yamada Level3
10000th digit is 4 Mellor-Yamada Level2
100000th digit is 0 coefficient of turbulent diffusion : 0.1 (2003/12/26)
not available for more than 20000
100000th digit is 1 coefficient of turbulent diffusion:0.2
100000th digit is 2 coefficient of turbulent diffusion:0.2 to 0.1 in PBL and 0.1 above PBL
100000th digit is 3 coefficient of turbulent diffusion:0.2 in PBL and 0.1 above PBL
mswsys(20) basic equation -1 anelastic (AE), hydrostatic  
0 anelastic (AE)  
1 Elastic (HI-VI)  
2 Elastic (HE-VI)  
mswsys(21) fall-out of rain drops and graupel (23 recommended) 0 Euler scheme (2005/04/11)
1 box-Lagrangian scheme for rain drops
10 same as 1, but only for graupel
11 same as 1, but for both rain drops and graupel
20 0 + cloud ice fallout by Euler scheme (2005/04/11)(2006/03/02)
21 1 + cloud ice fallout by Euler scheme
22 10 + cloud ice fallout by Euler scheme
23 11 + cloud ice fallout by Euler scheme
30 Rain drops, snow and graupel fallout by analytical scheme (2005/03/30)
31 30 + cloud ice fallout by analytical scheme (2005/04/11)
101 1 + vertical advection of rain drops by box-Lagrangian valid for mswsys(18)==-40,0,1
110 10 + vertical advection of graupel by box-lagrangian
111 11 + vertical advection of both rain drops and graupel by box-lagrangian
121 21 + vertical advection of rain drops by box-Lagrangian
122 22 + vertical advection of both rain drops and graupel by box-lagrangian
123 23 + vertical advection of both rain drops and graupel by box-lagrangian
mswsys(22) cumulus parametrization 0 w/o cumulus parametrization  
1 cloud microphysics w/moist convective adjustment (MCA), condensation in MCA added to cloud water  
2 same as 1, but condensation in MCA added to precipitation  
3 grid scale condensation w/MCA qc and qr not predicted
4 grid scale condensation w/o parametrization qc and qr not predicted
5 same as 1, but condensation in MCA added to rain and graupel  
11 same as 1, but w/Arakawa-Schubert scheme (AS)  
13 same as 3, but w/AS  
15 Grell scheme and cloud microphysics w/moist convective adjustment (MCA) (2006/03/23)(explanation by Ohmori)(explanation by Shinpo)
20 same as 1, but w/Kain-Fritsch scheme (KF)  
mswsys(23) BC for pressure 0 without damping  
1 Rayleigh damping used in upper layers  
2 Rayleigh damping used in upper and lateral boundary layers When using cyclic condition, this switch is ignored in lateral.
mswsys(24) mean pressure and buoyancy at model top -1 mean pressure fixed  
0 fixed reference atmosphere  
1 buoyancy adjusted  
mswsys(25) : lateral boundary relaxation for U, V, W, q and qv 0 w/o damping  
1 Rayleigh damping used When using cyclic condition, this switch is ignored.
mswsys(26) mass flux through lateral boundaries 0 without adjustment  
1 with adjustment preserving total mass  
2 with adjustment tracking mean pressure of outer model  
3 with adjustment for observed total mass disapproval using it in the present condition
+10 vapor diffusion is considered in continuity equation and evaporation is also considered in flux adjustment (additional to above) (2003.09.09)
+20 vapor diffusion is considered in continuity equation and evaporation is not considered in flux adjustment (additional to above)
+30 vapor diffusion and evaporation is not considered either in continuity equation and flux adjustment (additional to above)
mswsys(27) vertical spacing 0 stretched using DZL, ZDR, IZ1, IZ2 (given in NAMGRD)  
1 arbitrary setting for scalar levels Only for Idealized Experiment: the plane of the model is given by "Z(M)" of the card of input parameter
mswsys(28) order of accuracy for advection terms 0 centered, second order, flux form  
1 upstream, first order, flux form only for scalar, momentum is fourth order
2 second order centered, advective form only for scalar, momentum is fourth order
3 upstream, third order, advective form only for scalar, momentum is fourth order
4 centered fourth order, advective form only for scalar, momentum is fourth order
30 upstream, third order, flux form reference to section2.3
31 upstream, third order, flux form same as WRF
35 scalar is third , momentum is fourth and flux form (2003.08.11)
36 scalar is third , momentum is fourth and flux form same as WRF (2003.08.11)
40 centered, fourth order, flux form reference to section 2.3
41 centered, fourth order, flux form same as WRF
50 upstream, fifth order, flux form reference to section 2.3
51 upstream, fifth order, flux form same as WRF
55 scalar is fifth , momentum is fourth and flux form (2003.08.11)
56 scalar is fifth , momentum is fourth and flux form same as WRF (2003.08.11)
4040 horizontal fourth, vertical fourth order flux form advection (2006/03/22)
4000 horizontal second, vertical fourth order flux form advection
mswsys(29) evaporation in sub-grid scale 0 not considered  
1 considered using predicted cloud amount  
mswsys(30) modified advection scheme 0 without modification  
1 modified q, water substances (qv, qc, qr, qi, qs and qg) mass is not conserve, density is considered
2 modified U, V, W, q, water substances (qv, qc, qr, qi, qs and qg) application at a maximum (turbulent flow energy not corrected now). mass ... not conserve , density ... considered
3 same as 2, but only for upper limit mass is not conserve , density is considered
4 modified U, V, W mass is not conserve , density is considered
5 same as 4, but only for upper limit mass is not conserve , density is considered
+10 (opposite each as 1st digit is 1 and over) mass is not conserve, density is not considered (Aug. 13, 2003)
+20 (opposite each as 1st digit is 1 and over) mass is conserve, density is considered (Aug. 20 and 21, 2003)

NAMCLD: parameter for cloud physical process
name of parameter type of parameter meaning default value note
arsl_fcst integer Aerosol forecast
0: not forecast
1***: forecast the AerosolAx1 variable (mixing ratio)
2***: forecast the AerosolAx2 variable (mixing ratio / concentration)
0***: not forecast the AerosolA
*1**: forecast the AerosolBx1 variable (mixing ratio)
*2**: forecast the AerosolBx2 variable (mixing ratio / concentration)
*0**: not forecast the AerosolB
**1*: forecast the AerosolCx1 variable (mixing ratio)
**2*: forecast the AerosolCx2 variable (mixing ratio / concentration)
**0*: not forecast the AerosolC
***1: forecast the AerosolDx1 variable (mixing ratio)
***2: forecast the AerosolDx2 variable (mixing ratio / concentration)
***0: not forecast the AerosolD
0 if you set iemit over 1, arsl_fcst is automatically set arsl_fcst=***1 or ***2.
A,B,C,D represent a kind of aerosol respectively (under construction now)
iarsl_cld integer interaction between cloud and aerosol
0: no-interaction
1: activated cloud nuclei are calculated from diagnostic super saturation (arsl_fcst=1***)
2: activated cloud nuclei are calculated from look-up table (arsl_fcst=2***)
0 (under construction now)
iemit integer Aerosol are emitted from moving point source
0: not emit
1: emit dry-ice (arsl_fcst=***2)
2: emit liquid carbon-dioxide (arsl_fcst=***1)
3: emit AgI (arsl_fcst=***1)(not implemented)
 0 In case of dry-ice and liquid carbon-dioxide, it is reproduced that cloud ice nucleation by sublimation of themselves
iemit=3,13 are not available (but reserved)
Aerosol are emitted from a point source (not moving)
10: not emit
11: emit dry-ice (arsl_fcst=***2)
12: emit liquid carbon-dioxide (arsl_fcst=***1)
13: emit AgI (arsl_fcst=***1) (not implemented)
Aerosol are emitted from column-like source
20: not emit
21: emit from column source from surface to specified height(arsl_fcst=***2)
Aerosol are emitted from plane-like source
30: not emit
31: emit from plane-like source(arsl_fcst=***2)
Aerosol are emitted in volume
40: not emit
41: emit in volume(arsl_fcst=***2)
emit_rate real Aerosol emission rate
emitted mass rate (kg s-1) in case of moving point or point source,
emitted mixing ratio rate (kg kg-1 s-1) in case of column-like, plane-like or volume source,
0.0333
emit_st
emit_en		 real start/end time(hour) of emitting Aerosol emit_st=3.0
emit_en=6.0		 specify by forecast time
emit_xc(n)
emit_yc(n)
emit_zc(n)
(1<=n<=32)		 real grid points of emitted aerosol
0 or positive value : distance from leftbelow (m)
negative value : grid point number (-ix,-jy,-kz)
emit_xc(1)=57760.0
emit_yc(1)=32980.0
emit_zc(1)=2600.0
emit_xc(2:32)=1.0e+30
emit_yc(2:32)=1.0e+30
emit_zc(2:32)=1.0e+30		 vertical coordinate of emission is; 0 or positive value is real height (m), negative value is grid point in the model
In column-like emit, vertical points indicate upper point of emission
Setting 32 points are available,
but in moving point only 1 point is available (emit_xc(1),emit_yc(1),emit_zc(1))
emit_xc and emit_yc is a grid point of south-west corner in case of plane-like or volume source.
emit_zc is a grid point of lower level in case of volume source.
emit_xc_ne(n)
emit_yc_ne(n)
emit_zc_ne(n)
(1<=n<=32)		 real grid points of north-east corner of emitted aerosol
0 or positive value : distance from leftbelow (m)
negative value : grid point number (-ix,-jy,-kz)
emit_xc_ne=emit_xc
emit_yc_ne=emit_yc
emit_zc_ne=emit_zc		 emit_xc_ne and emit_yc_ne is available in case of plane-like or volume sorce.
emit_zc_ne is a grid point of upper level in case of volume sorce.
and the others are same as emit_xc, emit_yc,emit_zc
if overlaps are happened when plural sorces are set, the united region of themselves are used.
emit_vel real moving velocity of point source (m s-1) 100.0 available in moving point case
emit_rad real half of longer distance in the moving point track with like a shape of 8 (m) 8000.0 available in moving point case
emit_ang real rotation of y-direction in the moving point track with like a shape of 8 (degree in anti-clockwise) -45.0 available in moving point case
emit_asp real the ratio of longer and shorter radius in the moving point track with like a shape of 8 0.1 available in moving point case
inuclei integer the structure making ice crystal
1: method of "Murakami (1990)"
2: method of "MM5"
3: method of "Huffman and Vali (1973)"
4: method of "Meyers et al. (1992)"
5: same as "inuclei = 4", but probability making ice crystal is function of updraft.
6: method of "Murakami (1999)"
7: same as "inuclei = 6", but probability making ice crystal is function of updraft.
8: same as "inuclei = 4", but time variation of ice super-saturation is diagnosed from adiabatic cooling rate and growth rate of cloud water		 6 (Reference : NPD-report Bessatsu Vol.49 "Japan Meteorological Agency NonHydrostatics Model", pp52-76)
ipccnd integer saturation adjustment(additional heat is released heat from water vapor condensation when saturation ratio is change from water-saturation to ice-saturation with the moment of freezing of all of super-cooled water droplets under iso-pressure.
10: method of " MRI", non-added heat
11: method of "MRI", added head
12: same as "ipccnd = 10", but non-control of icy saturation
20: method of "MM5", non-added heat
21: method of "MM5", added heat
30: method of "Ferrier", non-added heat
31: method of "Ferrier", added heat
40: method of "Soong and Ogura", non-added heat
41: method of "Soong and Ogura", added heat
50: method of "Tao et al." , non-added heat
51: method of "Tao et al." , non-added heat
90: scheme for part of condensation, non-added heat
91: scheme for part of condensation, added heat		 10
irain integer expression of "velocity falling raindrop"
1: exponential function
2: polynomial expression		 2
icoleff integer factor of collision-acquisition between raindrop, snow, graupel and cloud droplet
1: NHM (cloud ice crystal and cloud ice crystal: 0.1,
rain and cloud water : function of "Stokes number",
snow and cloud water : function of Stokes number,
graupel and cloud water : function of Stokes number,
snow and cloud water: function of temperature,
graupel and cloud water : 0.1, snow and hail: function of temperature,
snow and snow : 0.1,
rain and graupel : 1.0,
rain and snow : 1.0)
2: MM5 (cloud ice crystal and aloud ice crystal : 0.1, rain and cloud water : 1.0, snow and cloud water : 1.0, graupel and cloud water : 1.0,
snow and cloud ice crystal : 1.0,
graupel and cloud ice crystal : 0.05,
snow and graupel : 0.001,
snow and snow : 0.1,
rain and graupel: 1.0,
rain and snow: 1.0)		 2
t_pccnr character method of transform from cloud water to rain
'ks': type of "Kessler"
'br': type of "Berry and Reinhardt"
'rc': type of "Richard and Chaumerliac"		 ks
fac_revap real control evaporation of rain at cloud physical process 1.0 -100 < mswsys(18) ≤ 1 .and. mswsys(18) ≠ -3
available below 1.0d0
0.0 means not considered evaporation
not applied for cumulus parametrization (2003/12/08-09)
fac_sevap real control evaporation of snow at cloud physical process 1.0 -100 < mswsys(18) ≤ 0
available below 1.0d0
0.0 means not considered evaporation
not applied for cumulus parametrization (2003/12/08-09)
fac_gevap real control evaporation of hail at cloud physical process 1.0
tm_scng real factor for transform from snow to hail 160.0 (MRI Tech. Rep. vol. 28 fig. B-11-3)2006/09/02
auts1 real parameter of auto-conversion from cloud ice crystal to snow 1.0e-3(s^-1) (Reference : NPD-report Bessatsu Vol.49 "Japan Meteorological Agency NonHydrostatics Model", pp52-76)
auts2 real parameter of auto-conversion from cloud ice crystal to snow 1.0e-4(kg/kg)
var_x real broadness of droplet spectrum for t_pccnr = 'rc' and 'br' (0.25 =< var_x <= 1.0) 1.0
val_micro integer option for selecting cloud physical process - mode. select from 1 to 3. 1: MRI - mode, 2: MM5(ver2) - mode, 3: the mode defined by user 3
↓Hereafter, available in val_micro = 3. And parameters of below are ignored in val_micro = 1 or 2
qn0cw_usr real density of cloud water 1.0e+8(m^-3)
pn0r_usr real intercept parameter of cloud ice 3.0e+11(m^-4)
pn0r_usr real intercept parameter of raindrop 8.0e+6(m^-4)
dnsr_usr real density of raindrop 1.0e+3(kg/m^3)
p1r_usr real parameter for calculating velocity of falling raindrop 842.0
p2r_usr real parameter for calculating velocity of falling raindrop 0.8
autr1_usr real parameter of auto-conversion from cloud water to rain water (conversion rate) 1.0e-3(s^-1)
autr2_usr real parameter of auto-conversion from cloud water to rain water (threshold of conversion) 1.0e-4(kg/kg)
dnsci_usr real density of cloud water 0.15e+3(kg/m^3)
eii_usr real factor of collision-collection between cloud ice crystals 0.1
dispvi_usr real variance of velocity falling cloud ice crystal 0.25
pn0s_usr real slope parameter of snow 1.8e+6(m^-4)
dnss_usr real density of snow 0.84e+2(kg/m^3)
p1s_usr real parameter for calculating velocity of falling snow 17.0
p2s_usr real parameter for calculating velocity of falling snow 0.5
rads0_usr real minimum radius of snow flake 0.75e-4(m)
esi_usr real factor of collision-collection between snow and cloud ice crystal 0.1
ess_usr real factor of collision-collection between snow 0.1
pn0g_usr real slope parameter of hail 1.1e+6
dnsg_usr real density of graupel 0.3e+3(kg/m^3)
p1g_usr real parameter for calculating velocity of falling graupel 1.24e+2
p2g_usr real parameter for calculating velocity of falling graupel 0.64
egi_usr real factor of collision-collection between graupel and cloud ice crystal 0.1
egiw_usr real factor of collision-collection between graupel and cloud ice crystal at moist developing regime 1.0
cegacs_usr real factor of collision-collection between snow and graupel 1.0e-3
cimnms_usr real minimum mass of cloud ice crystal 1.0e-12(kg)
cimxms_usr real maximum mass of cloud ice crystal 1.35e-10(kg)
swmnms_usr real minimum mass of snow 1.5e-10(kg)
swmxms_usr real maximum mass of snow 0.5e-3(kg)
grmnms_usr real minimum mass of graupel 5.5e-10(kg)
grmxms_usr real maximum mass of graupel 6.0e-3(kg)
cwmnms_usr real minimum mass of cloud water 4.19e-15(kg)
cwmxms_usr real maximum mass of cloud water 6.54e-11(kg)
rnmnms_usr real minimum mass of raindrop 2.86e-10(kg)
rnmxms_usr real maximum mass of raindrop 1.13e-4(kg)
dnscnb_usr real density of the AerosolA 1.769e+3(kg/m^3)  
dnsinb_usr real density of the AerosolB 1.769e+3(kg/m^3)
dnscn_usr real density of the AerosolC 1.769e+3(kg/m^3)
dnsin_usr real density of the AerosolD 1.769e+3(kg/m^3)
dnsdr_usr real density of the dry ice 1.56e+3(kg/m^3)
cnbmnms_usr real minimum mass of the AerosolA 9.26e-25(kg)
cnbmxms_usr real maximum mass of the AerosolA 9.26e-17(kg)
inbmnms_usr real minimum mass of the AerosolB 9.26e-25(kg)
inbmxms_usr real maximum mass of the AerosolB 9.26e-17(kg)
cnmnms_usr real minimum mass of the AerosolC 9.26e-25(kg)
cnmxms_usr real maximum mass of the AerosolC 9.26e-17(kg)
inmnms_usr real minimum mass of the AerosolD 9.26e-25(kg)
inmxms_usr real maximum mass of the AerosolD 9.26e-17(kg)
drmnms_usr real minimum mass of the dry ice pellet 8.16e-10(kg)
drmxms_usr real maximum mass of the dry ice pellet 5.22e-5(kg)

NAMCLD: parameter for cumulus parametrization "Kain-Fritsch"
name of parameter type of parameter meaning default value note
cudt real time interval of KF activation [min] 5.0 
warm_rain logical treatment of the condensed water substances. usage of ice-phase is decided by warm_rain and p_qi (Reference: NPD-report Bessatsu Vol.49 "Japan Meteorological Agency NonHydrostatics Model", pp87-88) .false..false. or .true.
p_qr integer 1 0 or 1
p_qi integer 1 0 or 1
p_qs integer 1 0 or 1
nu_strt integer starting level of the search for the cloud forming 1 
nu_width integer grid number in inshore areas for decreasing nu_strt by 1 when nu_strt on the land is set 4 for zs < 500 m and 5 for zs >= 500m. Here, nu_strt in the other areas over the sea is set 1. This setting is only valid for nu_width >= 1. 0 
del_abe real fraction of left CAPE after applied KF
It's means the rest of CAPE after convection beginning. KF is applied until "(1 - del_abe) * CAPE0" is removed from CAPE0 before applied KF. The larger "del_abe" is, the larger the rest of CAPE remains and the weaker cumulus parametrization effects.		 0.1 
dtlcl_fct real Increment for temperature in updraft is calculated by vertical wind at LCL.
Increment for temperature become small by the smaller "dtlcl_fct". Therefore, occurrence of cumulus convection is reduced.		 1.0 
dtrh_fct real Increment for temperature in updraft is calculated by relative humidity at LCL.
Increment for temperature become small by the smaller "dtrh_fct". Therefore, occurrence of cumulus convection is reduced.		 0.0 (2006.08.24)
rate_kf real conversion rate of condensation or precipitation by Ogura and Cho (1973) (cond = 0) 0.03 
depth_usl real search upper limit from the surface (in hPa) of cloud forming updraft 300.0 
cond integer conversion scheme from producted condensation hydrometeor to precipitation, cloud water or clod ice.
0: the way of Ogura and Cho (1973), 1: the type of Kessler		 0 (2003/09/18)
wide_kf logical treatment of environmental U, V, PT, QV in searching cloud base
.true.: average 9 grids, .false.: not average		 .false. (2004/01/14)
w_time integer W used in KF is averaged w_time steps in searching cloud base 16
w_wide_kf logical W used in KF is averaged w_time steps in searching cloud base
.true.: average 9 grids, .false.: not average		 .false.
w_vert_ave logical treatment of environmental W used to calculate DTLCL in searching cloud base
.true.: average 850 to 500 hPa, .false.: average with W at the LCL		 .false.
kf_thresh real threshold of Kessler type autoconversion [kg/kg] 1.0d-3 (2004/02/10)
cu_lifetime_min real minimum of life time of deep convection [sec] 1800.0 
shallow_lifetime real life time of shallow convection [sec] 2400.0 
wklcl_fct real the factor for multiplication of the temperature perturbation based on the updraft velocity 1.0d0 (2009/04/20) (2010/04/14)
kfrad_var integer radius (rad [m]) of an updraft at the LCL
0: according to the updraft velocity, rad = 1000 .. 2000
1: according to the updraft velocity and the grid spacing of the model (dx [m]),
rad = kfrad_fct * 1000 * dx / 25000 .. kfrad_fct * 2000 * dx / 25000, where kfrad_fct is a factor
2: according to the pressure at the LCL (plcl [Pa]),
rad = kfrad_min for plcl > kfrad_plcl_l,
rad = kfrad_fct * 1000 * dx / 25000 .. kfrad_fct * 2000 * dx / 25000 for plcl <= kfrad_plcl_h,
rad is interpolated between the values at kfrad_plcl_h and kfrad_plcl_l for kfrad_plcl_h < plcl <= kfrad_plcl_l. 		0
kfrad_fct real factor for the radius of updaft at the LCL if kfrad_var is 1 or 2 5.0d0
kfrad_min real minimum value of the radius [m] of updraft at the LCL if kfrad_var is 2 500.0d0
kfrad_plcl_h real radius of updraft is kfrad_fct * 1000 * dx / 25000 .. kfrad_fct * 2000 * dx / 25000 for the pressure of the LCL is smaller than kfrad_plcl_h [Pa] if kfrad_var is 2 800.0d2
kfrad_plcl_l real radius of updraft is kfrad_min for the pressure of the LCL is larger than kfrad_plcl_l [Pa] if kfrad_var is 2 950.0d2

NAMCLD: parameter for cumulus parametrization "Grell"
name of parameter type of parameter meaning default note
lambda_grell real percentage of entrainment of ascending current 1.0d-4 (1/m) (Mar. 23, 2006)

NAMCLD: list of "NAMELIST" parameter for moist convection adjustment
name of parameter type of parameter meaning default note
ice_adj integer 0: not consider ice phase by Tetens-expression
1: consider ice phase by Tetens-expression		 1  
dtmp_adj real Temperature treated as completely ice in Tetens-expression
(Treated as ice below "-dtmp_adj" centigrade, treated as water above 0 centigrade)		 15.0(deg)
dtcons_adj real time scale of convection (for grid spacing enough big)
instability is removed in this time scale		 300.(sec)
dxmax_adj,
dxmin_adj		 real A instability of a grid which has greater interval than dxmax_adj is removed dtcons_adj in time. A tendency is set to multiply tendency to dt/dtcons_adj. A instability of a grid which has smaller interval than dxmax_adj has tendency = 0. If a gird spacing is during dxmax_adj and dxmin_adj, tendency is divided by according to its size. dxmax_adj = 20000(m), dxmin_adj = 1000(m)
indxm_adj integer limitation for number of change between stability/instability layers 3
bcrit_adj,
bdrylp_adj		 real Using for decision of stability or instability of layers.
RH = "bdrylp_adj", then dry-adiabatic is threshold RH = "bcrit_adj", then moist-adiabatic is threshold RH = during "bdrylp_adj" and "bcrit_adj", then threshold is divide by rate between bdrylp_adj / bcrit_adj.		 bcrit_adj = 1.0, bdrylp = 0.5
drymp_adj,
wetmp_adj		 real adjustment parameters for decision of stability or instability of layers.
The threshold is that dry-adiabatic and moist-adiabatic are multiplied by "drymp_adj", "wetmp_adj"		 drymp_adj = 1.0,
wetmp_adj = 1.0
cbcal_adj integer set the limit of cloud base
0 : Cloud base lower limit is set to the layer just above the LCL of the lowest air-mass.
1 : Cloud base lower limit is set to the layer just above the minimum LCL of air-mass below zlcl_adj.
2 : Cloud base lower limit is set to the layer just above the PBL height.
3 : Cloud base lower limit is set to the layer +1 just above the PBL height.
0
zlcl_adj real the height of calculating LCL when cbcal_adj=1 2500.0 [m]
zmbtmx_adj real When the limit of cloudbase which calculates above mentioned method is smaller than zmbtmx_adj, the limit is replaced to zmbtmx_adj 250.0

NAMPAR: basic parameter for the model (for example : time step, and so on)
name of parameter type of parameter mean default note
itst integer step number of the start of time integration 1 restart for more than 2 cases
itend integer step number of the end of time integration    
itout integer step interval of output   not used
itchk integer step interval of monitoring in standard output (log file)    
dt real time step [sec] 30 adjust accordingly horizontal and vertical resolution of the model
it's usually set 3 times as big as horizontal resolution (km)
dt_base real (previously used) basic dt of a long term run, when using a temporal short dt in the restart run   Remarks: 1) itst and itend are these under the assumption of the short dt from the initial. 2) the short dt should be the divisor of the basic dt. 3) output intervals for the restart and the forecasts are accounted in the unit of the dt_base, NOT SHORT DT, so they should be the same as the previous run.
dx real grid spacing in x-direction [m]    
dy real grid spacing in y-direction [m]    
dz real grid spacing in z-direction [m]   accorded "dzr" for variable cases
ptrf real reference value of q 300.0 (usually) 300.0
presrf real reference pressure [Pa] 100000.0 disapproval except 100000.0
idfi integer whether switch for digital filter 0 0: not use, 1: use
output_node integer the setting about output node 0 0: no output node
1: one node parted for output as output node
n(>=2) : 1 node parted for output node, n - 1 nodes do nothing, other for calc. to avoid prime number of calc. nodes
save_memory integer the setting about quantity using memory and mfbd 1 1: normal
0: not deallocate the already allocated array, so that overhead of allocation is reduced. But it need more physical memory.
-1: add 0, boundary file (mfbd : connected to unit number 21) is kept on physical memory. IT need more and more physical memory.
Using save_memory <= 0, it is possible that the bugs are hiden. Therefore, in developping, save_memory = 1 is recommended.
use_sx_pauto integer use special code for NEC-SX 0 0: no use,
1: standard output is slightly accelerated on NEC-SX
it2utc integer display UTC and true sun time on standard output 0 0: not output
1: display
(07/10/26)
scm_mode integer Switch for use of Single Column Model(SCM) 0 0: Not SCM, 1: Run as SCM
exact_log integer Switch so that stdout of model coincide exactly with different number of PI process 0 0: Not concide, 1: Concide
sw_div_order integer Order of the discretization of the divergence term in the pressure equation and the acoustic filter 22 1st digit : Ordef of the discretization in the pressure equation, 10th digit : Order of the discretization in the acoustic filter
2 : 2nd order of discretization, 4: 4th order of discretization

NAMFIL: parameter for controlling input and output of model
.
name of parameter type of parameter meaning default note
out_zspairf_notbaseit0 integer switch to output ZS and PAIRF at it = 0 which is not the basetime. 1 0: Not output. 1: Output.
fake_type2_auto integer switch to change type2 automatically when validtime is former than basetime 1 0: Not changed automatically. 1: Changed automatically.
ibase_fake integer dislocate basetime by ibase_fake 0 if validtime is former than basetime, type2(1:2) is change 'AF' for output. (except fcst_const.nus)
(06/09/29)
output_ninfo integer whether input "INFO" record of "NuSDaS"
0: not output "INFO"
1: output "INFO"
 1  
imt_rf integer the input device number of reference data -1 imt_ini is used when imt_rf = -1
fname_rf character(20) the input file name of reference data blank the data is read from the input device number imt_rf when fname_rf is blank
imt_ini integer the input device number of initial number data 23  
fname_ini character(20) the input file name of initial number data blank the data is read from the input device number imt_ini when fname_ini is blank
imt_ex integer the input device number of boundary number data 22  
fname_ex character(20) the input file name of boundary number data blank the data is read from the input device number imt_ex when fname_ex is blank
imt_ptgrd integer the input device number of file of "PTGRD" 25  
fname_ptgrd character(20) the input file name of file "PTGRD" blank the data is read from the input device number imt_ptgrd when fname_ptgrd is blank
imt_hm integer the input device number of topographic file 29  
fname_hm character(20) the input file name of file topographic file blank the data is read from the input device number imt_hm when fname_hm is blank
imt_sst integer the input device number of sst (sea surface temperature) file 53  
fname_sst character(20) the input file name of file sst file blank the data is read from the input device number imt_sst when fname_sst is blank
imt_gsm9603_veg integer the input device number of MJ-SiB constants file 32  
fname_gsm9603_veg character(20) the input file name of MJ-SiB constants file blank the data is read from the input device number imt_gsm9603_veg when fname_gsm9603_veg is blank
imt_gsm9603_add integer the input device number of added MJ-SiB constants file -i if this file is used, the number is set to 38
fname_gsm9603_add character(20) the input file name of MJ-SiB constants file blank the data is read from the input device number imt_gsm9603_add when fname_gsm9603_add is blank
imt_glcc_sib integer the input device number of the SiB vegetation file from GLCC -1  
fname_glcc_sib character(20) the input file name of the SiB vegetation file from GLCC blank the data is read from the input device number imt_glcc_sib when fname_glcc_sib is blank
imt_sibini integer the input device number of the MJ-SiB initial file from the Offline MJ-SiB climatology or an Outer NHM with MJ-SiB -1  
fname_sibini character(40) the input file name of the MJ-SiB initial file from the Offline MJ-SiB climatology or an Outer NHM with MJ-SiB blank the data is read from the input device number imt_sibini when fname_sibini is blank
imt_ab integer the input device number of file which contains "a, b" of eta level -1  
fname_ab character(20) the input file name of file which contains "a, b" of eta level blank the data is read from the input device number imt_ab when fname_ab is blank
imt_band integer the input device number of file calculating radiation 30 for "mswsys(13)= from 6 to 9" cases
fname_band character(20) the input file name of file calculating radiation blank the data is read from the input device number imt_band when fname_band is blank
imt_band_gsm0103 integer the input device number of file calculating radiation 31 for "mswsys(13)=10 or 11"cases(2005/03/31)
fname_band_gsm0103 character(20) the input file name of file calculating radiation blank the data is read from the input device number imt_band_gsm0103 when fname_band_gsm0103 is blank
imt_band_2004 integer the input device number of file calculating radiation 33 for "mswsys(13)=12 or 13"cases(2006/04/28)
fname_band_2004 character(20) the input file name of file calculating radiation blank the data is read from the input device number imt_band_2004 when fname_band_2004 is blank
imt_ozon_2004 integer the input device number of monthly 3D ozone climate 34 for "mswsys(13)>= 6 and ozone == 1"cases(2006/04/28)
fname_ozon_2004 character(20) the input file name of monthly 3D ozone climate blank the data is read from the input device number imt_ozon_2004 when fname_ozon_2004 is blank
imt_aer_2d integer the input device number of monthly 2D aersol cliamte(optical depth) 37 for "mswsys(13)>= 10 and aersl == 2"cases(2008/06/11)
fname_aer_2d character(40) the input file name of monthly 2D aersol cliamte(optical depth) blank the data is read from the input device number imt_aer_2d when fname_aer_2d is blank
imt_bd integer the input device number of working file 21  
imt_ev integer the input or output device number of file of the proper function 19  
imt_time integer the output device number of file of elapse time 90  
fname_time character(20) the output file name of file of elapse time blank the elapse time profile is outputted to the device number imt_time when fname_time is blank
imt_mri integer the output device number by "MRI" format   not inputted for minus value cases
fname_mri character(20) the output file name of "MRI" format blank the data is outputted to the device number imt_mri when fname_mri is blank
istr_mri integer interval outputted by "MRI" format (stepping unit)    
imt_mri2 integer the output device number by "MRI" format 2   not inputted for minus value cases
istr_mri2 integer interval outputted by "MRI" format 2 (stepping unit)    
imt_mri3 integer the output device number by "MRI" format 3   not inputted for minus value cases
istr_mri3 integer interval outputted by "MRI" format 3 (stepping unit)    
imt_mri4 integer the output device number by "MRI" format 4   not inputted for minus value cases
istr_mri4 integer interval outputted by "MRI" format 4 (stepping unit)    
imt_mri_sib integer the output device number of the SIB monitor of MRI format   not inputted for minus value cases(2006/01/30)(Mar. 01, 2006)
fname_mri_sib character(20) the output file name of the SIB monitor for MRI format blank the data is outputted to the device number imt_mri_sib when fname_mri_sib is blank
istr_mri_sib integer interval outputted of the SIB monitor for MRI (stepping unit)   (2006/01/30)(Mar. 01, 2006)
imt_heat_monitor integer control outputting to the monitor of potential temperature tendency  not inputted for minus value cases(2006/11/01)
fname_heat_monitor character(20) the output file name of the monitor of potential temperature tendency blank the data is outputted to the device number imt_heat_monitor when fname_heat_monitor is blank
imt_cf integer control outputting the file of constants   *1
yrev_cf integer order to "y"direction about the file of constants   *2
type1_cf character a type1 of the file of constants for NuSDaS _NHMLMLY  
type2_cf character a type2 of the file of constants for NuSDaS CCSV  
type3_cf character a type3 of the file of constants for NuSDaS STD1  
member_cf character member of the file of constants for NuSDaS speace  
imt_sf integer control outputting the file about forecasting value on the surface   *1
istr_sf integer interval outputted for the file about forecasting value on the surface (stepping unit)    
istr_sf_st integer the step started the output file about forecasting value on the surface 0 (06/08/25)
istr_sf_en integer the step ended the output file about forecasting value on the surface itend
yrev_sf integer order of the y-direction's line of the earth's surface forecast data file   *2
type1_sf character a type1 of the file about forecasting value on the surface for NuSDaS _NHMLMLY  
type2_sf character a type2 of the file about forecasting value on the surface for NuSDaS FCSV  
type3_sf character a type3 of the file about forecasting value on the surface for NuSDaS STD1  
member_sf character member of forecasting value on the surface for NuSDaS space  
imt_pb integer control outputting the file about forecasting value on the isobaric surface   *1
istr_pb integer interval outputted for the file about forecasting value on the isobaric surface (stepping unit)    
istr_pb_st integer the step started the output file about forecasting value on the isobaric surface 0  
istr_pb_en integer the step ended the output file about forecasting value on the isobaric surface itend  
yrev_pb integer order of the y-direction's line of the isobaric surface forecast data file   *2
type1_pb character a type1 of the file about forecasting value on the isobaric surface _NHMLMPP  
type2_pb character a type2 of the file about forecasting value on the isobaric surface FCSV  
type3_pb character a type3 of the file about forecasting value on the isobaric surface STD1  
member_pb character member of the file about forecasting value on the isobaric surface space  
imt_zb integer control outputting the file about forecasting value on the model plane   *1
istr_zb integer interval outputted for the file about forecasting value on the model plane (stepping unit)    
istr_zb_st integer the step started the output file about forecasting value on the model 0  
istr_zb_en integer the step ended the output file about forecasting value on the model itend  
yrev_zb integer order of the y-direction's line of the model's surface forecast data file   *2
type1_zb character a type1 of the file of forecasting value on the "NuSDaS"-model plane _NHMLMZS  
type2_zb character a type2 of the file of forecasting value on the "NuSDaS"-model plane FCSV  
type3_zb character a type3 of the file of forecasting value on the "NuSDaS"-model plane STD1  
member_zb character member of the file of forecasting value on the "NuSDaS"-model plane blank  
imt_gs integer control outputting the guess-file for QC   *1
istr_gs integer interval outputted for the guess-file for QC(stepping unit)    
istr_gs_st integer the step started the output guess-file for QC 0  
istr_gs_en integer the step ended the output guess-file for QC itend  
yrev_gs integer order of the y-direction's line of the guess file for GC   *2
type1_gs character a type1 of the file about guess for QC _NHMLMZS  
type2_gs character a type2 of the file about guess for QC GSSV  
type3_gs character a type3 of the file about guess for QC STD1  
member_gs character member of the file about guess for QC blank  
imt_ln integer control outputting about fcst_land.nus   (06/05/29)
imt_ln:*1
yrev_ln:*2
istr_ln integer interval outputted for fcst_land.nus(stepping unit)  
istr_ln_st integer the step started the output for fcst_land.nus 0
istr_ln_en integer the step ended the output for fcst_land.nus itend
yrev_ln integer >order of the y-direction's line of fcst_land.nus  
type1_ln character a type1 of "fcst_land.nus" _NHMLMLY
type2_ln character a type2 of "fcst_land.nus" FCSV
type3_ln character a type3 of "fcst_land.nus" LND1
member_ln character member of fcst_land.nus blank
imt_2m integer control outputting the "NuSDas/MRI" physical monitor(2-dimension)   *1
"MRI" output is valid for imt_2m=imt_mri*.
istr_2m integer interval outputted the "NuSDas" physical monitor(2-dimension) (stepping unit)   This setting is neglected for imt_2m=imt_mri*.
istr_2m_st integer the step started the output for "NuSDas" physical monitor(2-dimension) 0   
istr_2m_en integer the step ended the output for "NuSDas" physical monitor(2-dimension) itend  
yrev_2m integer order of the y-direction's line of "NuSDas" physical monitor(2-dimension)   *2
type1_2m character a type1 of 2D physics monitor of NuSDaS _NHMLMLY  
type2_2m character a type2 of 2D physics monitor of NuSDaS FCSV  
type3_2m character a type3 of 2D physics monitor of NuSDaS 2D_1  
member_2m character member of "NuSDaS"physical monitor(2-dimension) blank  
imt_3m integer control outputting the "NuSDas/MRI" physical monitor(3-dimension)   *1
"MRI" output is valid for imt_3m=imt_mri*.
istr_3m integer interval outputted the "NuSDas" physical monitor(3-dimension) (stepping unit)   This setting is neglected for imt_3m=imt_mri*.
istr_3m_st integer the step started the output for "NuSDas" physical monitor(3-dimension) 0  
istr_3m_en integer the step ended the output for "NuSDas" physical monitor(2-dimension) itend  
kdd_3m integer KDD for "MRI" 3-dimentional physical monitor -1 kdd_3m is automatically set using KDD*(1)-KDD*(20) for kdd_3m=-1.
yrev_3m integer order of the y-direction's line of "NuSDas" physical monitor(3-dimension)   *2
type1_3m character a type1 of 3D physics monitor of NuSDaS _NHMLMZS  
type2_3m character a type2 of 3D physics monitor of NuSDaS FCSV  
type3_3m character a type3 of 3D physics monitor of NuSDaS 3D_1  
member_3m character member of "NuSDas" physical monitor(3-dimension)NuSDaS blank  
imt_sb integer control outputting the "NuSDas" SiB monitor   (06/10/10)
imt_sb:*1
yrev_sb:*2
istr_sb integer output interval of NuSDaS SiB monitor (in time-step)  
istr_sb_st integer start step number for output of NuSDaS SiB monitor 0
istr_sb_en integer end step number for output of NuSDaS SiB monitor itend
yrev_sb integer order of the y-direction's line of "NuSDas" SiB monitor  
type1_sb character a typ1 of SiB monitor of NuSDaS _NHMLMLY
type2_sb character a type2 of SiB monitor of NuSDaS FCSV
type3_sb character a type3 of SiB monitor of NuSDaS SIB1
member_sb character member of NuSDaS SiB monitor blank
imt_eb integer output control for NuSDaS ETA-plane output   (07/04/28)
istr_eb integer step interval of NuSDaS ETA-plane output  
istr_eb_st integer start step number of NuSDaS ETA^plane output 0
istr_eb_en integer end step number of NuSDaS ETA-plane output itend
yrev_eb integer the order of y-direction's line of NuSDaS ETA-plane output  
type1_eb character type1 of NuSDaS ETA-plane output _NHMLMET
type2_eb character type2 of NuSDaS ETA-plane output FCSV
type3_eb character type3 of NuSDaS ETA-plane output STD1
member_eb character menber of NuSDaS ETA-plane output  
imt_ab integer unit number of record of coefficient "a" and "b" for ETA-plane
The file is fortran sequential. In 1st record, nz_eta + 1(4byte integer), aa(0:nz_eta)(single float), bb(0:nz_eta)(single float) are recorded. 		 
imt_dfi integer unit number of digital filter output    
pout real output levels on isobaric plane(hPa) 950., 900., 850., 800., 700., 500. available by 100th-layer
not use below 0.01hPa
elem_2m_usr character(6) output elements for physical monitor(2-dimension)   (06/12/25) (07/01/04)
*3
elem_3m_usr character(6) output elements for physical monitor(3-dimension)  
imt_zb_half integer     reserved
istr_zb_half integer     reserved
imt_zb_full integer     reserved
istr_zb_full integer     reserved
imt_rs_out integer output control of restart file -1 imt_rs_out ≤ 0 : do not output restart file
0 < imt_rs_out < 100 : Output by NuSDaS format
100 < imt_rs_out < 200 : Output by 8 byte binary format from rank 00 of MPI processes
200 < imt_rs_out < 300 : Output by 8 byte binary from each rank of MPI processes
300 < imt_rs_out < 400 : Output by NuSDaS format. Nusdas definition file are automatically created.
400 < imt_rs_out < 500 : Output with a text format suitable for SCM.
500 < imt_rs_out < 600 : Output by netCDF format. Installing the netCDF library is required to use this option. In order to activate the function in the model, you should add the following into your config file before compiling the model. 
  INC_NCDF=/usr/include/netcdf 
  LIB_NCDF="-L/usr/lib -lnetcdf"
yrev_rs integer the order of y-direction's line of NuSDaS restart file output  *2
type1_rs character type1 of NuSDaS restart output _NHMLMZS 
type2_rs character type2 of NuSDaS restart output GSSV 
type3_rs character type3 of NuSDaS restart output REST 
member_rs character member of NuSDaS restart output blank 
istr_rs_outinteger output interval for restart file 1  
imt_rs_in integer input control for restart -1 imt_rs_in ≤0 : do not read restart file
0 < imt_rs_out < 100 : Input by NuSDaS format
100 < imt_rs_out < 200 : Input by 8 byte binary format from rank 00 of MPI processes
200 < imt_rs_out < 300 : Input by 8 byte binary from each rank of MPI processes
300 < imt_rs_out < 400 : Input by NuSDaS format
sw_rest_packing integer control of packing of restart file when the file is in NuSDaS format 1 1: Output by original precisions of prognostic variables
2: packin by 2UPC
rest_bufsize_max_user integer Buffer size used in restart process (Unit:byte). Communication of the data and output unit is decided by this value. When the value is too small, it is automatically reset by model the number of grids. When the value is too large, it is adjusted to 2GB. When NuSDaS output is not used, same value of this parameter should be used before and after restart. 200000000  
sw_rest_it integer add indicator for restart file by output time-step number -1 >below 0: overwrite restart file,
above 0: write restart file in each files named by time-step.
rest_dname character(256) If you want to specify the output directory for restart file, use this option. You must prepare output directory before model run.   Not specified, output current directory on model running.
232 character is limit for directory name because it's added restart filename.(07/10/10)
nusdef_ver integer version number of NuSDaS for output of NuSDaS definition file 11 (07/10/02)
nusdef_path integer path of NuSDaS definition file
0: path = nwp_path_s (output in each validtime)
1: path = nwp_path_bs (output in 1basetime file except model-plane file)		 0
fake_step_as_min integer An option to fake 1 time step as 1 minute (because NuSDaS does not support output in second unit)
0: output normally
1: fake 1 time step as 1 minutes in output		 0
nx0out integer number of grids to clip output data (x direction) *4
ny0out integer number of grids to clip output data (y direction) *4
out_xsft integer index of start grid point to clip output data (x direction) *4
out_ysft integer index of start grid point to clip output data (y direction) *4

*1:
imt_XX<0, not output
0<imt_XX<100, output NuSDaS format (need to prepare NuSDaS definition file)
100<imt_XX<200, output 4byte binary file (fortran sequential) named by following name.
At the same time, output GrADS ctl in same name except filename extension (.ctl)
300<imt_XX<400, "mod(imt_XX, 100)" is set to NuSDaS Root Directory (NRD) and NuSDaS definition file is automatically generated on there.
and then output on NuSDaS format(07/10/02)
(In case of ibase_fake /= 0, "mod(imt_XX, 100) + 20" is set to NuSDaS Root Directory (NRD). And then NuSDaS definition file is generated for FT < 0.)
400<imt_XX<500, output with format for SCM, in which nz * time values of each element are outputted to separated files in ASCII.

*2:
In case of yrev_XX=0, array data is output from south to north in y-direction.
In case of yreb_XX=1, array data is output from north to south in y-direction.

*3:
In case of elem_2m_usr='_RSDB ', '_RLUT ', _RSDB and _RLUT are outputted to the 2-dimensional physical monitor In case of elem_3m_usr='_RADPT', '_RSHPT', _RADPT and _RSHPT are outputted to the 3-dimensional physical monitor if the registration(modification of the source code, update of make file dependences, preparation of NuSDaS definition file) of the output elements are finished in advance.
*4:
When clipping output data with NuSDaS format, the head character of type 1 is converted to 'C' automatically (such as type1 = 'CNHMLMLY'). NuSDaS definition file is automatically produced with the new type1.

NAMGRD: information for grid setting
name of parameter type of parameter meaning default note
nx0 integer grid number in x-direction    
ny0 integer grid number in y-direction    
nz0 integer grid number in z-direction    
ngm integer number of layer in soil 4 it's fixed now
dxl real grid spacing in x-direction around the west side   for variable horizontal grid spacing
dxr real grid spacing in x-direction around the east side   for variable horizontal grid spacing
ix1 integer the grid starting dx   for variable horizontal grid spacing
ix2 integer the grid ending dx   for variable horizontal grid spacing
dyl real grid spacing in x-direction around the south side   for variable horizontal grid spacing
dyr real grid spacing in x-direction around the north side   for variable horizontal grid spacing
iy1 integer the grid starting dy   for variable horizontal grid spacing
iy2 integer the grid starting dy   for variable horizontal grid spacing
dzl real layer spacing in z-direction in low level   for variable vertical layer spacing
dzr real layer spacing in z-direction in upper level   for variable vertical layer spacing
iz1 integer the grid starting dz   for variable vertical layer spacing
iz2 integer the grid ending dz   for variable vertical layer spacing
h_depth(3) integer width of halo region (1): x-direction, (2): y-direction 1 1 is only available now. (3) is reserved.
b_depth(3) integer width of physical boundary. (1): x-direction, (2): y-direction 1 1 is only available now. (3) is reserved.
x_num integer MPI process number in x-direction for 2-dimensional division of MPI 1 Not specified, using 1-dimensional division in y-direction. Specified negative value, 2-dimensional division is automatically selected(2006/08/18)
y_num integer MPI process number in y-direction for 2-dimensional division of MPI   (2003/11/07)
z_num integer MPI process number in z-direction for 3-dimensional division of MPI 1 1 is only available now.
vctrans_type integer type of translation function (hybrid vertical coordinate) 0 (1st comment)(2nd comment)
zl_vctrans real zl decide for translation function"f" 1000.d0
zh_vctrans real zh decide for translation function"f" 9000.d0
n_vctrans integer n decide for translation function"f" 1

NAMVAL: the rest of control parameters
name of parameter type of parameter meaning default note
ratioi real weighting parameter at inflow boundary (except wind) 1.0 set to 0.5-1.0
ratioo real same as RATIOI but at outflow boundary 1.0 set to 0.0-1.0
ratio2 real same as RATIOO (gravity wave) 1.0 set to 0.0-1.0
ruvni real same as RATIOI, but for wind   set to 0.5-1.0
ruvno real same as RATIOO, but for wind   set to 0.0-1.0
ruvn2 real same as RUNVNO (gravity wave)   set to 0.0-1.0
fnltr real factor of damping coefficient in boundary region 0.0 set to 0.0-1.0
idifx integer width of lateral boundary relaxation sponge layers 24  
difnl real coefficient for non-linear numerical damping -2400.0 DIFNL>0: t=DIFNL*DT, DIFNL<0: t=|DIFNL|
rtnl_uv real renormalization factor of U and V for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for U and V
rtnl_w real renormalization factor of W for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for W
rtnl_pt real renormalization factor of PT for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for PT
rtnl_qv real renormalization factor of QV for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for QV
rtnl_qc real renormalization factor of QC for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for QC
rtnl_qr real renormalization factor of QR for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QR
rtnl_qci real renormalization factor of QCI for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for QCI
rtnl_qs real renormalization factor of QS for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QS
rtnl_qg real renormalization factor of QG for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QG
rtnl_qh real renormalization factor of QH for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QH
rtnl_qncw real renormalization factor of QNCW for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QNCW
rtnl_qnr real renormalization factor of QNR for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QNR
rtnl_qnci real renormalization factor of QNCI for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QNCI
rtnl_qns real renormalization factor of QNS for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QNS
rtnl_qng real renormalization factor of QNG for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QNG
rtnl_qnh real renormalization factor of QNH for non-linear numerical dumping 0.0 0.0: no non-linear numerical dumping for QNH
rtnl_etrb real renormalization factor of ETRB for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for ETRB
rtnl_tsq real renormalization factor of TSQ for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for TSQ
rtnl_qsq real renormalization factor of QSQ for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for QSQ
rtnl_cov real renormalization factor of COV for non-linear numerical dumping 1.0 0.0: no non-linear numerical dumping for COV
dif2d real coefficient for 4th order linear numerical damping -1200.0 same branches available as DIFNL
rt2d_uv real renormalization factor of U and V for 4th order linear numerical dumping 1.0 0.0: no 4th order linear numerical dumping for U and V
rt2d_w real renormalization factor of W for 4th order linear numerical dumping 1.0 0.0: no 4th order linear numerical dumping for W
rt2d_pt real renormalization factor of PT for 4th order linear numerical dumping 1.0 0.0: no 4th order linear numerical dumping for PT
rt2d_qv real renormalization factor of QV for 4th order linear numerical dumping 1.0 00.0: no 4th order linear numerical dumping for QV
rt2d_qc real renormalization factor of QC for 4th order linear numerical dumping 1.0 0.0: no 4th order linear numerical dumping for QC
rt2d_qr real renormalization factor of QR for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QR
rt2d_qci real renormalization factor of QCI for 4th order linear numerical dumping 1.0 0.0: no 4th order linear numerical dumping for QCI
rt2d_qs real renormalization factor of QS for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QS
rt2d_qg real renormalization factor of QG for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QG
rt2d_qh real renormalization factor of QH for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QH
rt2d_qncw real renormalization factor of QNCW for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QNCW
rt2d_qnr real renormalization factor of QNR for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QNR
rt2d_qnci real renormalization factor of QNCI for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QNCI
rt2d_qns real renormalization factor of QNS for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QNS
rt2d_qng real renormalization factor of QNG for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QNG
rt2d_qnh real renormalization factor of QNH for 4th order linear numerical dumping 0.0 0.0: no 4th order linear numerical dumping for QNH
rt2d_etrb real renormalization factor of ETRB for 4th order linear numerical dumping 1.0 0.0: no 4th order linear numerical dumping for ETRB
rt2d_tsq real renormalization factor of TSQ for 4th order linear numerical dumping 1.0 0.0: no 4th order linear numerical dumping for TSQ
rt2d_qsq real renormalization factor of QSQ for 4th order linear numerical dumping 1.0 0.0: no 4th order linear numerical dumping for QSQ
rt2d_cov real renormalization factor of COV for 4th order linear numerical dumping 1.0 0.0: no 4th order linear numerical dumping for COV
diftg real coefficient for Targeted Moisture Diffusion 0.0 same branches available as DIFNL (2004/01/21)
astfc real coefficient for Asselin time filter 0.2  
stdlon real standard longitude 140.0 In some map projection, need to set appropriate value
stdlat real standard latitude 60.0 In some map projection, need to set appropriate value
kzdst real start index for upper Rayleigh damping layer 32 About nz0 * 0.8 is appropriate value.
It is appropriate that the depth of Rayleigh friction is 1.5 times as thick as a vertical wave length of inertial gravity wave.
kzden real end index for upper Rayleigh damping layer 40 equal to nz0
rldmpx real coefficient for lateral boundary relaxation -2400.0 "1/e-folding time" is "rldmpx" * dt in positive value, absolute value of "rldmpx" in negative value.
rldmpz real coefficient for upper boundary layer -2400.0 "1/e-folding time" is "rldmpz" * dt in positive value, absolute value of "rldmpz" in negative value.
rldmpo real coefficient for whole domain Rayleigh damping 0.0 not used usually
ptgrds real sea surface potential temperature[K] 288.3 available at mswsys(13)=1
ptgrdr real amplitude for diurnal change of ground potential temperature[K] 0.0 available at mswsys(13)=1
ptgrdl real ground surface potential temperature[K] 0.0  
itgrow integer end time step for wind / mountain grow initiation 0  
ubias real bias for u 0.0  
vbias real bias for v 0.0  
itsst integer start time step of elastic equation 0 Anelastic calculation until it<itsst
eover real coefficient for implicit treatment for HI-VI 0.5  
wthrt real threshold of upward motion (m/s) 2.0 diffusion is applied "omw(ix,jy,kz)>wthrt" (2004/01/21)
negative value: WTHRT is automatically set 2 ~ 20 [m/s] according to the horizontal resolution
krn_thrt real threshold of kuran number (m/s) 2.0 diffusion is applied "kuran number>wthrt"
vdv_wthrt_amp real amplitude to apply 1st order upwind scheme for vertical advection 2.0 1st order upwind scheme is used when "abs(omw(ix,jy,kz)) >wthrt * vadv_wthrt_amp"
0: means out of using 1st order upwind scheme for vertical advection
tkf_wthrt_amp real amplitude to increase the effect of K-F scheme 1.0 the effect of K-F scheme is increased when "abs(wdiag_sc(ix,jy,kz)) >wthrt * tkf_wthrt_amp"
0: means out of using the amplitude to increase the effect of K-F scheme
cfl_sound_max real safety coefficient for CFL conditions of acoustic wave 0.7 available at HE-VI (mswsys(20)=2).
Short time step is calculated using CFL_SOUND_MAX.
Not specified 0.7 is used. (2004/04/22)
sw_qsmodify integer qcorrection of qs, qr, qh 1 0: 0 at lateral boundaries, correction by chkmn0 for internal area.
1: correction by column at lateral boundaries, chkmn0 for internal area.
2: correction by column for all area
(2005/12/15)(2005/12/08)
divdmp_h real Divergence damping coefficient for horizontal direction 0.06 (2006/08/23)
divdmp_v real Divergence damping coefficient for vertical direction 0.05
tin_init real(the number of underground layers) the initial values of soil temperature for the standalone mode (including SCM)(K)
   If tin_init is not specifiled, the initial values are generated assuming that potential temperature at the lowest atmospheric layer is identical to the one at the surface and soil temperature is vertically uniform.

NAMNST: parameter for nesting
name of parameter type of parameter meaning default value note
KTSTO integer start time of nesting file 0 usually 0
KTENO integer end time of nesting file    
KTDTO integer time interval of nesting file    
DTRATIO real unit of nesting file [sec] 3600.0 time unit for KTSTO, KTENO
ALPHA real ratio of weighting parameter at variational calculus (VC) 0.5  
ITRMX integer maximum iteration number for successive over relaxation in VC 20000  
RLXCON real minimum to stop the iteration 1.0e-4  
OVERLX real coefficient for SOR(Successive Over Relaxation) 1.8  
IO_INIT_W integer whether contain W in the initial value 1 1 : include W, 0 : not include W
IO_INIT_QC integer whether contain QC in the initial value 0
IO_INIT_QCI integer whether contain QCI in the initial value 0
IO_INIT_QR integer whether contain QR in the initial value 0
IO_INIT_QS integer whether contain QS in the initial value 0
IO_INIT_QG integer whether contain QG in the initial value 0
IO_INIT_QH integer whether contain QH in the initial value 0
IO_INIT_QNCW integer whether contain QNCW in the initial value 0 1 : include QNCW, 0 : not include QNCW
(2005/12/23)(2005/12/26)
IO_INIT_QNCI integer whether contain QNCI in the initial value 0
IO_INIT_QNR integer whether contain QNR in the initial value 0
IO_INIT_QNS integer whether contain QNS in the initial value 0
IO_INIT_QNG integer whether contain QNG in the initial value 0
IO_INIT_QNH integer whether contain QNH in the initial value 0
IO_INIT_ETB integer
  • 1st digit: "eturb" is contained in initial file?
  • 10th digit: "tsq" (autocorrelation of potential temperature perturbation) is contained in initial file?
  • 100th digit: "qsq" (autocorrelation of water volume perturbation) is contained in initial file?
  • 1000th digit: "cov" (correlation between potential temperature and water volume perturbation) is contained in initial file?
 0 1 : contain, 0 : not contain
IO_INIT_W_G integer "W_G" is contained in initial file? 0
IO_INIT_PRS integer whether contain PRS_LOC in the initial value 0
IO_BNDRY_W integer whether contain W in the boundary value 1
IO_BNDRY_QC integer whether contain QC in the boundary value 0 -1: Use Qx=0 at inflow. Not use Qx at outflow , replace the inner value.
0: not contain, 1: contain, 2: contain, use at inflow and not use at outflow
(2005/12/23)(2005/12/26)(2006/07/23)(2006/08/03)
IO_BNDRY_QCI integer whether contain QCI in the boundary value 0
IO_BNDRY_QR integer whether contain QR in the boundary value 0
IO_BNDRY_QS integer whether contain QS in the boundary value 0
IO_BNDRY_QG integer whether contain QG in the boundary value 0
IO_BNDRY_QH integer whether contain QH in the boundary value 0
IO_BNDRY_QNCW integer whether contain QNCW in the boundary value 0
IO_BNDRY_QNCI integer whether contain QNCI in the boundary value 0
IO_BNDRY_QNR integer whether contain QNR in the boundary value 0
IO_BNDRY_QNS integer whether contain QNS in the boundary value 0
IO_BNDRY_QNG integer whether contain QNG in the boundary value 0
IO_BNDRY_QNH integer whether contain QNH in the boundary value 0
IO_BNDRY_ETB integer whether contain ETURB in the boundary value 0 1 : contain, 0 : not contain
IO_BNDRY_PRS integer whether contain PRS_LOC in the boundary value 0
SW_BNDRH integer Use relative humidity for water vapor nesting 0 0: not use (qv is liner-interpolated in time)
1: use (RH is liner-interpolated in time, and qv is calculated from the.)
(06/09/19)
SW_CHANGE_SST integer treatment of SST and underground temperature 0 0: fixed tin(:,:,1) at initial time
1: SST is varied by time variation of tin(:,:,1) between boundary files.
(06/10/04)
2: add "1", tin(:,:,4) over LAND is varied by time variation of tin(:,:,4) between boundary files.
(07/07/03)

NAMRAD: parameters for radiation
name of parameter type of parameter meaning default memo
DTRADS real time interval of calculation of short radiation [sec] 300.0 available at "mswsys(13) ≥ 6". 900.0 is used in routine.
thinout_x_rd,
thinout_y_rd		 integer constant for tool of reduced gird 2 (2003/10/01)
ZCLMIN real lower limit height for diagnosing cloud amounts [m] 200.0 available at LOW_CLD = 1 (2004/02/16)
LOW_CLD integer decision method of lower limit of diagnosing lower cloud amounts (MSWSYS(13)>=10 and CLD_RD=0 or MSWSYS(13)=6, 7))
0 : diagnose above 5th layer
1 : use ZCLMIN as lower limit
2 : ignore ZCLMIN, diagnose above lowest layer
 0 (2004/02/16)(2007/02/14)
LOW_CLD_QC integer decision method of lower limit of diagnosing lower cloud amounts (MSWSYS(13)>=10 and CLD_RD=1, 2, 3 .or. MSWSYS(13)=8, 9)
0 : diagnose above 5th layer
1 : use ZCLMIN as lower limit
2 : ignore ZCLMIN, diagnose above lowest layer
 2 (2007/02/14)
MONIT_LOW_CLD integer the relationship between KCLMIN and cloud amount in NUSDAS output
0 : output cloud amount limited by KCLMIN
not 0 : ignore KCLMIN, output cloud amount diagnosed above lowest layer in NUSDAS		 0 (2007/02/14)
CLD_RD integer selection of cloud amount and cloud water/ice for radiation
  • 0 : Cloud amounts is diagnosed from RH in using the same function at mswsys(13) = 6.
    And cloud water is diagnosed by the diagnosed cloud amounts.
    Diagnostic methods of cloud water are switched by "DIAG_CWC", the ratio between ice/water is decided by air temperature.
  • 1 : Predicted cloud water / ice are used in radiation calculation.
    Cloud amounts are diagnosed by Xu and Randall(1996)
  • 2 : Predicted cloud water / ice are used in radiation calculation.
    Cloud fraction in a grid is set to 1.0 or 0.0 by QC, QCI, CLDRD2_QC_TH, CLDRD2_QCI_TH.
    This is effective in cloud resolving simulation.
  • 3 : Calculated cloud amount and cloud water/ice by partial condensation scheme are used in radiation calculation.(see details 2006/12/06)
 0 (valid for mswsys(13) >= 10)(2005/03/31)
CLD_PHYS integer selection of cloud raction and cloud water/ice for monitor
  • 1 : Predicted cloud water / ice are used for cloud fraction monitor.
    Cloud fraction are diagnosed by the method of Xu and Randall(1996)
  • 2 : Predicted cloud water / ice are used for cloud fraction monitor.
    Cloud fraction in a grid is set to only 0.0 or 1.0 by QC, QCI, CLDRD2_QC_TH, CLDRD2_QCI_TH(all or nothing).
 1 (valid when cloud micro-physics is used)
DIAG_CWC integer diagnostic method of cloud water when cloud amounts are diagnosed
0 : Heymsfield(1977)
1 : Hack(1998)		 0 (valid for mswsys(13) >= 10) (2005/06/09)
CWCL0 real maximum value of cloud water in using Hack(1998) [g/m^3] 0.21
AERSL integer Utilization of aerosol climatological value (valid for mswsys(13) >= 10)
0 : do not use
1 : use
2 : use 2D aerosol monthly climate(optical depth)		 1 (valid for mswsys(13) >= 10) (2005/03/31)
RAISA integer Utilization of the method of Räisänen(1998)(Correction of cloud fraction by emissivity)
0 : do not use the method of Räisänen(1998)
1 : use the method of Räisänen(1998)
0 (valid for mswsys(13) = 12)
SW_SLIQ integer effective radius of cloud water diagnosis method
0 : determined by SLIQ0
1 : 10μm over land, 13μm over sea(by Kawamoto et al. 2001)
0 (valid for mswsys(13) = 10, 12)
SW_SICE integer effective radius of cloud ice diagnosis method
0 : diagnosed from temperature by the method of Ou and Liou(1995)(tuned for previous GSM)
1 : diagnosed from temperature by the method of Ou and Liou(1995)(from original paper)
2 : diagnosed from temperature and cloud ice content by the method of Wyser(1998)
0
PFACT real P value for Xu and Randall (1996) 0.25 (valid for mswsys(13) >= 10) (2005/03/31)
GAM real γ value for Xu and Randall (1996) 0.20
ALPHAZERO real α0 value for Xu and Randall (1996) 0
OZONE integer Utilization of ozone climatological value 0 recommended to use when mswsys(13)=12, 13 (2006/04/28)
VL_EXP integer switch parameter to expand of the vertical layer to calculate the radiation from until nz (default, nz) 0
SLIQ0 real effective radius of water cloud(4.2d0<=SLIQ0<=16.6d0) 15(µm) Using for optical properties of water cloud (2006/07/27)
CUT_FLSH_S integer cut sensible heat flux at sea surface 0 1 : cut latent or sensible heat flux at sea surface if CUT_FL[SL]H = 1(2007/03/02)
CUT_FLLH_S integer cut latent heat flux at sea surface 0
CO2PPM real concentration CO2 (ppm) 375.0 used only for long wave radiation if mswsys(13) >= 12 (2008/02/22)
PPMCH4 real concentration of CH4 (ppm) 1.75 valid for mswsys(13) >= 12 (2008/02/22)
PPMN2O real concentration of N2O (ppm) 0.28
PPMF11 real concentration of CFC-11 (ppm) 0.3d-3
PPMF12 real concentration of CFC-12 (ppm) 0.5d-3
PPMF22 real concentration of HCFC-22 (ppm) 0.2d-3
CLDRD2_QC_TH real Threshold to determine cloud fraction (1.0 or 0.0) on QC 1.d-6(kg/kg) (valid for (mswsys(13) >= 10 .and. CLD_RD = 2) .or. MSWSYS(13) = 8, 9)
CLDRD2_QCI_TH real Threshold to determine cloud fraction (1.0 or 0.0) on QCI 1.d-6(kg/kg) (valid for (mswsys(13) >= 10 .and. CLD_RD = 2) .or. MSWSYS(13) = 8, 9)
TMEAN integer On radiation calculation, time averaged input(TMEAN>=1) or not(TMEAN == 0) 0 valid for MSWSYS(13) >= 6
MONIT_OVERLAP integer overlap assumption used in diagnosis of output high, middle, low level, total cloud fraction
0 : maximum overlap
1 : random overlap
2 : maximum-random overlap
 0 (2008/03/16)

NAMPTG: parameter about the calculation of surface temperature, ground and boundary layers scheme
name of parameter type of parameter meaning default value note
DAY0 real start day   number of days from Jan. 1st
recommend using timecard, not recommend using here
GTIME0 real start time   number of hours from Jan. 1st
recommend using timecard, not recommend using here
ALBEDL real albedo of land 0.2 Using at stand alone run or nesting run without albedo in ptgrd file
Precede ALBED_USR(1:4).
ALBEDS real albedo of sea 0.1 Using at stand alone run or nesting run without albedo in ptgrd file
Precede ALBED_USR(1:4).
WETL real wetness of land 0.1 Using at stand alone run or nesting run without wetness in ptgrd file
Precede WET_USR(1:4).
WETS real wetness of sea 1.0 Using at stand alone run or nesting run without wetness in ptgrd file
Precede WET_USR(1:4).
g_emit real rate of outgoing from surface 1.0 (2005/02/14)
FKTG_USR(1:4) real heat diffusion coefficient at surface [m^2/s^4]
(1): land, (2): sea,
(3): snow, (4): ice 		7.0d-7(1), 1.3d-6(2), 7.0d-7(3), 7.0d-7(4) Using at stand alone run or nesting run without fktg in ptgrd file
SST is not predictable, so (2)sea is not used.
ROCTG_USR(1:4) real heat capacity of surface [J/K/m^3]
(1): land, (2): sea,
(3): snow, (4): ice 		2.3d6(1), 1.9d6(2), 2.3d6(3), 2.3d6(4) Using at stand alone run or nesting run without roctg in ptgrd file
SST is not predictable, so (2)sea is not used.
WET_USR(1:4) real wetness
(1): land, (2): sea,
(3): snow, (4): ice 		0.1d0(1), 1.0d0(2), 0.1d0(3), 0.1d0(4) Using at stand alone run or nesting run without wetness in ptgrd file
In case of that WETL, WETS are defined, WET_USR(1, 3, 4) = WETL, WET_USR(2) = WETS
ALBED_USR(1:4) real albedo
(1): land, (2): sea,
(3): snow, (4): ice 		0.6d0(1), 0.1d0(2), 0.6d0(3), 0.6d0(4) Using at stand alone run or nesting run without albedo in ptgrd file
In case of that ALBEDL, ALBEDS are defined, ALBED_USR(1, 3, 4) = ALBEDL, ALBED_USR(2) = ALBEDS
ROUGH_RATIO_USR(1:4) real Constant for calculation of heat roughness
(1): land, (2): sea,
(3): snow, (4): ice 		1.d0 / 7.4d0(1), 1.d0(2), 1.d0 / 7.4d0(3), 1.d0 / 7.4d0(4) Using at stand alone run or nesting run without z_0h in ptgrd file, z_0h os calculated by using them and z_0m.
USE_URBAN integer 0 : not use urban heat
1 : use urban heat		 0 Using at nesting run within urban heat data
URB_SCHEME integer selection of the urban scheme 0 ones place digit
  = 0 : do not use urban scheme
  = 1 : single layer urban canopy scheme*
          (treat buildings as one body)
  = 2 : single layer urban canopy scheme*
          (distinguish roof and walls)
tens place digit
  = 0 : not consider the latent heat fluxes
        from/to the buildings
  = 1 : consider the latent heat fluxes
        from/to the buildings

*another setting file for the urban scheme is needed
URB_CHTC integer setting of the heat transfer on building surfaces 1   = 1 : Deardorff(1978) Ch=0.01*(1+0.3/uaf)
  = 2 : Kusaka et al.(2001) h=4.19*uaf+6.8
  = 3 : Masson(2000) h=4.2*uaf+11.8
  = 4 : Hagishima et al.(2003)
                h_roof=3.96*uaf+6.42; h_wall=10.21*uaf+4.47
URB_W_G_LIMIT logical set or do not set the limitation of wetness for the ground inside the canyon .true.   = .true. : set the limit
              upper=1.5*climatic value
              lower=0.7*climatic value
  = .false. : do not set the limit
URB_W_MAX_R real maximum water content of the reservoir at the roof 1.0d-3 [m3/m2]
URB_W_MAX_W real maximum water content of the reservoir at the wall 1.0d-4 [m3/m2]
URB_W_MAX_G real muximum value of soil moisture content at the 1st layer W_G 0.4  
RI_B_MAX real maximum value of bulk Richardson number 1.d10 Using at strong stable condition (2005/07)
SW_ROUGH_ZERO_SHIFT integer switch parameter for zero plane displacement d
0 : d=0
1 : d=-z_0m
    (momentum roughness)		 0 Using large roughness (above height of lowest layer or the limit of diagnosed roughness is exceed 1[m]), have to be set "SW_ROUGH_ZERO_SHIFT"=1 (2005/07)
MY3_SGM_MAX_FCT real coefficient to determine the maximum value of σ of partial condensation scheme 0.25d0  
MY3_SGM_MIN_FCT real coefficient to determine the minimum value of σ of partial condensation scheme 0.07d0  
MY3_BUOY integer 0:switch parameter of buoyancy flux ... no consideration of buoyancy flux for calculation of MY3
1:switch parameter of buoyancy flux ... consideration of buoyancy flux for calculation of MY3 		1  
MY3_PCOND integer 0:switch parameter of partial condensation scheme ... no consideration of partial condensation scheme for calculation of MY3
1:switch parameter of partial condensation scheme ... consideration of partial condensation scheme for calculation of MY3		 1  
MY3_ITR_INI integer MLevel 2 diagnosis iteration number in the beginning 10  
MY3_ADV integer switch parameter of consideration of turbulence variables advection
0:No consideration of advection
1:Use advection scheme of mswsys (do not use flux correction transport)
2:same as 1 but using flux correction transport
3:1st order upwind advection (do not use flux correction transport).
4:same as 3 but using flux correction transport		 1 (06/09/19)
MY3_PDF integer switch parameter of probability distribution function for partial condensation scheme
0:normal distribution
1:uniform distribution		 0  
MY3_PCOND_RAD_QCTH integer threshold of cloud water and cloud ice using cloud water content of partial condensation in the radiation scheme 1.d-6  
MY3_TRP_NOCLD integer switch parameter to suppress the production of cloud by partial condensation scheme over tropopause
0:no control
1:control		 1 (06/11/15)
MY3_IMPROVED_LENGTH integer Use improved mixing length by Nakanishi and Niino
0: not use (use original MY mixing length)
1: use		 1 (07/09/11)
MY3_IMPROVED_CLOSURE integer Use improved closure constants by Nakanishi and Niino
0: not use (use original MY closure constants)
1: use		 1
DDF_IMPL integer Switch whether implicit scheme is used for time integration of TKE when ten-thousand's place of MSWSYS(19) equals to 0.
0:not use
1:use 		0
W_G_MAX_FCT real factor for maxium value of soil moisture content at the 1st layer W_G against climatic values 1.5d0  
W_G_MIN_FCT real factor for minium value of soil moisture content at the 1st layer W_G against climatic values 0.7d0  
W_2_MAX_FCT real factor for maxium value of soil moisture content at the 2nd layer W_2 against climatic values 1.5d0  
W_2_MIN_FCT real factor for minium value of soil moisture content at the 2nd layer W_2 against climatic values 0.7d0  

NAMKDD: parameter of the factor for outputting file on the model plane ( KDD(*) < 0 : output both time level KT and KTN. )
name of parameter outputted factor meaning of factor scale of array value setting
KDD(1) U momentum in x-direction nx0*ny0*nz0 0 not output
1 output only time level = KTN
KDD(2) V momentum in y-direction nx0*ny0*nz0 0 not output
1 output only time level = KTN
KDD(3) W momentum in z-direction nx0*ny0*nz0 0 not output
1 output only time level = KTN
KDD(4) PT potential temperature nx0*ny0*nz0 0 not output
1 output only time level = KTN
2 output TIN(under ground temperature)
KDD(5) QV mixing ratio of water vapor nx0*ny0*nz0 0 not output
1 output only time level = KTN
2 output W_G(water contents in soil)
KDD(6) QC mixing ratio of cloud water nx0*ny0*nz0 0 not output
1 output QC (time level=KTN)
2 output QNCW
3 output both QC and QNCW
KDD(7) QR mixing ratio of rain nx0*ny0*nz0 0 not output
1 output QR (time level = KTN)
2 output QNR
3 output both QR and QNR
KDD(8) ETURB,
PTSQ, QWSQ, PTQW		 energy of turbulence and variables of turbulence nx0*ny0*nz0 0 not output
1 output ETURB at time level = KTN
2 output ETURB, PTSQ, QWSQ, PTQW at time level = KTN
KDD(9) EDDYxx coefficient of eddy diffusion nx0*ny0*nz0 0 not output
1 output EDDYKM (coefficient of eddy diffusion of momentum)
2 output EDDYKH (coefficient of eddy diffusion of heat)
3 output both EDDYKM and EDDYKH
KDD(10) PRS deviation of pressure nx0*ny0*nz0 0 not output
1 output at time level = KTN
KDD(11) PPT production rate of potential temperature nx0*ny0*nz0 0 not output
1 output
2 add "1" to PQV
3 add "2" to PQV and PQCW
4 add "3" to PQVA (adjustment amount in PQV)
KDD(12) PQR Production term of rain nx0*ny0*nz0 0 not output
1 output
2 add "1" to PQCI(for cloud ice)
3 add "2" to PQNR, PQNCW and PQNCI(for rain, cloudwater and cloudice)
KDD(13) QCI mixing ratio of cloud ice nx0*ny0*nz0 0 not output
1 output QCI
2 output QNCI
3 output both QCI and QNCI
KDD(14) QS mixing ratio of snow nx0*ny0*nz0 0 not output
1 out put QS at time level = KTN
2 output QNS
3 output both QS and QNS
KDD(15) QG mixing ratio of graupel nx0*ny0*nz0 0 not output
1 output QG at time level = KTN
2 output QNG
3 output both QG and QNG
KDD(16) PQS production term of snow / graupel nx0*ny0*nz0 0 not output
1 output at time level = KTN
2 add "1" to PQG
3 add "2" to PQNS and PQNG
KDD(17) DNSG2 Density of air * G^(1/2) nx0*ny0*nz0 -1 output at IT=ITST
0 not output
1 output
KDD(18) TPPT(TPQV) accumulated PPT(PQV) in time nx0*ny0*nz0 0 not output
1 output TPPT
2 output TPQV
3 output both TPPT and TPQV
KDD(19) BUOYD buoyancy term nx0*ny0*nz0 0 not output
1 output BUOYD
KDD(20) QH mixing ratio / production term of hail nx0*ny0*nz0 0 not output
1 output QH
2 output both QH and QNH
3 output QH, QNH and PQH
4 output QH, QNH, PQH and PQNH
KDD(21) QCNB,QCN,
NCNB,NCN		 mixing ratio of aerosolA and C nx0*ny0*nz0 0 not output
1 output QCNB,QCN
2 output NCNB,NCN
3 output QCNB,NCNB,QCN,NCN
KDD(22) SOLAR,DLWB,
TSOLAR,TDLWB,
CLL,CLM,CLH,
OLR,
and
QINB,QIN,
NINB,NIN		 for MRI format
short / long radiation
(its accumulation), cloud mounts
and
mixing ratio of aerosolB and D		 nx0*ny0 0 not output
***1 output SOLAR,DLWB
***2 output SOLAR,DSWB,DLWB
**1* TSOLAR,TDLWB
**2* output TSOLAR,TDSWB,TDLWB
*1** output CLL,CLM,CLH
*2** output OLR
*3** output CLL,CLM,CLH,OLR
nx0*ny0*nz0 1*** output QINB,QIN
2*** output NINB,NIN
3*** output QINB,QIN,NINB,NIN
KDD(23) RADPT(TRADPT) variation of potential temperature by radiation
(its accumulation in time)		 nx0*ny0*nz0 0 not output
1 output RADPT
2 output TRADPT
3 output RADPT, TRADPT
KDD(24) RSTPTW potential temperature/water vapor flux on surface by turbulence
(their accumulations in time only for MRI format)		 nx0*ny0*2 0 not output
*1 output RSTPTW
(level 2: output RSTUW)
*2 output RSTPTW,RSTQVW
(level 2 of RSTQVW: output RSTVW)
nx0*ny0 1* output TRSTPTW
2* output TRSTPTW,TRSTQVW
KDD(25) for MRI format
2-dimensional output		 PSEA,SM*,PTGRD,
TS,TMAX,TMIN,
U10,V10,UVMAX,
SM*CLD,SM*CV		 nx0*ny0 -1 not output
0 =11111 output PSEA,SM*,PTGRD,TS,U10,V10,
UVMAX,SM*CLD,SM*CV
****1 output PSEA
***1* output SMQR,SMQI,SMQS,SMQG,SMQH
***2* output SMQR,SMSW=SMQI+SMQS+SMQG+SMQH
***3* output SMAL=SMQR+SMQI+SMQS+SMQG+SMQH
**1** output PTGRD,TS
**2** output TS
**3** output TMAX,TMIN
**4** output TS,TMAX,TMIN
**5** output PTGRD,TS,TMAX,TMIN
*1*** output U10,V10
*2*** output UVMAX
*3*** output U10,V10,UVMAX
1**** output SM*CLD,SM*CV
KDD(26) SIB monitor SIB monitor for MRI-format
(2006/01/30)(2006/03/01)		 nx0*ny0 0 not output
1 minimized output (latent/sensible heat, friction and radiation at surface)
2 add "1" to gird average variables
3 add "2" to predicted variables in SIB
KDD(27) IST left end for MRI-format 1~nx0-1 0 not specified(IST=1)
IST left end in X-direction
KDD(28) IEN right end for MRI-format 2~nx0 0 not specified(IST=nx0)
IST right end in X-direction
KDD(29) JST south end for MRI-format 1 - ny0-1 0 not specified(JST=1)
IST south end in Y-direction
KDD(30) JEN north end for MRI-format 2~ny0 0 not specified(JEN=ny0)
IST north end in Y-direction
KDD(31) ISTR_START start time step for MRI-format itst~itend 0 not specified(ISTR_START=itst)
IT start time step for MRI-format
KDD(32) ISTR_END end time step for MRI-format itst~itend 0 not specified(ISTR_END=itend)
IT end time step for MRI-format

NAMOUT: parameter for surface element diagnosis
name of parameter type of parameter meaning default value note
Z_0M_LIMIT real upper limit of roughness for surface element diagnosis, negative value for no limit [m] -99.0 (2005/07)
Z_0H_FACT real When Z_0M(momentum roughness) exceed upper limit at the time of surface diagnosis, Z_0H(heat roughness) is given Z_0M_LIMIT*Z_0H_FACT
Used at Z_0M_LIMIT>0 only.		 0.1d0
SW_SFC_DIAG integer surface flux scheme for surface element diagnosis
1st and 10th digit : same as mswsys(1), but Louis for sea corrected by sub-layer Dalton number and sub-layer Stanton number
100th digit : not used
1000th digit : same as mswsys(1)		 2044
PTV_FACT_MAX
PTV_FACT_MIN		 real Upper and lower limit of weight of the lowest layer to diagnose virtual potential temperature of 1.5 m height 1.d0
0.d0
QV_FACT_MAX
QV_FACT_MIN		 real Upper and lower limit of weight of the lowest layer to diagnose relative humidity of 1.5 m height 1.d0
0.d0
SFC_WIND_HEIGHT real diagnosed height [m] of surface wind 10.d0 (2006/03/15)
SFC_T_HEIGHT real diagnosed height [m] of surface temperature (and TTD) 1.5d0

NAMSIB: parameter for MRI/JMA-SiB
name of parameter type of parameter meaning default value
IFLAG_USE_NHM_WET integer IFLAG_USE_NHM_WET=1 => convert WET in PTGRD file (evaporation efficiency "bete") into MJ-SiB saturation ratio
IFLAG_USE_NHM_WET/=1 => specify soil saturation ratios of three soil layers with 0.5 (default) or any values 		1
W1_SIB real specify saturation ration of 1st soil layer with W1_SIB, with IFLAG_USE_NHM_WET/=1 0.5
W2_SIB real specify saturation ration of 2nd soil layer with W2_SIB, with IFLAG_USE_NHM_WET/=1 0.5
W3_SIB real specify saturation ration of 3rd soil layer with W3_SIB, with IFLAG_USE_NHM_WET/=1 0.5
EXPL_TIME real stabilize time integration for canopy temperature, sensible heat and/or latent heat under dt=30-90 sec, with constant "arpha" in Krank-Nicolson method modified into expl_time/dt, which is used in coupling of SiB prognostic variables with T and Qv in the lowest atmosphere.
More effective, less EXPL_TIME is used. Do not work with dt<30 		10.0
IFLAG_SIB_TRB_IMPLICIT integer 1: use implicit solver of vertical turbulent diffusion terms in SiB.
0: use explicit solver of vertical turbulent diffusion terms in SiB.
0
SIB_EDDYK_Z_MAX real maximum of vertical turbulent diffusive coefficients in SiB 1.d10
SIB_KZ_TOP integer top of the atmosphere in SiB NZ-1

NAMSCM: parameters for SCM(Single Column Model)(It is included only if scm_mode = 1.)
name of parameter type of parameter meaning default value
SCM_ZS real Height above sea level in SCM (m) 0.0
SCM_SKIND integer Surface kind in SCM
0: not set(It is determined by SCM_SL)
1: land
2: sea
3: snow on land
4: ice on sea 		0
SCM_SL real Ratio of sea and land in SCM(0--1)
It is not used when SCM_SKIND is set. 		1.0
SCM_Z0M real Roughness for momentum(m) in SCM0.0
SCM_LAT real Latitude in SCM (degree)30.0
SCM_LON real Longitude in SCM (degree)140.0
SCM_PSURF real Initial surface pressure in SCM (Pa)101300.0
SCM_GEO_FORCE integer Switch for geostorophic wind forcing in SCM
0: OFF
1: forced with vertically and temporally constant values specified by SCM_UG and SCM_VG
2: forced with values varing vertically and temporally, specified by SCM_UG_INPUT and SCM_VG_INPUT
0
SCM_UG real forced geostrophic zonal wind (m/s)
used only when SCM_GEO_FORCE=1 		0.0
SCM_VG real forced geostrophic meridional wind (m/s)
used only when SCM_GEO_FORCE=1 		0.0
SCM_DT_UVG_INPUT real Time interval to give the time series of the geostrophic winds (s)
used only when SCM_GEO_FORCE=2 		0.0
SCM_UG_INPUT real(the number of times x the number of vertical layers)
max 105000 		Time series of the geostrophic zonal wind (m/s)
used only when SCM_GEO_FORCE=2 		 
SCM_VG_INPUT real(the number of times x the number of vertical layers)
max 105000 		Time series of the geostrophic meridional wind (m/s)
used only when SCM_GEO_FORCE=2 		 
SCM_UINC_FORCE integer Switch for increment forcing of zonal wind velocity in SCM
0:OFF, 1:ON 		0
SCM_VINC_FORCE integer Switch for increment forcing of meridional wind velocity in SCM
0:OFF, 1:ON 		0
SCM_TINC_FORCE integer Switch for increment forcing of temperature in SCM
0:OFF
1:ON, specified in temperature in SCM_TINC_INPUT
2:ON, specified in potential temperature in SCM_TINC_INPUT
0
SCM_QINC_FORCE integer Switch for increment forcing of specific humidity in SCM
0:OFF, 1:ON 		0
SCM_DT_INC_INPUT real Time interval to give time series of increment forcing (wind velocity, temperature and specific humidity) (s)
0.0
SCM_UINC_INPUT real(the number of times x the number of vertical layers)
max 105000 		Time series of vertical profiles forcing increments of zonal wind velocity (m/s/s)
 
SCM_VINC_INPUT real(the number of times x the number of vertical layers)
max 105000 		Time series of vertical profiles forcing increments of meridional wind velocity (m/s/s)
 
SCM_TINC_INPUT real(the number of times x the number of vertical layers)
max 105000 		Time series of vertical profiles forcing increments of temperature or potential temperature (K/s)
 
SCM_QINC_INPUT real(the number of times x the number of vertical layers)
max 105000 		Time series of vertical profiles forcing increments of specific humidity (kg/kg/s)
 
SCM_WFORC_FORCE integer Switch for vertical velocity forcing to vertically advect heat and moisture in SCM
0:OFF
1:ON, specified in m/s in SCM_WFORC_INPUT
2:ON, specified in Pa/s in SCM_WFORC_INPUT
0
SCM_DT_WFORC_INPUT real Time interval to give time series of forced vertical velocity (s)
0.0
SCM_WFORC_INPUT real(the number of times x the number of vertical layers)
max 105000 		Time series of vertical profiles forcing vertical velocity (m/s or Pa/s)
Note that these vertical velocities shoild be assigned at the full levels.
 
SCM_TIN_FORCE integer Switch for surface temperature forcing 0: OFF, 1: ON
0
SCM_DT_TIN_INPUT real Time width of time series of surface temperature for forcing (sec) 0
SCM_TIN_INPUT real(max 500) Time Series of forced surface temperature (K)
An initial value should be stored in the first item of the array and its time series are given with the time width of SCM_DT_TIN_INPUT. 		 
SCM_FLUX_H_FORCE integer Switch for surface sensible heat flux forcing in SCM
0:OFF, 1:ON 		0
SCM_FLUX_E_FORCE integer Switch for surface latent heat flux forcing in SCM
0:OFF, 1:ON 		0
SCM_FLUX_M_FORCE integer Switch for surface latent momentum flux forcing in SCM
0:OFF
1:use the forcing specified by time series -1: friction velocity is calculated assuming the Monin-Obukov similarity and momentum fluxes are derived with the friction velocity (It is valid only when forcing of sensible and latent heat fluxes are turned on.) 		0
SCM_DT_FLUX_INPUT real Time interval to give time series of forced surface fluxes (sensible heat, latent heat, momentum)(s)
0.0
SCM_FLUX_H_INPUT real(max 1500) Time series of forced surface sensible heat flux (W/m2)
 
SCM_FLUX_E_INPUT real(max 1500) Time series of forced surface latent heat flux (W/m2)
 
SCM_FLUX_U_INPUT real(max 1500) Time series of forced surface momentum flux of zonal element (kg/ms2)
should be specified by -rho*(u'w') 		 
SCM_FLUX_V_INPUT real(max 1500) Time series of forced surface momentum flux of meridional element (kg/ms2)
should be specified by -rho*(v'w') 		 

NAMRCM: parameters for RCM (Regional Climate Model)
name of parametertype of parametermeaningdefault value
use_sbcinteger1: Use SBC method(Kida et al. 1991)
0: Not use SBC method		0
sbc_x_wavenumrealHow small scale information in inner model to be replaced.
>0 :max. mode numbers in X-direction
<0 : min. wavelength (m)
=0 error (to be aborted)		-1000.0d3
sbc_y_wavenumrealHow small scale information in inner model to be replaced.
>0 :max. mode numbers in Y-direction
<0 : min. wavelength (m)
=0 error (to be aborted)		-1000.0d3
sbc_hminrealSBC only works above sbc_hmin(m). The SBC intensity increase linerly from sbc_hmin to sbc_hmax.
should be >0		7000.0d0
sbc_hmaxreallowest altitude (m) for SBC working at SBC_MAXRATIO
should be ≥sbc_hmin		7000.0d0
sbc_intervalreal>0: SBC works evety sbc_interval (sec)
<0: SBC works every abs(sbc_interval) steps.		1200.0d0
sbc_maxratiorealshould be 0≤ sbc_maxratio ≤11.0d0