| 1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 2 | !! NEMO/OPA : 1D configuration based on Kato-Philipps () exp. |
|---|
| 3 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 4 | |
|---|
| 5 | !!====================================================================== |
|---|
| 6 | !! *** Domain & Run management namelists *** !! |
|---|
| 7 | !! !! |
|---|
| 8 | !! namrun parameters of the run |
|---|
| 9 | !! namdom space and time domain |
|---|
| 10 | !! namcfg parameters of the configuration (default: user defined GYRE) |
|---|
| 11 | !! namwad Wetting and drying (default NO) |
|---|
| 12 | !! namtsd data: temperature & salinity (default NO) |
|---|
| 13 | !! namcrs coarsened grid (for outputs and/or TOP) (ln_crs =T) |
|---|
| 14 | !! namc1d 1D configuration options ("key_c1d") |
|---|
| 15 | !! namc1d_dyndmp 1D newtonian damping applied on currents ("key_c1d") |
|---|
| 16 | !! namc1d_uvd 1D data (currents) ("key_c1d") |
|---|
| 17 | !!====================================================================== |
|---|
| 18 | ! |
|---|
| 19 | !----------------------------------------------------------------------- |
|---|
| 20 | &namrun ! parameters of the run |
|---|
| 21 | !----------------------------------------------------------------------- |
|---|
| 22 | cn_exp = "PAPA" ! experience name |
|---|
| 23 | nn_it000 = 1 ! first time step |
|---|
| 24 | nn_itend = 87600 ! last time step |
|---|
| 25 | nn_date0 = 20100615 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) |
|---|
| 26 | nn_leapy = 1 ! Leap year calendar (1) or not (0) |
|---|
| 27 | nn_stock = 87600 ! frequency of creation of a restart file (modulo referenced to 1) |
|---|
| 28 | nn_write = 240 ! frequency of write in the output file (modulo referenced to nn_it000) |
|---|
| 29 | / |
|---|
| 30 | !----------------------------------------------------------------------- |
|---|
| 31 | &namdom ! space and time domain (bathymetry, mesh, timestep) |
|---|
| 32 | !----------------------------------------------------------------------- |
|---|
| 33 | nn_bathy = 0 ! compute (=0) or read (=1) the bathymetry file |
|---|
| 34 | rn_bathy = 4198. ! value of the bathymetry. if (=0) bottom flat at jpkm1 |
|---|
| 35 | ! |
|---|
| 36 | rn_rdt = 360. ! time step for the dynamics |
|---|
| 37 | ppglam0 = -150.0 ! longitude of first raw and column T-point (jphgr_msh = 1) |
|---|
| 38 | ppgphi0 = 50.0 ! latitude of first raw and column T-point (jphgr_msh = 1) |
|---|
| 39 | ppe1_deg = 0.1 ! zonal grid-spacing (degrees) |
|---|
| 40 | ppe2_deg = 0.1 ! meridional grid-spacing (degrees) |
|---|
| 41 | ppe1_m = 999999.0 ! zonal grid-spacing (degrees) |
|---|
| 42 | ppe2_m = 999999.0 ! meridional grid-spacing (degrees) |
|---|
| 43 | ppsur = -3958.951371276829 ! ORCA r4, r2 and r05 coefficients |
|---|
| 44 | ppa0 = 103.9530096000000 ! (default coefficients) |
|---|
| 45 | ppa1 = 2.415951269000000 ! |
|---|
| 46 | ppkth = 15.35101370000000 ! |
|---|
| 47 | ppacr = 7.0 ! |
|---|
| 48 | ppdzmin = 999999.0 ! Minimum vertical spacing |
|---|
| 49 | pphmax = 999999.0 ! Maximum depth |
|---|
| 50 | ldbletanh = .TRUE. ! Use/do not use double tanf function for vertical coordinates |
|---|
| 51 | ppa2 = 100.7609285000000 ! Double tanh function parameters |
|---|
| 52 | ppkth2 = 48.02989372000000 ! |
|---|
| 53 | ppacr2 = 13.00000000000 ! |
|---|
| 54 | / |
|---|
| 55 | !----------------------------------------------------------------------- |
|---|
| 56 | &namcfg ! parameters of the configuration |
|---|
| 57 | !----------------------------------------------------------------------- |
|---|
| 58 | cp_cfg = "papa" ! name of the configuration |
|---|
| 59 | jp_cfg = 1 ! resolution of the configuration |
|---|
| 60 | ! jpidta = 3 ! 1st lateral dimension ( >= jpi ) = 30*jp_cfg+2 |
|---|
| 61 | ! jpjdta = 3 ! 2nd " " ( >= jpj ) = 20*jp_cfg+2 |
|---|
| 62 | ! jpkdta = 75 ! number of levels ( >= jpk ) |
|---|
| 63 | ! jpiglo = 3 ! 1st dimension of global domain --> i = jpidta |
|---|
| 64 | ! jpjglo = 3 ! 2nd - - --> j = jpjdta |
|---|
| 65 | jperio = 0 ! lateral cond. type (between 0 and 6) |
|---|
| 66 | / |
|---|
| 67 | !----------------------------------------------------------------------- |
|---|
| 68 | &namzgr ! vertical coordinate |
|---|
| 69 | !----------------------------------------------------------------------- |
|---|
| 70 | ln_zps = .true. ! z-coordinate - partial steps |
|---|
| 71 | / |
|---|
| 72 | !----------------------------------------------------------------------- |
|---|
| 73 | &namzgr_sco ! s-coordinate or hybrid z-s-coordinate |
|---|
| 74 | !----------------------------------------------------------------------- |
|---|
| 75 | / |
|---|
| 76 | !----------------------------------------------------------------------- |
|---|
| 77 | &namc1d ! 1D configuration options ("key_c1d" default: PAPA station) |
|---|
| 78 | !----------------------------------------------------------------------- |
|---|
| 79 | rn_lat1d = 50. ! Column latitude (default at PAPA station) |
|---|
| 80 | rn_lon1d = -145. ! Column longitude (default at PAPA station) |
|---|
| 81 | ln_c1d_locpt= .true. ! Localization of 1D config in a grid (T) or independant point (F) |
|---|
| 82 | / |
|---|
| 83 | !----------------------------------------------------------------------- |
|---|
| 84 | &namc1d_uvd ! data: U & V currents ("key_c1d" default: NO) |
|---|
| 85 | !----------------------------------------------------------------------- |
|---|
| 86 | / |
|---|
| 87 | !----------------------------------------------------------------------- |
|---|
| 88 | &namc1d_dyndmp ! U & V newtonian damping ("key_c1d" default: NO) |
|---|
| 89 | !----------------------------------------------------------------------- |
|---|
| 90 | / |
|---|
| 91 | !----------------------------------------------------------------------- |
|---|
| 92 | &namtsd ! Temperature & Salinity Data (default: NO) |
|---|
| 93 | !----------------------------------------------------------------------- |
|---|
| 94 | ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
|---|
| 95 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
|---|
| 96 | sn_tem = 'init_PAPASTATION' , 24 ,'votemper', .false. , .true., 'daily' , '' , ' ' , '' |
|---|
| 97 | sn_sal = 'init_PAPASTATION' , 24 ,'vosaline', .false. , .true., 'daily' , '' , ' ' , '' |
|---|
| 98 | ! |
|---|
| 99 | ln_tsd_init = .true. ! Initialisation of ocean T & S with T &S input data (T) or not (F) |
|---|
| 100 | ln_tsd_dmp = .false. ! damping of ocean T & S toward T &S input data (T) or not (F) |
|---|
| 101 | cn_dir = './' ! root directory for the location of the runoff files |
|---|
| 102 | / |
|---|
| 103 | |
|---|
| 104 | !!====================================================================== |
|---|
| 105 | !! *** Surface Boundary Condition namelists *** !! |
|---|
| 106 | !!====================================================================== |
|---|
| 107 | ! |
|---|
| 108 | !----------------------------------------------------------------------- |
|---|
| 109 | &namsbc ! Surface Boundary Condition manager |
|---|
| 110 | !----------------------------------------------------------------------- |
|---|
| 111 | nn_fsbc = 1 ! frequency of the calling of SBC module |
|---|
| 112 | ! Type of air-sea fluxes |
|---|
| 113 | ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk ) |
|---|
| 114 | ! Misc. options of sbc : |
|---|
| 115 | ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr) |
|---|
| 116 | / |
|---|
| 117 | !----------------------------------------------------------------------- |
|---|
| 118 | &namsbc_flx ! surface boundary condition : flux formulation |
|---|
| 119 | !----------------------------------------------------------------------- |
|---|
| 120 | / |
|---|
| 121 | !----------------------------------------------------------------------- |
|---|
| 122 | &namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk = T) |
|---|
| 123 | !----------------------------------------------------------------------- |
|---|
| 124 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
|---|
| 125 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
|---|
| 126 | sn_wndi = 'forcing_PAPASTATION_1h' , 1 , 'wndwe' , .true. , .false. , 'yearly' , '' , '', '' |
|---|
| 127 | sn_wndj = 'forcing_PAPASTATION_1h' , 1 , 'wndsn' , .true. , .false. , 'yearly' , '' , '', '' |
|---|
| 128 | sn_qsr = 'forcing_PAPASTATION_1h' , 1 , 'qsr' , .false. , .false. , 'yearly' , '' , '', '' |
|---|
| 129 | sn_qlw = 'forcing_PAPASTATION_1h' , 1 , 'qlw' , .false. , .false. , 'yearly' , '' , '', '' |
|---|
| 130 | sn_tair = 'forcing_PAPASTATION_1h' , 1 , 'tair' , .true. , .false. , 'yearly' , '' , '', '' |
|---|
| 131 | sn_humi = 'forcing_PAPASTATION_1h' , 1 , 'humi' , .true. , .false. , 'yearly' , '' , '', '' |
|---|
| 132 | sn_prec = 'forcing_PAPASTATION_1h' , 1 , 'prec' , .false. , .false. , 'yearly' , '' , '', '' |
|---|
| 133 | sn_snow = 'forcing_PAPASTATION_1h' , 1 , 'snow' , .false. , .false. , 'yearly' , '' , '', '' |
|---|
| 134 | sn_slp = 'slp.15JUNE2009_fill' , 6 , 'SLP', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' |
|---|
| 135 | |
|---|
| 136 | ! ! bulk algorithm : |
|---|
| 137 | ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008) |
|---|
| 138 | ln_COARE_3p0= .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) |
|---|
| 139 | ln_COARE_3p5= .false. ! "COARE 3.5" algorithm (Edson et al. 2013) |
|---|
| 140 | ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31) |
|---|
| 141 | |
|---|
| 142 | rn_zqt = 2. ! Air temperature & humidity reference height (m) |
|---|
| 143 | / |
|---|
| 144 | !----------------------------------------------------------------------- |
|---|
| 145 | &namtra_qsr ! penetrative solar radiation (ln_traqsr =T) |
|---|
| 146 | !----------------------------------------------------------------------- |
|---|
| 147 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
|---|
| 148 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
|---|
| 149 | sn_chl ='chlorophyll_PAPASTATION', -1 , 'CHLA' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 150 | |
|---|
| 151 | ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration |
|---|
| 152 | |
|---|
| 153 | nn_chldta = 1 ! RGB : Chl data (=1) or cst value (=0) |
|---|
| 154 | cn_dir = './' ! root directory for the location of the runoff files |
|---|
| 155 | / |
|---|
| 156 | !----------------------------------------------------------------------- |
|---|
| 157 | &namdrg ! top/bottom drag coefficient (default: NO selection) |
|---|
| 158 | !----------------------------------------------------------------------- |
|---|
| 159 | ln_non_lin = .true. ! non-linear drag: Cd = Cd0 |U| |
|---|
| 160 | / |
|---|
| 161 | |
|---|
| 162 | !!====================================================================== |
|---|
| 163 | !! Tracer (T & S) namelists |
|---|
| 164 | !!====================================================================== |
|---|
| 165 | ! |
|---|
| 166 | !----------------------------------------------------------------------- |
|---|
| 167 | &nameos ! ocean physical parameters (default: NO selection) |
|---|
| 168 | !----------------------------------------------------------------------- |
|---|
| 169 | ln_eos80 = .true. ! = Use EOS80 equation of state |
|---|
| 170 | / |
|---|
| 171 | !----------------------------------------------------------------------- |
|---|
| 172 | &namtra_adv ! advection scheme for tracer (default: NO selection) |
|---|
| 173 | !----------------------------------------------------------------------- |
|---|
| 174 | ln_traadv_NONE= .true. ! No tracer advection |
|---|
| 175 | / |
|---|
| 176 | !----------------------------------------------------------------------- |
|---|
| 177 | &namtra_ldf ! lateral diffusion scheme for tracers |
|---|
| 178 | !----------------------------------------------------------------------- |
|---|
| 179 | ln_traldf_NONE= .true. ! No operator (no explicit diffusion) |
|---|
| 180 | / |
|---|
| 181 | |
|---|
| 182 | !!====================================================================== |
|---|
| 183 | !! *** Dynamics namelists *** !! |
|---|
| 184 | !!====================================================================== |
|---|
| 185 | ! |
|---|
| 186 | !----------------------------------------------------------------------- |
|---|
| 187 | &nam_vvl ! vertical coordinate options (default: z-star) |
|---|
| 188 | !----------------------------------------------------------------------- |
|---|
| 189 | / |
|---|
| 190 | !----------------------------------------------------------------------- |
|---|
| 191 | &namdyn_adv ! formulation of the momentum advection (default: selection) |
|---|
| 192 | !----------------------------------------------------------------------- |
|---|
| 193 | ln_dynadv_NONE= .true. ! linear dynamics (no momentum advection) |
|---|
| 194 | / |
|---|
| 195 | !----------------------------------------------------------------------- |
|---|
| 196 | &namdyn_vor ! Vorticity / Coriolis scheme (default: NO selection) |
|---|
| 197 | !----------------------------------------------------------------------- |
|---|
| 198 | ln_dynvor_ene = .true. ! energy conserving scheme |
|---|
| 199 | / |
|---|
| 200 | !----------------------------------------------------------------------- |
|---|
| 201 | &namdyn_hpg ! Hydrostatic pressure gradient option |
|---|
| 202 | !----------------------------------------------------------------------- |
|---|
| 203 | / |
|---|
| 204 | !----------------------------------------------------------------------- |
|---|
| 205 | &namdyn_spg ! surface pressure gradient (default: NO) |
|---|
| 206 | !----------------------------------------------------------------------- |
|---|
| 207 | / |
|---|
| 208 | !----------------------------------------------------------------------- |
|---|
| 209 | &namdyn_ldf ! lateral diffusion on momentum |
|---|
| 210 | !----------------------------------------------------------------------- |
|---|
| 211 | ln_dynldf_NONE= .true. ! No operator (no explicit diffusion) |
|---|
| 212 | / |
|---|
| 213 | |
|---|
| 214 | !!====================================================================== |
|---|
| 215 | !! vertical physics namelists !! |
|---|
| 216 | !!====================================================================== |
|---|
| 217 | ! |
|---|
| 218 | !----------------------------------------------------------------------- |
|---|
| 219 | &namzdf ! vertical physics (default: NO selection) |
|---|
| 220 | !----------------------------------------------------------------------- |
|---|
| 221 | ! ! type of vertical closure |
|---|
| 222 | ln_zdfcst = .false. ! constant mixing |
|---|
| 223 | ln_zdfric = .false. ! local Richardson dependent formulation (T => fill namzdf_ric) |
|---|
| 224 | ln_zdftke = .false. ! Turbulent Kinetic Energy closure (T => fill namzdf_tke) |
|---|
| 225 | ln_zdfgls = .true. ! Generic Length Scale closure (T => fill namzdf_gls) |
|---|
| 226 | ln_zdfosm = .false. ! OSMOSIS BL closure (T => fill namzdf_osm) |
|---|
| 227 | ! |
|---|
| 228 | ! ! convection |
|---|
| 229 | ln_zdfevd = .false. ! enhanced vertical diffusion |
|---|
| 230 | nn_evdm = 0 ! apply on tracer (=0) or on tracer and momentum (=1) |
|---|
| 231 | rn_evd = 100. ! mixing coefficient [m2/s] |
|---|
| 232 | ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm |
|---|
| 233 | nn_npc = 1 ! frequency of application of npc |
|---|
| 234 | nn_npcp = 365 ! npc control print frequency |
|---|
| 235 | ! |
|---|
| 236 | ln_zdfddm = .false. ! double diffusive mixing |
|---|
| 237 | rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity) |
|---|
| 238 | rn_hsbfr = 1.6 ! heat/salt buoyancy flux ratio |
|---|
| 239 | ! |
|---|
| 240 | ! ! gravity wave-driven vertical mixing |
|---|
| 241 | ln_zdfiwm = .false. ! internal wave-induced mixing (T => fill namzdf_iwm) |
|---|
| 242 | ln_zdfswm = .false. ! surface wave-induced mixing (T => ln_wave=ln_sdw=T ) |
|---|
| 243 | ! |
|---|
| 244 | ! ! coefficients |
|---|
| 245 | rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) |
|---|
| 246 | rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) |
|---|
| 247 | nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) |
|---|
| 248 | nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) |
|---|
| 249 | / |
|---|
| 250 | !----------------------------------------------------------------------- |
|---|
| 251 | &namzdf_ric ! richardson number dependent vertical diffusion (ln_zdfric =T) |
|---|
| 252 | !----------------------------------------------------------------------- |
|---|
| 253 | / |
|---|
| 254 | !----------------------------------------------------------------------- |
|---|
| 255 | &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T) |
|---|
| 256 | !----------------------------------------------------------------------- |
|---|
| 257 | / |
|---|
| 258 | !----------------------------------------------------------------------- |
|---|
| 259 | &namzdf_gls ! GLS vertical diffusion (ln_zdfgls =T) |
|---|
| 260 | !----------------------------------------------------------------------- |
|---|
| 261 | / |
|---|
| 262 | !----------------------------------------------------------------------- |
|---|
| 263 | &namzdf_osm ! OSM vertical diffusion (ln_zdfosm =T) |
|---|
| 264 | !----------------------------------------------------------------------- |
|---|
| 265 | / |
|---|
| 266 | !----------------------------------------------------------------------- |
|---|
| 267 | &namzdf_iwm ! internal wave-driven mixing parameterization (ln_zdfiwm =T) |
|---|
| 268 | !----------------------------------------------------------------------- |
|---|
| 269 | / |
|---|
