# control files for creofs, which is read in by shell script 

export DBASE_MET_NOW=NAM4	
export DBASE_MET_FOR=NAM4
export DBASE_WL_NOW=NCOM
export DBASE_WL_FOR=NCOM
export DBASE_TS_NOW=NCOM
export DBASE_TS_FOR=NCOM

export OCEAN_MODEL=SELFE
export LEN_FORECAST=48
export IGRD_MET=0
export IGRD_OBC=1
export BASE_DATE=2011010100
export TIME_START=2011090100
export MINLON=-128.5
export MINLAT=38.50
export MAXLON=-121.0 
export MAXLAT=50.25
export SCALE_HFLUX=1.0 
export CREATE_TIDEFORCING=1
########################################################
##  static input file name, do not include path name
########################################################
export GRIDFILE=nos.creofs.hgrid.gr3
export GRIDFILE_LL=nos.creofs.hgrid.ll
export HC_FILE_OBC=nos.creofs.HC.nc 
export HC_FILE_OFS=nos.creofs.HC.nc 
export RIVER_CTL_FILE=nos.creofs.river.ctl
export RIVER_CLIM_FILE=nos.ofs.river.clim.usgs.nc
export OBC_CTL_FILE=nos.creofs.obc.ctl
export OBC_CLIM_FILE=nos.ofs.clim.WOA05.nc
export STA_OUT_CTL=nos.creofs.station.in
export STA_NETCDF_CTL=nos.creofs.station.info
export VGRID_CTL=nos.creofs.vgrid.in
export RUNTIME_CTL=nos.creofs.run_control.param.in
export Nudging_weight=nos.creofs.t_nudge.gr3
export RUNTIME_MET_CTL=nos.creofs.sflux.inputs.ctl
export RUNTIME_COMBINE_RST=nos.creofs.nowcast.combine.hotstart.in
export RUNTIME_COMBINE_NETCDF=nos.creofs.combine.netcdf.field.in
export RUNTIME_COMBINE_NETCDF_STA=nos.creofs.combine.netcdf.station.in
export CORRECTION_STATION_CTL=nos.creofs.wl.correction.ctl
export WL_OFFSET_OLD=nos.creofs.wl.correction.old
########################################################
# parameters for SELFE RUN
########################################################
export ne_global=142684
export np_global=74061
export ns_global=216748
export nvrt=54
export NNODE=74061
export NELE=142684
export KBm=54
# Notes to change TIME STEP
# 1. DELT_MODEL is equal to EXTSTEP_SECONDS
# 2. 3600/DELT_MODEL has to be an integer (for station and field output time intervals)
# 3. Change DELT=DELT_MODEL/3600 in river control file
# 4. Change DELT=DELT_MODEL in OBC control file (in seconds)
export DELT_MODEL=120      
export EXTSTEP_SECONDS=120.0
export NDTFAST=20
export NRST=21600
export NSTA=360
export NSTATION=97
export NFLT=3600
export NHIS=3600
export NAVG=3600
export CPP_LON_VALUE=-124
export CPP_LAT_VALUE=45
export IHORCON_VALUE=0
export IHDIF_VALUE=0
export IDRAG_VALUE=2
export BFRIC_VALUE=0
export IHHAT_VALUE=1 
export INUNFL_VALUE=0
# STEP_NU_VALUE is time interval (in seconds) of T & S nudging 
export STEP_NU_VALUE=10800.0   
export MIN_DEPTH=0.01
export NWS_VALUE=2
export NRAMPWIND_VALUE=1
export DRAMPWIND_VALUE=1.0
export IWINDOFF_VALUE=0
export RDRG=3.0d-03
export RDRG2=0.005d0
export Zob=0.0005d0
export VISC2=0.0d0
export VISC4=0.0d0
export AKT_BAK="5.0d-6 5.0d-6 !m2/s"                   
export AKV_BAK="5.0d-5   !m2/s"                         
export AKK_BAK="5.0d-6   !m2/s"                        
export AKP_BAK="5.0d-6   !m2/s "                       
export DCRIT="0.10d0     !m"                 
export DSTART=151.0d0
export TIDE_START=0.0d0
export N_PROC=6
export TOTAL_TASKS=384
export NTILE_I=12
export NTILE_J=16
export NVAR=9
# #############################################################
# GLOSSARY
# #############################################################
# GRIDFILE    :ocean model grid netCDF file including lon, lat, depth, etc.
# DBASE       :Name of NCEP atmospheric operational products, e.g. NAM, GFS, RTMA, NDFD, etc.
# DBASE_MET_NOW : Data source Name of NCEP atmospheric operational products for Nowcast run.
# DBASE_MET_FOR : Data source Name of NCEP atmospheric operational products for Forecast run.
# DBASE_WL_NOW  : Data source Name of water level open boundary conditions for Nowcast run.
# DBASE_WL_FOR  : Data source Name of water level open boundary conditions for Forecast run.
# DBASE_TS_NOW  : Data source Name of T & S open boundary conditions for Nowcast run.
# DBASE_TS_FOR  : Data source Name of T & S open boundary conditions for Forecast run.
# OCEAN_MODEL :Name of Hydrodynamic Ocean Model, e.g. ROMS, FVCOM, SELFE, etc.
# LEN_FORECAST:Forecast length of OFS forecast cycle.
# IGRD_MET    :spatial interpolation method for atmospheric forcing fields
#           =0:on native grid of NCEP products with wind rotated to earth coordinates
#	    =1:on ocean model grid (rotated to local coordinates) interpolated using remesh routine.
#	    =2:on ocean model grid (rotated to local coordinates) interpolated using bicubic routine.
#	    =3:on ocean model grid (rotated to local coordinates) interpolated using bilinear routine.
#           =4:on ocean model grid (rotated to local coordinates) interpolated using nature neighbors routine.
# IGRD_OBC    :spatial interpolation method for ocean open boundary forcing fields
# BASE_DATE   :base date for the OFS time system, e.g. YYYYMMDDHH (2008010100)
# TIME_START  :forecast start time/current time, e.g. 2008110600
# MINLON      :longitude of lower left/southwest corner to cover the OFS domain
# MINLAT      :latitude of lower left /southwest corner to cover the OFS domain
# MAXLON      :longitude of upper right/northeast corner to cover the OFS domain
# MAXLAT      :latitude of  upper right/northeast corner to cover the OFS domain
# THETA_S     :S-coordinate surface control parameter, [0 < theta_s < 20].
# THETA_B     :S-coordinate bottom  control parameter, [0 < theta_b < 1].
# TCLINE      :Width (m) of surface or bottom boundary layer in which
#             :higher vertical resolution is required during stretching.
# SCALE_HFLUX :scaling factor (fraction) of surface heat flux (net short-wave and downward
#              long-wave radiation). if =1.0, no adjustment to atmospheric products.  
# CREATE_TIDEFORCING : > 0, generate tidal forcing file
# HC_FILE_ADCIRC     : ADCIRC EC2001 harmonic constant file 
# HC_FILE_ROMS     : Tidal forcing file of ROMS (contains tide constituents of WL, ubar, and vbar) 
# EL_HC_CORRECTION   : > 0, correction elevation harmonics with user provided data
# FILE_EL_HC_CORRECTION : file name contains elevation harmonics for correction                
# RIVER_CTL_FILE  : File name contains river attributes (Xpos, Epos, Flag, River name,etc.)
# OBC_CTL_FILE  : Control file name for generating open boundary conditions (WL, T and S).
# IM          :GRID Number of I-direction RHO-points, it is xi_rho for ROMS
# JM          :GRID Number of J-direction RHO-points, it is eta_rho for ROMS
# DELT_ROMS   :Time-Step size in seconds.  If 3D configuration, DT is the
#              size of baroclinic time-step.  If only 2D configuration, DT
#              is the size of the barotropic time-step.
#  NDTFAST     Number of barotropic time-steps between each baroclinic time
#              step. If only 2D configuration, NDTFAST should be unity since
#              there is not need to splitting time-stepping.
# KBm         :Number of vertical levels at temperature points of OFS
#  NRST        Number of time-steps between writing of re-start fields.
#
#  NSTA        Number of time-steps between writing data into stations file.
#              Station data is written at all levels.
#
#  NFLT        Number of time-steps between writing data into floats file.
#  NHIS        Number of time-steps between writing fields into history file.
#
#  RDRG2       Quadratic bottom drag coefficient.
#
#  Zob         Bottom roughness (m).
#  AKT_BAK     Background vertical mixing coefficient (m2/s) for active
#              (NAT) and inert (NPT) tracer variables.
#  AKV_BAK     Background vertical mixing coefficient (m2/s) for momentum.
#
#  AKK_BAK     Background vertical mixing coefficient (m2/s) for turbulent
#              kinetic energy.
#
#  AKP_BAK     Background vertical mixing coefficient (m2/s) for turbulent
#              generic statistical field, "psi".
#
#  TKENU2      Lateral, harmonic, constant, mixing coefficient (m2/s) for
#              turbulent closure variables.
#
#  TKENU4      Lateral, biharmonic, constant mixing coefficient (m4/s) for
#              turbulent closure variables.
#  DCRIT       Minimum depth (m) for wetting and drying.
#  DSTART      Time stamp assigned to model initialization (days).  Usually
#              a Calendar linear coordinate, like modified Julian Day.  For
#              Example:
#  TIDE_START  Reference time origin for tidal forcing (days). This is the
#              time used when processing input tidal model data. It is needed
#              in routine "set_tides" to compute the correct phase lag with
#              respect ROMS/TOMS initialization time.
# N_PROC       Number of computer processors to run parallel ocean model
# TOTAL_TASKS  Total tasks to be run
#  NTILE_I     Number of domain partitions in the I-direction (XI-coordinate).
#              It must be equal or greater than one.
#
#  NTILE_J     Number of domain partitions in the J-direction (ETA-coordinate).
#              It must be equal or greater than one.