Changeset 1113 for trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist
- Timestamp:
- 2008-06-14T16:58:05+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist
r1108 r1113 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/OPA : 1 - run manager (namrun , nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl)2 !! NEMO/OPA : 1 - run manager (namrun) 3 3 !! namelists 2 - Domain (nam_zgr, nam_zgr_sco, namdom) 4 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 6 !! 4 - lateral boundary (namlbc, namcla, namobc, nam bdy, namtide)6 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, namtide) 7 7 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 8 8 !! 6 - Tracer (nameos, nam_traadv, nam_traldf, namtdp) … … 10 10 !! 8 - Verical physics (namzdf, namnpc, namric, namtke, namkpp, namddm) 11 11 !! 9 - diagnostics (namtrd, namgap, namspr, namflo, namptr) 12 !! 9 - miscellaneous (namsol )12 !! 9 - miscellaneous (namsol, nam_mpp, nam_mpp_dyndist, namctl) 13 13 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 ! CAUTION: some scripts does not support CAPITALs for logical definition 15 ! ******* use .true. or .false. and NOT .TRUE. or .FALSE. 14 ! CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE. 16 15 17 16 !!====================================================================== … … 19 18 !!====================================================================== 20 19 !! namrun parameters of the run 21 !! nam_mpp Massively Parallel Processing 22 !! nam_mpp_dyndist 23 !! nam_ctl Control prints & Benchmark 24 !!====================================================================== 25 26 !----------------------------------------------------------------------- 27 ! namrun parameters of the run 28 !----------------------------------------------------------------------- 29 ! no job number 30 ! cexper experience name for vairmer format 31 ! ln_rstart boolean term for restart (true or false) 32 ! nrstdt restart control = 0 restart, do not control nit000 in the restart file. 33 ! ! = 1 restart, control nit000 in the restart file. Do not 34 ! ! use the date in the restart file (use ndate0 in namelist) 35 ! ! = 2 restart, control nit000 in the restart file, use the date 36 ! ! in the restart file. ndate0 in the namelist is ignored. 37 ! nit000 number of the first time step 38 ! nitend number of the last time step 39 ! ndate0 initial calendar date aammjj 40 ! nleapy Leap year calendar (0/1) 41 ! ninist initial state output flag (0/1) 42 ! nstock frequency of restart file 43 ! nwrite frequency of OUTPUT file 44 ! ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 45 &namrun 46 no = 0 47 cexper = "ORCA2" 48 ln_rstart = .false. 49 nrstdt = 0 50 nit000 = 1 51 nitend = 5475 52 ndate0 = 010101 53 nleapy = 0 54 ninist = 0 55 nstock = 5475 56 nwrite = 5475 57 ln_dimgnnn = .false. 58 / 59 !----------------------------------------------------------------------- 60 ! nam_mpp Massively Parallel Processing 61 !----------------------------------------------------------------------- 62 ! c_mpi_send mpi send/recieve type 63 ! = 'S' : standard blocking send 64 ! = 'B' : buffer blocking send 65 ! = 'I' : immediate non-blocking send 66 &nam_mpp 67 c_mpi_send = 'S' 68 / 69 !----------------------------------------------------------------------- 70 ! nam_mpp_dyndist Massively Parallel Distribution ("key_agrif" && "key_mpp_dyndist") 71 !----------------------------------------------------------------------- 72 ! jpni number of processors following i 73 ! jpnj number of processors following j 74 ! jpnij number of local domains 75 &nam_mpp_dyndist 76 jpni = 1 77 jpnj = 1 78 jpnij = 1 79 / 80 !----------------------------------------------------------------------- 81 ! nam_ctl Control prints & Benchmark 82 !----------------------------------------------------------------------- 83 ! ln_ctl trends control print (expensive!) 84 ! nprint level of print (0 no print) 85 ! nictls start i indice to make the control SUM (very usefull to compare mono- 86 ! nictle end i indice to make the control SUM (-versus multi processor runs) 87 ! njctls start j indice to make the control SUM (very usefull to compare mono- 88 ! njctle end j indice to make the control SUM (-versus multi processor runs) 89 ! nisplt number of processors following i 90 ! njsplt number of processors following j 91 ! nbench Bench parameter (0/1): CAUTION it must be zero except for bench 92 ! for which we don't care about physical meaning of the results 93 ! nbit_cmp bit comparison mode parameter (0/1): enables bit comparison between 94 ! single and multiple processor runs. 95 &namctl 96 ln_ctl = .false. 97 nprint = 0 98 nictls = 0 99 nictle = 0 100 njctls = 0 101 njctle = 0 102 isplt = 1 103 jsplt = 1 104 nbench = 0 105 nbit_cmp = 0 106 / 107 20 !!====================================================================== 21 22 !----------------------------------------------------------------------- 23 &namrun ! parameters of the run 24 !----------------------------------------------------------------------- 25 no = 0 ! job number 26 cexper = "ORCA2" ! experience name 27 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 28 nrstdt = 0 ! restart control = 0 nit000 is not compared to the restart file value 29 ! = 1 use ndate0 in namelist (not the value in the restart file) 30 ! = 2 calendar parameters read in the restart file 31 nit000 = 1 ! first time step 32 nitend = 5475 ! last time step 33 ndate0 = 010101 ! initial calendar date yymmdd (used if nrstdt=1) 34 nleapy = 0 ! Leap year calendar (1) or not (0) 35 ninist = 0 ! output the initial state (1) or not (0) 36 nstock = 5475 ! frequency of creation of a restart file 37 nwrite = 5475 ! frequency of write in the output file 38 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) 39 / 108 40 !!====================================================================== 109 41 !! *** Domain namelists *** … … 115 47 116 48 !----------------------------------------------------------------------- 117 ! nam_zgr vertical coordinate 118 !----------------------------------------------------------------------- 119 ! ln_zco z-coordinate - full steps (T/F) ("key_zco" may also be defined) 120 ! ln_zps z-coordinate - partial steps (T/F) 121 ! ln_sco s- or hybrid z-s-coordinate (T/F) 122 &nam_zgr 123 ln_zco = .false. 124 ln_zps = .true. 125 ln_sco = .false. 126 / 127 !----------------------------------------------------------------------- 128 ! nam_zgr_sco s-coordinate or hybrid z-s-coordinate 129 !----------------------------------------------------------------------- 130 ! sbot_min minimum depth of s-bottom surface (>0) (m) 131 ! sbot_max maximum depth of s-bottom surface (= ocean depth) (>0) (m) 132 ! theta surface control parameter (0<=theta<=20) 133 ! thetb bottom control parameter (0<=thetb<= 1) 134 ! r_max maximum cut-off r-value allowed (0<r_max<1) 135 &nam_zgr_sco 136 sbot_min = 300. 137 sbot_max = 5250. 138 theta = 6.0 139 thetb = 0.75 140 r_max = 0.15 141 / 142 !----------------------------------------------------------------------- 143 ! namdom space and time domain (bathymetry, mesh, timestep) 144 !----------------------------------------------------------------------- 145 ! ntopo = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 146 ! e3zps_min the thickness of the partial step is set larger than the 147 ! e3zps_rat the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 148 ! nmsh =1 create a mesh file (coordinates, scale factors, masks) 149 ! nacc acceleration of convergence method = 0, not used, rdt = rdttra 150 ! = 1, used, rdt < rdttra(k) 151 ! atfp asselin time filter parameter 152 ! rdt time step for the dynamics (and tracer if nacc=0) 153 ! rdtmin minimum time step on tracers 154 ! rdtmax maximum time step on tracers 155 ! rdth depth variation of tracer time step 156 ! rdtbt barotropic time step (for the time splitting algorithm) 157 ! nclosea =0 remove the closed sea from the global domain (orca configuration) 158 ! =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 159 &namdom 160 ntopo = 1 161 e3zps_min = 5. 162 e3zps_rat = 0.1 163 nmsh = 0 164 nacc = 0 165 atfp = 0.1 166 rdt = 5760. 167 rdtmin = 5760. 168 rdtmax = 5760. 169 rdth = 800. 170 rdtbt = 90. 171 nclosea = 0 172 / 173 49 &nam_zgr ! vertical coordinate 50 !----------------------------------------------------------------------- 51 ln_zco = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) 52 ln_zps = .true. ! z-coordinate - partial steps (T/F) 53 ln_sco = .false. ! s- or hybrid z-s-coordinate (T/F) 54 / 55 !----------------------------------------------------------------------- 56 &nam_zgr_sco ! s-coordinate or hybrid z-s-coordinate 57 !----------------------------------------------------------------------- 58 sbot_min = 300. ! minimum depth of s-bottom surface (>0) (m) 59 sbot_max = 5250. ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 60 theta = 6.0 ! surface control parameter (0<=theta<=20) 61 thetb = 0.75 ! bottom control parameter (0<=thetb<= 1) 62 r_max = 0.15 ! maximum cut-off r-value allowed (0<r_max<1) 63 / 64 !----------------------------------------------------------------------- 65 &namdom ! space and time domain (bathymetry, mesh, timestep) 66 !----------------------------------------------------------------------- 67 ntopo = 1 ! compute (=0) or read(=1) the bathymetry file 68 e3zps_min = 5. ! the thickness of the partial step is set larger than the minimum 69 e3zps_rat = 0.1 ! of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 70 nmsh = 0 ! create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 71 nacc = 0 ! =1 acceleration of convergence method used, rdt < rdttra(k) 72 ! =0, no acceleration, rdt = rdttra 73 atfp = 0.1 ! asselin time filter parameter 74 rdt = 5760. ! time step for the dynamics (and tracer if nacc=0) 75 rdtmin = 5760. ! minimum time step on tracers (used if nacc=1) 76 rdtmax = 5760. ! maximum time step on tracers (used if nacc=1) 77 rdth = 800. ! depth variation of tracer time step (used if nacc=1) 78 rdtbt = 90. ! barotropic time step (for the split explicit algorithm) ("key_dynspg_ts") 79 nclosea = 0 ! = 0 no closed sea in the model domain 80 ! = 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 81 / 174 82 !!====================================================================== 175 83 !! *** Surface Boundary Condition namelists *** … … 188 96 189 97 !----------------------------------------------------------------------- 190 ! namsbc Surface Boundary Condition (surface module) 191 !----------------------------------------------------------------------- 192 ! nn_fsbc frequency of surface boundary condition computation 193 ! (= the frequency of sea-ice model call) 194 ! ln_ana analytical formulation (fill namsbc_ana ) 195 ! ln_flx flux formulation (fill namsbc_flx ) 196 ! ln_blk_clio CLIO bulk formulation (fill namsbc_core) 197 ! ln_blk_core CORE bulk formulation (fill namsbc_clio) 198 ! ln_cpl Coupled formulation (fill namsbc_cpl ) 199 ! nn_ice =0 no ice boundary condition , 200 ! =1 observed ice-cover , 201 ! =2 ice-model used ("key_lim3" or "key_lim2) 202 ! nn_ico_cpl ice-ocean coupling : =0 LIM-3 old case 203 ! =1 stresses computed using now ocean velocity 204 ! =2 combination of 0 and 1 cases 205 ! ln_dm2dc Daily mean to Diurnal Cycle short wave (qsr) 206 ! ln_ssr Sea Surface Restoring on T and/or S (fill namsbc_ssr) 207 ! nn_fwb FreshWater Budget: =0 no control , 208 ! =1 annual global mean of e-p-r set to zero , 209 ! =2 global mean of e-p-r set to zero at each nn_fsbc time step 210 &namsbc 211 nn_fsbc = 5 212 ln_ana = .false. 213 ln_flx = .false. 214 ln_blk_clio = .true. 215 ln_blk_core = .false. 216 ln_cpl = .false. 217 nn_ice = 2 218 nn_ico_cpl = 0 219 ln_dm2dc = .false. 220 ln_rnf = .false. 221 ln_ssr = .false. 222 nn_fwb = 0 223 / 224 !----------------------------------------------------------------------- 225 ! namsbc_ana analytical surface boundary condition 226 !----------------------------------------------------------------------- 227 ! nn_tau000 gently increase the stress over the first ntau_rst time-steps 228 ! rn_utau0 uniform value for the i-stress 229 ! rn_vtau0 uniform value for the j-stress 230 ! rn_q0 uniform value for the total heat flux 231 ! rn_qsr0 uniform value for the solar radiation 232 ! rn_emp0 uniform value for the freswater budget (E-P) 233 &namsbc_ana 234 nn_tau000 = 0 235 rn_utau0 = 0.5 236 rn_vtau0 = 0.e0 237 rn_q0 = 0.e0 238 rn_qsr0 = 0.e0 239 rn_emp0 = 0.e0 240 / 241 !----------------------------------------------------------------------- 242 ! namsbc_flx surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 243 !----------------------------------------------------------------------- 244 ! cn_dir directory in which the model is executed 98 &namsbc ! Surface Boundary Condition (surface module) 99 !----------------------------------------------------------------------- 100 nn_fsbc = 5 ! frequency of surface boundary condition computation 101 ! (= the frequency of sea-ice model call) 102 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 103 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 104 ln_blk_clio = .true. ! CLIO bulk formulation (T => fill namsbc_clio) 105 ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) 106 ln_cpl = .false. ! Coupled formulation (T => fill namsbc_cpl ) 107 nn_ice = 2 ! =0 no ice boundary condition , 108 ! =1 use observed ice-cover , 109 ! =2 ice-model used ("key_lim3" or "key_lim2) 110 nn_ico_cpl = 0 ! ice-ocean coupling : =0 each nn_fsbc 111 ! =1 stresses recomputed each ocean time step ("key_lim3" only) 112 ! =2 combination of 0 and 1 cases ("key_lim3" only) 113 ln_dm2dc = .false. ! daily mean to diurnal cycle short wave (qsr) 114 ln_rnf = .false. ! runoffs (T => fill namsbc_ssr) 115 ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 116 nn_fwb = 0 ! FreshWater Budget: =0 unchecked , 117 ! =1 annual global mean of e-p-r set to zero , 118 ! =2 global mean of e-p-r set to zero at each nn_fsbc time step 119 / 120 !----------------------------------------------------------------------- 121 &namsbc_ana ! analytical surface boundary condition 122 !----------------------------------------------------------------------- 123 nn_tau000 = 0 ! gently increase the stress over the first ntau_rst time-steps 124 rn_utau0 = 0.5 ! uniform value for the i-stress 125 rn_vtau0 = 0.e0 ! uniform value for the j-stress 126 rn_q0 = 0.e0 ! uniform value for the total heat flux 127 rn_qsr0 = 0.e0 ! uniform value for the solar radiation 128 rn_emp0 = 0.e0 ! uniform value for the freswater budget (E-P) 129 / 130 !----------------------------------------------------------------------- 131 &namsbc_flx ! surface boundary condition : flux formulation 132 !----------------------------------------------------------------------- 133 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! starting ! 134 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 135 sn_utau = 'utau.nc' , 24. , 'utau' , .false. , 0 , 0 136 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .false. , 0 , 0 137 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .false. , 0 , 0 138 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .false. , 0 , 0 139 sn_emp = 'emp.nc' , 24. , 'emp' , .false. , 0 , 0 245 140 ! 246 ! THE ORDER OF THE FILES MATTER: 247 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 248 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 249 &namsbc_flx 250 cn_dir = './' 251 sn_utau = 'utau.nc' , 24. , 'utau' , .false. , 0 , 0 252 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .false. , 0 , 0 253 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .false. , 0 , 0 254 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .false. , 0 , 0 255 sn_emp = 'emp.nc' , 24. , 'emp' , .false. , 0 , 0 141 cn_dir = './' ! root directory for the location of the flux files 256 142 / 257 143 !----------------------------------------------------------------------- 258 ! namsbc_clio CLIO bulk formulea 259 !----------------------------------------------------------------------- 260 ! cn_dir directory in which the model is executed 144 &namsbc_clio ! namsbc_clio CLIO bulk formulea 145 !----------------------------------------------------------------------- 146 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! starting ! 147 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 148 sn_utau = 'taux_1m' , -12. , 'sozotaux' , .false. , 1 , 0 149 sn_vtau = 'tauy_1m' , -12. , 'sometauy' , .false. , 1 , 0 150 sn_wndm = 'flx' , -12. , 'socliowi' , .false. , 1 , 0 151 sn_tair = 'flx' , -12. , 'socliot2' , .false. , 1 , 0 152 sn_humi = 'flx' , -12. , 'socliohu' , .false. , 1 , 0 153 sn_ccov = 'flx' , -12. , 'socliocl' , .false. , 1 , 0 154 sn_prec = 'flx' , -12. , 'socliopl' , .false. , 1 , 0 261 155 ! 262 ! THE ORDER OF THE FILES MATTER: 263 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 264 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 265 &namsbc_clio 266 cn_dir = './' 267 sn_utau = 'taux_1m' , -12. , 'sozotaux', .false. , 1 , 0 268 sn_vtau = 'tauy_1m' , -12. , 'sometauy', .false. , 1 , 0 269 sn_wndm = 'flx' , -12. , 'socliowi', .false. , 1 , 0 270 sn_tair = 'flx' , -12. , 'socliot2', .false. , 1 , 0 271 sn_humi = 'flx' , -12. , 'socliohu', .false. , 1 , 0 272 sn_ccov = 'flx' , -12. , 'socliocl', .false. , 1 , 0 273 sn_prec = 'flx' , -12. , 'socliopl', .false. , 1 , 0 274 / 275 !----------------------------------------------------------------------- 276 ! namsbc_core CORE bulk formulea 277 !----------------------------------------------------------------------- 278 ! cn_dir directory in which the model is executed 279 ! ln_2m logical flag to use air temp. and hum referenced at 2m instead 10m 280 ! alpha_precip multiplication factor for precipitation (total & snow) 156 cn_dir = './' ! root directory for the location of the bulk files are 157 / 158 !----------------------------------------------------------------------- 159 &namsbc_core ! namsbc_core CORE bulk formulea 160 !----------------------------------------------------------------------- 161 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! starting ! 162 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 163 sn_wndi = 'u10' , 24. , 'U_10_MOD' , .false. , 1 , 0 164 sn_wndj = 'v10' , 24. , 'V_10_MOD' , .false. , 1 , 0 165 sn_qsr = 'rad' , 24. , 'SWDN_MOD' , .false. , 1 , 0 166 sn_qlw = 'rad' , 24. , 'LWDN_MOD' , .false. , 1 , 0 167 sn_tair = 't10' , 24. , 'T_10_MOD' , .false. , 1 , 0 168 sn_humi = 'q10' , 24. , 'Q_10_MOD' , .false. , 1 , 0 169 sn_prec = 'precip' , -12. , 'RAIN' , .false. , 1 , 0 170 sn_snow = 'precip' , -12. , 'SNOW' , .false. , 1 , 0 281 171 ! 282 ! THE ORDER OF THE FILES MATTER: 283 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 284 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 285 &namsbc_core 286 cn_dir = './' 287 ln_2m = .false. 288 alpha_precip = 1. 289 sn_wndi = 'u10' , 24. , 'U_10_MOD', .false. , 1 , 0 290 sn_wndj = 'v10' , 24. , 'V_10_MOD', .false. , 1 , 0 291 sn_qsr = 'rad' , 24. , 'SWDN_MOD', .false. , 1 , 0 292 sn_qlw = 'rad' , 24. , 'LWDN_MOD', .false. , 1 , 0 293 sn_tair = 't10' , 24. , 'T_10_MOD', .false. , 1 , 0 294 sn_humi = 'q10' , 24. , 'Q_10_MOD', .false. , 1 , 0 295 sn_prec = 'precip' , -12. , 'RAIN' , .false. , 1 , 0 296 sn_snow = 'precip' , -12. , 'SNOW' , .false. , 1 , 0 297 / 298 !----------------------------------------------------------------------- 299 ! namsbc_cpl coupled ocean/atmosphere model ("key_coupled") 300 !----------------------------------------------------------------------- 301 &namsbc_cpl 302 / 303 !----------------------------------------------------------------------- 304 ! namqsr penetrative solar radiation 305 !----------------------------------------------------------------------- 306 ! ln_traqsr : penetrative solar radiation (T) or not (F) (Default=T) 307 ! rabs fraction of qsr associated with xsi1 308 ! xsi1 first depth of extinction 309 ! xsi2 second depth of extinction 310 &namqsr 311 ln_traqsr = .true. 312 rabs = 0.58 313 xsi1 = 0.35 314 xsi2 = 23.0 315 / 316 !----------------------------------------------------------------------- 317 ! namsbc_rnf runoffs namelist surface boundary condition 318 !----------------------------------------------------------------------- 319 !cn_dir Root directory for location of ssr files 320 !nn_runoff =0 no, 1 runoff, 2 runoff+river mouth ups adv 321 !rn_hrnf runoffs, depth over which enhanced vertical mixing is used 322 !rn_avt_rnf runoffs, value of the additional vertical mixing coef. [m2/s] 323 !sn_rnf informations about the runoff file to be read 324 !sn_cnf informations about the runoff mouth file to be read 325 ! THE ORDER OF THE FILES MATTER: 326 ! ! file name ! frequency (hours) ! variable ! time interpolation ! 327 ! ! ! (if <0 months) ! name ! (logical) ! 328 &namsbc_rnf 329 cn_dir = './' 330 nn_runoff = 2 331 rn_hrnf = 0.e0 332 rn_avt_rnf = 1.e-3 333 sn_rnf = 'runoff_1m_nomask.nc' , -12. , 'sorunoff', .true. , 1 , 0 334 sn_cnf = 'runoff_1m_nomask.nc' , 0. , 'socoefr' , .false. , 1 , 0 335 / 336 !----------------------------------------------------------------------- 337 ! namsbc_ssr surface boundary condition : sea surface restoring 338 !----------------------------------------------------------------------- 339 ! cn_dir directory in which the model is executed 172 cn_dir = './' ! root directory for the location of the bulk files 173 ln_2m = .false. ! air temperature and humidity referenced at 2m (T) instead 10m (F) 174 alpha_precip= 1. ! multiplicative factor for precipitation (total & snow) 175 / 176 !----------------------------------------------------------------------- 177 &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") 178 !----------------------------------------------------------------------- 179 / 180 !----------------------------------------------------------------------- 181 &namqsr ! penetrative solar radiation 182 !----------------------------------------------------------------------- 183 ln_traqsr = .true. ! penetrative solar radiation (T) or not (F) 184 rabs = 0.58 ! fraction of qsr associated with xsi1 185 xsi1 = 0.35 ! first depth of extinction 186 xsi2 = 23.0 ! second depth of extinction 187 / 188 !----------------------------------------------------------------------- 189 &namsbc_rnf ! runoffs namelist surface boundary condition 190 !----------------------------------------------------------------------- 191 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! starting ! 192 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 193 sn_rnf = 'runoff_1m_nomask.nc', -12. , 'sorunoff' , .true. , 1 , 0 194 sn_cnf = 'runoff_1m_nomask.nc', 0. , 'socoefr' , .false. , 1 , 0 195 ! 196 cn_dir = './' ! root directory for the location of the runoff files 197 nn_runoff = 2 ! no runoff (0), runoff (1), runoff+river mouth ups adv with cen2 (2) 198 rn_hrnf = 0.e0 ! depth over which enhanced vertical mixing is used 199 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] 200 / 201 !----------------------------------------------------------------------- 202 &namsbc_ssr ! surface boundary condition : sea surface restoring 203 !----------------------------------------------------------------------- 204 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! starting ! 205 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 206 sn_sst = 'sst_data.nc' , 24. , 'sst' , .false. , 0 , 0 207 sn_sss = 'sss_data.nc' , -12. , 'sss' , .true. , 0 , 0 208 340 209 ! nn_sstr =0/1 add a retroaction term in the surface heat flux 341 210 ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux … … 344 213 ! dqdt magnitude of the retroaction on temperature [W/m2/K] 345 214 ! deds magnitude of the damping on salinity [mm/day] 346 ! THE ORDER OF THE FILES MATTER: 347 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 348 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 349 &namsbc_ssr 350 cn_dir = './' 351 nn_sstr = 0 352 nn_sssr = 0 353 dqdt = -40. 354 deds = 27.7 355 sn_sst = 'sst_data.nc' , 24. , 'sst' , .false. , 0 , 0 356 sn_sss = 'sss_data.nc' , -12. , 'sss' , .true. , 0 , 0 215 ! 216 cn_dir = './' ! root directory for the location of the runoff files 217 nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0) 218 nn_sssr = 0 ! add a damping term in the surface freshwater flux (=1) or not (=0) 219 dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 220 deds = 27.7 ! magnitude of the damping on salinity [mm/day/psu] 357 221 / 358 222 !----------------------------------------------------------------------- 359 ! namalb albedo parameters 360 !----------------------------------------------------------------------- 361 ! cgren correction of the snow or ice albedo to take into account 362 ! albice albedo of melting ice in the arctic and antarctic 363 ! alphd coefficients for linear interpolation used to compute albedo 364 ! between two extremes values (Pyane, 1972) 365 ! alphc " " 366 ! alphdi " " 367 &namalb 368 cgren = 0.06 369 albice = 0.5 370 alphd = 0.80 371 alphc = 0.65 372 alphdi = 0.72 373 / 374 223 &namalb ! albedo parameters 224 !----------------------------------------------------------------------- 225 cgren = 0.06 ! correction of the snow or ice albedo to take into account the 226 albice = 0.5 ! albedo of melting ice in the arctic and antarctic 227 alphd = 0.80 ! coefficients for linear interpolation used to 228 alphc = 0.65 ! compute albedo between two extremes values 229 alphdi = 0.72 ! (Pyane, 1972) 230 / 375 231 !!====================================================================== 376 232 !! *** Lateral boundary condition *** … … 379 235 !! namcla cross land advection 380 236 !! namobc open boundaries parameters ("key_obc") 237 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") 381 238 !! nambdy Unstructured open boundaries ("key_bdy") 382 239 !! namtide Tidal forcing at open boundaries ("key_bdy_tides") … … 384 241 385 242 !----------------------------------------------------------------------- 386 ! namlbc lateral momentum boundary condition 387 !----------------------------------------------------------------------- 388 ! shlat shlat = 0 , free slip ; 0 < shlat < 2 , partial slip 389 ! shlat = 2 , no slip ; 2 < shlat , strong slip 390 &namlbc 391 shlat = 2. 392 / 393 !----------------------------------------------------------------------- 394 ! namcla cross land advection 395 !----------------------------------------------------------------------- 396 ! n_cla advection between 2 ocean pts separates by land 397 &namcla 398 n_cla = 1 399 / 400 !----------------------------------------------------------------------- 401 ! namobc open boundaries parameters (#ifdef key_obc) 402 !----------------------------------------------------------------------- 403 ! nobc_dta = 0 the obc data are equal to the initial state 404 ! = 1 the obc data are read in 'obc .dta' files 405 ! rdpeob time relaxation (days) for the east open boundary 406 ! rdpwob time relaxation (days) for the west open boundary 407 ! rdpnob time relaxation (days) for the north open boundary 408 ! rdpsob time relaxation (days) for the south open boundary 409 ! zbsic1 barotropic stream function on isolated coastline 1 410 ! zbsic2 barotropic stream function on isolated coastline 2 411 ! zbsic3 barotropic stream function on isolated coastline 3 412 ! ln_obc_clim climatological obc data files (default T) 413 ! ln_vol_cst total volume conserved 414 &namobc 415 nobc_dta = 0 416 rdpein = 1. 417 rdpwin = 1. 418 rdpnin = 30. 419 rdpsin = 1. 420 rdpeob = 1500. 421 rdpwob = 15. 422 rdpnob = 150. 423 rdpsob = 15. 424 zbsic1 = 140.e+6 425 zbsic2 = 1.e+6 426 zbsic3 = 0. 427 ln_obc_clim = .true. 428 ln_vol_cst = .false. 429 / 430 !----------------------------------------------------------------------- 431 ! nambdy unstructured open boundaries parameters ("key_bdy") 432 !----------------------------------------------------------------------- 433 ! filbdy_mask Name of mask file (at T-points) 434 ! filbdy_data_T Name of data file (at T-points) 435 ! filbdy_data_U Name of data file (at U-points) 436 ! filbdy_data_V Name of data file (at V-points) 437 ! ln_bdy_clim bdy files must contain 1 or 12 time dumps and be cyclic. 438 ! ln_bdy_vol Total volume correction (see volbdy parameter) 439 ! ln_bdy_fla Flather boundary conditions 440 ! nbdy_dta 0 the bdy data are equal to the initial state 441 ! 1 the bdy data are read in 'bdydata .nc' files 442 ! nb_rimwidth width of the relaxation zone 443 ! volbdy 0 the total water flux across open boundaries is zero 444 ! 1 the total volume of the system is conserved 445 &nambdy 446 filbdy_mask = 'bdymask_grid_T.nc' 447 filbdy_data_T = 'bdydata_grid_T.nc' 448 filbdy_data_U = 'bdydata_grid_U.nc' 449 filbdy_data_V = 'bdydata_grid_V.nc' 450 ln_bdy_clim = .false. 451 ln_bdy_vol = .false. 452 ln_bdy_fla = .true. 453 nbdy_dta = 1 454 nb_rimwidth = 1 455 volbdy = 1 456 / 457 !----------------------------------------------------------------------- 458 ! namtide tidal forcing at unstructured boundaries ("key_bdy_tides") 459 !----------------------------------------------------------------------- 460 ! filtide File name root of tidal forcing files 461 ! tide_cpt Names of tidal components used 462 ! tide_speed Phase speeds of tidal components (deg/hour) 463 ! ln_tide_date If true, adjust tidal harmonics for start date of run. 464 &namtide 465 filtide = 'bdytide_' 466 tide_cpt = 'M2','S1' 467 tide_speed = 28.984106, 15.000001 468 ln_tide_date = .true. 469 / 470 243 &namlbc ! lateral momentum boundary condition 244 !----------------------------------------------------------------------- 245 shlat = 2. ! shlat = 0 : free slip 246 ! 0 < shlat < 2 : partial slip 247 ! shlat = 2 : no slip 248 ! 2 < shlat : strong slip 249 / 250 !----------------------------------------------------------------------- 251 &namcla ! cross land advection 252 !----------------------------------------------------------------------- 253 n_cla = 1 ! advection between 2 ocean pts separates by land 254 / 255 !----------------------------------------------------------------------- 256 &namobc ! open boundaries parameters ("key_obc") 257 !----------------------------------------------------------------------- 258 nobc_dta = 0 ! = 0 the obc data are equal to the initial state 259 ! = 1 the obc data are read in 'obc.dta' files 260 rdpein = 1. ! ??? 261 rdpwin = 1. ! ??? 262 rdpnin = 30. ! ??? 263 rdpsin = 1. ! ??? 264 rdpeob = 1500. ! time relaxation (days) for the east open boundary 265 rdpwob = 15. ! " " " west " 266 rdpnob = 150. ! " " " north " 267 rdpsob = 15. ! " " " south " 268 zbsic1 = 140.e+6 ! barotropic stream function on first isolated coastline 269 zbsic2 = 1.e+6 ! " " second " 270 zbsic3 = 0. ! " " thrid " 271 ln_obc_clim= .true. ! climatological obc data files (T) or not (F) 272 ln_vol_cst = .false. ! impose the total volume conservation (T) or not (F) 273 / 274 !----------------------------------------------------------------------- 275 &namagrif ! ("key_agrif") 276 !----------------------------------------------------------------------- 277 nbclineupdate = 3 ! baroclinic update frequency 278 ln_spc_dyn = .true. ! use 0 as special value for dynamics 279 visc_tra = 2880. ! viscosity coeeficient for tracers sponge layer 280 visc_dyn = 2880. ! viscosity coeeficient for dynamics sponge layer 281 / 282 !----------------------------------------------------------------------- 283 &nambdy ! unstructured open boundaries parameters ("key_bdy") 284 !----------------------------------------------------------------------- 285 filbdy_mask = 'bdymask_grid_T.nc' ! name of mask file (at T-points) 286 filbdy_data_T = 'bdydata_grid_T.nc' ! name of data file (at T-points) 287 filbdy_data_U = 'bdydata_grid_U.nc' ! name of data file (at U-points) 288 filbdy_data_V = 'bdydata_grid_V.nc' ! name of data file (at V-points) 289 ln_bdy_clim = .false. ! must contain 1 or 12 time dumps and be cyclic 290 ln_bdy_vol = .false. ! total volume correction (see volbdy parameter) 291 ln_bdy_fla = .true. ! flather boundary conditions 292 nbdy_dta = 1 ! = 0, bdy data are equal to the initial state 293 ! = 1, bdy data are read in 'bdydata .nc' files 294 nb_rimwidth = 1 ! width of the relaxation zone 295 volbdy = 1 ! = 0, the total water flux across open boundaries is zero 296 ! = 1, the total volume of the system is conserved 297 / 298 !----------------------------------------------------------------------- 299 &namtide ! tidal forcing at unstructured boundaries ("key_bdy_tides") 300 !----------------------------------------------------------------------- 301 filtide = 'bdytide_' ! file name root of tidal forcing files 302 tide_cpt = 'M2','S1' ! names of tidal components used 303 tide_speed = 28.984106, 15.000001 ! phase speeds of tidal components (deg/hour) 304 ln_tide_date = .true. ! adjust tidal harmonics for start date of run 305 / 471 306 !!====================================================================== 472 307 !! *** Bottom boundary condition *** … … 476 311 !! nambbl bottom boundary layer scheme ("key_trabbl_dif","key_trabbl_adv") 477 312 !!====================================================================== 478 !----------------------------------------------------------------------- 479 ! nambfr bottom friction 480 !----------------------------------------------------------------------- 481 ! nbotfr type of bottom friction : = 0 : no slip ; = 1 : linear friction 482 ! = 3 : free slip ; = 2 : non linear friction 483 ! bfri1 bottom drag coefficient (linear case) 484 ! bfri2 bottom drag coefficient (non linear case) 485 ! bfeb2 bottom turbulent kinetic energy (m2/s2) 486 &nambfr 487 nbotfr = 1 488 bfri1 = 4.e-4 489 bfri2 = 1.e-3 490 bfeb2 = 2.5e-3 491 / 492 !----------------------------------------------------------------------- 493 ! nambbc bottom temperature boundary condition 494 !----------------------------------------------------------------------- 495 ! ngeo_flux = 0 no geothermal heat flux 496 ! = 1 constant geothermal heat flux 497 ! = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 498 ! ngeo_flux_const Constant value of geothermal heat flux (W/m2) 499 &nambbc 500 ngeo_flux = 2 501 ngeo_flux_const = 86.4e-3 502 / 503 !----------------------------------------------------------------------- 504 ! nambbl bottom boundary layer scheme 505 !----------------------------------------------------------------------- 506 ! ! diffusive bbl ("key_trabbl") 507 ! ! advective bbl ("key_trabbl_adv") 508 ! atrbbl mixing coefficient of the bottom boundary layer scheme (m2/s) 509 &nambbl 510 atrbbl = 10000. 511 / 512 313 314 !----------------------------------------------------------------------- 315 &nambfr ! bottom friction 316 !----------------------------------------------------------------------- 317 nbotfr = 1 ! type of bottom friction : = 0 : no slip, = 2 : nonlinear friction 318 ! = 3 : free slip, = 1 : linear friction 319 bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 320 bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) 321 bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2) 322 / 323 !----------------------------------------------------------------------- 324 &nambbc ! bottom temperature boundary condition 325 !----------------------------------------------------------------------- 326 ngeo_flux = 2 ! geothermal heat flux = 0 no flux considered 327 ! = 1 constant flux 328 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 329 ngeo_flux_const = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 330 / 331 !----------------------------------------------------------------------- 332 &nambbl ! bottom boundary layer scheme 333 !----------------------------------------------------------------------- 334 ! ! diffusive bbl ("key_trabbl") 335 ! ! advective bbl ("key_trabbl_adv") 336 atrbbl = 10000. ! lateral mixing coefficient in the bbl [m2/s] 337 / 513 338 !!====================================================================== 514 339 !! Tracer (T & S ) namelists … … 521 346 522 347 !----------------------------------------------------------------------- 523 ! nameos ocean physical parameters 524 !----------------------------------------------------------------------- 525 ! neos type of equation of state and Brunt-Vaisala frequency 526 ! ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 527 ! ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 528 ! ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 529 ! ralpha thermal expension coefficient (neos= 1 or 2) 530 ! rbeta saline expension coefficient (neos= 2) 531 &nameos 532 neos = 0 533 ralpha = 2.e-4 534 rbeta = 0.001 535 / 536 !----------------------------------------------------------------------- 537 ! nam_traadv advection scheme for tracer 538 !----------------------------------------------------------------------- 539 ! ln_traadv_cen2 2nd order centered scheme (default T) 540 ! ln_traadv_tvd TVD scheme (default F) 541 ! ln_traadv_muscl MUSCL scheme (default F) 542 ! ln_traadv_muscl2 MUSCL scheme + cen2 at boundaries (default F) 543 ! ln_traadv_ubs UBS scheme (default F) 544 ! ln_traadv_qck QUICKEST scheme (default F) 545 &nam_traadv 546 ln_traadv_cen2 = .true. 547 ln_traadv_tvd = .false. 548 ln_traadv_muscl = .false. 549 ln_traadv_muscl2 = .false. 550 ln_traadv_ubs = .false. 551 ln_traadv_qck = .false. 552 / 553 !----------------------------------------------------------------------- 554 ! nam_traldf lateral diffusion scheme for tracer 555 !----------------------------------------------------------------------- 556 ! Type of the operator : 557 ! ln_traldf_lap laplacian operator (default T) 558 ! ln_traldf_bilap bilaplacian operator (default F) 559 ! Direction of action : 560 ! ln_traldf_level iso-level (default F) 561 ! ln_traldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" when ln_sco=T) 562 ! ln_traldf_iso iso-neutral (default T) (require "key_ldfslp") 563 ! Coefficient 564 ! aht0 horizontal eddy diffusivity for tracers (m2/s) 565 ! ahtb0 background eddy diffusivity for isopycnal diffusion (m2/s) 566 ! aeiv0 eddy induced velocity coefficient (m2/s) (require "key_traldf_eiv") 567 &nam_traldf 568 ln_traldf_lap = .true. 569 ln_traldf_bilap = .false. 570 ln_traldf_level = .false. 571 ln_traldf_hor = .false. 572 ln_traldf_iso = .true. 573 aht0 = 2000. 574 ahtb0 = 0. 575 aeiv0 = 2000. 576 / 577 !----------------------------------------------------------------------- 578 ! namtdp tracer newtonian damping ('key_tradmp') 579 !----------------------------------------------------------------------- 580 ! ndmp type of damping in temperature and salinity 581 ! ! ='latitude', damping poleward of 'ndmp' degrees and function 582 ! ! of the distance-to-coast. Red and Med Seas as ndmp=-1 583 ! ! =-1 damping only in Med and Red Seas 584 ! ndmpf =1 create a damping.coeff NetCDF file (the 3D damping array) 585 ! nmldmp type of damping in the mixed layer 586 ! ! =0 damping throughout the water column 587 ! ! =1 no damping in the mixed layer defined by avt >5cm2/s 588 ! ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 589 ! sdmp surface time scale for internal damping (days) 590 ! bdmp bottom time scale for internal damping (days) 591 ! hdmp depth of transition between sdmp and bdmp (meters) 592 &namtdp 593 ndmp = -1 594 ndmpf = 1 595 nmldmp = 1 596 sdmp = 50. 597 bdmp = 360. 598 hdmp = 800. 599 / 600 348 &nameos ! ocean physical parameters 349 !----------------------------------------------------------------------- 350 neos = 0 ! type of equation of state and Brunt-Vaisala frequency 351 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 352 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 353 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 354 ralpha = 2.e-4 ! thermal expension coefficient (neos= 1 or 2) 355 rbeta = 0.001 ! saline expension coefficient (neos= 2) 356 / 357 !----------------------------------------------------------------------- 358 &nam_traadv ! advection scheme for tracer 359 !----------------------------------------------------------------------- 360 ln_traadv_cen2 = .true. ! 2nd order centered scheme 361 ln_traadv_tvd = .false. ! TVD scheme 362 ln_traadv_muscl = .false. ! MUSCL scheme 363 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 364 ln_traadv_ubs = .false. ! UBS scheme 365 / 366 !----------------------------------------------------------------------- 367 &nam_traldf ! lateral diffusion scheme for tracer 368 !----------------------------------------------------------------------- 369 ! ! Type of the operator : 370 ln_traldf_lap = .true. ! laplacian operator 371 ln_traldf_bilap = .false. ! bilaplacian operator 372 ! Direction of action : 373 ln_traldf_level = .false. ! iso-level 374 ln_traldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 375 ln_traldf_iso = .true. ! iso-neutral (require "key_ldfslp") 376 ! ! Coefficient 377 aht0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 378 ahtb0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 379 aeiv0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 380 / 381 !----------------------------------------------------------------------- 382 &namtdp ! tracer newtonian damping ('key_tradmp') 383 !----------------------------------------------------------------------- 384 ndmp = -1 ! type of damping in temperature and salinity 385 ! ='latitude', damping poleward of 'ndmp' degrees and function 386 ! of the distance-to-coast. Red and Med Seas as ndmp=-1 387 ! =-1 damping only in Med and Red Seas 388 ndmpf = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) 389 nmldmp = 1 ! type of damping: =0 damping throughout the water column 390 ! =1 no damping in the mixed layer defined by avt >5cm2/s ) 391 ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 392 sdmp = 50. ! surface time scale for internal damping (days) 393 bdmp = 360. ! bottom time scale for internal damping (days) 394 hdmp = 800. ! depth of transition between sdmp and bdmp (meters) 395 / 601 396 !!====================================================================== 602 397 !! *** Dynamics namelists *** … … 611 406 612 407 !----------------------------------------------------------------------- 613 ! nam_dynadv formulation of the momentum advection 614 !----------------------------------------------------------------------- 615 ! ln_dynadv_vec vector form (T) or flux form (F) (default T) 616 ! ln_dynadv_cen2 flux form - 2nd order centered scheme (default T) 617 ! ln_dynadv_ubs flux form - 3rd order UBS scheme (default F) 618 &nam_dynadv 619 ln_dynadv_vec = .true. 620 ln_dynadv_cen2 = .false. 621 ln_dynadv_ubs = .false. 408 &nam_dynadv ! formulation of the momentum advection 409 !----------------------------------------------------------------------- 410 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) 411 ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme 412 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 622 413 / 623 414 !----------------------------------------------------------------------- 624 ! nam_dynvor option of physics/algorithm (not control by CPP keys) 625 !----------------------------------------------------------------------- 626 ! ln_dynvor_ens vorticity trends: enstrophy conserving scheme (default T) 627 ! ln_dynvor_ene " " : energy conserving scheme (default F) 628 ! ln_dynvor_mix " " : mixed scheme (default F) 629 ! ln_dynvor_een " " : energy & enstrophy scheme (default F) 630 &nam_dynvor 631 ln_dynvor_ene = .false. 632 ln_dynvor_ens = .true. 633 ln_dynvor_mix = .false. 634 ln_dynvor_een = .false. 635 / 636 !----------------------------------------------------------------------- 637 ! nam_dynhpg Hydrostatic pressure gradient option 638 !----------------------------------------------------------------------- 639 ! type of pressure gradient scheme (choose one only!) 640 ! ln_hpg_zco z-coordinate - full steps (default T) 641 ! ln_hpg_zps z-coordinate - partial steps (interpolation) 642 ! ln_hpg_sco s-coordinate (standard jacobian formulation) 643 ! ln_hpg_hel s-coordinate (helsinki modification) 644 ! ln_hpg_wdj s-coordinate (weighted density jacobian) 645 ! ln_hpg_djc s-coordinate (Density Jacobian with Cubic polynomial) 646 ! ln_hpg_rot s-coordinate (ROTated axes scheme) 647 ! parameters 648 ! gamm weighting coefficient (wdj scheme) 649 &nam_dynhpg 650 ln_hpg_zco = .false. 651 ln_hpg_zps = .true. 652 ln_hpg_sco = .false. 653 ln_hpg_hel = .false. 654 ln_hpg_wdj = .false. 655 ln_hpg_djc = .false. 656 ln_hpg_rot = .false. 657 gamm = 0.e0 658 / 659 !----------------------------------------------------------------------- 660 ! namflg algorithm flags (algorithm not control by CPP keys) 661 !----------------------------------------------------------------------- 662 ! ln_dynhpg_imp hydrostatic pressure gradient: semi-implicit time scheme (T) 663 ! centered time scheme (F) 664 ! nn_dynhpg_rst add dynhpg implicit variables in restart ot not (1/0) 665 &namflg 666 ln_dynhpg_imp = .false. 667 nn_dynhpg_rst = 0 668 / 669 !----------------------------------------------------------------------- 670 ! nam_dynspg surface pressure gradient (CPP key only) 671 !----------------------------------------------------------------------- 672 ! ! explicit free surface ("key_dynspg_exp") 673 ! ! filtered free surface ("key_dynspg_flt") 674 ! ! split-explicit free surface ("key_dynspg_ts") 675 ! ! rigid-lid ("key_dynspg_rl") 676 677 !----------------------------------------------------------------------- 678 ! nam_dynldf lateral diffusion on momentum 679 !----------------------------------------------------------------------- 680 ! Type of the operator : 681 ! ln_dynldf_lap laplacian operator (default T) 682 ! ln_dynldf_bilap bilaplacian operator (default F) 683 ! Direction of action : 684 ! ln_dynldf_level iso-level (default F) 685 ! ln_dynldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" if ln_sco=T) 686 ! ln_dynldf_iso iso-neutral (default T) (require "key_ldfslp") 687 ! Coefficient 688 ! ahm0 horizontal eddy viscosity for the dynamics (m2/s) 689 ! ahmb0 background eddy viscosity for isopycnal diffusion (m2/s) 690 &nam_dynldf 691 ln_dynldf_lap = .true. 692 ln_dynldf_bilap = .false. 693 ln_dynldf_level = .false. 694 ln_dynldf_hor = .true. 695 ln_dynldf_iso = .false. 696 ahm0 = 40000. 697 ahmb0 = 0. 698 / 699 415 &nam_dynvor ! option of physics/algorithm (not control by CPP keys) 416 !----------------------------------------------------------------------- 417 ln_dynvor_ene = .false. ! enstrophy conserving scheme 418 ln_dynvor_ens = .true. ! energy conserving scheme 419 ln_dynvor_mix = .false. ! mixed scheme 420 ln_dynvor_een = .false. ! energy & enstrophy scheme 421 / 422 !----------------------------------------------------------------------- 423 &nam_dynhpg ! Hydrostatic pressure gradient option 424 !----------------------------------------------------------------------- 425 ln_hpg_zco = .false. ! z-coordinate - full steps 426 ln_hpg_zps = .true. ! z-coordinate - partial steps (interpolation) 427 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) 428 ln_hpg_hel = .false. ! s-coordinate (helsinki modification) 429 ln_hpg_wdj = .false. ! s-coordinate (weighted density jacobian) 430 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 431 ln_hpg_rot = .false. ! s-coordinate (ROTated axes scheme) 432 gamm = 0.e0 ! weighting coefficient (wdj scheme) 433 / 434 !----------------------------------------------------------------------- 435 &namflg ! algorithm flags (algorithm not control by CPP keys) 436 !----------------------------------------------------------------------- 437 ln_dynhpg_imp = .false. ! hydrostatic pressure gradient: semi-implicit time scheme (T) 438 ! centered time scheme (F) 439 nn_dynhpg_rst = 0 ! add dynhpg implicit variables in restart ot not (1/0) 440 / 441 !----------------------------------------------------------------------- 442 !nam_dynspg ! surface pressure gradient (CPP key only) 443 !----------------------------------------------------------------------- 444 ! ! explicit free surface ("key_dynspg_exp") 445 ! ! filtered free surface ("key_dynspg_flt") 446 ! ! split-explicit free surface ("key_dynspg_ts") 447 ! ! rigid-lid ("key_dynspg_rl") 448 449 !----------------------------------------------------------------------- 450 &nam_dynldf ! lateral diffusion on momentum 451 !----------------------------------------------------------------------- 452 ! ! Type of the operator : 453 ln_dynldf_lap = .true. ! laplacian operator 454 ln_dynldf_bilap = .false. ! bilaplacian operator 455 ! ! Direction of action : 456 ln_dynldf_level = .false. ! iso-level 457 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 458 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") 459 ! Coefficient 460 ahm0 = 40000. ! horizontal eddy viscosity [m2/s] 461 ahmb0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 462 / 700 463 !!====================================================================== 701 464 !! Tracers & Dynamics vertical physics namelists … … 710 473 711 474 !----------------------------------------------------------------------- 712 ! namzdf vertical physics 713 !----------------------------------------------------------------------- 714 ! ln_zdfevd enhanced vertical diffusion (default T) 715 ! ln_zdfnpc Non-Penetrative Convection (default T) 716 ! avm0 vertical eddy viscosity for the dynamic (m2/s) (also background Kz if not "key_zdfcst") 717 ! avt0 vertical eddy diffusivity for tracers (m2/s) (also background Kz if not "key_zdfcst") 718 ! avevd vertical coefficient for enhanced diffusion scheme (m2/s) 719 ! n_evdm = 0 apply enhanced mixing on tracer only 720 ! = 1 apply enhanced mixing on both tracer and momentum 721 ! ln_zdfexp vertical physics: (=T) time splitting (T) (Default=F) 722 ! (=F) euler backward (F) 723 ! n_zdfexp number of sub-timestep for time splitting scheme 724 &namzdf 725 ln_zdfevd = .true. 726 ln_zdfnpc = .false. 727 avm0 = 1.2e-4 728 avt0 = 1.2e-5 729 avevd = 100. 730 n_evdm = 1 731 ln_zdfexp = .false. 732 n_zdfexp = 3 733 / 734 !----------------------------------------------------------------------- 735 ! namnpc non penetrative convection 736 !----------------------------------------------------------------------- 737 ! nnpc1 non penetrative convective scheme frequency 738 ! nnpc2 non penetrative convective scheme print frequency 739 &namnpc 740 nnpc1 = 1 741 nnpc2 = 365 742 / 743 !----------------------------------------------------------------------- 744 ! namric richardson number dependent vertical diffusion ("key_zdfric" ) 745 !----------------------------------------------------------------------- 746 ! avmri maximum value of the vertical viscosity 747 ! alp coefficient of the parameterization 748 ! nric coefficient of the parameterization 749 &namric 750 avmri = 100.e-4 751 alp = 5. 752 nric = 2 753 / 754 !----------------------------------------------------------------------- 755 ! namtke turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 756 !----------------------------------------------------------------------- 757 ! ln_rstke flag to restart with tke from a run without tke (default F) 758 ! rn_ediff coef. to compute vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 759 ! rn_ediss coef. of the Kolmogoroff dissipation 760 ! rn_ebb coef. of the surface input of tke 761 ! rn_efave coef. to applied to the tke diffusion ( avtke=rn_efave*avm ) 762 ! rn_emin minimum value of tke (m^2/s^2) 763 ! rn_emin0 surface minimum value of tke (m^2/s^2) 764 ! nn_itke number of restart iterative loops 765 ! rn_cri critic richardson number = 2/9 = 0.22222222 (hard coded) 766 ! nn_mxl length used = 0 bounded by the distance to surface and bottom 767 ! ! = 1 bounded by the local vertical scale factor 768 ! ! = 2 first vertical derivative of mixing length bounded by 1 769 ! ! = 3 same criteria as case 2 but applied in a different way 770 ! nn_pdl Prandtl number = 0 no vertical prandtl number (avt=avm) 771 ! ! = 1 prandtl number function of richarson number (avt=nn_pdl*avm) 772 ! nn_ave = horizontal averaged (=1) or not (=0) of avt (default =1) 773 ! nn_avb = 0 cst background avt0, avm0 / =1 profile used on avtb 774 ! ln_mxl0 mixing length scale surface value as function of wind stress or not 775 ! rn_lmin interior buoyancy lenght scale minimum value 776 ! rn_lmin0 surface buoyancy lenght scale minimum value 777 ! nn_etau to add tke induced by wind = 0 no add tke induced by wind 778 ! ! = 1 add tke induced by wind 779 ! ! = 2 add tke induced by wind only at the mxl base 780 ! nn_htau TKE penetration profile type 781 ! rn_efr fraction of TKE surface value which penetrates inside the thermocline 782 ! ln_lc flag to take into account Langmuir Circulation (LC) 783 ! rn_lc coef to compute vertical velocity of LC 784 ! nn_havtb horizontal shape or not for avtb (=0/1) 785 &namtke 786 ln_rstke = .false. 787 rn_ediff = 0.1 788 rn_ediss = 0.7 789 rn_ebb = 3.75 790 rn_efave = 1. 791 rn_emin = 1.e-6 792 rn_emin0 = 1.e-4 793 nn_itke = 50 794 nn_mxl = 2 795 nn_pdl = 1 796 nn_avb = 0 797 nn_ave = 1 798 ln_mxl0 = .false. 799 rn_lmin = 0.4 800 rn_lmin0 = 0.4 801 nn_etau = 0 802 nn_htau = 2 803 rn_efr = 0.05 804 ln_lc = .false. 805 rn_lc = 0.15 806 nn_havtb = 0 807 / 808 !----------------------------------------------------------------------- 809 ! namkpp K-Profile Parameterization dependent vertical mixing ("key_zdfkpp" ) 810 !----------------------------------------------------------------------- 811 ! ln_kpprimix shear instability mixing (default T) 812 ! difmiw constant internal wave viscosity (m2/s) 813 ! difsiw constant internal wave diffusivity (m2/s) 814 ! Riinfty local Richardson Number limit for shear instability 815 ! difri maximum shear mixing at Rig = 0 (m2/s) 816 ! bvsqcon Brunt-Vaisala squared (1/s**2) for maximum convection 817 ! difcon maximum mixing in interior convection (m2/s) 818 ! nave = 0/1 flag for horizontal average on avt, avmu, avmv 819 ! navb = 0/1 flag for constant or profile background avt 820 &namkpp 821 ln_kpprimix = .true. 822 difmiw = 1.e-04 823 difsiw = 0.1e-04 824 Riinfty = 0.8 825 difri = 0.0050 826 bvsqcon = -0.01e-07 827 difcon = 1. 828 navb = 0 829 nave = 1 830 / 831 !----------------------------------------------------------------------- 832 ! namddm double diffusive mixing parameterization ("key_zdfddm") 833 !----------------------------------------------------------------------- 834 ! avts maximum avs for dd mixing 835 ! hsbfr heat/salt buoyancy flux ratio 836 &namddm 837 avts = 1.e-4 838 hsbfr = 1.6 839 / 840 475 &namzdf ! vertical physics 476 !----------------------------------------------------------------------- 477 avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") 478 avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") 479 ln_zdfnpc = .false. ! convection: Non-Penetrative algorithm (T) or not (F) 480 ln_zdfevd = .true. ! convection: enhanced vertical diffusion (T) or not (F) 481 avevd = 100. ! vertical coefficient for enhanced diffusion scheme [m2/s] 482 n_evdm = 1 ! enhanced mixing apply on tracer (=0) or on tracer and momentum (=1) 483 ln_zdfexp = .false. ! split explicit (T) or implicit (F) time stepping 484 n_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T 485 / 486 !----------------------------------------------------------------------- 487 &namnpc ! non penetrative convection 488 !----------------------------------------------------------------------- 489 nnpc1 = 1 ! non penetrative convective scheme computation frequency 490 nnpc2 = 365 ! non penetrative convective scheme print frequency 491 / 492 !----------------------------------------------------------------------- 493 &namric ! richardson number dependent vertical diffusion ("key_zdfric" ) 494 !----------------------------------------------------------------------- 495 avmri = 100.e-4 ! maximum value of the vertical viscosity 496 alp = 5. ! coefficient of the parameterization 497 nric = 2 ! coefficient of the parameterization 498 / 499 !----------------------------------------------------------------------- 500 &namtke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 501 !----------------------------------------------------------------------- 502 ln_rstke = .false. ! restart with tke from a run without tke (T) or not (F) 503 nn_itke = 50 ! number of iterative loops if ln_rstke=T 504 rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 505 rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation 506 rn_ebb = 3.75 ! coef. of the surface input of tke 507 rn_efave = 1. ! boost of the tke diffusion ( avtke=rn_efave*avm ) 508 rn_emin = 1.e-6 ! minimum value of tke [m2/s2] 509 rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2] 510 nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom 511 ! = 1 bounded by the local vertical scale factor 512 ! = 2 first vertical derivative of mixing length bounded by 1 513 ! = 3 same criteria as case 2 but applied in a different way 514 nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 515 nn_avb = 0 ! profile for constant background used on avt & avm (=1) or not (=0) 516 nn_ave = 1 ! horizontal averaged on avt (=1) or not (=0) 517 ln_mxl0 = .false. ! mixing length scale surface value as function of wind stress (T) or not (F) 518 rn_lmin = 0.4 ! interior buoyancy lenght scale minimum value 519 rn_lmin0 = 0.4 ! surface buoyancy lenght scale minimum value 520 nn_etau = 0 ! exponentially deceasing penetration of tke due to internal & intertial waves 521 ! = 0 no penetration ( O(2 km) resolution) 522 ! = 1 additional tke source 523 ! = 2 additional tke source applied only at the base of the mixed layer 524 nn_htau = 2 ! type of exponential decrease of tke penetration 525 ! = 0 constant 10 m length scale 526 ! = 1 ??? 527 ! = 2 ??? 528 rn_efr = 0.05 ! fraction of surface tke value which penetrates inside the ocean 529 ln_lc = .false. ! Langmuir cell effect 530 rn_lc = 0.15 ! coef. associated to Langmuir cells 531 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) 532 / 533 !------------------------------------------------------------------------ 534 &namkpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionnally: 535 ! "key_kppcustom" or "key_kpplktb") 536 !------------------------------------------------------------------------ 537 ln_kpprimix = .true. ! shear instability mixing 538 difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s] 539 difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s] 540 Riinfty = 0.8 ! local Richardson Number limit for shear instability 541 difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s] 542 bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 543 difcon = 1. ! maximum mixing in interior convection [m2/s] 544 navb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv 545 nave = 1 ! constant (=0) or profile (=1) background on avt 546 / 547 !----------------------------------------------------------------------- 548 &namddm ! double diffusive mixing parameterization ("key_zdfddm") 549 !----------------------------------------------------------------------- 550 avts = 1.e-4 ! maximum avs (vertical mixing on salinity) 551 hsbfr = 1.6 ! heat/salt buoyancy flux ratio 552 / 841 553 !!====================================================================== 842 554 !! *** Miscelaneous namelists *** 843 555 !!====================================================================== 844 !! namsol elliptic solver / island / free surface 845 !!====================================================================== 846 847 !----------------------------------------------------------------------- 848 ! namsol elliptic solver / island / free surface 849 !----------------------------------------------------------------------- 850 ! nsolv elliptic solver (=1 preconditioned conjugate gradient: pcg) 851 ! (=2 successive-over-relaxation: sor) 852 ! (=3 FETI currently it does not work! ("key_feti") 853 ! nsol_arp absolute/relative (0/1) precision convergence test 854 ! nmin minimum of iterations for the SOR solver 855 ! nmax maximum of iterations for the SOR solver 856 ! nmod frequency of test for the SOR solver 857 ! eps absolute precision of the solver 858 ! resmax absolute precision for the SOR solver 859 ! sor optimal coefficient for SOR solver 860 ! epsisl absolute precision on stream function solver 861 ! nmisl maximum pcg iterations for island 862 ! rnu strength of the additional force used in free surface b.c. 863 &namsol 864 nsolv = 1 865 nsol_arp = 0 866 nmin = 300 867 nmax = 800 868 nmod = 10 869 eps = 1.E-6 870 resmax = 1.E-10 871 sor = 1.92 872 epsisl = 1.e-10 873 nmisl = 4000 874 rnu = 1. 556 !! nam_mpp Massively Parallel Processing ("key_mpp_mpi) 557 !! nam_mpp_dyndist Massively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist") 558 !! namctl Control prints & Benchmark 559 !! namsol elliptic solver / island / free surface 560 !!====================================================================== 561 562 !----------------------------------------------------------------------- 563 &namsol ! elliptic solver / island / free surface 564 !----------------------------------------------------------------------- 565 nsolv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) 566 ! =2 successive-over-relaxation (sor) 567 ! =3 FETI (fet) ("key_feti") 568 ! =4 sor with extra outer halo 569 nsol_arp = 0 ! absolute/relative (0/1) precision convergence test 570 nmin = 300 ! minimum of iterations for the SOR solver 571 nmax = 800 ! maximum of iterations for the SOR solver 572 nmod = 10 ! frequency of test for the SOR solver 573 eps = 1.e-6 ! absolute precision of the solver 574 resmax = 1.e-10 ! absolute precision for the SOR solver 575 sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 576 epsisl = 1.e-10 ! absolute precision on stream function solver 577 nmisl = 4000 ! maximum pcg iterations for island ("key_islands") 578 rnu = 1. ! strength of the additional force used in filtered free surface 579 / 580 !----------------------------------------------------------------------- 581 &nam_mpp ! Massively Parallel Processing ("key_mpp_mpi) 582 !----------------------------------------------------------------------- 583 c_mpi_send = 'S' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, 584 ! buffer blocking send or immediate non-blocking sends, resp. 585 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 586 / 587 !----------------------------------------------------------------------- 588 &nam_mpp_dyndist ! Massively Parallel Distribution ("key_agrif" && "key_mpp_dyndist") 589 !----------------------------------------------------------------------- 590 jpni = 1 ! jpni number of processors following i 591 jpnj = 1 ! jpnj number of processors following j 592 jpnij = 1 ! jpnij number of local domains 593 / 594 !----------------------------------------------------------------------- 595 &namctl ! Control prints & Benchmark 596 !----------------------------------------------------------------------- 597 ln_ctl = .false. ! trends control print (expensive!) 598 nprint = 0 ! level of print (0 no extra print) 599 nictls = 0 ! start i indice of control sum (use to compare mono versus 600 nictle = 0 ! end i indice of control sum multi processor runs 601 njctls = 0 ! start j indice of control over a subdomain) 602 njctle = 0 ! end j indice of control 603 isplt = 1 ! number of processors in i-direction 604 jsplt = 1 ! number of processors in j-direction 605 nbench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 606 ! (no physical validity of the results) 607 nbit_cmp = 0 ! bit comparison mode (1/0): CAUTION use zero except for test 608 ! of comparison between single and multiple processor runs 875 609 / 876 610 !!====================================================================== … … 885 619 886 620 !----------------------------------------------------------------------- 887 ! namtrd diagnostics on dynamics and/or tracer trends("key_trddyn" and/or "key_trdtra")621 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 888 622 ! or mixed-layer trends ('key_trdmld') 889 623 ! or barotropic vorticity ("key_trdvor") 890 624 !----------------------------------------------------------------------- 891 ! ntrd time step frequency dynamics and tracers trends 892 ! nctls control surface type in mixed-layer trends (0,1 or n<jpk) 893 ! ln_trdmld_restart restart for ML diagnostics 894 ! ucf unit conversion factor (=1 -> /seconds | =86400. -> /day) 895 ! ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 896 &namtrd 897 ntrd = 365 898 nctls = 0 899 ln_trdmld_restart = .false. 900 ucf = 1. 901 ln_trdmld_instant = .false. 902 / 903 !----------------------------------------------------------------------- 904 ! namgap level mean model-data gap ('key_diagap') 905 !----------------------------------------------------------------------- 906 ! ngap time-step frequency of model-data gap computation 907 ! nprg time-step frequency of gap print in model output 908 &namgap 909 ngap = 15 910 nprg = 10 911 / 912 !----------------------------------------------------------------------- 913 ! namspr surface pressure diagnostic 914 !----------------------------------------------------------------------- 915 ! nmaxp maximum of iterations for the solver 916 ! epsp absolute precision of the solver 917 ! niterp number of iteration done by the solver 918 &namspr 919 nmaxp = 1000 920 epsp = 1.e-3 921 niterp = 400 922 / 923 !----------------------------------------------------------------------- 924 ! namflo float parameters ("key_float") 925 !----------------------------------------------------------------------- 926 ! ln_rstflo boolean term for float restart (true or false) 927 ! nwritefl frequency of float output file 928 ! nstockfl frequency of float restart file 929 ! ln_argo Argo type floats (stay at the surface each 10 days) 930 ! ln_flork4 = T trajectories computed with a 4th order Runge-Kutta 931 ! = F (default) computed with Blanke' scheme 932 &namflo 933 ln_rstflo = .false. 934 nwritefl = 75 935 nstockfl = 5475 936 ln_argo = .false. 937 ln_flork4 = .false. 938 / 939 !----------------------------------------------------------------------- 940 ! namptr Poleward Transport Diagnostic 941 !----------------------------------------------------------------------- 942 ! ln_diaptr logical flag for Poleward transport computation 943 ! ln_subbas logical flag for Atlantic/Pacific/Indian basins computation 944 ! need input basins mask file named "subbasins.nc" 945 ! nf_ptr Frequency of computation 946 &namptr 947 ln_diaptr = .false. 948 ln_subbas = .false. 949 nf_ptr = 15 950 / 625 ntrd = 365 ! time step frequency dynamics and tracers trends 626 nctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 627 ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 628 ln_trdmld_restart = .false. ! restart for ML diagnostics 629 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 630 / 631 !----------------------------------------------------------------------- 632 &namgap ! level mean model-data gap ('key_diagap') 633 !----------------------------------------------------------------------- 634 ngap = 15 ! time-step frequency of model-data gap computation 635 nprg = 10 ! time-step frequency of gap print in model output 636 / 637 !----------------------------------------------------------------------- 638 &namspr ! surface pressure diagnostic 639 !----------------------------------------------------------------------- 640 nmaxp = 1000 ! maximum of iterations for the solver 641 epsp = 1.e-3 ! absolute precision of the solver 642 niterp = 400 ! number of iteration done by the solver 643 / 644 !----------------------------------------------------------------------- 645 &namflo ! float parameters ("key_float") 646 !----------------------------------------------------------------------- 647 ln_rstflo = .false. ! float restart (T) or not (F) 648 nwritefl = 75 ! frequency of writing in float output file 649 nstockfl = 5475 ! frequency of creation of the float restart file 650 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 651 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 652 ! or computed with Blanke' scheme (F) 653 / 654 !----------------------------------------------------------------------- 655 &namptr ! Poleward Transport Diagnostic 656 !----------------------------------------------------------------------- 657 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 658 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 659 ! (orca configuration only, need input basins mask file named "subbasins.nc" 660 nf_ptr = 15 ! Frequency of ptr computation [time step] 661 /
Note: See TracChangeset
for help on using the changeset viewer.