
###############################
##### Grid Specifications #####
###############################

im: 291                     # I gridpoints
jm: 361                     # J gridpoints
lm:  35                     # Number of atmospheric layers

tph0d: 0.                   # Central geographic latitude of grid (degrees); NOT NEEDED FOR GRID-ASSOCIATED NESTS.
tlm0d: 0.                   # Central geographic longitude of grid (degrees, positive east); NOT NEEDED FOR GRID-ASSOCIATED NESTS.

wbd: 0.                     # Grid's western boundary (rotated degrees); NOT NEEDED FOR GRID-ASSOCIATED NESTS.
sbd: 0.                     # Grid's southern boundary (rotated degrees); NOT NEEDED FOR GRID-ASSOCIATED NESTS.


##############################
##### Grid Decomposition #####
##############################

inpes:  06                  # Number of compute tasks in the I direction
jnpes:  07                  # Number of compute tasks in the J direction


###############################
#### Specify the I/O tasks ####
###############################

quilting:              true    # Do you want asynchronous quilting/history writes?
write_groups:          1
write_tasks_per_group: 1


##########################################
##### Fundamental Timestep (seconds) #####
##########################################

dt_int: 10                  # Integer seconds
dt_num:  0                  # Numerator of fractional second
dt_den:  1                  # Denominator of fractional second

filt_dt_int: 10                  # Integer seconds
filt_dt_num:  0                  # Numerator of fractional second
filt_dt_den:  1                  # Denominator of fractional second

###############################
##### Forecast Start Time #####
###############################

start_year:   2013
start_month:  01
start_day:    20
start_hour:   00
start_minute: 0
start_second: 0


#############################################################
#####  Forecast Length / Restarting / Output Frequency  #####
#############################################################

nhours_fcst:      3        # Length of the forecast (hours)
minutes_history: 60        # Frequency of history output (hours)
nhours_dfini:      0

restart:         false    # True--> A restarted run
minutes_restart: 720        # Frequency of restart output (hours)

avg_max_length: 3600        # Time in seconds when averaged fields are reset

rst_out_00:          false  # True--> Write 00h history in restarted run
write_last_restart:  true   # True--> Write last restart file
freerun:             true   # True--> Original
                            # False--> 5X divergence dumping in NDAS forec.

#########################
##### General modes #####
#########################

global:  false              # True--> Global ; False--> Regional
hydro:   false              # True--> Hydrostatic ;  False--> Nonhydrostatic
adiabatic: false            # True--> Adiabatic ;  False--> Diabatic
oper:      false


###########################
#####  Printouts      #####
###########################

print_all:     false        # Print all statements to err file
print_diag:    false        # Print diagnostics (CALL FIELD_STATS)
print_esmf:    true         # Print ESMF return signals
print_output:  false        # Print info on records written into history/restart file
print_timing:  true         # Print timing info (CPU time)


####################
##### Dynamics #####
####################

secadv:  true               # True--> 2nd order advection ; False--> 4th order advection

smag2:   0.4                # Smagorinsky constant for 2nd order diffusion

codamp:  9.0                # Divergence damping constant
wcor:    0.18               # Divergence correction factor

idtadt:  2                  # Number of adjustment timesteps between tracer advection calls (normally 1)

num_tracers_chem:  0        # Number of specified chem/aerosol "tracer" scalars

lnsh:  1                    # Number of boundary blending rows for H points
lnsv:  1                    # Number of boundary blending rows for V points


###################################
##### Digital Filtering Options ###
###################################

filter_method: 3            # 0--->none ; 1--->DFL ; 2---> DDFI ; 3--->TDFI

nsecs_dfl: 1800             # HALF foward filter window (s)          (DFL)
nsecs_bckddfi: 1800         # HALF backward integration duration (s) (DDFI)
nsecs_fwdddfi: 1800         # HALF foward filter window (s)          (DDFI)
nsecs_bcktdfi: 1800         # HALF backward filter window (s)        (TDFI)
nsecs_fwdtdfi: 1800         # HALF forward filter window (s)         (TDFI)

radar_init: 0

#############################
#####  Global Summations ####
#############################

use_allreduce:  false       # False--> use mpi send/recv for global sums
                            # True--> use mpi_allreduce for global sums

##################################
### Precipitation Assimilation####
##################################

pcpflg:  false              # True--> assimilation on ; False--> assimilation off
pcphr: 6
write_prec_adj: false       # True--> create baseline prec. files ; False--> regular run


###########################
### Gravity Wave Drag  ####
###########################

gwdflg:  true               # True--> GWD on ; False--> GWD off

cdmb:       0.1             # tunable parameter for adjusting mountain blocking
cleffamp:   1.00            # tunable parameter, depends on resolution
sigfac:     0.00            # tunable parameter
factop:     0.50            # de-acceleration limiter
rlolev:   500.00            # critical pressure level (check units)
dpmin:      5.00            # minimum thickness of the reference layer (check units)


##################################
##### Physics Scheme Options #####
##################################

#---------------------------------
shortwave: rrtm             # Radiation schemes
longwave:  rrtm
                                # gfdl--> GFDL (shortwave and longwave)
                                # gsfc--> Goddard shortwave
                                # dudh--> Dudhia shortwave
                                # rrtm--> RRTM longwave

np3d: 3                     # cloud properties for radiation
                                # 3--> ferrier's microphysics cloud scheme (only stratiform cloud)
                                # 4--> zhao/carr/sundqvist microphysics cloud (now available in the NMMB)
                                # 5--> NAM stratiform + convective cloud optical depth and fraction (GFDL type convective cloud)
                                # 8--> Thompson microphysics and special Thompson rrtm SW radiation

cldfraction: default        # cloud fraction scheme for radiation
                                # default  --> according to the choice of np3d
                                # thompson --> thompson

co2tf:  1                       # co2tf--> Read (0) or generate internally (1)
                                #          the GFDL CO2 transmission functions

#---------------------------------
convection:  bmj            # Convection scheme
                                # bmj --> Betts-Miller-Janjic
                                # bmj_dev--> BMJ (development)
                                # kf  --> Kain-Fritsch
                                # sas --> Simplified Arakawa-Schubert
                                # gd  --> Grell-Devenyi
                                # none--> No convection

#---------------------------------
microphysics:  fer_hires    # Microphysics scheme
                                # fer --> Ferrier
                                # fer_hires--> Ferrier operational
                                # kes --> Kessler
                                # lin --> Lin
                                # tho --> Thompson
                                # wsm3--> WSM3

rhgrd: 0.98                 # Onset of condensation
spec_adv:  .false.          # Individual cloud species advected
lmprate:  .false.           # .true.: 4D array called MPRATES containing
                            # 3D microphysics source/sink terms
                            # .false.: One 3D array called MPRATES
                            # with zero values

#---------------------------------
turbulence:  myj            # Turbulence schemes
                                # myj--> Mellor-Yamada-Janjic

#---------------------------------

sfc_layer:  myj             # Surface layer schemes
                                # myj--> Mellor-Yamada-Jamnjic
                                # mm5--> MM5 sfc layer

#---------------------------------

land_surface:  noah         # Land surface schemes
                                # noah--> NOAH scheme

ucmcall: 0                  # Use (1) or do not use (0) the Urban Canopy Model
ivegsrc: 1                  # Vegetation map identifier, 0--> USGS, 1--> IGBP

#---------------------------------

gfs:  false                 # Select entire GFS physics suite


######################################
##### Shallow Convection Switches ####
######################################

fres:     1.00               # resolution factor for dsp's (default)
fr:       1.00               # land factor for dsp's (dafault)
fsl:      0.85               # reduction factor for "slow" dsp's over land (dafault)
fss:      0.85               # reduction factor for "slow" dsp's over water (dafault)

#fres:     0.80               # resolution factor for dsp's
#fr:       1.00               # land factor for dsp's
#fsl:      0.85               # reduction factor for "slow" dsp's over land
#fss:      0.85               # reduction factor for "slow" dsp's over water

entrain:    .false.         # entrainment
newall:     .false.         # new cloud used at all shallow points
newswap:    .false.         # new clouds at swap shallow points
newupup:    .false.         # new cloud used for both heat and moisture up shallow pts.
nodeep:     .false.         # all deep convection diverted to shallow swap algorythm


#############################
##### Physics Timesteps #####
#############################

nrads:    360               # Number of dynamics timesteps between calls to shortwave
nradl:    360               # Number of dynamics timesteps  between calls to longwave

nphs:      18               # Number of dynamics timesteps between calls to landsurface and turbulence

nprecip:   18               # Number of dynamics timesteps between calls to convection and microphysics

####
#### NOTE:  User must set nhrs_udef to .TRUE. and set the emptying frequencies (nhrs_*) to the desired values
####        or else all accumulations will automatically be emptied hourly.
####

nhrs_udef:  true            # User defined when fluxes are emptied.
                            # False--> Auto ; True--> User defined:

nhrs_prec:  3               # Frequency in hours between times when precip arrays are emptied
nhrs_heat:  3               # Frequency in hours between times when heating arrays are emptied
nhrs_clod:  3               # Frequency in hours between times when cloud arrays are emptied
nhrs_rdlw:  3               # Frequency in hours between times when LW radiation arrays are emptied
nhrs_rdsw:  3               # Frequency in hours between times when SW radiation arrays are emptied
nhrs_srfc:  3               # Frequency in hours between times when sfc evap/flux arrays are emptied


#################################################
####  Write History/Restart Specifications  ####
#################################################

nemsio_input: false

write_hst_bin:           .TRUE.      # True--> Write history files in binary format
write_hst_nemsio:        .TRUE.      # True--> Write history files in NEMSIO format
write_rst_bin:           .TRUE.      # True--> Write restart files in binary format
write_rst_nemsio:        .TRUE.      # True--> Write restart files in NEMSIO format

write_nemsioctl:         .FALSE.     # True--> Yes write ctl for nemsio run history files
write_donefileflag:      .FALSE.
write_fsyncflag:         .FALSE.     # True--> Yes write fsync
write_dopost:        false  # True--> Yes run post on quilt
post_gribversion:   'grib1' # grib version for post output

hst_name_base:           'nmmb_hst'
rst_name_base:           'nmmb_rst'


##################################
####  Nesting Specifications  ####
##################################

my_domain_id:       3
my_parent_id:       2
n_children:         0
i_parent_start:   247
j_parent_start:    82
parent_child_space_ratio:  3

my_domain_moves:  false
ntrack: 00                    # the number of nphys steps between tracker calls

input_ready:  .FALSE.

nest_mode:  '1-way'          # If there is nesting, select '1-way' or '2-way' interaction between parents/children

#######################
####  Timing Info  ####
#######################

nhours_clocktime:   1        # Frequency in hours between clocktime diagnostic prints.
npe_print:         55        # The MPI task that will provide the clocktimes.


#######################
####  Hurricane    ####
#######################

run_tc:             false    # true=use hwrf phys

sas_pgcon:          0.55     # convectively forced pressure gradient factor,default=0.55
sas_shal_pgcon:     -1       # convectively forced pressure gradient factor sas shallow conv,
                             # '-1' means sas_pgcon is used
sas_shalconv:       1        # 1=enable shallow conv,better with gfspblhur scheme
sas_mass_flux:      9.0e9    # mass flux limit,default=9e9
sas_mommix:         1.0      # SAS momentum mixing coef
var_ric:            1.0      # 1=use variable critical Ri# in gfspblhur
coef_ric_l:         0.16     # Regression coef for land Ric,default=0.16
coef_ric_s:         0.16     # Regression coef for sea Ric,default=0.16
alpha:              0.7      # adjustment coef for K in gfspblhur
sfenth:             0.0      # GFDL surface-layer enhancement coef
disheat:            true     # true='consider diss heating'

