Skip to content
Snippets Groups Projects
Commit 4c42515b authored by Tomas Lovato's avatar Tomas Lovato
Browse files

namelist update for V4.2.1

parent 0a4fa8f2
No related branches found
No related tags found
No related merge requests found
!------------------------------------------------------------------------------
&namdyn_adv ! Ice advection
!------------------------------------------------------------------------------
ln_adv_Pra = .true. ! Advection scheme (Prather)
ln_adv_UMx = .false. ! Advection scheme (Ultimate-Macho)
nn_UMx = 5 ! order of the scheme for UMx (1-5 ; 20=centered 2nd order)
/
!------------------------------------------------------------------------------
&namsbc ! Ice surface boundary conditions
!------------------------------------------------------------------------------
rn_cio = 5.0e-03 ! ice-ocean drag coefficient (-)
nn_snwfra = 2 ! calculate the fraction of ice covered by snow (for zdf and albedo)
! = 0 fraction = 1 (if snow) or 0 (if no snow)
! = 1 fraction = 1-exp(-0.2*rhos*hsnw) [MetO formulation]
! = 2 fraction = hsnw / (hsnw+0.02) [CICE formulation]
rn_snwblow = 0.66 ! mesure of snow blowing into the leads
! = 1 => no snow blowing, < 1 => some snow blowing
nn_flxdist = -1 ! Redistribute heat flux over ice categories
! =-1 Do nothing (needs N(cat) fluxes)
! = 0 Average N(cat) fluxes then apply the average over the N(cat) ice
! = 1 Average N(cat) fluxes then redistribute over the N(cat) ice using T-ice and albedo sensitivity
! = 2 Redistribute a single flux over categories
ln_cndflx = .false. ! Use conduction flux as surface boundary conditions (i.e. for Jules coupling)
ln_cndemulate = .false. ! emulate conduction flux (if not provided in the inputs)
nn_qtrice = 0 ! Solar flux transmitted thru the surface scattering layer:
! = 0 Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow)
! = 1 Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities)
/
! Variable setting
ctrcnm (jp_age) = 'Age'
ctrcln (jp_age) = 'Sea water age since surface contact'
ctrcun (jp_age) = 'year'
ln_trc_ini(jp_age) = .false.
ln_trc_sbc(jp_age) = .false.
ln_trc_cbc(jp_age) = .false.
ln_trc_obc(jp_age) = .false.
......@@ -12,8 +12,10 @@
! ! diagnostics:
ln_bergdia = .true. ! Calculate budgets
nn_verbose_level = 0 ! Turn on more verbose output if level > 0
nn_verbose_write = 15 ! Timesteps between verbose messages
nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage
!
! nn_verbose_write and nn_sample_rate need to be a multiple of nn_fsbc
nn_verbose_write = 16 ! Timesteps between verbose messages
nn_sample_rate = 16 ! Timesteps between sampling for trajectory storage
!
! ! iceberg setting:
! ! Initial mass required for an iceberg of each class
......
!-----------------------------------------------------------------------
&namdta_dyn ! offline ocean input files (OFF_SRC only)
!-----------------------------------------------------------------------
cn_dir = './' ! root directory for the ocean data location
!___________!_________________________!___________________!___________!_____________!________!___________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' !
sn_tem = 'dyna_grid_T' , 120. , 'votemper' , .true. , .true. , 'yearly'
sn_sal = 'dyna_grid_T' , 120. , 'vosaline' , .true. , .true. , 'yearly'
sn_mld = 'dyna_grid_T' , 120. , 'somixhgt' , .true. , .true. , 'yearly'
sn_emp = 'dyna_grid_T' , 120. , 'sowaflup' , .true. , .true. , 'yearly'
sn_fmf = 'dyna_grid_T' , 120. , 'iowaflup' , .true. , .true. , 'yearly'
sn_ice = 'dyna_grid_T' , 120. , 'soicecov' , .true. , .true. , 'yearly'
sn_qsr = 'dyna_grid_T' , 120. , 'soshfldo' , .true. , .true. , 'yearly'
sn_wnd = 'dyna_grid_T' , 120. , 'sowindsp' , .true. , .true. , 'yearly'
sn_uwd = 'dyna_grid_U' , 120. , 'uocetr_eff', .true. , .true. , 'yearly'
sn_vwd = 'dyna_grid_V' , 120. , 'vocetr_eff', .true. , .true. , 'yearly'
sn_wwd = 'dyna_grid_W' , 120. , 'wocetr_eff', .true. , .true. , 'yearly'
sn_avt = 'dyna_grid_W' , 120. , 'voddmavs' , .true. , .true. , 'yearly'
sn_ubl = 'dyna_grid_U' , 120. , 'sobblcox' , .true. , .true. , 'yearly'
sn_vbl = 'dyna_grid_V' , 120. , 'sobblcoy' , .true. , .true. , 'yearly'
/
!-----------------------------------------------------------------------
&namdta_dyn ! offline ocean input files (OFF_SRC only)
!-----------------------------------------------------------------------
ln_dynrnf = .true. ! runoffs option enabled (T) or not (F)
ln_dynrnf_depth = .false. ! runoffs is spread in vertical (T) or not (F)
!
cn_dir = './' ! root directory for the ocean data location
!___________!_________________________!___________________!_____________!_____________!________!___________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' !
sn_tem = 'dyna_grid_T' , 120. , 'thetao' , .true. , .true. , 'yearly'
sn_sal = 'dyna_grid_T' , 120. , 'so' , .true. , .true. , 'yearly'
sn_div = 'dyna_grid_T' , 120. , 'hdivtr' , .true. , .true. , 'yearly'
sn_mld = 'dyna_grid_T' , 120. , 'mldr10_1' , .true. , .true. , 'yearly'
sn_emp = 'dyna_grid_T' , 120. , 'wfo' , .true. , .true. , 'yearly'
sn_empb = 'dyna_grid_T' , 120. , 'wfob' , .true. , .true. , 'yearly'
sn_fmf = 'dyna_grid_T' , 120. , 'fmmflx' , .true. , .true. , 'yearly'
sn_rnf = 'dyna_grid_T' , 120. , 'runoffs' , .true. , .true. , 'yearly'
sn_ice = 'dyna_grid_T' , 120. , 'siconc' , .true. , .true. , 'yearly'
sn_qsr = 'dyna_grid_T' , 120. , 'rsntds' , .true. , .true. , 'yearly'
sn_wnd = 'dyna_grid_T' , 120. , 'windsp' , .true. , .true. , 'yearly'
sn_uwd = 'dyna_grid_U' , 120. , 'uocetr_eff', .true. , .true. , 'yearly'
sn_vwd = 'dyna_grid_V' , 120. , 'vocetr_eff', .true. , .true. , 'yearly'
sn_wwd = 'dyna_grid_W' , 120. , 'wocetr_eff', .true. , .true. , 'yearly'
sn_avt = 'dyna_grid_W' , 120. , 'difvsolog' , .true. , .true. , 'yearly'
sn_ubl = 'dyna_grid_U' , 120. , 'ahu_bbl' , .true. , .true. , 'yearly'
sn_vbl = 'dyna_grid_V' , 120. , 'ahv_bbl' , .true. , .true. , 'yearly'
/
:
:-----------------------------------------------------------------------------
: Atmospheric Boundary Layer preprocessing tool
:-----------------------------------------------------------------------------
:
&nml_dom
jpka = 50, ! ABL vertical levels number
hmax = 2000., ! ABL last level altitude
theta_s = 2., ! vertical grid stretching parameters
hc = 100., !
ln_impose_z1 = .true., ! force ABL first level altitude
z1 = 10., ! ABL firt level imposed altitude [m]
/
&nml_opt
ptemp_method = 3 , ! potential temperature computation method
ln_slp_smth = .true. , ! smooth slp and ghw at high latitudes only
ln_drw_smth = .false. , ! smooth after drowning
ln_slp_log = .false. , ! read log(slp)
ln_read_zsurf = .false. , ! read surface geopotential
ln_hpg_frc = .true. , ! compute horizontal pressure gradient
ln_geo_wnd = .false. , ! compute geostrophic wind
ln_c1d = .false. , ! 1D case
ln_read_mask = .true. , ! read mask file
ln_lsm_land = .false. , ! inverse land & sea masks
ln_perio_latbc = .true. , ! periodic lateral boundary conditions
/
&nml_fld
cn_dir = '',
mask_var = 'LSM',
file_m = 'MASK.nc',
file_u = 'U3D.nc',
file_v = 'V3D.nc',
file_t = 'T3D.nc',
file_q = 'Q3D.nc',
file_p = 'P2D.nc',
file_z = 'Z2D.nc',
file_geos = 'UVG_OUT.nc',
file_hpg = 'HPG_OUT.nc',
/
&nml_out
grd_file = 'dom_cfg_abl_L50Z10.nc',
abl_file = 'ABL_L50Z10_OUT.nc',
drwn_file = 'ABL_DRWN_L50Z10_OUT.nc',
var_name = '',
/
&nml_c1d
iloc = 283,
jloc = 52,
/
!-----------------------------------------------------------------------
&namisf ! Top boundary layer (ISF) (default: OFF)
!-----------------------------------------------------------------------
!
! ---------------- ice shelf melt formulation -------------------------------
!
ln_isf = .true. ! activate ice shelf module
!
! ---------------- cavities opened -------------------------------
!
ln_isfcav_mlt = .false. ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc)
cn_isfcav_mlt = '3eq' ! ice shelf melting formulation (spe/2eq/3eq/oasis)
! ! spe = fwfisf is read from a forcing field
! ! 2eq = ISOMIP like: 2 equations formulation (Hunter et al., 2006 for a short description)
! ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description
rn_htbl = 30. ! thickness of the top boundary layer (Losh et al. 2008)
! ! 0 => thickness of the tbl = thickness of the first wet cell
!
!* 'spe' and 'oasis' case
!---------------------------------------------------------------------------------------------
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' !
sn_isfcav_fwf = 'isfmlt_cav', -12. , 'fwflisf' , .false. , .true. , 'yearly'
!
! ---------------- cavities parametrised -------------------------------
!
ln_isfpar_mlt = .true. ! ice shelf melting parametrised
cn_isfpar_mlt = 'spe' ! ice shelf melting parametrisation (spe/bg03/oasis)
! ! spe = fwfisf is read from a forcing field
!
!---------------------------------------------------------------------------------------------
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' !
sn_isfpar_zmax = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' , -12 ,'sodepmax_isf' , .false. , .true. , 'yearly'
sn_isfpar_zmin = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' , -12 ,'sodepmin_isf' , .false. , .true. , 'yearly'
!* 'spe' and 'oasis' case
sn_isfpar_fwf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' , -12 ,'sornfisf', .false. , .true. , 'yearly'
!----------------------------------------------------------------------------------
&namldf_eke ! GEOMETRIC param. (total EKE equation) (nn_aei_ijk_t = 32)
!----------------------------------------------------------------------------------
rn_ekedis = 100. ! dissipation time scale of EKE [days]
nn_eke_dis = 0 ! dissipation option
! ! = 0 constant in space
! ! =-20 read in geom_diss_2D.nc file
rn_geom = 0.04 ! geometric parameterization master coefficient (>0 & <1)
rn_eke_init = 1.e-1 ! initial total EKE value
rn_eke_min = 1.e+0 ! background value of total EKE
rn_ross_min = 4.e+3 ! tapering of aeiv based on min Rossby radius [m]
! ! set to zero to not taper it
rn_eke_lap = 500. ! Laplacian diffusion coefficient of EKE
! ! this is in all options below, so set it to zero and nothing is done
rn_aeiv_min = 1.e+1 ! minimum bound of eiv coefficient
rn_aeiv_max = 1.5e+4 ! maximum bound of eiv coefficient
rn_SFmin = 1.0 ! minimum bound of Structure Function
rn_SFmax = 1.0 ! maximum bound of Structure Function
nn_eke_opt = 1 ! options for terms to include in EKE budget
! ! = 0 PE->EKE conversion, dissipation only
! ! = 1 as 0 but with advection
! ! = 2 as 1 but with additional KE->EKE conversion
! ! for testing purposes:
! ! = 88 only advection by depth-averaged flow
! ! = 99 only Laplacian diffusion
ln_adv_wav = .false. ! include advection at long Rossby speed
ln_beta_plane = .false. ! beta plane option for computing long Rossby speed (default: sphere option)
/
!-----------------------------------------------------------------------
&namtrc_opt ! light availability in the water column
!-----------------------------------------------------------------------
! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !
sn_par = 'par.orca' , 24 , 'fr_par' , .true. , .true. , 'yearly' , '' , '' , ''
cn_dir = './' ! root directory for the location of the dynamical files
ln_varpar = .true. ! Read PAR from file
parlux = 0.43 ! Fraction of shortwave as PAR
light_loc = 'center' ! Light location in the water cell ('center', 'integral')
/
......@@ -3,7 +3,6 @@
!-----------------------------------------------------------------------
pislopen = 2. ! P-I slope
pisloped = 2. ! P-I slope for diatoms
xadap = 0. ! Adaptation factor to low light
excretn = 0.05 ! excretion ratio of phytoplankton
excretd = 0.05 ! excretion ratio of diatoms
bresp = 0.033 ! Basal respiration rate
......
......@@ -7,7 +7,6 @@
excretn = 0.05 ! excretion ratio of phytoplankton
excretp = 0.05 ! excretion ratio of picophytoplankton
excretd = 0.05 ! excretion ratio of diatoms
xadap = 0. ! Adaptation factor to low light
bresp = 0.02 ! Basal respiration rate
thetannm = 0.3 ! Maximum Chl/N in nanophytoplankton
thetanpm = 0.3 ! Maximum Chl/N in picophytoplankton
......
!-----------------------------------------------------------------------
&namsbc_rnf ! runoffs (ln_rnf =T)
!-----------------------------------------------------------------------
ln_rnf_mouth = .false. ! specific treatment at rivers mouths
rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used (ln_rnf_mouth=T)
rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T)
rn_rfact = 1.e0 ! multiplicative factor for runoff
ln_rnf_icb = .true. ! read iceberg flux
cn_dir = './' ! root directory for the location of the runoff files
!---------------------------------------------------------------------------------------------
! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' !
sn_rnf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc', -1 , 'sorunoff', .true. , .true. , 'yearly'
sn_i_rnf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc', -1 , 'Icb_flux', .true. , .true. , 'yearly'
sn_cnf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc', 0 , 'socoefr' , .false. , .true. , 'yearly'
sn_s_rnf = 'runoffs' , 24. , 'rosaline', .true. , .true. , 'yearly'
sn_t_rnf = 'runoffs' , 24. , 'rotemper', .true. , .true. , 'yearly'
sn_dep_rnf = 'runoffs' , 0. , 'rodepth' , .false. , .true. , 'yearly'
......@@ -11,7 +11,6 @@
! ! = 20 F(i,j) =ldf_c2d
! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation
! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d
! ! = 32 F(i,j,t) = GEOMETRIC parameterization (=> fill namldf_eke)
! ! time invariant coefficients: aei0 = 1/2 Ue*Le
rn_Ue = 0.02 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
rn_Le = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10)
......
......@@ -19,7 +19,7 @@
!
jp_dia3d = 0 ! Number of 3D diagnostic variables
jp_dia2d = 0 ! Number of 2D diagnostic variables
!_____________!___________!_________________________________________!____________!________________!
! ! name ! title of the field ! units ! init from file !
! sn_tracer(1) = 'tracer ', 'Tracer Concentration ', ' - ' , .false.
/
!-----------------------------------------------------------------------
&namtrc_ais ! Representation of Antarctic Ice Sheet tracers supply
!-----------------------------------------------------------------------
rn_trafac(3) = 4.476e-07 ! ( 0.5e-3 / 55.85 * 0.05 )
!
nn_ais_tr = 1 ! tracer concentration in iceberg and ice shelf
! = 0 (null concentrations)
! = 1 prescribed concentrations
rn_icbdep = 120. ! Mean underwater depth of iceberg (m)
!----------------------------------------------------------------------
&namtrc_bc ! data for boundary conditions
!----------------------------------------------------------------------
! Surface and coastal input data must be specified individually for each tracer.
! Lateral open boundary inputs are specified for each BDY segment and related inputfiles
! must contain data for the variables with active open boundary (set through &namtrc).
! By default, variable names of BDY inputfiles are the same as in &namtrc, but using
! cn_tronam it is possible to specify alternative variable names to match with inputfiles.
! ! file name ! freq ! variable ! time interp. ! clim !'yearly' ! weights ! rotation | land/sea
! sn_trcobc(1) = 'bdy_dta_trc_North_y1980', -1 , 'dummy' , .false. , .true. ,'yearly' , '', '', ''
!-----------------------------------------------------------------------
cn_dir_sbc = './' ! root directory for the location of SURFACE data files
cn_dir_cbc = './' ! root directory for the location of COASTAL data files
......@@ -7,4 +15,5 @@
ln_rnf_ctl = .false. ! Remove runoff dilution on tracers with absent river load
rn_sbc_time = 86400. ! Time scaling factor for SBC data (seconds in a day)
rn_cbc_time = 86400. ! Time scaling factor for CBC data (seconds in a day)
! cn_tronam(1) = 'var1' ! Tracer-name to variable-name translation
/
!----------------------------------------------------------------------
&namtrc_bc ! data for boundary conditions
!-----------------------------------------------------------------------
! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' !
sn_trcsbc(2) = 'dust.orca.new' , -1 , 'dustfer' , .true. , .true. , 'yearly'
sn_trcsbc(3) = 'ndeposition.orca', -12 , 'ndep' , .false. , .true. , 'yearly'
rn_trsfac(2) = 6.266e-04 ! Multiplicative factor
rn_trsfac(3) = 5.4464e-01 !
!
sn_trccbc(1) = 'river.orca' , -12 , 'riverdic' , .true. , .true. , 'yearly'
sn_trccbc(2) = 'river.orca' , -12 , 'riverdfe' , .true. , .true. , 'yearly'
sn_trccbc(3) = 'river.orca' , -12 , 'riverdin' , .true. , .true. , 'yearly'
rn_trcfac(1) = 1.0 ! Multiplicative factor
rn_trcfac(2) = 1.0 !
rn_trcfac(3) = 1.0 !
rn_cbc_time = 3.1536e+7 ! Time scaling factor for CBC data (seconds in a year)
......@@ -8,4 +8,6 @@
! = 0 NO damping of tracers at open boudaries
! = 1 Only for tracers forced with external data
! = 2 Damping applied to all tracers
ln_zintobc = .false. ! T if a vertical interpolation is required. Variables gdep[t] and e3[t] must exist in the file
! automatically defined to T if the number of vertical levels in bdy dta /= jpk
/
!-----------------------------------------------------------------------
&namtrc ! tracers definition
!-----------------------------------------------------------------------
jp_bgc = 24
!
ln_pisces = .true.
ln_my_trc = .false.
ln_age = .false.
ln_cfc11 = .false.
ln_cfc12 = .false.
ln_c14 = .false.
!
!
ln_trcdta = .true. ! Initialisation from data input file (T) or not (F)
ln_trcdmp = .false. ! add a damping termn (T) or not (F)
ln_trcdmp_clo = .false. ! damping term (T) or not (F) on closed seas
ln_trcbc = .true. ! Surface, Lateral or Open Boundaries conditions
ln_trcais = .true. ! Antarctic Ice Sheet nutrient supply
! ! ! ! ! !
! ! name ! title of the field ! units ! init ! sbc ! cbc ! obc ! ais
sn_tracer(1) = 'DIC ' , 'Dissolved inorganic Concentration ', 'mol-C/L' , .true. , .false., .true. , .false. , .false.
sn_tracer(2) = 'Fer ' , 'Dissolved Iron Concentration ', 'mol-C/L' , .true. , .true. , .true. , .false. , .true.
sn_tracer(3) = 'NO3 ' , 'Nitrates Concentration ', 'mol-C/L' , .true. , .true. , .true. , .false. , .false.
sn_tracer(4) = 'PHY ' , 'Nanophytoplankton Concentration ', 'mol-C/L' , .false. , .false., .false., .false. , .false.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment