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