# CONFIGURATION TEMPLATE (EMBM + GOLDSTEIN + GOLDSTEIN SEA-ICE + ATCHEM + BIOGEM)


# *******************************************************************
# COMPILER OPTIONS
# *******************************************************************
# Set floating point precision
# NOTE: currently BIGOEM will ony work in conjunction the IGCM compiled if everything is real*8
GENIEDP=TRUE
IGCMATMOSDP=TRUE
# Set build options (un-comment to compile and run with more comprehsneive de-bugging)
#BUILD=DEBUG


# *******************************************************************
# GENIE COMPONENT SELECTION
# *******************************************************************
ma_flag_ebatmos=.TRUE.
ma_flag_goldsteinocean=.TRUE.
ma_flag_goldsteinseaice=.TRUE.
ma_flag_igcmatmos=.FALSE.
ma_flag_fixedocean=.FALSE.
ma_flag_fixedseaice=.FALSE.
ma_flag_fixedicesheet=.FALSE.
ma_flag_fixedchem=.FALSE.
ma_flag_biogem=.TRUE.
ma_flag_atchem=.TRUE.


# *******************************************************************
# GRID CONFIGURATION
# *******************************************************************
# Set ATCHEM grid the same as GOLDSTEIN grid (but one vertical level)
GENIENXOPTS='-DGENIENX=36'
GENIENYOPTS='-DGENIENY=36'
GENIENLOPTS='-DGENIENL=1'
# Grid resolution
GENIENXOPTS='$(DEFINE)GENIENX=36'
GENIENYOPTS='$(DEFINE)GENIENY=36'
GOLDSTEINNLONSOPTS='$(DEFINE)GOLDSTEINNLONS=36'
GOLDSTEINNLATSOPTS='$(DEFINE)GOLDSTEINNLATS=36'
GOLDSTEINNLEVSOPTS='$(DEFINE)GOLDSTEINNLEVS=8'
# Specify the observational data and topography: EMBM
ea_topo='p0650a'
ea_taux_u="p0650a_taux_u.interp"
ea_tauy_u="p0650a_tauy_u.interp"
ea_taux_v="p0650a_taux_v.interp"
ea_tauy_v="p0650a_tauy_v.interp"
ea_adv_u="p0650a_uncep.silo"
ea_adv_v="p0650a_vncep.silo"
# Specify the observational data and topography: GOLDSTEIN
go_topo='p0650a'
# Specify the observational data and topography: GOLDSTEIN sea-ice
gs_topo='p0650a'
#assumed longitudinal offset of the grid
gm_par_grid_lon_offset=-180.0


# *******************************************************************
# TRACER CONFIGURATION
# *******************************************************************
#
# Set number of tracers
GOLDSTEINNTRACSOPTS='$(DEFINE)GOLDSTEINNTRACS=14'
# Set selected tracers
gm_atm_select_3=.true.
gm_atm_select_4=.true.
gm_atm_select_6=.true.
gm_atm_select_16=.true.
gm_ocn_select_3=.true.
gm_ocn_select_4=.true.
gm_ocn_select_8=.true.
gm_ocn_select_10=.true.
gm_ocn_select_12=.true.
gm_ocn_select_15=.true.
gm_ocn_select_16=.true.
gm_ocn_select_20=.true.
gm_ocn_select_35=.true.
gm_ocn_select_38=.true.
gm_ocn_select_40=.true.
gm_ocn_select_50=.true.
gm_sed_select_3=.true.
gm_sed_select_4=.true.
gm_sed_select_8=.true.
gm_sed_select_14=.true.
gm_sed_select_15=.true.
gm_sed_select_22=.true.
gm_sed_select_32=.true.
gm_sed_select_33=.true.
gm_sed_select_34=.true.
gm_sed_select_36=.true.
# Set tracer initial values
ac_atm_init_3=278.0E-06
ac_atm_init_4=-6.5
ac_atm_init_6=0.2095
bg_ocn_init_3=2.244E-03
bg_ocn_init_4=0.4
bg_ocn_init_8=2.159E-06
bg_ocn_init_10=1.696E-04
bg_ocn_init_12=2.363E-03
bg_ocn_init_35=1.025E-02
bg_ocn_init_38=2.916E-02
bg_ocn_init_40=0.0
bg_ocn_init_50=5.282E-02


# *******************************************************************
# OUTPUT CONFIGURATION
# *******************************************************************
# Set experiment ID
EXPID=cgenie_eb_go_gs_ac_bg.p0650a.NONE
# this is to only write ocean-grid data in genie-main
ma_write_flag_atm=.false.
ma_write_flag_sic=.false.
# this is to turn graphics off
ma_lgraphics=.false.


# *******************************************************************
# TIME CONTROL AND TIME-STEPPING
# *******************************************************************
# NOTE: For igcm where dt is usually 1 hour, use koverall_total=8640 for a 1 year run.
#       For c-goldstein, where dt is usually 365.25/500 days, use koverall_total=500 for a 1 year run.
# NOTE: 720=1 month of igcm (timestep=1 hour) to 2500000=5000 years of c-goldstein
#       BIOGEM run length = ma_koverall_total/500
#
# 1. TIME_STEP DURATION
# c-goldstein; ma_genie_timestep = 365.25*24.0/(5*96) * 3600.0 (GOLDSTEIN year length)
ma_genie_timestep=65745.0
#
# 2. RELATIVE EMBM/GOLDSTEIN/SEA-ICE CALLING FREQUENCY
ma_ksic_loop=5
ma_kocn_loop=5
# 3. set reduced ocean:biogeochem time-stepping gearing
ma_conv_kocn_katchem=2
ma_conv_kocn_kbiogem=2


# *******************************************************************
# PHYSICAL CLIMATE CALIBRATION
# *******************************************************************
# this is to set the climate model goin parameters as per Ridgwell et al. [2007]
# temp0
go_10=0.0
# temp1
go_11=0.0
# rel
go_12=0.9000000
# scl_tau
go_13=1.931692
# ocean diffusivites iso (or horiz)
go_14=4488.812
# ocean diffusivites dia (or vert)
go_15=2.7242320E-05
# inverse minimum drag in days
go_16=2.939598
# scl_tau (should be same as GOLDSTEIN's value)
ea_11=1.931692
# atm. diff. amp. for T
ea_12=4667550.0
# atm. diff. amp. for q
ea_13=1097107.0
# dist'n width
ea_14=1.083292
# dist'n slope
ea_15=6.3252553E-02
# atm. advection factor for T_z
ea_16=0.1121861
# atm. advection factor for T_m
ea_17=0.0000000E+00
# atm. advection factor for q_z
ea_18=0.2269510
# atm. advection factor for q_m
ea_19=0.2269510
# extra1a
ea_25=0.0
# extra1b
ea_26=0.0
# extra1c
ea_27=0.0
# sea-ice eddy diffusivity
gs_11=6200.000


# *******************************************************************
# CLIMATE ADJUSTMENT
# *******************************************************************
# set no seasonal cycle
ea_dosc=.false.
go_dosc=.false.
gs_dosc=.false.


# *******************************************************************
# USER-APPENDED OPTIONS FOLLOW ...
# *******************************************************************