Changeset 2061 for CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ORCA2
- Timestamp:
- 04/26/13 13:01:38 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ORCA2
r970 r2061 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 2 !! NEMO/OPA : 1 - run manager (namrun) 3 !! namelists 2 - Domain (namzgr, namzgr_sco, namdom )3 !! namelists 2 - Domain (namzgr, namzgr_sco, namdom, namtsd) 4 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namsbc_alb) 5 !! namsbc_cpl, namtra_qsr, namsbc_rnf, 6 !! namsbc_apr, namsbc_ssr, namsbc_alb) 6 7 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 7 8 !! 5 - bottom boundary (nambfr, nambbc, nambbl) … … 9 10 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 10 11 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 11 !! 9 - diagnostics (namtrd, namgap, namspr, namflo, namptr) 12 !! 9 - miscellaneous (namsol, nammpp, nammpp_dyndist, namctl) 12 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namptr, namhsb) 13 !! 10 - miscellaneous (namsol, nammpp, nammpp_dyndist, namctl) 14 !! 11 - Obs & Assim (namobs, nam_asminc) 13 15 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 ! CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE.15 16 16 17 !!====================================================================== 17 18 !! *** Run management namelists *** 18 19 !!====================================================================== 19 !! namrun 20 !!====================================================================== 21 20 !! namrun parameters of the run 21 !!====================================================================== 22 ! 22 23 !----------------------------------------------------------------------- 23 24 &namrun ! parameters of the run 24 25 !----------------------------------------------------------------------- 25 nn_no = 0 ! job number 26 cn_exp = "ORCA2" ! AUTO - experience name 26 nn_no = 0 ! job number (no more used...) 27 cn_exp = "ORCA2" ! AUTO - experience name 27 28 nn_it000 = 1 ! AUTO - first time step 28 nn_itend = 5475 ! AUTO - last time step (std 5475)29 nn_date0 = 010101 ! AUTO - initial calendar date yymmdd (used if nn_rstctl=1)29 nn_itend = 5475 ! AUTO - last time step (std 5475) 30 nn_date0 = 010101 ! AUTO - date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 30 31 nn_leapy = 0 ! AUTO - Leap year calendar (1) or not (0) 32 ln_rstart = .false. ! AUTO - start from rest (F) or from a restart file (T) 33 nn_rstctl = 0 ! AUTO - restart control => activated only if ln_rstart = T 34 ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 35 ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 36 ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart 37 cn_ocerst_in = "restartopa" ! suffix of ocean restart name (input) 38 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 31 39 nn_istate = 0 ! output the initial state (1) or not (0) 32 40 nn_stock = 5475 ! AUTO - frequency of creation of a restart file (modulo referenced to 1) 33 41 nn_write = 5475 ! AUTO - frequency of write in the output file (modulo referenced to nn_it000) 34 42 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) 35 ln_mskland = .true. 43 ln_mskland = .true. ! mask land points in NetCDF outputs (costly: + ~15%) 36 44 ln_clobber = .false. ! clobber (overwrite) an existing file 37 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (working only with iom_nf90 routines) 38 ln_rstart = .false. ! AUTO - start from rest (F) or from a restart file (T) 39 nn_rstctl = 0 ! AUTO - restart control = 0 nn_it000 is not compared to the restart file value 40 ! = 1 use nn_date0 in namelist (not the value in the restart file) 41 ! = 2 calendar parameters read in the restart file 42 cn_ocerst_in = "restartopa"! suffix of ocean restart name (input) 43 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 44 / 45 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 46 / 47 45 48 !!====================================================================== 46 49 !! *** Domain namelists *** … … 49 52 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 50 53 !! namdom space and time domain (bathymetry, mesh, timestep) 51 !!====================================================================== 52 54 !! namtsd data: temperature & salinity 55 !!====================================================================== 56 ! 53 57 !----------------------------------------------------------------------- 54 58 &namzgr ! vertical coordinate … … 61 65 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate 62 66 !----------------------------------------------------------------------- 63 rn_sbot_min = 300. ! minimum depth of s-bottom surface (>0) (m) 64 rn_sbot_max = 5250. ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 67 ln_s_sh94 = .true. ! Song & Haidvogel 1994 hybrid S-sigma (T)| 68 ln_s_sf12 = .false. ! Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 69 ln_sigcrit = .false. ! use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 70 ! stretching coefficients for all functions 71 rn_sbot_min = 10.0 ! minimum depth of s-bottom surface (>0) (m) 72 rn_sbot_max = 7000.0 ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 73 rn_hc = 150.0 ! critical depth for transition to stretched coordinates 74 !!!!!!! Envelop bathymetry 75 rn_rmax = 0.3 ! maximum cut-off r-value allowed (0<r_max<1) 76 !!!!!!! SH94 stretching coefficients (ln_s_sh94 = .true.) 65 77 rn_theta = 6.0 ! surface control parameter (0<=theta<=20) 66 rn_thetb = 0.75 ! bottom control parameter (0<=thetb<= 1) 67 rn_rmax = 0.15 ! maximum cut-off r-value allowed (0<r_max<1) 68 ln_s_sigma = .false. ! hybrid s-sigma coordinates 69 rn_bb = 0.8 ! stretching with s-sigma 70 rn_hc = 150.0 ! critical depth with s-sigma 78 rn_bb = 0.8 ! stretching with SH94 s-sigma 79 !!!!!!! SF12 stretching coefficient (ln_s_sf12 = .true.) 80 rn_alpha = 4.4 ! stretching with SF12 s-sigma 81 rn_efold = 0.0 ! efold length scale for transition to stretched coord 82 rn_zs = 1.0 ! depth of surface grid box 83 ! bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b 84 rn_zb_a = 0.024 ! bathymetry scaling factor for calculating Zb 85 rn_zb_b = -0.2 ! offset for calculating Zb 86 !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 87 rn_thetb = 1.0 ! bottom control parameter (0<=thetb<= 1) 71 88 / 72 89 !----------------------------------------------------------------------- 73 90 &namdom ! space and time domain (bathymetry, mesh, timestep) 74 91 !----------------------------------------------------------------------- 75 nn_bathy = 1 ! compute (=0) or read(=1) the bathymetry file 76 nn_closea = 1 ! closed seas and lakes are removed (=0) or kept (=1) from the ORCA domain 77 nn_msh = 0 ! AUTO create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 78 rn_e3zps_min= 20. ! the thickness of the partial step is set larger than the minimum 79 rn_e3zps_rat= 0.1 ! of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 92 nn_bathy = 1 ! compute (=0) or read (=1) the bathymetry file 93 nn_closea = 0 ! remove (=0) or keep (=1) closed seas and lakes (ORCA) 94 nn_msh = 0 ! AUTO - create (=1) a mesh file or not (=0) 95 rn_hmin = -3. ! min depth of the ocean (>0) or min number of ocean level (<0) 96 rn_e3zps_min= 20. ! partial step thickness is set larger than the minimum of 97 rn_e3zps_rat= 0.1 ! rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 80 98 ! 81 rn_rdt = 5760. ! time step for the dynamics (and tracer if n acc=0) ==> 576082 nn_baro = 64 ! number of barotropic time step (for the split explicit algorithm)("key_dynspg_ts")99 rn_rdt = 5760. ! time step for the dynamics (and tracer if nn_acc=0) 100 nn_baro = 64 ! number of barotropic time step ("key_dynspg_ts") 83 101 rn_atfp = 0.1 ! asselin time filter parameter 84 102 nn_acc = 0 ! acceleration of convergence : =1 used, rdt < rdttra(k) 85 103 ! =0, not used, rdt = rdttra 86 rn_rdtmin = 28800. ! minimum time step on tracers (used if nacc=1) 87 rn_rdtmax = 28800. ! maximum time step on tracers (used if nacc=1) 88 rn_rdth = 800. ! depth variation of tracer time step (used if nacc=1) 104 rn_rdtmin = 28800. ! minimum time step on tracers (used if nn_acc=1) 105 rn_rdtmax = 28800. ! maximum time step on tracers (used if nn_acc=1) 106 rn_rdth = 800. ! depth variation of tracer time step (used if nn_acc=1) 107 / 108 !----------------------------------------------------------------------- 109 &namtsd ! data : Temperature & Salinity 110 !----------------------------------------------------------------------- 111 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim !'yearly' or ! weights ! rotation ! 112 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 113 sn_tem = 'data_1m_potential_temperature_nomask', -1,'votemper', .true. , .true., 'yearly' , ' ' , ' ' 114 sn_sal = 'data_1m_salinity_nomask' , -1,'vosaline', .true. , .true., 'yearly' , '' , ' ' 115 ! 116 cn_dir = './' ! root directory for the location of the runoff files 117 ln_tsd_init = .true. ! Initialisation of ocean T & S with T &S input data (T) or not (F) 118 ln_tsd_tradmp = .true. ! damping of ocean T & S toward T &S input data (T) or not (F) 89 119 / 90 120 !!====================================================================== 91 121 !! *** Surface Boundary Condition namelists *** 92 122 !!====================================================================== 93 !! namsbc surface boundary condition 94 !! namsbc_ana analytical formulation 95 !! namsbc_flx flux formulation 96 !! namsbc_clio CLIO bulk formulea formulation 97 !! namsbc_core CORE bulk formulea formulation 98 !! namsbc_cpl CouPLed formulation ("key_coupled") 99 !! namtra_qsr penetrative solar radiation 100 !! namsbc_rnf river runoffs 101 !! namsbc_ssr sea surface restoring term (for T and/or S) 102 !! namsbc_alb albedo parameters 103 !!====================================================================== 104 123 !! namsbc surface boundary condition 124 !! namsbc_ana analytical formulation 125 !! namsbc_flx flux formulation 126 !! namsbc_clio CLIO bulk formulae formulation 127 !! namsbc_core CORE bulk formulae formulation 128 !! namsbc_mfs MFS bulk formulae formulation 129 !! namsbc_cpl CouPLed formulation ("key_coupled") 130 !! namtra_qsr penetrative solar radiation 131 !! namsbc_rnf river runoffs 132 !! namsbc_apr Atmospheric Pressure 133 !! namsbc_ssr sea surface restoring term (for T and/or S) 134 !! namsbc_alb albedo parameters 135 !!====================================================================== 136 ! 105 137 !----------------------------------------------------------------------- 106 138 &namsbc ! Surface Boundary Condition (surface module) 107 139 !----------------------------------------------------------------------- 108 nn_fsbc = 5 ! frequency of surface boundary condition computation 109 ! (= the frequency of sea-ice model call) 110 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 111 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 112 ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio) 113 ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) 114 ln_cpl = .true. ! Coupled formulation (T => fill namsbc_cpl ) 140 nn_fsbc = 5 ! frequency of surface boundary condition computation 141 ! (also = the frequency of sea-ice model call) 142 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 143 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 144 ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio) 145 ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) 146 ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs ) 147 ln_cpl = .true. ! Coupled formulation (T => fill namsbc_cpl ) 148 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr ) 115 149 nn_ice = 2 ! =0 no ice boundary condition , 116 150 ! =1 use observed ice-cover , 117 ! =2 ice-model used ("key_lim3" or "key_lim2) 118 nn_ico_cpl = 0 ! ice-ocean coupling : =0 each nn_fsbc 119 ! =1 stresses recomputed each ocean time step ("key_lim3" only) 120 ! =2 combination of 0 and 1 cases ("key_lim3" only) 121 ln_dm2dc = .false. ! daily mean to diurnal cycle short wave (qsr) 122 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 123 ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 124 nn_fwb = 0 ! FreshWater Budget: =0 unchecked 125 ! =1 global mean of e-p-r set to zero at each time step 126 ! =2 annual global mean of e-p-r set to zero 127 ! =3 global emp set to zero and spread out over erp area 151 ! =2 ice-model used ("key_lim3" or "key_lim2) 152 nn_ice_embd = 0 ! =0 levitating ice (no mass exchange, concentration/dilution effect) 153 ! =1 levitating ice with mass and salt exchange but no presure effect 154 ! =2 embedded sea-ice (full salt and mass exchanges and pressure) 155 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 156 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 157 ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 158 nn_fwb = 0 ! FreshWater Budget: =0 unchecked 159 ! =1 global mean of e-p-r set to zero at each time step 160 ! =2 annual global mean of e-p-r set to zero 161 ! =3 global emp set to zero and spread out over erp area 162 ln_wave = .false. ! Activate coupling with wave (either Stokes Drift or Drag coefficient, or both) (T => fill namsbc_wave) 163 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => fill namsbc_wave) 164 ln_sdw = .false. ! Computation of 3D stokes drift (T => fill namsbc_wave) 128 165 / 129 166 !----------------------------------------------------------------------- … … 140 177 &namsbc_flx ! surface boundary condition : flux formulation 141 178 !----------------------------------------------------------------------- 142 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim! 'yearly'/ ! weights ! rotation !143 ! ! ! (if <0 months) ! name ! (logical) ! (T/F)! 'monthly' ! filename ! pairing !144 sn_utau = 'utau' , 24 , 'utau' , .false. , .false., 'yearly' , '' , ''145 sn_vtau = 'vtau' , 24 , 'vtau' , .false. , .false., 'yearly' , '' , ''146 sn_qtot = 'qtot' , 24 , 'qtot' , .false. , .false., 'yearly' , '' , ''147 sn_qsr = 'qsr' , 24 , 'qsr' , .false. , .false., 'yearly' , '' , ''148 sn_emp = 'emp' , 24 , 'emp' , .false. , .false., 'yearly' , '' , ''149 ! 179 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 180 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 181 sn_utau = 'utau' , 24 , 'utau' , .false. , .false., 'yearly' , '' , '' 182 sn_vtau = 'vtau' , 24 , 'vtau' , .false. , .false., 'yearly' , '' , '' 183 sn_qtot = 'qtot' , 24 , 'qtot' , .false. , .false., 'yearly' , '' , '' 184 sn_qsr = 'qsr' , 24 , 'qsr' , .false. , .false., 'yearly' , '' , '' 185 sn_emp = 'emp' , 24 , 'emp' , .false. , .false., 'yearly' , '' , '' 186 150 187 cn_dir = './' ! root directory for the location of the flux files 151 / 152 !----------------------------------------------------------------------- 153 &namsbc_clio ! namsbc_clio CLIO bulk formul ea154 !----------------------------------------------------------------------- 155 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim! 'yearly'/ ! weights ! rotation !156 ! ! ! (if <0 months) ! name ! (logical) ! (T/F)! 'monthly' ! filename ! pairing !157 sn_utau = 'taux_1m' , -1 , 'sozotaux' , .true. , .true., 'yearly' , '' , ''158 sn_vtau = 'tauy_1m' , -1 , 'sometauy' , .true. , .true., 'yearly' , '' , ''159 sn_wndm = 'flx' , -1 , 'socliowi' , .true. , .true., 'yearly' , '' , ''160 sn_tair = 'flx' , -1 , 'socliot2' , .true. , .true., 'yearly' , '' , ''161 sn_humi = 'flx' , -1 , 'socliohu' , .true. , .true., 'yearly' , '' , ''162 sn_ccov = 'flx' , -1 , 'socliocl' , .false. , .true., 'yearly' , '' , ''163 sn_prec = 'flx' , -1 , 'socliopl' , .false. , .true., 'yearly' , '' , ''164 ! 188 / 189 !----------------------------------------------------------------------- 190 &namsbc_clio ! namsbc_clio CLIO bulk formulae 191 !----------------------------------------------------------------------- 192 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 193 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 194 sn_utau = 'taux_1m' , -1 , 'sozotaux', .true. , .true. , 'yearly' , '' , '' 195 sn_vtau = 'tauy_1m' , -1 , 'sometauy', .true. , .true. , 'yearly' , '' , '' 196 sn_wndm = 'flx' , -1 , 'socliowi', .true. , .true. , 'yearly' , '' , '' 197 sn_tair = 'flx' , -1 , 'socliot2', .true. , .true. , 'yearly' , '' , '' 198 sn_humi = 'flx' , -1 , 'socliohu', .true. , .true. , 'yearly' , '' , '' 199 sn_ccov = 'flx' , -1 , 'socliocl', .false. , .true. , 'yearly' , '' , '' 200 sn_prec = 'flx' , -1 , 'socliopl', .false. , .true. , 'yearly' , '' , '' 201 165 202 cn_dir = './' ! root directory for the location of the bulk files are 166 203 / 167 204 !----------------------------------------------------------------------- 168 &namsbc_core ! namsbc_core CORE bulk formul ea169 !----------------------------------------------------------------------- 170 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights! rotation !171 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename! pairing !172 sn_wndi = 'u 10_core' , -1 , 'u10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'U1'173 sn_wndj = 'v 10_core' , -1 , 'v10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'V1'174 sn_qsr = ' qsw_core' , -1 , 'swdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''175 sn_qlw = ' qlw_core' , -1 , 'lwdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''176 sn_tair = 't 2_core' , -1 , 't2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''177 sn_humi = 'q 2_core' , -1 , 'q2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''178 sn_prec = ' precip_core', -1 , 'precip' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''179 sn_snow = ' snow_core' , -1 , 'snow' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''180 sn_tdif = 'taudif_core' , 24 , 'taudif' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''181 ! 205 &namsbc_core ! namsbc_core CORE bulk formulae 206 !----------------------------------------------------------------------- 207 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 208 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 209 sn_wndi = 'u_10.15JUNE2009_orca2' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , '' , 'Uwnd' 210 sn_wndj = 'v_10.15JUNE2009_orca2' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , '' , 'Vwnd' 211 sn_qsr = 'ncar_rad.15JUNE2009_orca2' , 24 , 'SWDN_MOD', .false. , .true. , 'yearly' , '' , '' 212 sn_qlw = 'ncar_rad.15JUNE2009_orca2' , 24 , 'LWDN_MOD', .false. , .true. , 'yearly' , '' , '' 213 sn_tair = 't_10.15JUNE2009_orca2' , 6 , 'T_10_MOD', .false. , .true. , 'yearly' , '' , '' 214 sn_humi = 'q_10.15JUNE2009_orca2' , 6 , 'Q_10_MOD', .false. , .true. , 'yearly' , '' , '' 215 sn_prec = 'ncar_precip.15JUNE2009_orca2', -1 , 'PRC_MOD1', .false. , .true. , 'yearly' , '' , '' 216 sn_snow = 'ncar_precip.15JUNE2009_orca2', -1 , 'SNOW' , .false. , .true. , 'yearly' , '' , '' 217 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , '' , '' 218 182 219 cn_dir = './' ! root directory for the location of the bulk files 183 ln_2m = . true.! air temperature and humidity referenced at 2m (T) instead 10m (F)184 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data ?220 ln_2m = .false. ! air temperature and humidity referenced at 2m (T) instead 10m (F) 221 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 185 222 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 186 223 / 187 224 !----------------------------------------------------------------------- 188 &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") 189 !----------------------------------------------------------------------- 190 ! send 191 cn_snd_temperature= 'weighted oce and ice' ! 'oce only' 'weighted oce and ice' 'mixed oce-ice' 192 cn_snd_albedo = 'weighted ice' ! 'none' 'weighted ice' 'mixed oce-ice' 193 cn_snd_thickness = 'none' ! 'none' 'weighted ice and snow' 194 cn_snd_crt_nature = 'mixed oce-ice' ! 'none' 'oce only' 'weighted oce and ice' 'mixed oce-ice' 195 cn_snd_crt_refere = 'cartesian' ! 'spherical' 'cartesian' 196 cn_snd_crt_orient = 'eastward-northward' ! 'eastward-northward' or 'local grid' 197 cn_snd_crt_grid = 'T' ! 'T' 198 ! receive 199 cn_rcv_w10m = 'coupled' ! 'none' 'coupled' 200 cn_rcv_taumod = 'none' ! 'none' 'coupled' 201 cn_rcv_tau_nature = 'mixed oce-ice' ! 'oce only' 'oce and ice' 'mixed oce-ice' 202 cn_rcv_tau_refere = 'cartesian' ! 'spherical' 'cartesian' 203 cn_rcv_tau_orient = 'eastward-northward' ! 'eastward-northward' or 'local grid' 204 cn_rcv_tau_grid = 'U,V' ! 'T' 'U,V' 'U,V,F' 'U,V,I' 'T,F' 'T,I' 'T,U,V' 205 cn_rcv_dqnsdt = 'coupled' ! 'none' 'coupled' 206 cn_rcv_qsr = 'conservative' ! 'conservative' 'oce and ice' 'mixed oce-ice' 207 cn_rcv_qns = 'conservative' ! 'conservative' 'oce and ice' 'mixed oce-ice' 208 cn_rcv_emp = 'conservative' ! 'conservative' 'oce and ice' 'mixed oce-ice' 209 cn_rcv_rnf = 'coupled' ! 'coupled' 'climato' 'mixed' 210 cn_rcv_cal = 'coupled' ! 'none' 'coupled' 211 / 212 !----------------------------------------------------------------------- 213 &namsbc_cpl_co2 ! coupled ocean/biogeo/atmosphere model ("key_cpl_carbon_cycle") 214 !----------------------------------------------------------------------- 215 cn_snd_co2 = 'coupled' ! send : 'none' 'coupled' 216 cn_rcv_co2 = 'coupled' ! receive : 'none' 'coupled' 225 &namsbc_mfs ! namsbc_mfs MFS bulk formulae 226 !----------------------------------------------------------------------- 227 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 228 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 229 sn_wndi = 'ecmwf' , 6 , 'u10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' 230 sn_wndj = 'ecmwf' , 6 , 'v10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' 231 sn_clc = 'ecmwf' , 6 , 'clc' , .true. , .false. , 'daily' ,'bilinear.nc', '' 232 sn_msl = 'ecmwf' , 6 , 'msl' , .true. , .false. , 'daily' ,'bicubic.nc' , '' 233 sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false. , 'daily' ,'bicubic.nc' , '' 234 sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false. , 'daily' ,'bilinear.nc', '' 235 sn_prec = 'ecmwf' , 6 , 'precip' , .true. , .true. , 'daily' ,'bicubic.nc' , '' 236 237 cn_dir = './ECMWF/' ! root directory for the location of the bulk files 238 / 239 !----------------------------------------------------------------------- 240 &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") 241 !----------------------------------------------------------------------- 242 ! ! description ! multiple ! vector ! vector ! vector ! 243 ! ! ! categories ! reference ! orientation ! grids ! 244 ! send 245 sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , '' 246 sn_snd_alb = 'weighted ice' , 'no' , '' , '' , '' 247 sn_snd_thick = 'none' , 'no' , '' , '' , '' 248 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T' 249 sn_snd_co2 = 'coupled' , 'no' , '' , '' , '' 250 ! receive 251 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' 252 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' 253 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' 254 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' 255 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' 256 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' 257 sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' 258 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' 259 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 260 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 217 261 / 218 262 !----------------------------------------------------------------------- 219 263 &namtra_qsr ! penetrative solar radiation 220 264 !----------------------------------------------------------------------- 221 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim! 'yearly'/ ! weights ! rotation !222 ! ! ! (if <0 months) ! name ! (logical) ! (T/F)! 'monthly' ! filename ! pairing !223 sn_chl = 'chlorophyll', -1 , 'CHLA' , .true. , .true., 'yearly' , '' , ''224 265 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 266 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 267 sn_chl ='chlorophyll', -1 , 'CHLA' , .true. , .true. , 'yearly' , '' , '' 268 225 269 cn_dir = './' ! root directory for the location of the runoff files 226 270 ln_traqsr = .true. ! Light penetration (T) or not (F) … … 228 272 ln_qsr_2bd = .false. ! 2 bands light penetration 229 273 ln_qsr_bio = .true. ! bio-model light penetration 230 nn_chldta = 0 ! RGB : Chl 2D data (=1), 3D data (=2) or cst value (=0)274 nn_chldta = 1 ! RGB : Chl data (=1) or cst value (=0) 231 275 rn_abs = 0.58 ! RGB & 2 bands: fraction of light (rn_si1) 232 276 rn_si0 = 0.35 ! RGB & 2 bands: shortess depth of extinction 233 277 rn_si1 = 23.0 ! 2 bands: longest depth of extinction 234 rn_si2 = 62.0 ! 3 bands: longest depth of extinction (for blue waveband & 0.01 mg/m2 Chl)235 278 / 236 279 !----------------------------------------------------------------------- 237 280 &namsbc_rnf ! runoffs namelist surface boundary condition 238 281 !----------------------------------------------------------------------- 239 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 240 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 241 sn_rnf = 'runoff_1m_nomask' , -1 , 'sorunoff' , .true. , .true. , 'yearly' , '' , '' 242 sn_cnf = 'runoff_1m_nomask' , 0 , 'socoefr' , .false. , .true. , 'yearly' , '' , '' 243 282 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 283 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 284 sn_rnf = 'runoff_core_monthly', -1 , 'sorunoff', .true. , .true. , 'yearly' , '' , '' 285 sn_cnf = 'runoff_core_monthly', 0 , 'socoefr0', .false. , .true. , 'yearly' , '' , '' 286 sn_s_rnf = 'runoffs' , 24 , 'rosaline', .true. , .true. , 'yearly' , '' , '' 287 sn_t_rnf = 'runoffs' , 24 , 'rotemper', .true. , .true. , 'yearly' , '' , '' 288 sn_dep_rnf = 'runoffs' , 0 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' 289 244 290 cn_dir = './' ! root directory for the location of the runoff files 245 ln_rnf_emp = .true.! runoffs included into precipitation field (T) or into a file (F)246 ln_rnf_mouth = .true.! specific treatment at rivers mouths247 rn_hrnf = 15.e0! depth over which enhanced vertical mixing is used291 ln_rnf_emp = .false. ! runoffs included into precipitation field (T) or into a file (F) 292 ln_rnf_mouth = .true. ! specific treatment at rivers mouths 293 rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used 248 294 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] 249 295 rn_rfact = 1.e0 ! multiplicative factor for runoff 296 ln_rnf_depth = .false. ! read in depth information for runoff 297 ln_rnf_tem = .false. ! read in temperature information for runoff 298 ln_rnf_sal = .false. ! read in salinity information for runoff 299 / 300 !----------------------------------------------------------------------- 301 &namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk 302 !----------------------------------------------------------------------- 303 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 304 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 305 sn_apr = 'patm' , -1 ,'somslpre', .true. , .true. , 'yearly' , '' , '' 306 307 cn_dir = './' ! root directory for the location of the bulk files 308 rn_pref = 101000._wp ! reference atmospheric pressure [N/m2]/ 309 ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F) 310 ln_apr_obc = .false. ! inverse barometer added to OBC ssh data 250 311 / 251 312 !----------------------------------------------------------------------- 252 313 &namsbc_ssr ! surface boundary condition : sea surface restoring 253 314 !----------------------------------------------------------------------- 254 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim! 'yearly'/ ! weights ! rotation !255 ! ! ! (if <0 months) ! name ! (logical) ! (T/F)! 'monthly' ! filename ! pairing !256 sn_sst = 'sst_data' , 24 , 'sst' , .false. , .false., 'yearly' , '' , ''257 sn_sss = 'sss_data' , -1 , 'sss' , .true. , .true., 'yearly' , '' , ''258 315 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 316 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 317 sn_sst = 'sst_data' , 24 , 'sst' , .false. , .false., 'yearly' , '' , '' 318 sn_sss = 'sss_data' , -1 , 'sss' , .true. , .true. , 'yearly' , '' , '' 319 259 320 cn_dir = './' ! root directory for the location of the runoff files 260 321 nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0) 261 nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) 322 nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) 262 323 ! or to SSS only (=1) or no damping term (=0) 263 324 rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 264 rn_deds = -27.7 ! magnitude of the damping on salinity [mm/day/psu]325 rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day] 265 326 ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) 266 327 rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] 267 / 328 / 268 329 !----------------------------------------------------------------------- 269 330 &namsbc_alb ! albedo parameters 270 331 !----------------------------------------------------------------------- 271 rn_cloud = 0.0 ! cloud correction to snow and ice albedo272 rn_albice = 0.5 332 rn_cloud = 0.06 ! cloud correction to snow and ice albedo 333 rn_albice = 0.53 ! albedo of melting ice in the arctic and antarctic 273 334 rn_alphd = 0.80 ! coefficients for linear interpolation used to 274 rn_alphc = 0.65 ! compute albedo between two extremes values 335 rn_alphc = 0.65 ! compute albedo between two extremes values 275 336 rn_alphdi = 0.72 ! (Pyane, 1972) 337 / 338 !----------------------------------------------------------------------- 339 &namberg ! iceberg parameters 340 !----------------------------------------------------------------------- 341 ln_icebergs = .false. 342 ln_bergdia = .true. ! Calculate budgets 343 nn_verbose_level = 1 ! Turn on more verbose output if level > 0 344 nn_verbose_write = 15 ! Timesteps between verbose messages 345 nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage 346 ! Initial mass required for an iceberg of each class 347 rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 348 ! Proportion of calving mass to apportion to each class 349 rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 350 ! Ratio between effective and real iceberg mass (non-dim) 351 ! i.e. number of icebergs represented at a point 352 rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 353 ! thickness of newly calved bergs (m) 354 rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 355 rn_rho_bergs = 850. ! Density of icebergs 356 rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs 357 ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics 358 rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits 359 rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 360 ln_passive_mode = .false. ! iceberg - ocean decoupling 361 nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no) 362 ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 363 rn_test_box = 108.0, 116.0, -66.0, -58.0 364 rn_speed_limit = 0. ! CFL speed limit for a berg 365 366 ! filename ! freq (hours) ! variable ! time interp. ! clim !'yearly' or ! weights ! rotation ! 367 ! ! (<0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 368 sn_icb = 'calving' , -1 , 'calvingmask', .true. , .true., 'yearly' , ' ' , ' ' 369 370 cn_dir = './' 276 371 / 277 372 … … 282 377 !! namcla cross land advection 283 378 !! namobc open boundaries parameters ("key_obc") 284 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") 379 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") 285 380 !! nambdy Unstructured open boundaries ("key_bdy") 286 381 !! namtide Tidal forcing at open boundaries ("key_bdy_tides") 287 382 !!====================================================================== 288 383 ! 289 384 !----------------------------------------------------------------------- 290 385 &namlbc ! lateral momentum boundary condition … … 292 387 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 293 388 ! free slip ! partial slip ! no slip ! strong slip 389 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical eqs. 294 390 / 295 391 !----------------------------------------------------------------------- … … 301 397 &namobc ! open boundaries parameters ("key_obc") 302 398 !----------------------------------------------------------------------- 303 ln_obc_clim= .false. ! climatological obc data files (T) or not (F)304 ln_vol_cst= .true. ! impose the total volume conservation (T) or not (F)305 ln_obc_fla = .false. ! Flather open boundary condition306 nn_obcdta= 1 ! = 0 the obc data are equal to the initial state399 ln_obc_clim = .false. ! climatological obc data files (T) or not (F) 400 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) 401 ln_obc_fla = .false. ! Flather open boundary condition 402 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state 307 403 ! = 1 the obc data are read in 'obc.dta' files 308 cn_obcdta= 'annual' ! set to annual if obc datafile hold 1 year of data404 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data 309 405 ! set to monthly if obc datafile hold 1 month of data 310 rn_dpein= 1. ! damping time scale for inflow at east open boundary311 rn_dpwin= 1. ! - - - west - -312 rn_dpnin= 1. ! - - - north - -313 rn_dpsin= 1. ! - - - south - -314 rn_dpeob= 3000. ! time relaxation (days) for the east open boundary315 rn_dpwob= 15. ! - - - west - -316 rn_dpnob= 3000. ! - - - north - -317 rn_dpsob= 15. ! - - - south - -318 rn_volemp = 1.! = 0 the total volume change with the surface flux (E-P-R)406 rn_dpein = 1. ! damping time scale for inflow at east open boundary 407 rn_dpwin = 1. ! - - - west - - 408 rn_dpnin = 1. ! - - - north - - 409 rn_dpsin = 1. ! - - - south - - 410 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary 411 rn_dpwob = 15. ! - - - west - - 412 rn_dpnob = 3000. ! - - - north - - 413 rn_dpsob = 15. ! - - - south - - 414 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R) 319 415 ! = 1 the total volume remains constant 320 416 / … … 322 418 &namagrif ! AGRIF zoom ("key_agrif") 323 419 !----------------------------------------------------------------------- 324 nn_cln_update = 3 ! baroclinic update frequency 325 ln_spc_dyn = .true. ! use 0 as special value for dynamics 326 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [s] 327 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [s] 420 nn_cln_update = 3 ! baroclinic update frequency 421 ln_spc_dyn = .true. ! use 0 as special value for dynamics 422 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 423 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 424 / 425 !----------------------------------------------------------------------- 426 &nam_tide ! tide parameters (#ifdef key_tide) 427 !----------------------------------------------------------------------- 428 ln_tide_pot = .true. ! use tidal potential forcing 429 nb_harmo = 11 ! number of constituents used 430 clname(1) = 'M2' ! name of constituent 431 clname(2) = 'S2' 432 clname(3) = 'N2' 433 clname(4) = 'K1' 434 clname(5) = 'O1' 435 clname(6) = 'Q1' 436 clname(7) = 'M4' 437 clname(8) = 'K2' 438 clname(9) = 'P1' 439 clname(10) = 'Mf' 440 clname(11) = 'Mm' 328 441 / 329 442 !----------------------------------------------------------------------- 330 443 &nambdy ! unstructured open boundaries ("key_bdy") 331 444 !----------------------------------------------------------------------- 332 filbdy_mask = '' ! name of mask file (if ln_bdy_mask=.TRUE.) 333 filbdy_data_T = 'bdydata_grid_T.nc' ! name of data file (T-points) 334 filbdy_data_U = 'bdydata_grid_U.nc' ! name of data file (U-points) 335 filbdy_data_V = 'bdydata_grid_V.nc' ! name of data file (V-points) 336 ln_bdy_clim = .false. ! contain 1 (T) or 12 (F) time dumps and be cyclic 337 ln_bdy_vol = .true. ! total volume correction (see volbdy parameter) 338 ln_bdy_mask = .false. ! boundary mask from filbdy_mask (T) or boundaries are on edges of domain (F) 339 ln_bdy_tides = .true. ! Apply tidal harmonic forcing with Flather condition 340 ln_bdy_dyn_fla = .true. ! Apply Flather condition to velocities 341 ln_bdy_tra_frs = .false. ! Apply FRS condition to temperature and salinity 342 ln_bdy_dyn_frs = .false. ! Apply FRS condition to velocities 343 nbdy_dta = 1 ! = 0, bdy data are equal to the initial state 445 nb_bdy = 1 ! number of open boundary sets 446 ln_coords_file = .true. ! =T : read bdy coordinates from file 447 cn_coords_file = 'coordinates.bdy.nc' ! bdy coordinates files 448 ln_mask_file = .false. ! =T : read mask from file 449 cn_mask_file = '' ! name of mask file (if ln_mask_file=.TRUE.) 450 nn_dyn2d = 2 ! boundary conditions for barotropic fields 451 nn_dyn2d_dta = 3 ! = 0, bdy data are equal to the initial state 344 452 ! = 1, bdy data are read in 'bdydata .nc' files 345 nb_rimwidth = 9 ! width of the relaxation zone 346 volbdy = 0 ! = 0, the total water flux across open boundaries is zero 347 ! = 1, the total volume of the system is conserved 348 / 349 !----------------------------------------------------------------------- 350 &nambdy_tide ! tidal forcing at unstructured boundaries 351 !----------------------------------------------------------------------- 352 filtide = 'bdytide_' ! file name root of tidal forcing files 353 tide_cpt = 'M2','S1' ! names of tidal components used 354 tide_speed = 28.984106, 15.000001 ! phase speeds of tidal components (deg/hour) 355 ln_tide_date = .false. ! adjust tidal harmonics for start date of run 356 / 357 453 ! = 2, use tidal harmonic forcing data from files 454 ! = 3, use external data AND tidal harmonic forcing 455 nn_dyn3d = 0 ! boundary conditions for baroclinic velocities 456 nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state 457 ! = 1, bdy data are read in 'bdydata .nc' files 458 nn_tra = 1 ! boundary conditions for T and S 459 nn_tra_dta = 1 ! = 0, bdy data are equal to the initial state 460 ! = 1, bdy data are read in 'bdydata .nc' files 461 nn_rimwidth = 10 ! width of the relaxation zone 462 ln_vol = .false. ! total volume correction (see nn_volctl parameter) 463 nn_volctl = 1 ! = 0, the total water flux across open boundaries is zero 464 / 465 !----------------------------------------------------------------------- 466 &nambdy_dta ! open boundaries - external data ("key_bdy") 467 !----------------------------------------------------------------------- 468 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 469 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 470 bn_ssh = 'amm12_bdyT_u2d' , 24 , 'sossheig' , .true. , .false. , 'daily' , '' , '' 471 bn_u2d = 'amm12_bdyU_u2d' , 24 , 'vobtcrtx' , .true. , .false. , 'daily' , '' , '' 472 bn_v2d = 'amm12_bdyV_u2d' , 24 , 'vobtcrty' , .true. , .false. , 'daily' , '' , '' 473 bn_u3d = 'amm12_bdyU_u3d' , 24 , 'vozocrtx' , .true. , .false. , 'daily' , '' , '' 474 bn_v3d = 'amm12_bdyV_u3d' , 24 , 'vomecrty' , .true. , .false. , 'daily' , '' , '' 475 bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper' , .true. , .false. , 'daily' , '' , '' 476 bn_sal = 'amm12_bdyT_tra' , 24 , 'vosaline' , .true. , .false. , 'daily' , '' , '' 477 cn_dir = 'bdydta/' 478 ln_full_vel = .false. 479 / 480 !----------------------------------------------------------------------- 481 &nambdy_tide ! tidal forcing at open boundaries 482 !----------------------------------------------------------------------- 483 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 484 tide_cpt(1) ='Q1' ! names of tidal components used 485 tide_cpt(2) ='O1' ! names of tidal components used 486 tide_cpt(3) ='P1' ! names of tidal components used 487 tide_cpt(4) ='S1' ! names of tidal components used 488 tide_cpt(5) ='K1' ! names of tidal components used 489 tide_cpt(6) ='2N2' ! names of tidal components used 490 tide_cpt(7) ='MU2' ! names of tidal components used 491 tide_cpt(8) ='N2' ! names of tidal components used 492 tide_cpt(9) ='NU2' ! names of tidal components used 493 tide_cpt(10) ='M2' ! names of tidal components used 494 tide_cpt(11) ='L2' ! names of tidal components used 495 tide_cpt(12) ='T2' ! names of tidal components used 496 tide_cpt(13) ='S2' ! names of tidal components used 497 tide_cpt(14) ='K2' ! names of tidal components used 498 tide_cpt(15) ='M4' ! names of tidal components used 499 tide_speed(1) = 13.398661 ! phase speeds of tidal components (deg/hour) 500 tide_speed(2) = 13.943036 ! phase speeds of tidal components (deg/hour) 501 tide_speed(3) = 14.958932 ! phase speeds of tidal components (deg/hour) 502 tide_speed(4) = 15.000001 ! phase speeds of tidal components (deg/hour) 503 tide_speed(5) = 15.041069 ! phase speeds of tidal components (deg/hour) 504 tide_speed(6) = 27.895355 ! phase speeds of tidal components (deg/hour) 505 tide_speed(7) = 27.968210 ! phase speeds of tidal components (deg/hour) 506 tide_speed(8) = 28.439730 ! phase speeds of tidal components (deg/hour) 507 tide_speed(9) = 28.512585 ! phase speeds of tidal components (deg/hour) 508 tide_speed(10) = 28.984106 ! phase speeds of tidal components (deg/hour) 509 tide_speed(11) = 29.528479 ! phase speeds of tidal components (deg/hour) 510 tide_speed(12) = 29.958935 ! phase speeds of tidal components (deg/hour) 511 tide_speed(13) = 30.000002 ! phase speeds of tidal components (deg/hour) 512 tide_speed(14) = 30.082138 ! phase speeds of tidal components (deg/hour) 513 tide_speed(15) = 57.968212 ! phase speeds of tidal components (deg/hour) 514 ln_tide_date = .true. ! adjust tidal harmonics for start date of run 515 / 358 516 !!====================================================================== 359 517 !! *** Bottom boundary condition *** 360 518 !!====================================================================== 361 519 !! nambfr bottom friction 362 !! nambbc bottom temperature boundary condition ("key_trabbc")363 !! nambbl bottom boundary layer scheme ("key_trabbl _dif","key_trabbl_adv")364 !!====================================================================== 365 520 !! nambbc bottom temperature boundary condition 521 !! nambbl bottom boundary layer scheme ("key_trabbl") 522 !!====================================================================== 523 ! 366 524 !----------------------------------------------------------------------- 367 525 &nambfr ! bottom friction … … 371 529 rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 372 530 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) 373 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2) 374 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 375 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d = .true.) 531 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) 532 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 533 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T) 534 ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true) 376 535 / 377 536 !----------------------------------------------------------------------- 378 537 &nambbc ! bottom temperature boundary condition 379 538 !----------------------------------------------------------------------- 380 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 539 ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom 540 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 381 541 ! = 1 constant flux 382 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 542 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 383 543 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 384 544 / … … 386 546 &nambbl ! bottom boundary layer scheme 387 547 !----------------------------------------------------------------------- 388 ! ! diffusive bbl ("key_trabbl") 389 ! ! advective bbl ("key_trabbl_adv") 390 rn_ahtbbl = 10000. ! lateral mixing coefficient in the bbl [m2/s] 391 / 548 nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0) 549 nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0) 550 rn_ahtbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s] 551 rn_gambbl = 10. ! advective bbl coefficient [s] 552 / 553 392 554 !!====================================================================== 393 555 !! Tracer (T & S ) namelists … … 396 558 !! namtra_adv advection scheme 397 559 !! namtra_ldf lateral diffusion scheme 398 !! namtra_dmp T & S newtonian damping ("key_tradmp")399 !!====================================================================== 400 560 !! namtra_dmp T & S newtonian damping 561 !!====================================================================== 562 ! 401 563 !----------------------------------------------------------------------- 402 564 &nameos ! ocean physical parameters 403 565 !----------------------------------------------------------------------- 404 nn_eos = 0! type of equation of state and Brunt-Vaisala frequency566 nn_eos = 0 ! type of equation of state and Brunt-Vaisala frequency 405 567 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 406 568 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 407 569 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 408 rn_alpha = 2.e-4 ! thermal expension coefficient (neos= 1 or 2) 409 rn_beta = 0.001 ! saline expension coefficient (neos= 2) 410 / 411 !----------------------------------------------------------------------- 412 &namtra_adv ! advection scheme for tracer 413 !----------------------------------------------------------------------- 414 ln_traadv_cen2 = .false. ! 2nd order centered scheme 415 ln_traadv_tvd = .true. ! TVD scheme 416 ln_traadv_muscl = .false. ! MUSCL scheme 417 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 418 ln_traadv_ubs = .false. ! UBS scheme 419 / 420 !----------------------------------------------------------------------- 421 &namtra_ldf ! lateral diffusion scheme for tracer 422 !----------------------------------------------------------------------- 423 ! Type of the operator : 424 ln_traldf_lap = .true. ! laplacian operator 425 ln_traldf_bilap = .false. ! bilaplacian operator 426 ! Direction of action : 427 ln_traldf_level = .false. ! iso-level 428 ln_traldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 429 ln_traldf_iso = .true. ! iso-neutral (require "key_ldfslp") 430 ! Coefficient 431 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 432 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 433 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 434 / 435 !----------------------------------------------------------------------- 436 &namtra_dmp ! tracer: T & S newtonian damping ('key_tradmp') 437 !----------------------------------------------------------------------- 570 rn_alpha = 2.0e-4 ! thermal expension coefficient (nn_eos= 1 or 2) 571 rn_beta = 7.7e-4 ! saline expension coefficient (nn_eos= 2) 572 / 573 !----------------------------------------------------------------------- 574 &namtra_adv ! advection scheme for tracer 575 !----------------------------------------------------------------------- 576 ln_traadv_cen2 = .false. ! 2nd order centered scheme 577 ln_traadv_tvd = .true. ! TVD scheme 578 ln_traadv_muscl = .false. ! MUSCL scheme 579 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 580 ln_traadv_ubs = .false. ! UBS scheme 581 ln_traadv_qck = .false. ! QUICKEST scheme 582 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 583 / 584 !---------------------------------------------------------------------------------- 585 &namtra_ldf ! lateral diffusion scheme for tracers 586 !---------------------------------------------------------------------------------- 587 ! ! Operator type: 588 ln_traldf_lap = .true. ! laplacian operator 589 ln_traldf_bilap = .false. ! bilaplacian operator 590 ! ! Direction of action: 591 ln_traldf_level = .false. ! iso-level 592 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 593 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 594 ! ! Griffies parameters (all need "key_ldfslp") 595 ln_traldf_grif = .false. ! use griffies triads 596 ln_traldf_gdia = .false. ! output griffies eddy velocities 597 ln_triad_iso = .false. ! pure lateral mixing in ML 598 ln_botmix_grif = .false. ! lateral mixing on bottom 599 ! ! Coefficients 600 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 601 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 602 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 603 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 604 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 605 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 606 ! (normally=0; not used with Griffies) 607 / 608 !----------------------------------------------------------------------- 609 &namtra_dmp ! tracer: T & S newtonian damping 610 !----------------------------------------------------------------------- 611 ln_tradmp = .true. ! add a damping termn (T) or not (F) 438 612 nn_hdmp = -1 ! horizontal shape =-1, damping in Med and Red Seas only 439 613 ! =XX, damping poleward of XX degrees (XX>0) 440 614 ! + F(distance-to-coast) + Red and Med Seas 441 nn_zdmp = 1! vertical shape =0 damping throughout the water column615 nn_zdmp = 0 ! vertical shape =0 damping throughout the water column 442 616 ! =1 no damping in the mixing layer (kz criteria) 443 617 ! =2 no damping in the mixed layer (rho crieria) … … 445 619 rn_bot = 360. ! bottom time scale of damping [days] 446 620 rn_dep = 800. ! depth of transition between rn_surf and rn_bot [meters] 447 nn_file = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) 448 / 621 nn_file = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) 622 / 623 449 624 !!====================================================================== 450 625 !! *** Dynamics namelists *** … … 456 631 !! namdyn_ldf lateral diffusion scheme 457 632 !!====================================================================== 458 633 ! 459 634 !----------------------------------------------------------------------- 460 635 &namdyn_adv ! formulation of the momentum advection 461 636 !----------------------------------------------------------------------- 462 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) 637 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) 463 638 ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme 464 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 465 / 639 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 640 / 466 641 !----------------------------------------------------------------------- 467 642 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 468 643 !----------------------------------------------------------------------- 469 ln_dynvor_ene = .false. ! enstrophy conserving scheme 470 ln_dynvor_ens = .false. ! energy conserving scheme 471 ln_dynvor_mix = .false. ! mixed scheme 472 ln_dynvor_een = .true. ! energy & enstrophy scheme 644 ln_dynvor_ene = .false. ! enstrophy conserving scheme 645 ln_dynvor_ens = .false. ! energy conserving scheme 646 ln_dynvor_mix = .false. ! mixed scheme 647 ln_dynvor_een = .true. ! energy & enstrophy scheme 473 648 / 474 649 !----------------------------------------------------------------------- 475 650 &namdyn_hpg ! Hydrostatic pressure gradient option 476 651 !----------------------------------------------------------------------- 477 ln_hpg_zco = .false. ! z-coordinate - full steps 652 ln_hpg_zco = .false. ! z-coordinate - full steps 478 653 ln_hpg_zps = .true. ! z-coordinate - partial steps (interpolation) 479 654 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) 480 ln_hpg_hel = .false. ! s-coordinate (helsinki modification)481 ln_hpg_wdj = .false. ! s-coordinate (weighted density jacobian)482 655 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 483 ln_hpg_rot = .false. ! s-coordinate (ROTated axes scheme) 484 rn_gamma = 0.e0 ! weighting coefficient (wdj scheme) 656 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) 485 657 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 486 658 ! centered time scheme (F) 487 nn_dynhpg_rst = 0 ! =1 dynhpg restartable restart or not (=0)488 659 / 489 660 !----------------------------------------------------------------------- … … 497 668 &namdyn_ldf ! lateral diffusion on momentum 498 669 !----------------------------------------------------------------------- 499 ! Type of the operator :500 ln_dynldf_lap = .true. ! laplacian operator 501 ln_dynldf_bilap = .false. ! bilaplacian operator 502 ! Direction of action :503 ln_dynldf_level = .false. ! iso-level 670 ! ! Type of the operator : 671 ln_dynldf_lap = .true. ! laplacian operator 672 ln_dynldf_bilap = .false. ! bilaplacian operator 673 ! ! Direction of action : 674 ln_dynldf_level = .false. ! iso-level 504 675 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 505 676 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") 506 ! Coefficient 507 rn_ahm_0 = 40000. ! horizontal eddy viscosity [m2/s] 508 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 509 / 677 ! ! Coefficient 678 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 679 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 680 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 681 / 682 510 683 !!====================================================================== 511 684 !! Tracers & Dynamics vertical physics namelists 512 685 !!====================================================================== 513 !! 514 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric")515 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke")516 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp")517 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm")518 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx")519 !!====================================================================== 520 686 !! namzdf vertical physics 687 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 688 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 689 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp") 690 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 691 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") 692 !!====================================================================== 693 ! 521 694 !----------------------------------------------------------------------- 522 695 &namzdf ! vertical physics … … 529 702 nn_evdm = 0 ! evd apply on tracer (=0) or on tracer and momentum (=1) 530 703 rn_avevd = 100. ! evd mixing coefficient [m2/s] 531 ln_zdfnpc = .false. ! Non-Penetrative algorithm (T) or not (F)704 ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm (T) or not (F) 532 705 nn_npc = 1 ! frequency of application of npc 533 706 nn_npcp = 365 ! npc control print frequency … … 541 714 rn_alp = 5. ! coefficient of the parameterization 542 715 nn_ric = 2 ! coefficient of the parameterization 716 rn_ekmfc = 0.7 ! Factor in the Ekman depth Equation 717 rn_mldmin = 1.0 ! minimum allowable mixed-layer depth estimate (m) 718 rn_mldmax =1000.0 ! maximum allowable mixed-layer depth estimate (m) 719 rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer 720 rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer 721 ln_mldw = .true. ! Flag to use or not the mized layer depth param. 543 722 / 544 723 !----------------------------------------------------------------------- … … 547 726 rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 548 727 rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation 549 rn_ebb = 6 0. ! coef. of the surface input of tke728 rn_ebb = 67.83 ! coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T) 550 729 rn_emin = 1.e-6 ! minimum value of tke [m2/s2] 551 730 rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2] 552 rn_bshear = 1.e-20 ! background shear (>0)553 731 nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom 554 732 ! = 1 bounded by the local vertical scale factor 555 733 ! = 2 first vertical derivative of mixing length bounded by 1 556 ! = 3 same criteria as case 2 but applied in a different way734 ! = 3 as =2 with distinct disspipative an mixing length scale 557 735 nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 558 ln_mxl0 = .true. ! mixing length scale surface value as function of wind stress (T) or not (F) 559 rn_lmin = 0.001 ! interior buoyancy lenght scale minimum value 560 rn_lmin0 = 0.01 ! surface buoyancy lenght scale minimum value 561 nn_etau = 0 ! exponentially deceasing penetration of tke due to internal & intertial waves 562 ! = 0 no penetration ( O(2 km) resolution) 563 ! = 1 additional tke source (rn_efr * en) 564 ! = 2 additional tke source (rn_efr * en) applied only at the base of the mixed layer 565 ! = 3 additional tke source (HF contribution: mean of stress module - module of mean stress) 566 nn_htau = 1 ! type of exponential decrease of tke penetration 736 ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) 737 rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value 738 ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002) 739 rn_lc = 0.15 ! coef. associated to Langmuir cells 740 nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to internal & intertial waves 741 ! = 0 no penetration 742 ! = 1 add a tke source below the ML 743 ! = 2 add a tke source just at the base of the ML 744 ! = 3 as = 1 applied on HF part of the stress ("key_coupled") 745 rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 746 nn_htau = 1 ! type of exponential decrease of tke penetration below the ML 567 747 ! = 0 constant 10 m length scale 568 ! = 1 0.5m at the equator to 30m at high latitudes 569 ! = 2 30 meters constant depth penetration 570 ! otion used only if nn_etau = 1 or 2: 571 rn_efr = 0.05 ! fraction of surface tke value which penetrates inside the ocean 572 ! otion used only if nn_etau = 3: 573 rn_addhft = -1.e-3 ! add offset applied to the "mean of stress module - module of mean stress" (always kept > 0) 574 rn_sclhft = 1. ! scale factor applied to the "mean of stress module - module of mean stress" 575 ln_lc = .true. ! Langmuir cell parameterisation 576 rn_lc = 0.15 ! coef. associated to Langmuir cells 748 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 577 749 / 578 750 !------------------------------------------------------------------------ 579 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and option nally:751 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally: 580 752 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 581 ln_kpprimix = .true. ! shear instability mixing 753 ln_kpprimix = .true. ! shear instability mixing 582 754 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s] 583 755 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s] 584 756 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability 585 757 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s] 586 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 587 rn_difcon = 1. ! maximum mixing in interior convection [m2/s] 758 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 759 rn_difcon = 1. ! maximum mixing in interior convection [m2/s] 588 760 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv 589 761 nn_ave = 1 ! constant (=0) or profile (=1) background on avt 762 / 763 !----------------------------------------------------------------------- 764 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") 765 !----------------------------------------------------------------------- 766 rn_emin = 1.e-6 ! minimum value of e [m2/s2] 767 rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3] 768 ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988) 769 rn_clim_galp = 0.53 ! galperin limit 770 ln_crban = .true. ! Use Craig & Banner (1994) surface wave mixing parametrisation 771 ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case 772 rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux 773 rn_charn = 70000. ! Charnock constant for wb induced roughness length 774 nn_tkebc_surf = 1 ! surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 775 nn_tkebc_bot = 1 ! bottom tke condition (0/1=Dir/Neum) 776 nn_psibc_surf = 1 ! surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 777 nn_psibc_bot = 1 ! bottom psi condition (0/1=Dir/Neum) 778 nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 779 nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 590 780 / 591 781 !----------------------------------------------------------------------- … … 601 791 rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1) 602 792 rn_tfe = 0.333 ! tidal dissipation efficiency 603 rn_me = 0.2 ! mixing efficiency 604 ln_tmx_itf = . TRUE. ! ITF specific parameterisation793 rn_me = 0.2 ! mixing efficiency 794 ln_tmx_itf = .true. ! ITF specific parameterisation 605 795 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency 606 796 / 607 !!====================================================================== 608 !! *** Miscelaneous namelists *** 797 798 !!====================================================================== 799 !! *** Miscellaneous namelists *** 609 800 !!====================================================================== 610 801 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 611 !! nammpp_dyndist Massively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist")612 802 !! namctl Control prints & Benchmark 613 !! namsol elliptic solver / island / free surface 614 !!====================================================================== 615 616 !----------------------------------------------------------------------- 617 &namsol ! elliptic solver / island / free surface 803 !! namsol elliptic solver / island / free surface 804 !!====================================================================== 805 ! 806 !----------------------------------------------------------------------- 807 &namsol ! elliptic solver / island / free surface 618 808 !----------------------------------------------------------------------- 619 809 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) … … 633 823 ! buffer blocking send or immediate non-blocking sends, resp. 634 824 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 635 / 636 !----------------------------------------------------------------------- 637 &nammpp_dyndist ! Massively Parallel Distribution for AGRIF zoom ("key_agrif" && "key_mpp_dyndist") 638 !----------------------------------------------------------------------- 639 jpni = 1 ! jpni number of processors following i 640 jpnj = 1 ! jpnj number of processors following j 641 jpnij = 1 ! jpnij number of local domains 825 ln_nnogather= .false. ! activate code to avoid mpi_allgather use at the northfold 826 jpni = 0 ! jpni number of processors following i (set automatically if < 1) 827 jpnj = 0 ! jpnj number of processors following j (set automatically if < 1) 828 jpnij = 0 ! jpnij number of local domains (set automatically if < 1) 642 829 / 643 830 !----------------------------------------------------------------------- … … 654 841 nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 655 842 ! (no physical validity of the results) 656 nn_ bit_cmp = 0 ! bit comparison mode (1/0): CAUTION use zero except for test657 ! of comparison between single and multiple processor runs 658 / 659 660 !! ======================================================================661 !! *** Diagnostics namelists ***662 !! ======================================================================843 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 844 / 845 846 !!====================================================================== 847 !! *** Diagnostics namelists *** 848 !!====================================================================== 849 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") 663 850 !! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld") 664 !! namgap level mean model-data gap ("key_diagap")665 851 !! namflo float parameters ("key_float") 666 852 !! namptr Poleward Transport Diagnostics 667 !!====================================================================== 668 853 !! namhsb Heat and salt budgets 854 !!====================================================================== 855 ! 856 !----------------------------------------------------------------------- 857 &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") 858 !----------------------------------------------------------------------- 859 nn_nchunks_i= 4 ! number of chunks in i-dimension 860 nn_nchunks_j= 4 ! number of chunks in j-dimension 861 nn_nchunks_k= 31 ! number of chunks in k-dimension 862 ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 863 ! is optimal for postprocessing which works exclusively with horizontal slabs 864 ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression 865 ! (F) ignore chunking information and produce netcdf3-compatible files 866 / 669 867 !----------------------------------------------------------------------- 670 868 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 671 ! ! or mixed-layer trends or barotropic vorticity ( 'key_trdmld' or"key_trdvor")869 ! ! or mixed-layer trends or barotropic vorticity ("key_trdmld" or "key_trdvor") 672 870 !----------------------------------------------------------------------- 673 871 nn_trd = 365 ! time step frequency dynamics and tracers trends … … 680 878 / 681 879 !----------------------------------------------------------------------- 682 &namgap ! level mean model-data gap ('key_diagap')683 !-----------------------------------------------------------------------684 nn_gap = 15 ! time-step frequency of model-data gap computation685 nn_prg = 10 ! AUTO - time-step frequency of gap print in model output686 /687 !-----------------------------------------------------------------------688 880 &namflo ! float parameters ("key_float") 689 881 !----------------------------------------------------------------------- 690 ln_rstflo = .false. ! float restart (T) or not (F) 691 nn_writefl= 75 ! frequency of writing in float output file 692 nn_stockfl= 5475 ! frequency of creation of the float restart file 693 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 694 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 695 ! or computed with Blanke' scheme (F) 882 jpnfl = 1 ! total number of floats during the run 883 jpnnewflo = 0 ! number of floats for the restart 884 ln_rstflo = .false. ! float restart (T) or not (F) 885 nn_writefl = 75 ! frequency of writing in float output file 886 nn_stockfl = 5475 ! frequency of creation of the float restart file 887 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 888 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 889 ! or computed with Blanke' scheme (F) 890 ln_ariane = .true. ! Input with Ariane tool convention(T) 891 ln_flo_ascii = .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) 696 892 / 697 893 !----------------------------------------------------------------------- 698 894 &namptr ! Poleward Transport Diagnostic 699 895 !----------------------------------------------------------------------- 700 ln_diaptr = .true. 896 ln_diaptr = .true. ! Poleward heat and salt transport (T) or not (F) 701 897 ln_diaznl = .true. ! Add zonal means and meridional stream functions 702 ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not 898 ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not 703 899 ! (orca configuration only, need input basins mask file named "subbasins.nc" 704 900 ln_ptrcomp = .true. ! Add decomposition : overturning 705 nf_ptr = 1 ! Frequency of ptr computation [time step] 706 nf_ptr_wri = 15 ! AUTO - Frequency of ptr outputs 707 / 901 nn_fptr = 1 ! Frequency of ptr computation [time step] 902 nn_fwri = 15 ! Frequency of ptr outputs [time step] 903 / 904 !----------------------------------------------------------------------- 905 &namhsb ! Heat and salt budgets 906 !----------------------------------------------------------------------- 907 ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) 908 / 909 !----------------------------------------------------------------------- 910 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 911 !----------------------------------------------------------------------- 912 nit000_han = 1 ! First time step used for harmonic analysis 913 nitend_han = 75 ! Last time step used for harmonic analysis 914 nstep_han = 15 ! Time step frequency for harmonic analysis 915 tname(1) = 'M2' ! Name of tidal constituents 916 tname(2) = 'K1' 917 / 918 !----------------------------------------------------------------------- 919 &namdct ! transports through sections 920 !----------------------------------------------------------------------- 921 nn_dct = 15 ! time step frequency for transports computing 922 nn_dctwri = 15 ! time step frequency for transports writing 923 nn_secdebug = 112 ! 0 : no section to debug 924 ! -1 : debug all section 925 ! 0 < n : debug section number n 926 / 927 928 !!====================================================================== 929 !! *** Observation & Assimilation namelists *** 930 !!====================================================================== 931 !! namobs observation and model comparison ('key_diaobs') 932 !! nam_asminc assimilation increments ('key_asminc') 933 !!====================================================================== 934 ! 935 !----------------------------------------------------------------------- 936 &namobs ! observation usage switch ('key_diaobs') 937 !----------------------------------------------------------------------- 938 ln_t3d = .false. ! Logical switch for T profile observations 939 ln_s3d = .false. ! Logical switch for S profile observations 940 ln_ena = .false. ! Logical switch for ENACT insitu data set 941 ! ! ln_cor Logical switch for Coriolis insitu data set 942 ln_profb = .false. ! Logical switch for feedback insitu data set 943 ln_sla = .false. ! Logical switch for SLA observations 944 945 ln_sladt = .false. ! Logical switch for AVISO SLA data 946 947 ln_slafb = .false. ! Logical switch for feedback SLA data 948 ! ln_ssh Logical switch for SSH observations 949 950 ln_sst = .true. ! Logical switch for SST observations 951 ln_reysst = .true. ! ln_reysst Logical switch for Reynolds observations 952 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 953 954 ln_sstfb = .false. ! Logical switch for feedback SST data 955 ! ln_sss Logical switch for SSS observations 956 ! ln_seaice Logical switch for Sea Ice observations 957 ! ln_vel3d Logical switch for velocity observations 958 ! ln_velavcur Logical switch for velocity daily av. cur. 959 ! ln_velhrcur Logical switch for velocity high freq. cur. 960 ! ln_velavadcp Logical switch for velocity daily av. ADCP 961 ! ln_velhradcp Logical switch for velocity high freq. ADCP 962 ! ln_velfb Logical switch for feedback velocity data 963 ! ln_grid_global Global distribtion of observations 964 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 965 ! grid_search_file Grid search lookup file header 966 ! enactfiles ENACT input observation file names 967 ! coriofiles Coriolis input observation file name 968 ! ! profbfiles: Profile feedback input observation file name 969 profbfiles = 'profiles_01.nc' 970 ! ln_profb_enatim Enact feedback input time setting switch 971 ! slafilesact Active SLA input observation file name 972 ! slafilespas Passive SLA input observation file name 973 ! ! slafbfiles: Feedback SLA input observation file name 974 slafbfiles = 'sla_01.nc' 975 ! sstfiles GHRSST input observation file name 976 ! ! sstfbfiles: Feedback SST input observation file name 977 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 978 ! seaicefiles Sea Ice input observation file name 979 ! velavcurfiles Vel. cur. daily av. input file name 980 ! velhvcurfiles Vel. cur. high freq. input file name 981 ! velavadcpfiles Vel. ADCP daily av. input file name 982 ! velhvadcpfiles Vel. ADCP high freq. input file name 983 ! velfbfiles Vel. feedback input observation file name 984 ! dobsini Initial date in window YYYYMMDD.HHMMSS 985 ! dobsend Final date in window YYYYMMDD.HHMMSS 986 ! n1dint Type of vertical interpolation method 987 ! n2dint Type of horizontal interpolation method 988 ! ln_nea Rejection of observations near land switch 989 nmsshc = 0 ! MSSH correction scheme 990 ! mdtcorr MDT correction 991 ! mdtcutoff MDT cutoff for computed correction 992 ln_altbias = .false. ! Logical switch for alt bias 993 ln_ignmis = .true. ! Logical switch for ignoring missing files 994 ! endailyavtypes ENACT daily average types 995 ln_grid_global = .true. 996 ln_grid_search_lookup = .false. 997 / 998 !----------------------------------------------------------------------- 999 &nam_asminc ! assimilation increments ('key_asminc') 1000 !----------------------------------------------------------------------- 1001 ln_bkgwri = .false. ! Logical switch for writing out background state 1002 ln_trainc = .false. ! Logical switch for applying tracer increments 1003 ln_dyninc = .false. ! Logical switch for applying velocity increments 1004 ln_sshinc = .false. ! Logical switch for applying SSH increments 1005 ln_asmdin = .false. ! Logical switch for Direct Initialization (DI) 1006 ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU) 1007 nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1] 1008 nitdin = 0 ! Timestep of background for DI in [0,nitend-nit000-1] 1009 nitiaustr = 1 ! Timestep of start of IAU interval in [0,nitend-nit000-1] 1010 nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1] 1011 niaufn = 0 ! Type of IAU weighting function 1012 ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin 1013 salfixmin = -9999 ! Minimum salinity after applying the increments 1014 nn_divdmp = 0 ! Number of iterations of divergence damping operator 1015 / 1016 !----------------------------------------------------------------------- 1017 &namsbc_wave ! External fields from wave model 1018 !----------------------------------------------------------------------- 1019 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 1020 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 1021 sn_cdg = 'cdg_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' ,'' , '' 1022 sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' ,'' , '' 1023 sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' ,'' , '' 1024 sn_wn = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' ,'' , '' 1025 ! 1026 cn_dir_cdg = './' ! root directory for the location of drag coefficient files 1027 / 1028 !----------------------------------------------------------------------- 1029 &namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed) 1030 !----------------------------------------------------------------------- 1031 ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 1032 ln_neptsimp = .false. ! yes/no use simplified neptune 1033 1034 ln_smooth_neptvel = .false. ! yes/no smooth zunep, zvnep 1035 rn_tslse = 1.2e4 ! value of lengthscale L at the equator 1036 rn_tslsp = 3.0e3 ! value of lengthscale L at the pole 1037 ! Specify whether to ramp down the Neptune velocity in shallow 1038 ! water, and if so the depth range controlling such ramping down 1039 ln_neptramp = .true. ! ramp down Neptune velocity in shallow water 1040 rn_htrmin = 100.0 ! min. depth of transition range 1041 rn_htrmax = 200.0 ! max. depth of transition range 1042 /
Note: See TracChangeset
for help on using the changeset viewer.