Changeset 1602 for trunk/CONFIG/ORCA2_LIM/EXP00
- Timestamp:
- 2009-08-11T13:09:18+02:00 (15 years ago)
- Location:
- trunk/CONFIG/ORCA2_LIM/EXP00
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist
r1556 r1602 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 2 !! NEMO/OPA : 1 - run manager (namrun) 3 !! namelists 2 - Domain (nam _zgr, nam_zgr_sco, namdom)3 !! namelists 2 - Domain (namzgr, namzgr_sco, namdom) 4 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, nam alb)6 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nam tide)5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namsbc_alb) 6 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 7 7 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 8 !! 6 - Tracer (nameos, nam _traadv, nam_traldf, namtdp)9 !! 7 - dynamics (nam _dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf)10 !! 8 - Verical physics (nam _zdf, nam_npc, nam_ric, nam_tke, nam_kpp, nam_ddm, nam_tmx)11 !! 9 - diagnostics (namtrd, namgap, nam flo, namptr)12 !! 9 - miscellaneous (namsol, nam _mpp, nam_mpp_dyndist, namctl)8 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_dmp) 9 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 10 !! 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) 13 13 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 14 ! CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE. … … 23 23 &namrun ! parameters of the run 24 24 !----------------------------------------------------------------------- 25 no = 0 ! job number 26 cexper = "Agulhas" ! experience name 25 nn_no = 0 ! job number 26 cn_exp = "Agulhas"! experience name 27 nn_it000 = 1 ! first time step 28 nn_itend = 10950 ! last time step (std 5475) 29 nn_date0 = 010101 ! initial calendar date yymmdd (used if nrstdt=1) 30 nn_leapy = 0 ! Leap year calendar (1) or not (0) 31 nn_istate = 0 ! output the initial state (1) or not (0) 32 nn_stock = 10950 ! frequency of creation of a restart file (modulo referenced to 1) 33 nn_write = 10950 ! frequency of write in the output file (modulo referenced to nit000) 34 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) 35 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) 36 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. ! start from rest (F) or from a restart file (T) 39 nn_rstctl = 0 ! restart control = 0 nit000 is not compared to the restart file value 40 ! = 1 use ndate0 in namelist (not the value in the restart file) 41 ! = 2 calendar parameters read in the restart file 27 42 cn_ocerst_in = "restart" ! suffix of ocean restart name (input) 28 43 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 29 ln_rstart = .false. ! start from rest (F) or from a restart file (T)30 nrstdt = 0 ! restart control = 0 nit000 is not compared to the restart file value31 ! = 1 use ndate0 in namelist (not the value in the restart file)32 ! = 2 calendar parameters read in the restart file33 nit000 = 1 ! first time step34 nitend = 10950 ! last time step35 ndate0 = 010101 ! initial calendar date yymmdd (used if nrstdt=1)36 nleapy = 0 ! Leap year calendar (1) or not (0)37 ninist = 0 ! output the initial state (1) or not (0)38 nstock = 10950 ! frequency of creation of a restart file (modulo referenced to 1)39 nwrite = 10950 ! frequency of write in the output file (modulo referenced to nit000)40 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T)41 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%)42 ln_clobber = .false. ! clobber (overwrite) an existing file43 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (working only with iom_nf90 routines)44 44 / 45 45 !!====================================================================== 46 46 !! *** Domain namelists *** 47 47 !!====================================================================== 48 !! nam _zgr vertical coordinate49 !! nam _zgr_sco s-coordinate or hybrid z-s-coordinate50 !! namdom 51 !!====================================================================== 52 53 !----------------------------------------------------------------------- 54 &nam _zgr! vertical coordinate48 !! namzgr vertical coordinate 49 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 50 !! namdom space and time domain (bathymetry, mesh, timestep) 51 !!====================================================================== 52 53 !----------------------------------------------------------------------- 54 &namzgr ! vertical coordinate 55 55 !----------------------------------------------------------------------- 56 56 ln_zco = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) … … 59 59 / 60 60 !----------------------------------------------------------------------- 61 &nam _zgr_sco! s-coordinate or hybrid z-s-coordinate62 !----------------------------------------------------------------------- 63 sbot_min= 300. ! minimum depth of s-bottom surface (>0) (m)64 sbot_max= 5250. ! maximum depth of s-bottom surface (= ocean depth) (>0) (m)65 theta= 6.0 ! surface control parameter (0<=theta<=20)66 thetb= 0.75 ! bottom control parameter (0<=thetb<= 1)67 r _max= 0.15 ! maximum cut-off r-value allowed (0<r_max<1)61 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate 62 !----------------------------------------------------------------------- 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) 65 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 68 ln_s_sigma = .false. ! hybrid s-sigma coordinates 69 bb= 0.8 ! stretching with s-sigma70 hc = 150.0 ! critical depth with s-sigma69 rn_bb = 0.8 ! stretching with s-sigma 70 rn_hc = 150.0 ! critical depth with s-sigma 71 71 / 72 72 !----------------------------------------------------------------------- 73 73 &namdom ! space and time domain (bathymetry, mesh, timestep) 74 74 !----------------------------------------------------------------------- 75 ntopo = 1 ! compute (=0) or read(=1) the bathymetry file 76 e3zps_min = 5. ! the thickness of the partial step is set larger than the minimum 77 e3zps_rat = 0.1 ! of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 78 nmsh = 0 ! create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 79 nacc = 0 ! =1 acceleration of convergence method used, rdt < rdttra(k) 80 ! =0, no acceleration, rdt = rdttra 81 atfp = 0.1 ! asselin time filter parameter 82 rdt = 2880. ! time step for the dynamics (and tracer if nacc=0) 83 rdtmin = 2880. ! minimum time step on tracers (used if nacc=1) 84 rdtmax = 2880. ! maximum time step on tracers (used if nacc=1) 85 rdth = 800. ! depth variation of tracer time step (used if nacc=1) 75 nn_bathy = 1 ! compute (=0) or read(=1) the bathymetry file 76 nn_closea = 0 ! closed seas and lakes are removed (=0) or kept (=1) from the ORCA domain 77 nn_msh = 0 ! create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 78 rn_e3zps_min= 5. ! 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) 80 ! 81 rn_rdt = 2880. ! time step for the dynamics (and tracer if nacc=0) ==> 5760 86 82 nn_baro = 64 ! number of barotropic time step (for the split explicit algorithm) ("key_dynspg_ts") 87 nclosea = 0 ! = 0 no closed sea in the model domain 88 ! = 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 83 rn_atfp = 0.1 ! asselin time filter parameter 84 nn_acc = 0 ! acceleration of convergence : =1 used, rdt < rdttra(k) 85 ! =0, not used, rdt = rdttra 86 rn_rdtmin = 2880. ! minimum time step on tracers (used if nacc=1) 87 rn_rdtmax = 2880. ! maximum time step on tracers (used if nacc=1) 88 rn_rdth = 800. ! depth variation of tracer time step (used if nacc=1) 89 89 / 90 90 !!====================================================================== … … 97 97 !! namsbc_core CORE bulk formulea formulation 98 98 !! namsbc_cpl CouPLed formulation ("key_coupled") 99 !! nam qsrpenetrative solar radiation99 !! namtra_qsr penetrative solar radiation 100 100 !! namsbc_rnf river runoffs 101 101 !! namsbc_ssr sea surface restoring term (for T and/or S) 102 !! nam albalbedo parameters102 !! namsbc_alb albedo parameters 103 103 !!====================================================================== 104 104 … … 122 122 ln_rnf = .false. ! runoffs (T => fill namsbc_rnf) 123 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 nn_fsbc time step ,124 nn_fwb = 0 ! FreshWater Budget: =0 unchecked 125 ! =1 global mean of e-p-r set to zero at each time step 126 126 ! =2 annual global mean of e-p-r set to zero 127 127 / … … 139 139 &namsbc_flx ! surface boundary condition : flux formulation 140 140 !----------------------------------------------------------------------- 141 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !142 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !143 sn_utau = 'utau' , 24. , 'utau' , .false. , .false. ,'yearly' , '' , ''144 sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false. ,'yearly' , '' , ''145 sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false. ,'yearly' , '' , ''146 sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false. ,'yearly' , '' , ''147 sn_emp = 'emp' , 24. , 'emp' , .false. , .false. ,'yearly' , '' , ''141 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 142 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 143 sn_utau = 'utau' , 24. , 'utau' , .false. , .false. , 'yearly' , '' , '' 144 sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false. , 'yearly' , '' , '' 145 sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false. , 'yearly' , '' , '' 146 sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false. , 'yearly' , '' , '' 147 sn_emp = 'emp' , 24. , 'emp' , .false. , .false. , 'yearly' , '' , '' 148 148 ! 149 149 cn_dir = './' ! root directory for the location of the flux files … … 152 152 &namsbc_clio ! namsbc_clio CLIO bulk formulea 153 153 !----------------------------------------------------------------------- 154 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !155 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !156 sn_utau = 'taux_1m' , -1. , 'sozotaux' , .false. , .true. ,'yearly' , '' , ''157 sn_vtau = 'tauy_1m' , -1. , 'sometauy' , .false. , .true. ,'yearly' , '' , ''158 sn_wndm = 'flx' , -1. , 'socliowi' , .false. , .true. ,'yearly' , '' , ''159 sn_tair = 'flx' , -1. , 'socliot2' , .false. , .true. ,'yearly' , '' , ''160 sn_humi = 'flx' , -1. , 'socliohu' , .false. , .true. ,'yearly' , '' , ''161 sn_ccov = 'flx' , -1. , 'socliocl' , .false. , .true. ,'yearly' , '' , ''162 sn_prec = 'flx' , -1. , 'socliopl' , .false. , .true. ,'yearly' , '' , ''154 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 155 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 156 sn_utau = 'taux_1m' , -1. , 'sozotaux' , .true. , .true. , 'yearly' , '' , '' 157 sn_vtau = 'tauy_1m' , -1. , 'sometauy' , .true. , .true. , 'yearly' , '' , '' 158 sn_wndm = 'flx' , -1. , 'socliowi' , .true. , .true. , 'yearly' , '' , '' 159 sn_tair = 'flx' , -1. , 'socliot2' , .true. , .true. , 'yearly' , '' , '' 160 sn_humi = 'flx' , -1. , 'socliohu' , .true. , .true. , 'yearly' , '' , '' 161 sn_ccov = 'flx' , -1. , 'socliocl' , .false. , .true. , 'yearly' , '' , '' 162 sn_prec = 'flx' , -1. , 'socliopl' , .false. , .true. , 'yearly' , '' , '' 163 163 ! 164 164 cn_dir = './' ! root directory for the location of the bulk files are … … 167 167 &namsbc_core ! namsbc_core CORE bulk formulea 168 168 !----------------------------------------------------------------------- 169 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! weights! rotation !170 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename! pairing !171 sn_wndi = 'u10_core' , -1. , 'u10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' ,'U1'172 sn_wndj = 'v10_core' , -1. , 'v10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' ,'V1'173 sn_qsr = 'qsw_core' , -1. , 'swdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''174 sn_qlw = 'qlw_core' , -1. , 'lwdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''175 sn_tair = 't2_core' , -1. , 't2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''176 sn_humi = 'q2_core' , -1. , 'q2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''177 sn_prec = 'precip_core' , -1. , 'precip' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''178 sn_snow = 'snow_core' , -1. , 'snow' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''169 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 170 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 171 sn_wndi = 'u10_core' , -1. , 'u10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'U1' 172 sn_wndj = 'v10_core' , -1. , 'v10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'V1' 173 sn_qsr = 'qsw_core' , -1. , 'swdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 174 sn_qlw = 'qlw_core' , -1. , 'lwdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 175 sn_tair = 't2_core' , -1. , 't2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 176 sn_humi = 'q2_core' , -1. , 'q2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 177 sn_prec = 'precip_core', -1. , 'precip' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 178 sn_snow = 'snow_core' , -1. , 'snow' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 179 179 ! 180 180 cn_dir = './' ! root directory for the location of the bulk files 181 181 ln_2m = .true. ! air temperature and humidity referenced at 2m (T) instead 10m (F) 182 alpha_precip= 1. ! multiplicative factor for precipitation (total & snow)182 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 183 183 / 184 184 !----------------------------------------------------------------------- … … 187 187 / 188 188 !----------------------------------------------------------------------- 189 &nam qsr! penetrative solar radiation190 !----------------------------------------------------------------------- 191 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !192 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !193 sn_chl = 'chlorophyll' , -1. , 'CHLA' , .true. , .true. ,'yearly' , '' , ''194 ! 189 &namtra_qsr ! penetrative solar radiation 190 !----------------------------------------------------------------------- 191 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 192 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 193 sn_chl = 'chlorophyll', -1. , 'CHLA' , .true. , .true. , 'yearly' , '' , '' 194 195 195 cn_dir = './' ! root directory for the location of the runoff files 196 196 ln_traqsr = .true. ! Light penetration (T) or not (F) … … 207 207 &namsbc_rnf ! runoffs namelist surface boundary condition 208 208 !----------------------------------------------------------------------- 209 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !210 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !211 sn_rnf = 'runoff_1m_nomask' , -1. , 'sorunoff' , .true. , .true. ,'yearly' , '' , ''212 sn_cnf = 'runoff_1m_nomask' , 0. , 'socoefr' , .false. , .true. ,'yearly' , '' , ''213 ! 209 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 210 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 211 sn_rnf = 'runoff_1m_nomask' , -1. , 'sorunoff' , .true. , .true. , 'yearly' , '' , '' 212 sn_cnf = 'runoff_1m_nomask' , 0. , 'socoefr' , .false. , .true. , 'yearly' , '' , '' 213 214 214 cn_dir = './' ! root directory for the location of the runoff files 215 215 ln_rnf_emp = .false. ! runoffs included into precipitation field (T) or into a file (F) … … 217 217 rn_hrnf = 0.e0 ! depth over which enhanced vertical mixing is used 218 218 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] 219 rn_ mul_rnf= 1.e0 ! multiplicative factor for runoff219 rn_rfact = 1.e0 ! multiplicative factor for runoff 220 220 / 221 221 !----------------------------------------------------------------------- 222 222 &namsbc_ssr ! surface boundary condition : sea surface restoring 223 223 !----------------------------------------------------------------------- 224 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !225 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !226 sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false. ,'yearly' , '' , ''227 sn_sss = 'sss_data' , -1. , 'sss' , .true. , .false. ,'yearly' , '' , ''228 ! 224 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 225 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 226 sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false. , 'yearly' , '' , '' 227 sn_sss = 'sss_data' , -1. , 'sss' , .true. , .true. , 'yearly' , '' , '' 228 229 229 cn_dir = './' ! root directory for the location of the runoff files 230 230 nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0) 231 231 nn_sssr = 0 ! add a damping term in the surface freshwater flux (=1) or not (=0) 232 dqdt= -40. ! magnitude of the retroaction on temperature [W/m2/K]233 deds= -27.7 ! magnitude of the damping on salinity [mm/day/psu]232 rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 233 rn_deds = -27.7 ! magnitude of the damping on salinity [mm/day/psu] 234 234 ln_sssr_bnd = .false. ! flag to bound erp term (associated with nn_sssr=2) 235 235 rn_sssr_bnd = 0.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] 236 236 / 237 237 !----------------------------------------------------------------------- 238 &namalb ! albedo parameters 239 !----------------------------------------------------------------------- 240 cgren = 0.06 ! correction of the snow or ice albedo to take into account the 241 albice = 0.53 ! albedo of melting ice in the arctic and antarctic 242 alphd = 0.80 ! coefficients for linear interpolation used to 243 alphc = 0.65 ! compute albedo between two extremes values 244 alphdi = 0.72 ! (Pyane, 1972) 245 / 238 &namsbc_alb ! albedo parameters 239 !----------------------------------------------------------------------- 240 rn_cloud = 0.06 ! cloud correction to snow and ice albedo 241 rn_albice = 0.53 ! albedo of melting ice in the arctic and antarctic 242 rn_alphd = 0.80 ! coefficients for linear interpolation used to 243 rn_alphc = 0.65 ! compute albedo between two extremes values 244 rn_alphdi = 0.72 ! (Pyane, 1972) 245 / 246 246 247 !!====================================================================== 247 248 !! *** Lateral boundary condition *** … … 258 259 &namlbc ! lateral momentum boundary condition 259 260 !----------------------------------------------------------------------- 260 shlat = 2. ! shlat = 0 : free slip 261 ! 0 < shlat < 2 : partial slip 262 ! shlat = 2 : no slip 263 ! 2 < shlat : strong slip 261 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 262 ! free slip ! partial slip ! no slip ! strong slip 264 263 / 265 264 !----------------------------------------------------------------------- 266 265 &namcla ! cross land advection 267 266 !----------------------------------------------------------------------- 268 n _cla= 0 ! advection between 2 ocean pts separates by land267 nn_cla = 0 ! advection between 2 ocean pts separates by land 269 268 / 270 269 !----------------------------------------------------------------------- 271 270 &namobc ! open boundaries parameters ("key_obc") 272 271 !----------------------------------------------------------------------- 273 nobc_dta = 0 ! = 0 the obc data are equal to the initial state 272 ln_obc_clim= .false. ! climatological obc data files (T) or not (F) 273 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) 274 ln_obc_fla = .false. ! Flather open boundary condition 275 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state 274 276 ! = 1 the obc data are read in 'obc.dta' files 275 rdpein = 1. ! ??? 276 rdpwin = 1. ! ??? 277 rdpnin = 30. ! ??? 278 rdpsin = 1. ! ??? 279 rdpeob = 1500. ! time relaxation (days) for the east open boundary 280 rdpwob = 15. ! " " " west " 281 rdpnob = 150. ! " " " north " 282 rdpsob = 15. ! " " " south " 283 zbsic1 = 140.e+6 ! barotropic stream function on first isolated coastline 284 zbsic2 = 1.e+6 ! " " second " 285 zbsic3 = 0. ! " " thrid " 286 ln_obc_clim= .true. ! climatological obc data files (T) or not (F) 287 ln_vol_cst = .false. ! impose the total volume conservation (T) or not (F) 288 / 289 !----------------------------------------------------------------------- 290 &namagrif ! ("key_agrif") 291 !----------------------------------------------------------------------- 292 nbclineupdate = 3 ! baroclinic update frequency 277 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data 278 ! set to monthly if obc datafile hold 1 month of data 279 rn_dpein = 1. ! damping time scale for inflow at east open boundary 280 rn_dpwin = 1. ! - - - west - - 281 rn_dpnin = 1. ! - - - north - - 282 rn_dpsin = 1. ! - - - south - - 283 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary 284 rn_dpwob = 15. ! - - - west - - 285 rn_dpnob = 3000. ! - - - north - - 286 rn_dpsob = 15. ! - - - south - - 287 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R) 288 ! = 1 the total volume remains constant 289 / 290 !----------------------------------------------------------------------- 291 &namagrif ! AGRIF zoom ("key_agrif") 292 !----------------------------------------------------------------------- 293 nn_cln_update = 3 ! baroclinic update frequency 293 294 ln_spc_dyn = .true. ! use 0 as special value for dynamics 294 visc_tra = 2880. ! viscosity coeeficient for tracers sponge layer295 visc_dyn = 2880. ! viscosity coeeficient for dynamics sponge layer296 / 297 !----------------------------------------------------------------------- 298 &nambdy ! unstructured open boundaries parameters("key_bdy")295 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [s] 296 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [s] 297 / 298 !----------------------------------------------------------------------- 299 &nambdy ! unstructured open boundaries ("key_bdy") 299 300 !----------------------------------------------------------------------- 300 301 filbdy_mask = '' ! name of mask file (if ln_bdy_mask=.TRUE.) … … 316 317 / 317 318 !----------------------------------------------------------------------- 318 &nam tide! tidal forcing at unstructured boundaries319 &nambdy_tide ! tidal forcing at unstructured boundaries 319 320 !----------------------------------------------------------------------- 320 321 filtide = 'bdytide_' ! file name root of tidal forcing files … … 323 324 ln_tide_date = .false. ! adjust tidal harmonics for start date of run 324 325 / 326 325 327 !!====================================================================== 326 328 !! *** Bottom boundary condition *** … … 334 336 &nambfr ! bottom friction 335 337 !----------------------------------------------------------------------- 336 n botfr = 1 ! type of bottom friction : = 0 : no slip, = 2 : nonlinear friction338 nn_bfr = 1 ! type of bottom friction : = 0 : no slip, = 2 : nonlinear friction 337 339 ! = 3 : free slip, = 1 : linear friction 338 bfri1= 4.e-4 ! bottom drag coefficient (linear case)339 bfri2= 1.e-3 ! bottom drag coefficient (non linear case)340 bfeb2= 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2)340 rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 341 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) 342 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2) 341 343 / 342 344 !----------------------------------------------------------------------- 343 345 &nambbc ! bottom temperature boundary condition 344 346 !----------------------------------------------------------------------- 345 n geo_flux = 2 ! geothermal heat flux = 0 no flux considered346 ! 347 ! 348 ngeo_flux_const = 86.4e-3! Constant value of geothermal heat flux [W/m2]347 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 348 ! = 1 constant flux 349 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 350 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 349 351 / 350 352 !----------------------------------------------------------------------- … … 353 355 ! ! diffusive bbl ("key_trabbl") 354 356 ! ! advective bbl ("key_trabbl_adv") 355 atrbbl= 10000. ! lateral mixing coefficient in the bbl [m2/s]357 rn_ahtbbl = 10000. ! lateral mixing coefficient in the bbl [m2/s] 356 358 / 357 359 !!====================================================================== … … 359 361 !!====================================================================== 360 362 !! nameos equation of state 361 !! nam _traadv advection scheme362 !! nam _traldf lateral diffusion scheme363 !! namt dp tracer newtonian damping("key_tradmp")363 !! namtra_adv advection scheme 364 !! namtra_ldf lateral diffusion scheme 365 !! namtra_dmp T & S newtonian damping ("key_tradmp") 364 366 !!====================================================================== 365 367 … … 367 369 &nameos ! ocean physical parameters 368 370 !----------------------------------------------------------------------- 369 n eos= 0 ! type of equation of state and Brunt-Vaisala frequency371 nn_eos = 0 ! type of equation of state and Brunt-Vaisala frequency 370 372 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 371 373 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 372 374 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 373 r alpha= 2.e-4 ! thermal expension coefficient (neos= 1 or 2)374 r beta= 0.001 ! saline expension coefficient (neos= 2)375 / 376 !----------------------------------------------------------------------- 377 &nam _traadv ! advection scheme for tracer375 rn_alpha = 2.e-4 ! thermal expension coefficient (neos= 1 or 2) 376 rn_beta = 0.001 ! saline expension coefficient (neos= 2) 377 / 378 !----------------------------------------------------------------------- 379 &namtra_adv ! advection scheme for tracer 378 380 !----------------------------------------------------------------------- 379 381 ln_traadv_cen2 = .false. ! 2nd order centered scheme … … 384 386 / 385 387 !----------------------------------------------------------------------- 386 &nam _traldf ! lateral diffusion scheme for tracer387 !----------------------------------------------------------------------- 388 !! Type of the operator :388 &namtra_ldf ! lateral diffusion scheme for tracer 389 !----------------------------------------------------------------------- 390 ! Type of the operator : 389 391 ln_traldf_lap = .true. ! laplacian operator 390 392 ln_traldf_bilap = .false. ! bilaplacian operator 391 393 ! Direction of action : 392 394 ln_traldf_level = .false. ! iso-level 393 395 ln_traldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 394 396 ln_traldf_iso = .true. ! iso-neutral (require "key_ldfslp") 395 ! ! Coefficient 396 aht0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] 397 ahtb0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 398 aeiv0 = 0. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 399 / 400 !----------------------------------------------------------------------- 401 &namtdp ! tracer newtonian damping ('key_tradmp') 402 !----------------------------------------------------------------------- 403 ndmp = 90 ! type of damping in temperature and salinity 404 ! ='latitude', damping poleward of 'ndmp' degrees and function 405 ! of the distance-to-coast. Red and Med Seas as ndmp=-1 406 ! =-1 damping only in Med and Red Seas 407 ndmpf = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) 408 nmldmp = 1 ! type of damping: =0 damping throughout the water column 409 ! =1 no damping in the mixed layer defined by avt >5cm2/s ) 410 ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 411 sdmp = 50. ! surface time scale for internal damping (days) 412 bdmp = 360. ! bottom time scale for internal damping (days) 413 hdmp = 800. ! depth of transition between sdmp and bdmp (meters) 397 ! Coefficient 398 rn_aht_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] 399 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 400 rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 401 / 402 !----------------------------------------------------------------------- 403 &namtra_dmp ! tracer: T & S newtonian damping ('key_tradmp') 404 !----------------------------------------------------------------------- 405 nn_hdmp = 90 ! horizontal shape =-1, damping in Med and Red Seas only 406 ! =XX, damping poleward of XX degrees (XX>0) 407 ! + F(distance-to-coast) + Red and Med Seas 408 nn_zdmp = 1 ! vertical shape =0 damping throughout the water column 409 ! =1 no damping in the mixing layer (kz criteria) 410 ! =2 no damping in the mixed layer (rho crieria) 411 rn_surf = 50. ! surface time scale of damping [days] 412 rn_bot = 360. ! bottom time scale of damping [days] 413 rn_dep = 800. ! depth of transition between rn_surf and rn_bot [meters] 414 nn_file = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) 414 415 / 415 416 !!====================================================================== 416 417 !! *** Dynamics namelists *** 417 418 !!====================================================================== 418 !! nam_dynadv formulation of the momentum advection 419 !! nam_dynvor advection scheme 420 !! nam_dynhpg hydrostatic pressure gradient 421 !! namflg hydrostatic pressure gradient time stepping 422 !! nam_dynspg surface pressure gradient (CPP key only) 423 !! nam_dynldf lateral diffusion scheme 424 !!====================================================================== 425 426 !----------------------------------------------------------------------- 427 &nam_dynadv ! formulation of the momentum advection 419 !! namdyn_adv formulation of the momentum advection 420 !! namdyn_vor advection scheme 421 !! namdyn_hpg hydrostatic pressure gradient 422 !! namdyn_spg surface pressure gradient (CPP key only) 423 !! namdyn_ldf lateral diffusion scheme 424 !!====================================================================== 425 426 !----------------------------------------------------------------------- 427 &namdyn_adv ! formulation of the momentum advection 428 428 !----------------------------------------------------------------------- 429 429 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) … … 432 432 / 433 433 !----------------------------------------------------------------------- 434 &nam _dynvor ! option of physics/algorithm (not control by CPP keys)434 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 435 435 !----------------------------------------------------------------------- 436 436 ln_dynvor_ene = .false. ! enstrophy conserving scheme 437 437 ln_dynvor_ens = .false. ! energy conserving scheme 438 438 ln_dynvor_mix = .false. ! mixed scheme 439 ln_dynvor_een = .true. 440 / 441 !----------------------------------------------------------------------- 442 &nam _dynhpg ! Hydrostatic pressure gradient option439 ln_dynvor_een = .true. ! energy & enstrophy scheme 440 / 441 !----------------------------------------------------------------------- 442 &namdyn_hpg ! Hydrostatic pressure gradient option 443 443 !----------------------------------------------------------------------- 444 444 ln_hpg_zco = .false. ! z-coordinate - full steps … … 449 449 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 450 450 ln_hpg_rot = .false. ! s-coordinate (ROTated axes scheme) 451 gamm = 0.e0 ! weighting coefficient (wdj scheme) 452 / 453 !----------------------------------------------------------------------- 454 &namflg ! algorithm flags (algorithm not control by CPP keys) 455 !----------------------------------------------------------------------- 456 ln_dynhpg_imp = .false. ! hydrostatic pressure gradient: semi-implicit time scheme (T) 457 ! centered time scheme (F) 458 nn_dynhpg_rst = 0 ! add dynhpg implicit variables in restart ot not (1/0) 459 / 460 !----------------------------------------------------------------------- 461 !nam_dynspg ! surface pressure gradient (CPP key only) 451 rn_gamma = 0.e0 ! weighting coefficient (wdj scheme) 452 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 453 ! centered time scheme (F) 454 nn_dynhpg_rst = 0 ! =1 dynhpg restartable restart or not (=0) 455 / 456 !----------------------------------------------------------------------- 457 !namdyn_spg ! surface pressure gradient (CPP key only) 462 458 !----------------------------------------------------------------------- 463 459 ! ! explicit free surface ("key_dynspg_exp") … … 466 462 467 463 !----------------------------------------------------------------------- 468 &nam _dynldf ! lateral diffusion on momentum469 !----------------------------------------------------------------------- 470 !! Type of the operator :471 ln_dynldf_lap = .false. ! 472 ln_dynldf_bilap = .true. ! 473 !! Direction of action :474 ln_dynldf_level = .false. ! 475 ln_dynldf_hor = .true. ! horizontal (geopotential)(require "key_ldfslp" in s-coord.)476 ln_dynldf_iso = .false. ! iso-neutral(require "key_ldfslp")477 478 ahm0 = -8.5e+11 !horizontal eddy viscosity [m2/s]479 ahmb0 = 0. !background eddy viscosity for ldf_iso [m2/s]464 &namdyn_ldf ! lateral diffusion on momentum 465 !----------------------------------------------------------------------- 466 ! Type of the operator : 467 ln_dynldf_lap = .false. ! laplacian operator 468 ln_dynldf_bilap = .true. ! bilaplacian operator 469 ! Direction of action : 470 ln_dynldf_level = .false. ! iso-level 471 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 472 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") 473 ! Coefficient 474 rn_ahm_0 = -8.5e+11 ! horizontal eddy viscosity [m2/s] 475 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 480 476 / 481 477 !!====================================================================== 482 478 !! Tracers & Dynamics vertical physics namelists 483 479 !!====================================================================== 484 !! nam_zdf vertical physics 485 !! nam_npc non penetrative convection 486 !! nam_ric richardson number dependent vertical mixing ("key_zdfric" ) 487 !! nam_tke TKE dependent vertical mixing ("key_zdftke" ) 488 !! nam_kpp KPP dependent vertical mixing ("key_zdfkpp" ) 489 !! nam_ddm double diffusive mixing parameterization ("key_zdfddm" ) 490 !! nam_tmx tidal mixing parameterization ("key_zdftmx" ) 491 !!====================================================================== 492 493 !----------------------------------------------------------------------- 494 &nam_zdf ! vertical physics 495 !----------------------------------------------------------------------- 496 ! vertical eddy coef. or their background values 480 !! namzdf vertical physics 481 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric" ) 482 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke" ) 483 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp" ) 484 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm" ) 485 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx" ) 486 !!====================================================================== 487 488 !----------------------------------------------------------------------- 489 &namzdf ! vertical physics 490 !----------------------------------------------------------------------- 497 491 rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") 498 492 rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") 499 493 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 500 494 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) 501 ln_zdfevd = .true. ! convection:enhanced vertical diffusion (evd) (T) or not (F)502 nn_evdm = 0 ! 503 rn_avevd = 100. ! 504 ln_zdfnpc = .false. ! convection:Non-Penetrative algorithm (T) or not (F)505 nn_npc = 1 ! 506 nn_npcp = 365 ! 495 ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F) 496 nn_evdm = 0 ! evd apply on tracer (=0) or on tracer and momentum (=1) 497 rn_avevd = 100. ! evd mixing coefficient [m2/s] 498 ln_zdfnpc = .false. ! Non-Penetrative algorithm (T) or not (F) 499 nn_npc = 1 ! frequency of application of npc 500 nn_npcp = 365 ! npc control print frequency 507 501 ln_zdfexp = .false. ! time-stepping: split-explicit (T) or implicit (F) time stepping 508 502 nn_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T 509 503 / 510 504 !----------------------------------------------------------------------- 511 &nam _ric! richardson number dependent vertical diffusion ("key_zdfric" )512 !----------------------------------------------------------------------- 513 rn_avmri = 100.e-4 ! avm = rn_avmri / ( 1 + rn_alp * Ri**nn_ric )514 rn_alp = 5. ! avt = avm / ( 1 + rn_alp * Ri )515 nn_ric = 2 ! 516 / 517 !----------------------------------------------------------------------- 518 &nam _tke! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke")505 &namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" ) 506 !----------------------------------------------------------------------- 507 rn_avmri = 100.e-4 ! maximum value of the vertical viscosity 508 rn_alp = 5. ! coefficient of the parameterization 509 nn_ric = 2 ! coefficient of the parameterization 510 / 511 !----------------------------------------------------------------------- 512 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 519 513 !----------------------------------------------------------------------- 520 514 rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) … … 543 537 rn_lc = 0.15 ! coef. associated to Langmuir cells 544 538 / 545 !------------------------------------------------------------------------ 546 &nam _kpp! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionnally:539 !------------------------------------------------------------------------ 540 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionnally: 547 541 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 548 542 ln_kpprimix = .true. ! shear instability mixing … … 555 549 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv 556 550 nn_ave = 1 ! constant (=0) or profile (=1) background on avt 557 / 558 !----------------------------------------------------------------------- 559 &nam _ddm ! double diffusive mixing parameterization("key_zdfddm")551 / 552 !----------------------------------------------------------------------- 553 &namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm") 560 554 !----------------------------------------------------------------------- 561 555 rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity) … … 563 557 / 564 558 !----------------------------------------------------------------------- 565 &nam_tmx ! tidal mixing parameterization ("key_zdftmx") 566 !----------------------------------------------------------------------- 567 rn_htmx = 500. ! vertical decay scale for turbulence (meters) 568 rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1) 569 rn_tfe = 0.333 ! tidal dissipation efficiency 570 rn_me = 0.2 ! mixing efficiency 571 ln_tmx_itf = .FALSE. ! ITF specific parameterisation 572 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency 573 / 574 559 &namzdf_tmx ! tidal mixing parameterization ("key_zdftmx") 560 !----------------------------------------------------------------------- 561 rn_htmx = 500. ! vertical decay scale for turbulence (meters) 562 rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1) 563 rn_tfe = 0.333 ! tidal dissipation efficiency 564 rn_me = 0.2 ! mixing efficiency 565 ln_tmx_itf = .FALSE. ! ITF specific parameterisation 566 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency 567 / 575 568 !!====================================================================== 576 569 !! *** Miscelaneous namelists *** 577 570 !!====================================================================== 578 !! nam _mppMassively Parallel Processing ("key_mpp_mpi)579 !! nam _mpp_dyndistMassively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist")571 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 572 !! nammpp_dyndist Massively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist") 580 573 !! namctl Control prints & Benchmark 581 574 !! namsol elliptic solver / island / free surface … … 585 578 &namsol ! elliptic solver / island / free surface 586 579 !----------------------------------------------------------------------- 587 n solv= 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg)580 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) 588 581 ! =2 successive-over-relaxation (sor) 589 nsol_arp = 0 ! absolute/relative (0/1) precision convergence test 590 nmin = 300 ! minimum of iterations for the SOR solver 591 nmax = 800 ! maximum of iterations for the SOR solver 592 nmod = 10 ! frequency of test for the SOR solver 593 eps = 1.e-6 ! absolute precision of the solver 594 resmax = 1.e-10 ! absolute precision for the SOR solver 595 sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 596 rnu = 1. ! strength of the additional force used in filtered free surface 597 / 598 !----------------------------------------------------------------------- 599 &nam_mpp ! Massively Parallel Processing ("key_mpp_mpi) 600 !----------------------------------------------------------------------- 601 c_mpi_send = 'S' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, 582 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test 583 rn_eps = 1.e-6 ! absolute precision of the solver 584 nn_nmin = 300 ! minimum of iterations for the SOR solver 585 nn_nmax = 800 ! maximum of iterations for the SOR solver 586 nn_nmod = 10 ! frequency of test for the SOR solver 587 rn_resmax = 1.e-10 ! absolute precision for the SOR solver 588 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 589 / 590 !----------------------------------------------------------------------- 591 &nammpp ! Massively Parallel Processing ("key_mpp_mpi) 592 !----------------------------------------------------------------------- 593 cn_mpi_send = 'S' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, 602 594 ! buffer blocking send or immediate non-blocking sends, resp. 603 nn_buffer = 0! size in bytes of exported buffer ('B' case), 0 no exportation604 / 605 !----------------------------------------------------------------------- 606 &nam _mpp_dyndist ! Massively Parallel Distribution("key_agrif" && "key_mpp_dyndist")607 !----------------------------------------------------------------------- 608 jpni = 1! jpni number of processors following i609 jpnj = 1! jpnj number of processors following j610 jpnij = 1! jpnij number of local domains611 / 612 !----------------------------------------------------------------------- 613 &namctl ! Control prints & Benchmark614 !----------------------------------------------------------------------- 615 ln_ctl = .false.! trends control print (expensive!)616 n print = 0! level of print (0 no extra print)617 n ictls = 0! start i indice of control sum (use to compare mono versus618 n ictle = 0! end i indice of control sum multi processor runs619 n jctls = 0! start j indice of control over a subdomain)620 n jctle = 0 ! end j indice of control621 isplt = 1! number of processors in i-direction622 jsplt = 1! number of processors in j-direction623 n bench = 0! Bench mode (1/0): CAUTION use zero except for bench595 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 596 / 597 !----------------------------------------------------------------------- 598 &nammpp_dyndist ! Massively Parallel Distribution for AGRIF zoom ("key_agrif" && "key_mpp_dyndist") 599 !----------------------------------------------------------------------- 600 jpni = 1 ! jpni number of processors following i 601 jpnj = 1 ! jpnj number of processors following j 602 jpnij = 1 ! jpnij number of local domains 603 / 604 !----------------------------------------------------------------------- 605 &namctl ! Control prints & Benchmark 606 !----------------------------------------------------------------------- 607 ln_ctl = .false. ! trends control print (expensive!) 608 nn_print = 0 ! level of print (0 no extra print) 609 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 610 nn_ictle = 0 ! end i indice of control sum multi processor runs 611 nn_jctls = 0 ! start j indice of control over a subdomain) 612 nn_jctle = 0 ! end j indice of control 613 nn_isplt = 1 ! number of processors in i-direction 614 nn_jsplt = 1 ! number of processors in j-direction 615 nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 624 616 ! (no physical validity of the results) 625 n bit_cmp = 0! bit comparison mode (1/0): CAUTION use zero except for test617 nn_bit_cmp = 0 ! bit comparison mode (1/0): CAUTION use zero except for test 626 618 ! of comparison between single and multiple processor runs 627 619 / 620 628 621 !!====================================================================== 629 622 !! *** Diagnostics namelists *** … … 636 629 637 630 !----------------------------------------------------------------------- 638 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 639 ! or mixed-layer trends ('key_trdmld') 640 ! or barotropic vorticity ("key_trdvor") 641 !----------------------------------------------------------------------- 642 ntrd = 365 ! time step frequency dynamics and tracers trends 643 nctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 644 ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 645 cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 646 cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 647 ln_trdmld_restart = .false. ! restart for ML diagnostics 648 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 631 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 632 ! ! or mixed-layer trends or barotropic vorticity ('key_trdmld' or "key_trdvor") 633 !----------------------------------------------------------------------- 634 nn_trd = 365 ! time step frequency dynamics and tracers trends 635 nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 636 rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 637 cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 638 cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 639 ln_trdmld_restart = .false. ! restart for ML diagnostics 640 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 649 641 / 650 642 !----------------------------------------------------------------------- 651 643 &namgap ! level mean model-data gap ('key_diagap') 652 644 !----------------------------------------------------------------------- 653 n gap= 15 ! time-step frequency of model-data gap computation654 n prg= 10 ! time-step frequency of gap print in model output645 nn_gap = 15 ! time-step frequency of model-data gap computation 646 nn_prg = 10 ! time-step frequency of gap print in model output 655 647 / 656 648 !----------------------------------------------------------------------- … … 658 650 !----------------------------------------------------------------------- 659 651 ln_rstflo = .false. ! float restart (T) or not (F) 660 n writefl= 75 ! frequency of writing in float output file661 n stockfl= 5475 ! frequency of creation of the float restart file652 nn_writefl= 75 ! frequency of writing in float output file 653 nn_stockfl= 5475 ! frequency of creation of the float restart file 662 654 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 663 655 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) … … 667 659 &namptr ! Poleward Transport Diagnostic 668 660 !----------------------------------------------------------------------- 669 ln_diaptr = .false. 670 ln_diaznl = .false. 671 ln_subbas = .false. 661 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 662 ln_diaznl = .false. ! Add zonal means and meridional stream functions 663 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 672 664 ! (orca configuration only, need input basins mask file named "subbasins.nc" 673 665 nf_ptr = 1 ! Frequency of ptr computation [time step] -
trunk/CONFIG/ORCA2_LIM/EXP00/namelist
r1556 r1602 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 2 !! NEMO/OPA : 1 - run manager (namrun) 3 !! namelists 2 - Domain (nam _zgr, nam_zgr_sco, namdom)3 !! namelists 2 - Domain (namzgr, namzgr_sco, namdom) 4 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, nam alb)6 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nam tide)5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namsbc_alb) 6 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 7 7 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 8 !! 6 - Tracer (nameos, nam _traadv, nam_traldf, namtdp)9 !! 7 - dynamics (nam _dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf)10 !! 8 - Verical physics (nam _zdf, nam_npc, nam_ric, nam_tke, nam_kpp, nam_ddm, nam_tmx)11 !! 9 - diagnostics (namtrd, namgap, nam flo, namptr)12 !! 9 - miscellaneous (namsol, nam _mpp, nam_mpp_dyndist, namctl)8 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_dmp) 9 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 10 !! 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) 13 13 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 14 ! CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE. … … 23 23 &namrun ! parameters of the run 24 24 !----------------------------------------------------------------------- 25 no = 0 ! job number 26 cexper = "ORCA2" ! experience name 25 nn_no = 0 ! job number 26 cn_exp = "ORCA2" ! experience name 27 nn_it000 = 1 ! first time step 28 nn_itend = 315 ! last time step (std 5475) 29 nn_date0 = 010101 ! initial calendar date yymmdd (used if nrstdt=1) 30 nn_leapy = 0 ! Leap year calendar (1) or not (0) 31 nn_istate = 0 ! output the initial state (1) or not (0) 32 nn_stock = 5475 ! frequency of creation of a restart file (modulo referenced to 1) 33 nn_write = 5475 ! frequency of write in the output file (modulo referenced to nit000) 34 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) 35 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) 36 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. ! start from rest (F) or from a restart file (T) 39 nn_rstctl = 0 ! restart control = 0 nit000 is not compared to the restart file value 40 ! = 1 use ndate0 in namelist (not the value in the restart file) 41 ! = 2 calendar parameters read in the restart file 27 42 cn_ocerst_in = "restart" ! suffix of ocean restart name (input) 28 43 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 29 ln_rstart = .false. ! start from rest (F) or from a restart file (T)30 nrstdt = 0 ! restart control = 0 nit000 is not compared to the restart file value31 ! = 1 use ndate0 in namelist (not the value in the restart file)32 ! = 2 calendar parameters read in the restart file33 nit000 = 1 ! first time step34 nitend = 5475 ! last time step35 ndate0 = 010101 ! initial calendar date yymmdd (used if nrstdt=1)36 nleapy = 0 ! Leap year calendar (1) or not (0)37 ninist = 0 ! output the initial state (1) or not (0)38 nstock = 5475 ! frequency of creation of a restart file (modulo referenced to 1)39 nwrite = 5475 ! frequency of write in the output file (modulo referenced to nit000)40 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T)41 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%)42 ln_clobber = .false. ! clobber (overwrite) an existing file43 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (working only with iom_nf90 routines)44 44 / 45 45 !!====================================================================== 46 46 !! *** Domain namelists *** 47 47 !!====================================================================== 48 !! nam _zgr vertical coordinate49 !! nam _zgr_sco s-coordinate or hybrid z-s-coordinate50 !! namdom 51 !!====================================================================== 52 53 !----------------------------------------------------------------------- 54 &nam _zgr! vertical coordinate48 !! namzgr vertical coordinate 49 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 50 !! namdom space and time domain (bathymetry, mesh, timestep) 51 !!====================================================================== 52 53 !----------------------------------------------------------------------- 54 &namzgr ! vertical coordinate 55 55 !----------------------------------------------------------------------- 56 56 ln_zco = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) … … 59 59 / 60 60 !----------------------------------------------------------------------- 61 &nam _zgr_sco! s-coordinate or hybrid z-s-coordinate62 !----------------------------------------------------------------------- 63 sbot_min= 300. ! minimum depth of s-bottom surface (>0) (m)64 sbot_max= 5250. ! maximum depth of s-bottom surface (= ocean depth) (>0) (m)65 theta= 6.0 ! surface control parameter (0<=theta<=20)66 thetb= 0.75 ! bottom control parameter (0<=thetb<= 1)67 r _max= 0.15 ! maximum cut-off r-value allowed (0<r_max<1)61 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate 62 !----------------------------------------------------------------------- 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) 65 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 68 ln_s_sigma = .false. ! hybrid s-sigma coordinates 69 bb= 0.8 ! stretching with s-sigma70 hc= 150.0 ! critical depth with s-sigma69 rn_bb = 0.8 ! stretching with s-sigma 70 rn_hc = 150.0 ! critical depth with s-sigma 71 71 / 72 72 !----------------------------------------------------------------------- 73 73 &namdom ! space and time domain (bathymetry, mesh, timestep) 74 74 !----------------------------------------------------------------------- 75 ntopo = 1 ! compute (=0) or read(=1) the bathymetry file 76 e3zps_min = 5. ! the thickness of the partial step is set larger than the minimum 77 e3zps_rat = 0.1 ! of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 78 nmsh = 0 ! create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 79 nacc = 0 ! =1 acceleration of convergence method used, rdt < rdttra(k) 80 ! =0, no acceleration, rdt = rdttra 81 atfp = 0.1 ! asselin time filter parameter 82 rdt = 5760. ! time step for the dynamics (and tracer if nacc=0) 83 rdtmin = 5760. ! minimum time step on tracers (used if nacc=1) 84 rdtmax = 5760. ! maximum time step on tracers (used if nacc=1) 85 rdth = 800. ! depth variation of tracer time step (used if nacc=1) 75 nn_bathy = 1 ! compute (=0) or read(=1) the bathymetry file 76 nn_closea = 0 ! closed seas and lakes are removed (=0) or kept (=1) from the ORCA domain 77 nn_msh = 0 ! create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 78 rn_e3zps_min= 5. ! 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) 80 ! 81 rn_rdt = 5760. ! time step for the dynamics (and tracer if nacc=0) ==> 5760 86 82 nn_baro = 64 ! number of barotropic time step (for the split explicit algorithm) ("key_dynspg_ts") 87 nclosea = 0 ! = 0 no closed sea in the model domain 88 ! = 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 83 rn_atfp = 0.1 ! asselin time filter parameter 84 nn_acc = 0 ! acceleration of convergence : =1 used, rdt < rdttra(k) 85 ! =0, not used, rdt = rdttra 86 rn_rdtmin = 5760. ! minimum time step on tracers (used if nacc=1) 87 rn_rdtmax = 5760. ! maximum time step on tracers (used if nacc=1) 88 rn_rdth = 800. ! depth variation of tracer time step (used if nacc=1) 89 89 / 90 90 !!====================================================================== … … 97 97 !! namsbc_core CORE bulk formulea formulation 98 98 !! namsbc_cpl CouPLed formulation ("key_coupled") 99 !! nam qsrpenetrative solar radiation99 !! namtra_qsr penetrative solar radiation 100 100 !! namsbc_rnf river runoffs 101 101 !! namsbc_ssr sea surface restoring term (for T and/or S) 102 !! nam albalbedo parameters102 !! namsbc_alb albedo parameters 103 103 !!====================================================================== 104 104 … … 122 122 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 123 123 ln_ssr = .true. ! 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 nn_fsbc time step ,124 nn_fwb = 0 ! FreshWater Budget: =0 unchecked 125 ! =1 global mean of e-p-r set to zero at each time step 126 126 ! =2 annual global mean of e-p-r set to zero 127 127 / … … 139 139 &namsbc_flx ! surface boundary condition : flux formulation 140 140 !----------------------------------------------------------------------- 141 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !142 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !143 sn_utau = 'utau' , 24. , 'utau' , .false. , .false. ,'yearly' , '' , ''144 sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false. ,'yearly' , '' , ''145 sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false. ,'yearly' , '' , ''146 sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false. ,'yearly' , '' , ''147 sn_emp = 'emp' , 24. , 'emp' , .false. , .false. ,'yearly' , '' , ''141 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 142 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 143 sn_utau = 'utau' , 24. , 'utau' , .false. , .false. , 'yearly' , '' , '' 144 sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false. , 'yearly' , '' , '' 145 sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false. , 'yearly' , '' , '' 146 sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false. , 'yearly' , '' , '' 147 sn_emp = 'emp' , 24. , 'emp' , .false. , .false. , 'yearly' , '' , '' 148 148 ! 149 149 cn_dir = './' ! root directory for the location of the flux files … … 152 152 &namsbc_clio ! namsbc_clio CLIO bulk formulea 153 153 !----------------------------------------------------------------------- 154 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !155 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !156 sn_utau = 'taux_1m' , -1. , 'sozotaux' , .true. , .true. ,'yearly' , '' , ''157 sn_vtau = 'tauy_1m' , -1. , 'sometauy' , .true. , .true. ,'yearly' , '' , ''158 sn_wndm = 'flx' , -1. , 'socliowi' , .true. , .true. ,'yearly' , '' , ''159 sn_tair = 'flx' , -1. , 'socliot2' , .true. , .true. ,'yearly' , '' , ''160 sn_humi = 'flx' , -1. , 'socliohu' , .true. , .true. ,'yearly' , '' , ''161 sn_ccov = 'flx' , -1. , 'socliocl' , .false. , .true. ,'yearly' , '' , ''162 sn_prec = 'flx' , -1. , 'socliopl' , .false. , .true. ,'yearly' , '' , ''154 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 155 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 156 sn_utau = 'taux_1m' , -1. , 'sozotaux' , .true. , .true. , 'yearly' , '' , '' 157 sn_vtau = 'tauy_1m' , -1. , 'sometauy' , .true. , .true. , 'yearly' , '' , '' 158 sn_wndm = 'flx' , -1. , 'socliowi' , .true. , .true. , 'yearly' , '' , '' 159 sn_tair = 'flx' , -1. , 'socliot2' , .true. , .true. , 'yearly' , '' , '' 160 sn_humi = 'flx' , -1. , 'socliohu' , .true. , .true. , 'yearly' , '' , '' 161 sn_ccov = 'flx' , -1. , 'socliocl' , .false. , .true. , 'yearly' , '' , '' 162 sn_prec = 'flx' , -1. , 'socliopl' , .false. , .true. , 'yearly' , '' , '' 163 163 ! 164 164 cn_dir = './' ! root directory for the location of the bulk files are … … 167 167 &namsbc_core ! namsbc_core CORE bulk formulea 168 168 !----------------------------------------------------------------------- 169 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! weights! rotation !170 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename! pairing !171 sn_wndi = 'u10_core' , -1. , 'u10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' ,'U1'172 sn_wndj = 'v10_core' , -1. , 'v10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' ,'V1'173 sn_qsr = 'qsw_core' , -1. , 'swdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''174 sn_qlw = 'qlw_core' , -1. , 'lwdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''175 sn_tair = 't2_core' , -1. , 't2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''176 sn_humi = 'q2_core' , -1. , 'q2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''177 sn_prec = 'precip_core' , -1. , 'precip' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''178 sn_snow = 'snow_core' , -1. , 'snow' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc',''169 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 170 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 171 sn_wndi = 'u10_core' , -1. , 'u10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'U1' 172 sn_wndj = 'v10_core' , -1. , 'v10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'V1' 173 sn_qsr = 'qsw_core' , -1. , 'swdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 174 sn_qlw = 'qlw_core' , -1. , 'lwdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 175 sn_tair = 't2_core' , -1. , 't2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 176 sn_humi = 'q2_core' , -1. , 'q2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 177 sn_prec = 'precip_core', -1. , 'precip' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 178 sn_snow = 'snow_core' , -1. , 'snow' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', '' 179 179 ! 180 180 cn_dir = './' ! root directory for the location of the bulk files 181 181 ln_2m = .true. ! air temperature and humidity referenced at 2m (T) instead 10m (F) 182 alpha_precip= 1. ! multiplicative factor for precipitation (total & snow)182 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 183 183 / 184 184 !----------------------------------------------------------------------- 185 185 &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") 186 186 !----------------------------------------------------------------------- 187 ! SEND 187 ! send 188 188 cn_snd_temperature= 'weighted oce and ice' ! 'oce only' 'weighted oce and ice' 'mixed oce-ice' 189 189 cn_snd_albedo = 'weighted ice' ! 'none' 'weighted ice' 'mixed oce-ice' … … 193 193 cn_snd_crt_orient = 'eastward-northward' ! 'eastward-northward' or 'local grid' 194 194 cn_snd_crt_grid = 'T' ! 'T' 195 ! RECEIVE 195 ! receive 196 196 cn_rcv_w10m = 'coupled' ! 'none' 'coupled' 197 197 cn_rcv_tau_nature = 'oce only' ! 'oce only' 'oce and ice' 'mixed oce-ice' … … 209 209 &namsbc_cpl_co2 ! coupled ocean/biogeo/atmosphere model ("key_cpl_carbon_cycle") 210 210 !----------------------------------------------------------------------- 211 ! SEND 212 cn_snd_co2 = 'coupled' ! 'none' 'coupled' 213 ! RECEIVE 214 cn_rcv_co2 = 'coupled' ! 'none' 'coupled' 215 / 216 !----------------------------------------------------------------------- 217 &namqsr ! penetrative solar radiation 218 !----------------------------------------------------------------------- 219 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! weights ! rotation ! 220 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 221 sn_chl = 'chlorophyll' , -1. , 'CHLA' , .true. , .true. , 'yearly' , '' , '' 222 ! 211 cn_snd_co2 = 'coupled' ! send : 'none' 'coupled' 212 cn_rcv_co2 = 'coupled' ! receive : 'none' 'coupled' 213 / 214 !----------------------------------------------------------------------- 215 &namtra_qsr ! penetrative solar radiation 216 !----------------------------------------------------------------------- 217 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 218 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 219 sn_chl = 'chlorophyll', -1. , 'CHLA' , .true. , .true. , 'yearly' , '' , '' 220 223 221 cn_dir = './' ! root directory for the location of the runoff files 224 222 ln_traqsr = .true. ! Light penetration (T) or not (F) 225 ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration223 ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration 226 224 ln_qsr_2bd = .false. ! 2 bands light penetration 227 225 ln_qsr_bio = .false. ! bio-model light penetration … … 235 233 &namsbc_rnf ! runoffs namelist surface boundary condition 236 234 !----------------------------------------------------------------------- 237 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !238 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !239 sn_rnf = 'runoff_1m_nomask' , -1. , 'sorunoff' , .true. , .true. ,'yearly' , '' , ''240 sn_cnf = 'runoff_1m_nomask' , 0. , 'socoefr' , .false. , .true. ,'yearly' , '' , ''241 ! 235 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 236 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 237 sn_rnf = 'runoff_1m_nomask' , -1. , 'sorunoff' , .true. , .true. , 'yearly' , '' , '' 238 sn_cnf = 'runoff_1m_nomask' , 0. , 'socoefr' , .false. , .true. , 'yearly' , '' , '' 239 242 240 cn_dir = './' ! root directory for the location of the runoff files 243 241 ln_rnf_emp = .false. ! runoffs included into precipitation field (T) or into a file (F) … … 245 243 rn_hrnf = 0.e0 ! depth over which enhanced vertical mixing is used 246 244 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] 247 rn_ mul_rnf= 1.e0 ! multiplicative factor for runoff245 rn_rfact = 1.e0 ! multiplicative factor for runoff 248 246 / 249 247 !----------------------------------------------------------------------- 250 248 &namsbc_ssr ! surface boundary condition : sea surface restoring 251 249 !----------------------------------------------------------------------- 252 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or! weights ! rotation !253 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly'! filename ! pairing !254 sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false. , 'yearly', '' , ''255 sn_sss = 'sss_data' , -1. , 'sss' , .true. , .true. , 'yearly', '' , ''256 !250 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 251 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 252 sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false. , 'yearly' , '' , '' 253 sn_sss = 'sss_data' , -1. , 'sss' , .true. , .true. , 'yearly' , '' , '' 254 257 255 cn_dir = './' ! root directory for the location of the runoff files 258 256 nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0) 259 257 nn_sssr = 1 ! add a damping term in the surface freshwater flux (=1) or not (=0) 260 dqdt= -40. ! magnitude of the retroaction on temperature [W/m2/K]261 deds= -27.7 ! magnitude of the damping on salinity [mm/day/psu]258 rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 259 rn_deds = -27.7 ! magnitude of the damping on salinity [mm/day/psu] 262 260 ln_sssr_bnd = .false. ! flag to bound erp term (associated with nn_sssr=2) 263 261 rn_sssr_bnd = 0.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] 264 262 / 265 263 !----------------------------------------------------------------------- 266 &namalb ! albedo parameters 267 !----------------------------------------------------------------------- 268 cgren = 0.06 ! correction of the snow or ice albedo to take into account the 269 albice = 0.53 ! albedo of melting ice in the arctic and antarctic 270 alphd = 0.80 ! coefficients for linear interpolation used to 271 alphc = 0.65 ! compute albedo between two extremes values 272 alphdi = 0.72 ! (Pyane, 1972) 273 / 264 &namsbc_alb ! albedo parameters 265 !----------------------------------------------------------------------- 266 rn_cloud = 0.06 ! cloud correction to snow and ice albedo 267 rn_albice = 0.53 ! albedo of melting ice in the arctic and antarctic 268 rn_alphd = 0.80 ! coefficients for linear interpolation used to 269 rn_alphc = 0.65 ! compute albedo between two extremes values 270 rn_alphdi = 0.72 ! (Pyane, 1972) 271 / 272 274 273 !!====================================================================== 275 274 !! *** Lateral boundary condition *** … … 286 285 &namlbc ! lateral momentum boundary condition 287 286 !----------------------------------------------------------------------- 288 shlat = 2. ! shlat = 0 : free slip 289 ! 0 < shlat < 2 : partial slip 290 ! shlat = 2 : no slip 291 ! 2 < shlat : strong slip 287 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 288 ! free slip ! partial slip ! no slip ! strong slip 292 289 / 293 290 !----------------------------------------------------------------------- 294 291 &namcla ! cross land advection 295 292 !----------------------------------------------------------------------- 296 n _cla= 1 ! advection between 2 ocean pts separates by land293 nn_cla = 1 ! advection between 2 ocean pts separates by land 297 294 / 298 295 !----------------------------------------------------------------------- 299 296 &namobc ! open boundaries parameters ("key_obc") 300 297 !----------------------------------------------------------------------- 301 nobc_dta = 1 ! = 0 the obc data are equal to the initial state302 ! = 1 the obc data are read in 'obc.dta' files303 cffile = 'annual' ! set to annual if obc datafile hold 1 year of data304 ! set to monthly if obc datafile hold 1 month of data305 rdpein = 1. ! ???306 rdpwin = 1. ! ???307 rdpnin = 1. ! ???308 rdpsin = 1. ! ???309 rdpeob = 3000. ! time relaxation (days) for the east open boundary310 rdpwob = 15. ! " " " west "311 rdpnob = 3000. ! " " " north "312 rdpsob = 15. ! " " " south "313 zbsic1 = 140.e+6 ! barotropic stream function on first isolated coastline314 zbsic2 = 1.e+6 ! " " second "315 zbsic3 = 0. ! " " thrid "316 298 ln_obc_clim= .false. ! climatological obc data files (T) or not (F) 317 299 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) 318 / 319 !----------------------------------------------------------------------- 320 &namagrif ! ("key_agrif") 321 !----------------------------------------------------------------------- 322 nbclineupdate = 3 ! baroclinic update frequency 300 ln_obc_fla = .false. ! Flather open boundary condition 301 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state 302 ! = 1 the obc data are read in 'obc.dta' files 303 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data 304 ! set to monthly if obc datafile hold 1 month of data 305 rn_dpein = 1. ! damping time scale for inflow at east open boundary 306 rn_dpwin = 1. ! - - - west - - 307 rn_dpnin = 1. ! - - - north - - 308 rn_dpsin = 1. ! - - - south - - 309 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary 310 rn_dpwob = 15. ! - - - west - - 311 rn_dpnob = 3000. ! - - - north - - 312 rn_dpsob = 15. ! - - - south - - 313 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R) 314 ! = 1 the total volume remains constant 315 / 316 !----------------------------------------------------------------------- 317 &namagrif ! AGRIF zoom ("key_agrif") 318 !----------------------------------------------------------------------- 319 nn_cln_update = 3 ! baroclinic update frequency 323 320 ln_spc_dyn = .true. ! use 0 as special value for dynamics 324 visc_tra = 2880. ! viscosity coeeficient for tracers sponge layer325 visc_dyn = 2880. ! viscosity coeeficient for dynamics sponge layer326 / 327 !----------------------------------------------------------------------- 328 &nambdy ! unstructured open boundaries parameters("key_bdy")321 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [s] 322 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [s] 323 / 324 !----------------------------------------------------------------------- 325 &nambdy ! unstructured open boundaries ("key_bdy") 329 326 !----------------------------------------------------------------------- 330 327 filbdy_mask = '' ! name of mask file (if ln_bdy_mask=.TRUE.) … … 346 343 / 347 344 !----------------------------------------------------------------------- 348 &nam tide! tidal forcing at unstructured boundaries345 &nambdy_tide ! tidal forcing at unstructured boundaries 349 346 !----------------------------------------------------------------------- 350 347 filtide = 'bdytide_' ! file name root of tidal forcing files … … 353 350 ln_tide_date = .false. ! adjust tidal harmonics for start date of run 354 351 / 352 355 353 !!====================================================================== 356 354 !! *** Bottom boundary condition *** … … 364 362 &nambfr ! bottom friction 365 363 !----------------------------------------------------------------------- 366 n botfr = 1 ! type of bottom friction : = 0 : no slip, = 2 : nonlinear friction364 nn_bfr = 1 ! type of bottom friction : = 0 : no slip, = 2 : nonlinear friction 367 365 ! = 3 : free slip, = 1 : linear friction 368 bfri1= 4.e-4 ! bottom drag coefficient (linear case)369 bfri2= 1.e-3 ! bottom drag coefficient (non linear case)370 bfeb2= 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2)366 rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 367 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) 368 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2) 371 369 / 372 370 !----------------------------------------------------------------------- 373 371 &nambbc ! bottom temperature boundary condition 374 372 !----------------------------------------------------------------------- 375 n geo_flux = 2 ! geothermal heat flux = 0 no flux considered376 ! 377 ! 378 ngeo_flux_const = 86.4e-3! Constant value of geothermal heat flux [W/m2]373 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 374 ! = 1 constant flux 375 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 376 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 379 377 / 380 378 !----------------------------------------------------------------------- … … 383 381 ! ! diffusive bbl ("key_trabbl") 384 382 ! ! advective bbl ("key_trabbl_adv") 385 atrbbl= 10000. ! lateral mixing coefficient in the bbl [m2/s]383 rn_ahtbbl = 10000. ! lateral mixing coefficient in the bbl [m2/s] 386 384 / 387 385 !!====================================================================== … … 389 387 !!====================================================================== 390 388 !! nameos equation of state 391 !! nam _traadv advection scheme392 !! nam _traldf lateral diffusion scheme393 !! namt dp tracer newtonian damping("key_tradmp")389 !! namtra_adv advection scheme 390 !! namtra_ldf lateral diffusion scheme 391 !! namtra_dmp T & S newtonian damping ("key_tradmp") 394 392 !!====================================================================== 395 393 … … 397 395 &nameos ! ocean physical parameters 398 396 !----------------------------------------------------------------------- 399 n eos= 0 ! type of equation of state and Brunt-Vaisala frequency397 nn_eos = 0 ! type of equation of state and Brunt-Vaisala frequency 400 398 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 401 399 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 402 400 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 403 r alpha= 2.e-4 ! thermal expension coefficient (neos= 1 or 2)404 r beta= 0.001 ! saline expension coefficient (neos= 2)405 / 406 !----------------------------------------------------------------------- 407 &nam _traadv ! advection scheme for tracer401 rn_alpha = 2.e-4 ! thermal expension coefficient (neos= 1 or 2) 402 rn_beta = 0.001 ! saline expension coefficient (neos= 2) 403 / 404 !----------------------------------------------------------------------- 405 &namtra_adv ! advection scheme for tracer 408 406 !----------------------------------------------------------------------- 409 407 ln_traadv_cen2 = .true. ! 2nd order centered scheme … … 414 412 / 415 413 !----------------------------------------------------------------------- 416 &nam _traldf ! lateral diffusion scheme for tracer417 !----------------------------------------------------------------------- 418 !! Type of the operator :414 &namtra_ldf ! lateral diffusion scheme for tracer 415 !----------------------------------------------------------------------- 416 ! Type of the operator : 419 417 ln_traldf_lap = .true. ! laplacian operator 420 418 ln_traldf_bilap = .false. ! bilaplacian operator 421 419 ! Direction of action : 422 420 ln_traldf_level = .false. ! iso-level 423 421 ln_traldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 424 422 ln_traldf_iso = .true. ! iso-neutral (require "key_ldfslp") 425 ! ! Coefficient 426 aht0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 427 ahtb0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 428 aeiv0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 429 / 430 !----------------------------------------------------------------------- 431 &namtdp ! tracer newtonian damping ('key_tradmp') 432 !----------------------------------------------------------------------- 433 ndmp = -1 ! type of damping in temperature and salinity 434 ! ='latitude', damping poleward of 'ndmp' degrees and function 435 ! of the distance-to-coast. Red and Med Seas as ndmp=-1 436 ! =-1 damping only in Med and Red Seas 437 ndmpf = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) 438 nmldmp = 1 ! type of damping: =0 damping throughout the water column 439 ! =1 no damping in the mixed layer defined by avt >5cm2/s ) 440 ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 441 sdmp = 50. ! surface time scale for internal damping (days) 442 bdmp = 360. ! bottom time scale for internal damping (days) 443 hdmp = 800. ! depth of transition between sdmp and bdmp (meters) 423 ! Coefficient 424 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 425 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 426 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 427 / 428 !----------------------------------------------------------------------- 429 &namtra_dmp ! tracer: T & S newtonian damping ('key_tradmp') 430 !----------------------------------------------------------------------- 431 nn_hdmp = -1 ! horizontal shape =-1, damping in Med and Red Seas only 432 ! =XX, damping poleward of XX degrees (XX>0) 433 ! + F(distance-to-coast) + Red and Med Seas 434 nn_zdmp = 1 ! vertical shape =0 damping throughout the water column 435 ! =1 no damping in the mixing layer (kz criteria) 436 ! =2 no damping in the mixed layer (rho crieria) 437 rn_surf = 50. ! surface time scale of damping [days] 438 rn_bot = 360. ! bottom time scale of damping [days] 439 rn_dep = 800. ! depth of transition between rn_surf and rn_bot [meters] 440 nn_file = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) 444 441 / 445 442 !!====================================================================== 446 443 !! *** Dynamics namelists *** 447 444 !!====================================================================== 448 !! nam_dynadv formulation of the momentum advection 449 !! nam_dynvor advection scheme 450 !! nam_dynhpg hydrostatic pressure gradient 451 !! namflg hydrostatic pressure gradient time stepping 452 !! nam_dynspg surface pressure gradient (CPP key only) 453 !! nam_dynldf lateral diffusion scheme 454 !!====================================================================== 455 456 !----------------------------------------------------------------------- 457 &nam_dynadv ! formulation of the momentum advection 445 !! namdyn_adv formulation of the momentum advection 446 !! namdyn_vor advection scheme 447 !! namdyn_hpg hydrostatic pressure gradient 448 !! namdyn_spg surface pressure gradient (CPP key only) 449 !! namdyn_ldf lateral diffusion scheme 450 !!====================================================================== 451 452 !----------------------------------------------------------------------- 453 &namdyn_adv ! formulation of the momentum advection 458 454 !----------------------------------------------------------------------- 459 455 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) … … 462 458 / 463 459 !----------------------------------------------------------------------- 464 &nam _dynvor ! option of physics/algorithm (not control by CPP keys)460 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 465 461 !----------------------------------------------------------------------- 466 462 ln_dynvor_ene = .false. ! enstrophy conserving scheme … … 470 466 / 471 467 !----------------------------------------------------------------------- 472 &nam _dynhpg ! Hydrostatic pressure gradient option468 &namdyn_hpg ! Hydrostatic pressure gradient option 473 469 !----------------------------------------------------------------------- 474 470 ln_hpg_zco = .false. ! z-coordinate - full steps … … 479 475 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 480 476 ln_hpg_rot = .false. ! s-coordinate (ROTated axes scheme) 481 gamm = 0.e0 ! weighting coefficient (wdj scheme) 482 / 483 !----------------------------------------------------------------------- 484 &namflg ! algorithm flags (algorithm not control by CPP keys) 485 !----------------------------------------------------------------------- 486 ln_dynhpg_imp = .false. ! hydrostatic pressure gradient: semi-implicit time scheme (T) 487 ! centered time scheme (F) 488 nn_dynhpg_rst = 0 ! add dynhpg implicit variables in restart ot not (1/0) 489 / 490 !----------------------------------------------------------------------- 491 !nam_dynspg ! surface pressure gradient (CPP key only) 477 rn_gamma = 0.e0 ! weighting coefficient (wdj scheme) 478 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 479 ! centered time scheme (F) 480 nn_dynhpg_rst = 0 ! =1 dynhpg restartable restart or not (=0) 481 / 482 !----------------------------------------------------------------------- 483 !namdyn_spg ! surface pressure gradient (CPP key only) 492 484 !----------------------------------------------------------------------- 493 485 ! ! explicit free surface ("key_dynspg_exp") … … 496 488 497 489 !----------------------------------------------------------------------- 498 &nam _dynldf ! lateral diffusion on momentum499 !----------------------------------------------------------------------- 500 !! Type of the operator :501 ln_dynldf_lap = .true. ! 502 ln_dynldf_bilap = .false. ! 503 !! Direction of action :504 ln_dynldf_level = .false. ! 505 ln_dynldf_hor = .true. ! horizontal (geopotential)(require "key_ldfslp" in s-coord.)506 ln_dynldf_iso = .false. ! iso-neutral(require "key_ldfslp")507 508 ahm0 = 40000. !horizontal eddy viscosity [m2/s]509 ahmb0 = 0. !background eddy viscosity for ldf_iso [m2/s]490 &namdyn_ldf ! lateral diffusion on momentum 491 !----------------------------------------------------------------------- 492 ! Type of the operator : 493 ln_dynldf_lap = .true. ! laplacian operator 494 ln_dynldf_bilap = .false. ! bilaplacian operator 495 ! Direction of action : 496 ln_dynldf_level = .false. ! iso-level 497 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 498 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") 499 ! Coefficient 500 rn_ahm_0 = 40000. ! horizontal eddy viscosity [m2/s] 501 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 510 502 / 511 503 !!====================================================================== 512 504 !! Tracers & Dynamics vertical physics namelists 513 505 !!====================================================================== 514 !! nam_zdf vertical physics 515 !! nam_npc non penetrative convection 516 !! nam_ric richardson number dependent vertical mixing ("key_zdfric" ) 517 !! nam_tke TKE dependent vertical mixing ("key_zdftke" ) 518 !! nam_kpp KPP dependent vertical mixing ("key_zdfkpp" ) 519 !! nam_ddm double diffusive mixing parameterization ("key_zdfddm" ) 520 !! nam_tmx tidal mixing parameterization ("key_zdftmx" ) 521 !!====================================================================== 522 523 !----------------------------------------------------------------------- 524 &nam_zdf ! vertical physics 525 !----------------------------------------------------------------------- 526 ! vertical eddy coef. or their background values 506 !! namzdf vertical physics 507 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric" ) 508 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke" ) 509 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp" ) 510 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm" ) 511 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx" ) 512 !!====================================================================== 513 514 !----------------------------------------------------------------------- 515 &namzdf ! vertical physics 516 !----------------------------------------------------------------------- 527 517 rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") 528 518 rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") 529 519 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 530 520 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) 531 ln_zdfevd = .true. ! convection:enhanced vertical diffusion (evd) (T) or not (F)532 nn_evdm = 1 ! 533 rn_avevd = 100. ! 534 ln_zdfnpc = .false. ! convection:Non-Penetrative algorithm (T) or not (F)535 nn_npc = 1 ! 536 nn_npcp = 365 ! 521 ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F) 522 nn_evdm = 1 ! evd apply on tracer (=0) or on tracer and momentum (=1) 523 rn_avevd = 100. ! evd mixing coefficient [m2/s] 524 ln_zdfnpc = .false. ! Non-Penetrative algorithm (T) or not (F) 525 nn_npc = 1 ! frequency of application of npc 526 nn_npcp = 365 ! npc control print frequency 537 527 ln_zdfexp = .false. ! time-stepping: split-explicit (T) or implicit (F) time stepping 538 528 nn_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T 539 529 / 540 530 !----------------------------------------------------------------------- 541 &nam _ric! richardson number dependent vertical diffusion ("key_zdfric" )542 !----------------------------------------------------------------------- 543 rn_avmri = 100.e-4 ! avm = rn_avmri / ( 1 + rn_alp * Ri**nn_ric )544 rn_alp = 5. ! avt = avm / ( 1 + rn_alp * Ri )545 nn_ric = 2 ! 546 / 547 !----------------------------------------------------------------------- 548 &nam _tke! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke")531 &namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" ) 532 !----------------------------------------------------------------------- 533 rn_avmri = 100.e-4 ! maximum value of the vertical viscosity 534 rn_alp = 5. ! coefficient of the parameterization 535 nn_ric = 2 ! coefficient of the parameterization 536 / 537 !----------------------------------------------------------------------- 538 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 549 539 !----------------------------------------------------------------------- 550 540 rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) … … 573 563 rn_lc = 0.15 ! coef. associated to Langmuir cells 574 564 / 575 !------------------------------------------------------------------------ 576 &nam _kpp! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionnally:565 !------------------------------------------------------------------------ 566 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionnally: 577 567 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 578 568 ln_kpprimix = .true. ! shear instability mixing … … 585 575 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv 586 576 nn_ave = 1 ! constant (=0) or profile (=1) background on avt 587 / 588 !----------------------------------------------------------------------- 589 &nam _ddm ! double diffusive mixing parameterization("key_zdfddm")577 / 578 !----------------------------------------------------------------------- 579 &namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm") 590 580 !----------------------------------------------------------------------- 591 581 rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity) … … 593 583 / 594 584 !----------------------------------------------------------------------- 595 &nam_tmx ! tidal mixing parameterization ("key_zdftmx") 596 !----------------------------------------------------------------------- 597 rn_htmx = 500. ! vertical decay scale for turbulence (meters) 598 rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1) 599 rn_tfe = 0.333 ! tidal dissipation efficiency 600 rn_me = 0.2 ! mixing efficiency 601 ln_tmx_itf = .FALSE. ! ITF specific parameterisation 602 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency 603 / 604 585 &namzdf_tmx ! tidal mixing parameterization ("key_zdftmx") 586 !----------------------------------------------------------------------- 587 rn_htmx = 500. ! vertical decay scale for turbulence (meters) 588 rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1) 589 rn_tfe = 0.333 ! tidal dissipation efficiency 590 rn_me = 0.2 ! mixing efficiency 591 ln_tmx_itf = .FALSE. ! ITF specific parameterisation 592 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency 593 / 605 594 !!====================================================================== 606 595 !! *** Miscelaneous namelists *** 607 596 !!====================================================================== 608 !! nam _mppMassively Parallel Processing ("key_mpp_mpi)609 !! nam _mpp_dyndistMassively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist")597 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 598 !! nammpp_dyndist Massively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist") 610 599 !! namctl Control prints & Benchmark 611 600 !! namsol elliptic solver / island / free surface … … 615 604 &namsol ! elliptic solver / island / free surface 616 605 !----------------------------------------------------------------------- 617 n solv = 1! elliptic solver: =1 preconditioned conjugate gradient (pcg)606 nn_solv = 2 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) 618 607 ! =2 successive-over-relaxation (sor) 619 nsol_arp = 0 ! absolute/relative (0/1) precision convergence test 620 nmin = 300 ! minimum of iterations for the SOR solver 621 nmax = 800 ! maximum of iterations for the SOR solver 622 nmod = 10 ! frequency of test for the SOR solver 623 eps = 1.e-6 ! absolute precision of the solver 624 resmax = 1.e-10 ! absolute precision for the SOR solver 625 sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 626 rnu = 1. ! strength of the additional force used in filtered free surface 627 / 628 !----------------------------------------------------------------------- 629 &nam_mpp ! Massively Parallel Processing ("key_mpp_mpi) 630 !----------------------------------------------------------------------- 631 c_mpi_send = 'S' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, 608 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test 609 rn_eps = 1.e-6 ! absolute precision of the solver 610 nn_nmin = 300 ! minimum of iterations for the SOR solver 611 nn_nmax = 800 ! maximum of iterations for the SOR solver 612 nn_nmod = 10 ! frequency of test for the SOR solver 613 rn_resmax = 1.e-10 ! absolute precision for the SOR solver 614 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 615 / 616 !----------------------------------------------------------------------- 617 &nammpp ! Massively Parallel Processing ("key_mpp_mpi) 618 !----------------------------------------------------------------------- 619 cn_mpi_send = 'S' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, 632 620 ! buffer blocking send or immediate non-blocking sends, resp. 633 nn_buffer = 0! size in bytes of exported buffer ('B' case), 0 no exportation634 / 635 !----------------------------------------------------------------------- 636 &nam _mpp_dyndist ! Massively Parallel Distribution("key_agrif" && "key_mpp_dyndist")637 !----------------------------------------------------------------------- 638 jpni = 1! jpni number of processors following i639 jpnj = 1! jpnj number of processors following j640 jpnij = 1! jpnij number of local domains641 / 642 !----------------------------------------------------------------------- 643 &namctl ! Control prints & Benchmark644 !----------------------------------------------------------------------- 645 ln_ctl = .false.! trends control print (expensive!)646 n print = 0! level of print (0 no extra print)647 n ictls = 0! start i indice of control sum (use to compare mono versus648 n ictle = 0! end i indice of control sum multi processor runs649 n jctls = 0! start j indice of control over a subdomain)650 n jctle = 0 ! end j indice of control651 isplt = 1! number of processors in i-direction652 jsplt = 1! number of processors in j-direction653 n bench = 0! Bench mode (1/0): CAUTION use zero except for bench621 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 622 / 623 !----------------------------------------------------------------------- 624 &nammpp_dyndist ! Massively Parallel Distribution for AGRIF zoom ("key_agrif" && "key_mpp_dyndist") 625 !----------------------------------------------------------------------- 626 jpni = 1 ! jpni number of processors following i 627 jpnj = 1 ! jpnj number of processors following j 628 jpnij = 1 ! jpnij number of local domains 629 / 630 !----------------------------------------------------------------------- 631 &namctl ! Control prints & Benchmark 632 !----------------------------------------------------------------------- 633 ln_ctl = .false. ! trends control print (expensive!) 634 nn_print = 0 ! level of print (0 no extra print) 635 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 636 nn_ictle = 0 ! end i indice of control sum multi processor runs 637 nn_jctls = 0 ! start j indice of control over a subdomain) 638 nn_jctle = 0 ! end j indice of control 639 nn_isplt = 1 ! number of processors in i-direction 640 nn_jsplt = 1 ! number of processors in j-direction 641 nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 654 642 ! (no physical validity of the results) 655 n bit_cmp = 0! bit comparison mode (1/0): CAUTION use zero except for test643 nn_bit_cmp = 0 ! bit comparison mode (1/0): CAUTION use zero except for test 656 644 ! of comparison between single and multiple processor runs 657 645 / 646 658 647 !!====================================================================== 659 648 !! *** Diagnostics namelists *** … … 666 655 667 656 !----------------------------------------------------------------------- 668 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 669 ! or mixed-layer trends ('key_trdmld') 670 ! or barotropic vorticity ("key_trdvor") 671 !----------------------------------------------------------------------- 672 ntrd = 365 ! time step frequency dynamics and tracers trends 673 nctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 674 ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 675 cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 676 cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 677 ln_trdmld_restart = .false. ! restart for ML diagnostics 678 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 657 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 658 ! ! or mixed-layer trends or barotropic vorticity ('key_trdmld' or "key_trdvor") 659 !----------------------------------------------------------------------- 660 nn_trd = 365 ! time step frequency dynamics and tracers trends 661 nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 662 rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 663 cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 664 cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 665 ln_trdmld_restart = .false. ! restart for ML diagnostics 666 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 679 667 / 680 668 !----------------------------------------------------------------------- 681 669 &namgap ! level mean model-data gap ('key_diagap') 682 670 !----------------------------------------------------------------------- 683 n gap= 15 ! time-step frequency of model-data gap computation684 n prg= 10 ! time-step frequency of gap print in model output671 nn_gap = 15 ! time-step frequency of model-data gap computation 672 nn_prg = 10 ! time-step frequency of gap print in model output 685 673 / 686 674 !----------------------------------------------------------------------- … … 688 676 !----------------------------------------------------------------------- 689 677 ln_rstflo = .false. ! float restart (T) or not (F) 690 n writefl= 75 ! frequency of writing in float output file691 n stockfl= 5475 ! frequency of creation of the float restart file678 nn_writefl= 75 ! frequency of writing in float output file 679 nn_stockfl= 5475 ! frequency of creation of the float restart file 692 680 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 693 681 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) … … 697 685 &namptr ! Poleward Transport Diagnostic 698 686 !----------------------------------------------------------------------- 699 ln_diaptr = .false. 700 ln_diaznl = .false. 701 ln_subbas = .false. 687 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 688 ln_diaznl = .false. ! Add zonal means and meridional stream functions 689 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 702 690 ! (orca configuration only, need input basins mask file named "subbasins.nc" 703 691 nf_ptr = 1 ! Frequency of ptr computation [time step]
Note: See TracChangeset
for help on using the changeset viewer.