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