| 1 | !! NEMO/OCE : Reference namelist_ref !! |
|---|
| 2 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 3 | !! NEMO/OCE : 1 - Domain & run manager (namrun, namcfg, namdom, namzgr, namzgr_sco ) |
|---|
| 4 | !! 2 - diagnostics (namnc4) |
|---|
| 5 | !! 3 - miscellaneous (nammpp, namctl) |
|---|
| 6 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 7 | !----------------------------------------------------------------------- |
|---|
| 8 | &namrun ! parameters of the run |
|---|
| 9 | !----------------------------------------------------------------------- |
|---|
| 10 | nn_no = 0 ! Assimilation cycle index |
|---|
| 11 | cn_exp = "ORCA2" ! experience name |
|---|
| 12 | nn_it000 = 1 ! first time step |
|---|
| 13 | nn_itend = 5840 ! last time step (std 5840) |
|---|
| 14 | nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) |
|---|
| 15 | nn_time0 = 0 ! initial time of day in hhmm |
|---|
| 16 | nn_leapy = 0 ! Leap year calendar (1) or not (0) |
|---|
| 17 | ln_rstart = .false. ! start from rest (F) or from a restart file (T) |
|---|
| 18 | nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T |
|---|
| 19 | nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T |
|---|
| 20 | ! ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist |
|---|
| 21 | ! ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart |
|---|
| 22 | ! ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart |
|---|
| 23 | cn_ocerst_in = "restart" ! suffix of ocean restart name (input) |
|---|
| 24 | cn_ocerst_indir = "." ! directory from which to read input ocean restarts |
|---|
| 25 | cn_ocerst_out = "restart" ! suffix of ocean restart name (output) |
|---|
| 26 | cn_ocerst_outdir = "." ! directory in which to write output ocean restarts |
|---|
| 27 | ln_iscpl = .false. ! cavity evolution forcing or coupling to ice sheet model |
|---|
| 28 | nn_istate = 0 ! output the initial state (1) or not (0) |
|---|
| 29 | ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) |
|---|
| 30 | nn_stock = 5840 ! frequency of creation of a restart file (modulo referenced to 1) |
|---|
| 31 | nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written |
|---|
| 32 | nn_write = 5475 ! frequency of write in the output file (modulo referenced to nn_it000) |
|---|
| 33 | ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) |
|---|
| 34 | ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard |
|---|
| 35 | ln_clobber = .true. ! clobber (overwrite) an existing file |
|---|
| 36 | nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) |
|---|
| 37 | / |
|---|
| 38 | !----------------------------------------------------------------------- |
|---|
| 39 | &namdom ! space and time domain (bathymetry, mesh, timestep) |
|---|
| 40 | !----------------------------------------------------------------------- |
|---|
| 41 | nn_bathy = 1 ! compute analyticaly (=0) or read (=1) the bathymetry file |
|---|
| 42 | ! or compute (2) from external bathymetry |
|---|
| 43 | nn_interp = 1 ! type of interpolation (nn_bathy =2) |
|---|
| 44 | cn_topo = 'bathymetry_ORCA12_V3.3.nc' ! external topo file (nn_bathy =2) |
|---|
| 45 | cn_bath = 'Bathymetry' ! topo name in file (nn_bathy =2) |
|---|
| 46 | cn_lon = 'nav_lon' ! lon name in file (nn_bathy =2) |
|---|
| 47 | cn_lat = 'nav_lat' ! lat name in file (nn_bathy =2) |
|---|
| 48 | rn_bathy = 0. ! value of the bathymetry. if (=0) bottom flat at jpkm1 |
|---|
| 49 | nn_msh = 0 ! create (=1) a mesh file or not (=0) |
|---|
| 50 | rn_hmin = -3. ! min depth of the ocean (>0) or min number of ocean level (<0) |
|---|
| 51 | rn_e3zps_min= 20. ! partial step thickness is set larger than the minimum of |
|---|
| 52 | rn_e3zps_rat= 0.1 ! rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 |
|---|
| 53 | ! |
|---|
| 54 | rn_rdt = 5760. ! time step for the dynamics (and tracer if nn_acc=0) |
|---|
| 55 | rn_atfp = 0.1 ! asselin time filter parameter |
|---|
| 56 | ln_crs = .false. ! Logical switch for coarsening module |
|---|
| 57 | jphgr_msh = 0 ! type of horizontal mesh |
|---|
| 58 | ! = 0 curvilinear coordinate on the sphere read in coordinate.nc |
|---|
| 59 | ! = 1 geographical mesh on the sphere with regular grid-spacing |
|---|
| 60 | ! = 2 f-plane with regular grid-spacing |
|---|
| 61 | ! = 3 beta-plane with regular grid-spacing |
|---|
| 62 | ! = 4 Mercator grid with T/U point at the equator |
|---|
| 63 | ppglam0 = 0.0 ! longitude of first raw and column T-point (jphgr_msh = 1) |
|---|
| 64 | ppgphi0 = -35.0 ! latitude of first raw and column T-point (jphgr_msh = 1) |
|---|
| 65 | ppe1_deg = 1.0 ! zonal grid-spacing (degrees) |
|---|
| 66 | ppe2_deg = 0.5 ! meridional grid-spacing (degrees) |
|---|
| 67 | ppe1_m = 5000.0 ! zonal grid-spacing (degrees) |
|---|
| 68 | ppe2_m = 5000.0 ! meridional grid-spacing (degrees) |
|---|
| 69 | ppsur = -4762.96143546300 ! ORCA r4, r2 and r05 coefficients |
|---|
| 70 | ppa0 = 255.58049070440 ! (default coefficients) |
|---|
| 71 | ppa1 = 245.58132232490 ! |
|---|
| 72 | ppkth = 21.43336197938 ! |
|---|
| 73 | ppacr = 3.0 ! |
|---|
| 74 | ppdzmin = 10. ! Minimum vertical spacing |
|---|
| 75 | pphmax = 5000. ! Maximum depth |
|---|
| 76 | ldbletanh = .TRUE. ! Use/do not use double tanf function for vertical coordinates |
|---|
| 77 | ppa2 = 100.760928500000 ! Double tanh function parameters |
|---|
| 78 | ppkth2 = 48.029893720000 ! |
|---|
| 79 | ppacr2 = 13.000000000000 ! |
|---|
| 80 | // |
|---|
| 81 | !----------------------------------------------------------------------- |
|---|
| 82 | &namcfg ! parameters of the configuration |
|---|
| 83 | !----------------------------------------------------------------------- |
|---|
| 84 | ! |
|---|
| 85 | ln_e3_dep = .true. ! =T : e3=dk[depth] in discret sens. |
|---|
| 86 | ! ! ===>>> will become the only possibility in v4.0 |
|---|
| 87 | ! ! =F : e3 analytical derivative of depth function |
|---|
| 88 | ! ! only there for backward compatibility test with v3.6 |
|---|
| 89 | ! ! |
|---|
| 90 | cp_cfg = "orca" ! name of the configuration |
|---|
| 91 | jp_cfg = 2 ! resolution of the configuration |
|---|
| 92 | jpidta = 182 ! 1st lateral dimension ( >= jpi ) |
|---|
| 93 | jpjdta = 149 ! 2nd " " ( >= jpj ) |
|---|
| 94 | jpkdta = 31 ! number of levels ( >= jpk ) |
|---|
| 95 | jpiglo = 182 ! 1st dimension of global domain --> i =jpidta |
|---|
| 96 | jpjglo = 149 ! 2nd - - --> j =jpjdta |
|---|
| 97 | jperio = 4 ! lateral cond. type (between 0 and 6) |
|---|
| 98 | ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present |
|---|
| 99 | ! in netcdf input files, as the start j-row for reading |
|---|
| 100 | ln_domclo = .false. ! computation of closed sea masks (see namclo) |
|---|
| 101 | / |
|---|
| 102 | !----------------------------------------------------------------------- |
|---|
| 103 | &namzgr ! vertical coordinate (default: NO selection) |
|---|
| 104 | !----------------------------------------------------------------------- |
|---|
| 105 | ln_zco = .false. ! z-coordinate - full steps |
|---|
| 106 | ln_zps = .false. ! z-coordinate - partial steps |
|---|
| 107 | ln_sco = .false. ! s- or hybrid z-s-coordinate |
|---|
| 108 | ln_isfcav = .false. ! ice shelf cavity (T: see namzgr_isf) |
|---|
| 109 | ln_linssh = .false. ! linear free surface |
|---|
| 110 | / |
|---|
| 111 | !----------------------------------------------------------------------- |
|---|
| 112 | &namzgr_isf ! isf cavity geometry definition (default: OFF) |
|---|
| 113 | !----------------------------------------------------------------------- |
|---|
| 114 | rn_isfdep_min = 10. ! minimum isf draft tickness (if lower, isf draft set to this value) |
|---|
| 115 | rn_glhw_min = 1.e-3 ! minimum water column thickness to define the grounding line |
|---|
| 116 | rn_isfhw_min = 10 ! minimum water column thickness in the cavity once the grounding line defined. |
|---|
| 117 | ln_isfchannel = .false. ! remove channel (based on 2d mask build from isfdraft-bathy) |
|---|
| 118 | ln_isfconnect = .false. ! force connection under the ice shelf (based on 2d mask build from isfdraft-bathy) |
|---|
| 119 | nn_kisfmax = 999 ! limiter in level on the previous condition. (if change larger than this number, get back to value before we enforce the connection) |
|---|
| 120 | rn_zisfmax = 7000. ! limiter in m on the previous condition. (if change larger than this number, get back to value before we enforce the connection) |
|---|
| 121 | ln_isfcheminey = .false. ! close cheminey |
|---|
| 122 | ln_isfsubgl = .false. ! remove subglacial lake created by the remapping process |
|---|
| 123 | rn_isfsubgllon = 0.0 ! longitude of the seed to determine the open ocean |
|---|
| 124 | rn_isfsubgllat = 0.0 ! latitude of the seed to determine the open ocean |
|---|
| 125 | / |
|---|
| 126 | !----------------------------------------------------------------------- |
|---|
| 127 | &namzgr_sco ! s-coordinate or hybrid z-s-coordinate (default: OFF) |
|---|
| 128 | !----------------------------------------------------------------------- |
|---|
| 129 | ln_s_sh94 = .false. ! Song & Haidvogel 1994 hybrid S-sigma (T)| |
|---|
| 130 | ln_s_sf12 = .false. ! Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied |
|---|
| 131 | ln_sigcrit = .false. ! use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch |
|---|
| 132 | ! stretching coefficients for all functions |
|---|
| 133 | rn_sbot_min = 10.0 ! minimum depth of s-bottom surface (>0) (m) |
|---|
| 134 | rn_sbot_max = 7000.0 ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) |
|---|
| 135 | rn_hc = 150.0 ! critical depth for transition to stretched coordinates |
|---|
| 136 | !!!!!!! Envelop bathymetry |
|---|
| 137 | rn_rmax = 0.3 ! maximum cut-off r-value allowed (0<r_max<1) |
|---|
| 138 | !!!!!!! SH94 stretching coefficients (ln_s_sh94 = .true.) |
|---|
| 139 | rn_theta = 6.0 ! surface control parameter (0<=theta<=20) |
|---|
| 140 | rn_bb = 0.8 ! stretching with SH94 s-sigma |
|---|
| 141 | !!!!!!! SF12 stretching coefficient (ln_s_sf12 = .true.) |
|---|
| 142 | rn_alpha = 4.4 ! stretching with SF12 s-sigma |
|---|
| 143 | rn_efold = 0.0 ! efold length scale for transition to stretched coord |
|---|
| 144 | rn_zs = 1.0 ! depth of surface grid box |
|---|
| 145 | ! bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b |
|---|
| 146 | rn_zb_a = 0.024 ! bathymetry scaling factor for calculating Zb |
|---|
| 147 | rn_zb_b = -0.2 ! offset for calculating Zb |
|---|
| 148 | !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] |
|---|
| 149 | rn_thetb = 1.0 ! bottom control parameter (0<=thetb<= 1) |
|---|
| 150 | / |
|---|
| 151 | !----------------------------------------------------------------------- |
|---|
| 152 | &namclo ! (closed sea : need ln_domclo = .true. in namcfg) (default: OFF) |
|---|
| 153 | !----------------------------------------------------------------------- |
|---|
| 154 | rn_lon_opnsea = 0.0 ! longitude seed of open ocean |
|---|
| 155 | rn_lat_opnsea = 0.0 ! latitude seed of open ocean |
|---|
| 156 | nn_closea = 8 ! number of closed seas ( = 0; only the open_sea mask will be computed) |
|---|
| 157 | ! |
|---|
| 158 | ! name ! lon_src ! lat_src ! lon_trg ! lat_trg ! river mouth area ! correction scheme ! radius trg ! id trg |
|---|
| 159 | ! ! (degree)! (degree)! (degree)! (degree)! local/coast/global ! (glo/rnf/emp) ! (m) ! |
|---|
| 160 | ! North American lakes |
|---|
| 161 | sn_lake(1) = 'superior' , -86.57 , 47.30 , -66.49 , 50.45 , 'local' , 'rnf' , 550000.0 , 2 |
|---|
| 162 | sn_lake(2) = 'michigan' , -87.06 , 42.74 , -66.49 , 50.45 , 'local' , 'rnf' , 550000.0 , 2 |
|---|
| 163 | sn_lake(3) = 'huron' , -82.51 , 44.74 , -66.49 , 50.45 , 'local' , 'rnf' , 550000.0 , 2 |
|---|
| 164 | sn_lake(4) = 'erie' , -81.13 , 42.25 , -66.49 , 50.45 , 'local' , 'rnf' , 550000.0 , 2 |
|---|
| 165 | sn_lake(5) = 'ontario' , -77.72 , 43.62 , -66.49 , 50.45 , 'local' , 'rnf' , 550000.0 , 2 |
|---|
| 166 | ! African Lake |
|---|
| 167 | sn_lake(6) = 'victoria' , 32.93 , -1.08 , 30.44 , 31.37 , 'coast' , 'emp' , 100000.0 , 3 |
|---|
| 168 | ! Asian Lakes |
|---|
| 169 | sn_lake(7) = 'caspian' , 50.0 , 44.0 , 0.0 , 0.0 , 'global' , 'glo' , 0.0 , 1 |
|---|
| 170 | sn_lake(8) = 'aral' , 60.0 , 45.0 , 0.0 , 0.0 , 'global' , 'glo' , 0.0 , 1 |
|---|
| 171 | / |
|---|
| 172 | !----------------------------------------------------------------------- |
|---|
| 173 | &namlbc ! lateral momentum boundary condition (default: NO selection) |
|---|
| 174 | !----------------------------------------------------------------------- |
|---|
| 175 | ! ! free slip ! partial slip ! no slip ! strong slip |
|---|
| 176 | rn_shlat = 0 ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat |
|---|
| 177 | ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical Eqs. |
|---|
| 178 | / |
|---|
| 179 | !----------------------------------------------------------------------- |
|---|
| 180 | &namagrif ! AGRIF zoom ("key_agrif") |
|---|
| 181 | !----------------------------------------------------------------------- |
|---|
| 182 | ln_spc_dyn = .true. ! use 0 as special value for dynamics |
|---|
| 183 | rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] |
|---|
| 184 | rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] |
|---|
| 185 | ln_chk_bathy = .FALSE. ! |
|---|
| 186 | npt_connect = 2 |
|---|
| 187 | npt_copy = 2 |
|---|
| 188 | ln_bry_south = .TRUE. |
|---|
| 189 | / |
|---|
| 190 | !----------------------------------------------------------------------- |
|---|
| 191 | &nam_tide ! tide parameters ("key_tide") |
|---|
| 192 | !----------------------------------------------------------------------- |
|---|
| 193 | ln_tide_pot = .true. ! use tidal potential forcing |
|---|
| 194 | ln_tide_ramp= .false. ! |
|---|
| 195 | rdttideramp = 0. ! |
|---|
| 196 | clname(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg |
|---|
| 197 | / |
|---|
| 198 | !----------------------------------------------------------------------- |
|---|
| 199 | &nambdy ! unstructured open boundaries (default: OFF) |
|---|
| 200 | !----------------------------------------------------------------------- |
|---|
| 201 | ln_bdy = .false. ! Use unstructured open boundaries |
|---|
| 202 | nb_bdy = 0 ! number of open boundary sets |
|---|
| 203 | ln_coords_file = .true. ! =T : read bdy coordinates from file |
|---|
| 204 | cn_coords_file = 'coordinates.bdy.nc' ! bdy coordinates files |
|---|
| 205 | ln_mask_file = .false. ! =T : read mask from file |
|---|
| 206 | cn_mask_file = '' ! name of mask file (if ln_mask_file=.TRUE.) |
|---|
| 207 | cn_dyn2d = 'none' ! |
|---|
| 208 | nn_dyn2d_dta = 0 ! = 0, bdy data are equal to the initial state |
|---|
| 209 | ! ! = 1, bdy data are read in 'bdydata .nc' files |
|---|
| 210 | ! ! = 2, use tidal harmonic forcing data from files |
|---|
| 211 | ! ! = 3, use external data AND tidal harmonic forcing |
|---|
| 212 | cn_dyn3d = 'none' ! |
|---|
| 213 | nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state |
|---|
| 214 | ! ! = 1, bdy data are read in 'bdydata .nc' files |
|---|
| 215 | cn_tra = 'none' ! |
|---|
| 216 | nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state |
|---|
| 217 | ! ! = 1, bdy data are read in 'bdydata .nc' files |
|---|
| 218 | cn_ice = 'none' ! |
|---|
| 219 | nn_ice_dta = 0 ! = 0, bdy data are equal to the initial state |
|---|
| 220 | ! ! = 1, bdy data are read in 'bdydata .nc' files |
|---|
| 221 | rn_ice_tem = 270. ! si3 only: arbitrary temperature of incoming sea ice |
|---|
| 222 | rn_ice_sal = 10. ! si3 only: -- salinity -- |
|---|
| 223 | rn_ice_age = 30. ! si3 only: -- age -- |
|---|
| 224 | / |
|---|
| 225 | !----------------------------------------------------------------------- |
|---|
| 226 | &nammpp ! Massively Parallel Processing ("key_mpp_mpi") |
|---|
| 227 | !----------------------------------------------------------------------- |
|---|
| 228 | cn_mpi_send = 'I' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, |
|---|
| 229 | ! ! buffer blocking send or immediate non-blocking sends, resp. |
|---|
| 230 | nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation |
|---|
| 231 | ln_nnogather = .true. ! activate code to avoid mpi_allgather use at the northfold |
|---|
| 232 | jpni = 1 ! jpni number of processors following i (set automatically if < 1) |
|---|
| 233 | jpnj = 1 ! jpnj number of processors following j (set automatically if < 1) |
|---|
| 234 | / |
|---|
| 235 | !----------------------------------------------------------------------- |
|---|
| 236 | &namctl ! Control prints (default: OFF) |
|---|
| 237 | !----------------------------------------------------------------------- |
|---|
| 238 | ln_ctl = .FALSE. ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T |
|---|
| 239 | sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following |
|---|
| 240 | sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. |
|---|
| 241 | sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure |
|---|
| 242 | sn_cfctl%l_oceout = .FALSE. ! that all areas report. |
|---|
| 243 | sn_cfctl%l_layout = .FALSE. ! |
|---|
| 244 | sn_cfctl%l_mppout = .FALSE. ! |
|---|
| 245 | sn_cfctl%l_mpptop = .FALSE. ! |
|---|
| 246 | sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] |
|---|
| 247 | sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000] |
|---|
| 248 | sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1] |
|---|
| 249 | sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info |
|---|
| 250 | nn_print = 0 ! level of print (0 no extra print) |
|---|
| 251 | nn_ictls = 0 ! start i indice of control sum (use to compare mono versus |
|---|
| 252 | nn_ictle = 0 ! end i indice of control sum multi processor runs |
|---|
| 253 | nn_jctls = 0 ! start j indice of control over a subdomain) |
|---|
| 254 | nn_jctle = 0 ! end j indice of control |
|---|
| 255 | nn_isplt = 1 ! number of processors in i-direction |
|---|
| 256 | nn_jsplt = 1 ! number of processors in j-direction |
|---|
| 257 | ln_timing = .false. ! timing by routine write out in timing.output file |
|---|
| 258 | ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii |
|---|
| 259 | / |
|---|
