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