| 1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 2 | !! PISCES (key_pisces) reference namelist (see below for key_pisces_reduced) |
|---|
| 3 | !! 1 - air-sea exchange (nampisext) |
|---|
| 4 | !! 2 - biological parameters (nampisbio) |
|---|
| 5 | !! 3 - parameters for nutrient limitations (nampislim) |
|---|
| 6 | !! 4 - parameters for phytoplankton (nampisprod,nampismort) |
|---|
| 7 | !! 5 - parameters for zooplankton (nampismes,nampiszoo) |
|---|
| 8 | !! 6 - parameters for remineralization (nampisrem) |
|---|
| 9 | !! 7 - parameters for calcite chemistry (nampiscal) |
|---|
| 10 | !! 8 - parameters for inputs deposition (nampissed) |
|---|
| 11 | !! 9 - parameters for Kriest parameterization (nampiskrp, nampiskrs) |
|---|
| 12 | !! 10 - additional 2D/3D diagnostics (nampisdia) |
|---|
| 13 | !! 11 - Damping (nampisdmp) |
|---|
| 14 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 15 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 16 | &nampisext ! air-sea exchange |
|---|
| 17 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 18 | ln_co2int = .false. ! read atm pco2 from a file (T) or constant (F) |
|---|
| 19 | atcco2 = 280. ! Constant value atmospheric pCO2 - ln_co2int = F |
|---|
| 20 | clname = 'atcco2.txt' ! Name of atm pCO2 file - ln_co2int = T |
|---|
| 21 | nn_offset = 0 ! Offset model-data start year - ln_co2int = T |
|---|
| 22 | ! ! If your model year is iyy, nn_offset=(years(1)-iyy) |
|---|
| 23 | ! ! then the first atmospheric CO2 record read is at years(1) |
|---|
| 24 | / |
|---|
| 25 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 26 | &nampisatm ! Atmospheric prrssure |
|---|
| 27 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 28 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
|---|
| 29 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
|---|
| 30 | sn_patm = 'presatm' , -1 , 'patm' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 31 | sn_atmco2 = 'presatmco2_reg' , -1 , 'xco2' , .true. , .true. , 'yearly' , 'weights_AIRS_bilinear.nc' , '' , '' |
|---|
| 32 | cn_dir = './' ! root directory for the location of the dynamical files |
|---|
| 33 | ! |
|---|
| 34 | ln_presatm = .false. ! constant atmopsheric pressure (F) or from a file (T) |
|---|
| 35 | ln_presatmco2= .false. ! Read spatialized atm co2 files [ppm] if TRUE |
|---|
| 36 | / |
|---|
| 37 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 38 | &nampisbio ! biological parameters |
|---|
| 39 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 40 | nrdttrc = 3 ! time step frequency for biology |
|---|
| 41 | wsbio = 2. ! POC sinking speed |
|---|
| 42 | xkmort = 1.E-7 ! half saturation constant for mortality |
|---|
| 43 | ferat3 = 10.E-6 ! Fe/C in zooplankton |
|---|
| 44 | wsbio2 = 50. ! Big particles sinking speed |
|---|
| 45 | wsbio2max = 50. ! Big particles maximum sinking speed |
|---|
| 46 | wsbio2scale= 5000. ! Big particles length scale of sinking |
|---|
| 47 | niter1max = 1 ! Maximum number of iterations for POC |
|---|
| 48 | niter2max = 2 ! Maximum number of iterations for GOC |
|---|
| 49 | wfep = 0.2 ! FeP sinking speed |
|---|
| 50 | ldocp = 1.E-5 ! Phyto ligand production per unit doc |
|---|
| 51 | ldocz = 1.E-5 ! Zoo ligand production per unit doc |
|---|
| 52 | lthet = 0.5 ! Proportional loss of ligands due to Fe uptake |
|---|
| 53 | / |
|---|
| 54 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 55 | &nampislim ! parameters for nutrient limitations |
|---|
| 56 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 57 | concnno3 = 1.e-6 ! Nitrate half saturation of nanophytoplankton |
|---|
| 58 | concdno3 = 3.E-6 ! Nitrate half saturation for diatoms |
|---|
| 59 | concnnh4 = 1.E-7 ! NH4 half saturation for phyto |
|---|
| 60 | concdnh4 = 3.E-7 ! NH4 half saturation for diatoms |
|---|
| 61 | concnfer = 1.E-9 ! Iron half saturation for phyto |
|---|
| 62 | concdfer = 3.E-9 ! Iron half saturation for diatoms |
|---|
| 63 | concbfe = 1.E-11 ! Iron half-saturation for DOC remin. |
|---|
| 64 | concbnh4 = 2.E-8 ! NH4 half saturation for DOC remin. |
|---|
| 65 | concbno3 = 2.E-7 ! Nitrate half saturation for DOC remin. |
|---|
| 66 | xsizedia = 1.E-6 ! Minimum size criteria for diatoms |
|---|
| 67 | xsizephy = 1.E-6 ! Minimum size criteria for phyto |
|---|
| 68 | xsizern = 3.0 ! Size ratio for nanophytoplankton |
|---|
| 69 | xsizerd = 3.0 ! Size ratio for diatoms |
|---|
| 70 | xksi1 = 2.E-6 ! half saturation constant for Si uptake |
|---|
| 71 | xksi2 = 20E-6 ! half saturation constant for Si/C |
|---|
| 72 | xkdoc = 417.E-6 ! half-saturation constant of DOC remineralization |
|---|
| 73 | qnfelim = 7.E-6 ! Optimal quota of phyto |
|---|
| 74 | qdfelim = 7.E-6 ! Optimal quota of diatoms |
|---|
| 75 | caco3r = 0.3 ! mean rain ratio |
|---|
| 76 | / |
|---|
| 77 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 78 | &nampisopt ! parameters for optics |
|---|
| 79 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 80 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
|---|
| 81 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
|---|
| 82 | sn_par = 'par.orca' , 24 , 'fr_par' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 83 | cn_dir = './' ! root directory for the location of the dynamical files |
|---|
| 84 | ln_varpar = .true. ! boolean for PAR variable |
|---|
| 85 | parlux = 0.43 ! Fraction of shortwave as PAR |
|---|
| 86 | / |
|---|
| 87 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 88 | &nampisprod ! parameters for phytoplankton growth |
|---|
| 89 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 90 | pislope = 2. ! P-I slope |
|---|
| 91 | pislope2 = 2. ! P-I slope for diatoms |
|---|
| 92 | xadap = 0. ! Adaptation factor to low light |
|---|
| 93 | excret = 0.05 ! excretion ratio of phytoplankton |
|---|
| 94 | excret2 = 0.05 ! excretion ratio of diatoms |
|---|
| 95 | ln_newprod = .true. ! Enable new parame. of production (T/F) |
|---|
| 96 | bresp = 0.033 ! Basal respiration rate |
|---|
| 97 | chlcnm = 0.033 ! Maximum Chl/C in nanophytoplankton |
|---|
| 98 | chlcdm = 0.05 ! Maximum Chl/C in diatoms |
|---|
| 99 | chlcmin = 0.004 ! Minimum Chl/c in phytoplankton |
|---|
| 100 | fecnm = 40E-6 ! Maximum Fe/C in nanophytoplankton |
|---|
| 101 | fecdm = 40E-6 ! Maximum Fe/C in diatoms |
|---|
| 102 | grosip = 0.159 ! mean Si/C ratio |
|---|
| 103 | / |
|---|
| 104 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 105 | &nampismort ! parameters for phytoplankton sinks |
|---|
| 106 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 107 | wchl = 0.01 ! quadratic mortality of phytoplankton |
|---|
| 108 | wchld = 0.01 ! maximum quadratic mortality of diatoms |
|---|
| 109 | wchldm = 0.03 ! maximum quadratic mortality of diatoms |
|---|
| 110 | mprat = 0.01 ! phytoplankton mortality rate |
|---|
| 111 | mprat2 = 0.01 ! Diatoms mortality rate |
|---|
| 112 | / |
|---|
| 113 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 114 | &nampismes ! parameters for mesozooplankton |
|---|
| 115 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 116 | part2 = 0.75 ! part of calcite not dissolved in mesozoo guts |
|---|
| 117 | grazrat2 = 0.75 ! maximal mesozoo grazing rate |
|---|
| 118 | resrat2 = 0.005 ! exsudation rate of mesozooplankton |
|---|
| 119 | mzrat2 = 0.03 ! mesozooplankton mortality rate |
|---|
| 120 | xprefc = 1. ! mesozoo preference for diatoms |
|---|
| 121 | xprefp = 0.3 ! mesozoo preference for nanophyto. |
|---|
| 122 | xprefz = 1. ! mesozoo preference for microzoo. |
|---|
| 123 | xprefpoc = 0.3 ! mesozoo preference for poc |
|---|
| 124 | xthresh2zoo = 1E-8 ! zoo feeding threshold for mesozooplankton |
|---|
| 125 | xthresh2dia = 1E-8 ! diatoms feeding threshold for mesozooplankton |
|---|
| 126 | xthresh2phy = 1E-8 ! nanophyto feeding threshold for mesozooplankton |
|---|
| 127 | xthresh2poc = 1E-8 ! poc feeding threshold for mesozooplankton |
|---|
| 128 | xthresh2 = 3E-7 ! Food threshold for grazing |
|---|
| 129 | xkgraz2 = 20.E-6 ! half saturation constant for meso grazing |
|---|
| 130 | epsher2 = 0.35 ! Efficicency of Mesozoo growth |
|---|
| 131 | sigma2 = 0.6 ! Fraction of mesozoo excretion as DOM |
|---|
| 132 | unass2 = 0.3 ! non assimilated fraction of P by mesozoo |
|---|
| 133 | grazflux = 2.e3 ! flux-feeding rate |
|---|
| 134 | / |
|---|
| 135 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 136 | &nampiszoo ! parameters for microzooplankton |
|---|
| 137 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 138 | part = 0.5 ! part of calcite not dissolved in microzoo gutsa |
|---|
| 139 | grazrat = 3.0 ! maximal zoo grazing rate |
|---|
| 140 | resrat = 0.03 ! exsudation rate of zooplankton |
|---|
| 141 | mzrat = 0.004 ! zooplankton mortality rate |
|---|
| 142 | xpref2c = 0.1 ! Microzoo preference for POM |
|---|
| 143 | xpref2p = 1. ! Microzoo preference for Nanophyto |
|---|
| 144 | xpref2d = 0.5 ! Microzoo preference for Diatoms |
|---|
| 145 | xthreshdia = 1.E-8 ! Diatoms feeding threshold for microzooplankton |
|---|
| 146 | xthreshphy = 1.E-8 ! Nanophyto feeding threshold for microzooplankton |
|---|
| 147 | xthreshpoc = 1.E-8 ! POC feeding threshold for microzooplankton |
|---|
| 148 | xthresh = 3.E-7 ! Food threshold for feeding |
|---|
| 149 | xkgraz = 20.E-6 ! half sturation constant for grazing |
|---|
| 150 | epsher = 0.3 ! Efficiency of microzoo growth |
|---|
| 151 | sigma1 = 0.6 ! Fraction of microzoo excretion as DOM |
|---|
| 152 | unass = 0.3 ! non assimilated fraction of phyto by zoo |
|---|
| 153 | / |
|---|
| 154 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 155 | &nampisfer ! parameters for iron chemistry |
|---|
| 156 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 157 | ln_fechem = .false. ! complex iron chemistry ( T/F ) |
|---|
| 158 | ln_ligvar = .false. ! variable ligand concentration |
|---|
| 159 | ln_fecolloid = .false. ! variable colloidal fraction |
|---|
| 160 | xlam1 = 0.005 ! scavenging rate of Iron |
|---|
| 161 | xlamdust = 150.0 ! Scavenging rate of dust |
|---|
| 162 | ligand = 0.6E-9 ! Ligands concentration |
|---|
| 163 | kfep = 0. ! Nanoparticle formation rate constant |
|---|
| 164 | / |
|---|
| 165 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 166 | &nampisrem ! parameters for remineralization |
|---|
| 167 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 168 | xremik = 0.3 ! remineralization rate of DOC |
|---|
| 169 | nitrif = 0.05 ! NH4 nitrification rate |
|---|
| 170 | xsirem = 0.003 ! remineralization rate of Si |
|---|
| 171 | xsiremlab = 0.03 ! fast remineralization rate of Si |
|---|
| 172 | xsilab = 0.5 ! Fraction of labile biogenic silica |
|---|
| 173 | oxymin = 1.E-6 ! Half-saturation constant for anoxia |
|---|
| 174 | oxymin2 = 6.E-6 ! Minimum O2 concentration for oxic remin. |
|---|
| 175 | feratb = 10.E-6 ! Fe/C quota in bacteria |
|---|
| 176 | xkferb = 2.5E-10 ! Half-saturation constant for bacteria Fe/C |
|---|
| 177 | / |
|---|
| 178 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 179 | &nampispoc ! parameters for organic particles |
|---|
| 180 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 181 | xremip = 0.035 ! remineralisation rate of PON |
|---|
| 182 | jcpoc = 15 ! Number of lability classes |
|---|
| 183 | rshape = 1.0 ! Shape of the gamma function |
|---|
| 184 | / |
|---|
| 185 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 186 | &nampiscal ! parameters for Calcite chemistry |
|---|
| 187 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 188 | kdca = 6. ! calcite dissolution rate constant (1/time) |
|---|
| 189 | nca = 1. ! order of dissolution reaction (dimensionless) |
|---|
| 190 | / |
|---|
| 191 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 192 | &nampissbc ! parameters for inputs deposition |
|---|
| 193 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 194 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
|---|
| 195 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
|---|
| 196 | sn_dust = 'dust.orca' , -1 , 'dust' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 197 | sn_solub = 'solubility.orca' , -12 , 'solubility1' , .false. , .true. , 'yearly' , '' , '' , '' |
|---|
| 198 | sn_riverdic = 'river.orca' , 120 , 'riverdic' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 199 | sn_riverdoc = 'river.orca' , 120 , 'riverdoc' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 200 | sn_riverdin = 'river.orca' , 120 , 'riverdin' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 201 | sn_riverdon = 'river.orca' , 120 , 'riverdon' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 202 | sn_riverdip = 'river.orca' , 120 , 'riverdip' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 203 | sn_riverdop = 'river.orca' , 120 , 'riverdop' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 204 | sn_riverdsi = 'river.orca' , 120 , 'riverdsi' , .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 205 | sn_ndepo = 'ndeposition.orca', -12 , 'ndep' , .false. , .true. , 'yearly' , '' , '' , '' |
|---|
| 206 | sn_ironsed = 'bathy.orca' , -12 , 'bathy' , .false. , .true. , 'yearly' , '' , '' , '' |
|---|
| 207 | sn_hydrofe = 'hydrofe.orca' , -12 , 'epsdb' , .false. , .true. , 'yearly' , '' , '' , '' |
|---|
| 208 | ! |
|---|
| 209 | cn_dir = './' ! root directory for the location of the dynamical files |
|---|
| 210 | ln_dust = .true. ! boolean for dust input from the atmosphere |
|---|
| 211 | ln_solub = .true. ! boolean for variable solubility of atm. Iron |
|---|
| 212 | ln_river = .true. ! boolean for river input of nutrients |
|---|
| 213 | ln_ndepo = .true. ! boolean for atmospheric deposition of N |
|---|
| 214 | ln_ironsed = .true. ! boolean for Fe input from sediments |
|---|
| 215 | ln_ironice = .true. ! boolean for Fe input from sea ice |
|---|
| 216 | ln_hydrofe = .false. ! boolean for from hydrothermal vents |
|---|
| 217 | sedfeinput = 2.e-9 ! Coastal release of Iron |
|---|
| 218 | dustsolub = 0.02 ! Solubility of the dusta |
|---|
| 219 | mfrac = 0.035 ! Fe mineral fraction of dust |
|---|
| 220 | wdust = 2.0 ! Dust sinking speed |
|---|
| 221 | icefeinput = 15.e-9 ! Iron concentration in sea ice |
|---|
| 222 | nitrfix = 1.e-7 ! Nitrogen fixation rate |
|---|
| 223 | diazolight = 50. ! Diazotrophs sensitivity to light (W/m2) |
|---|
| 224 | concfediaz = 1.e-10 ! Diazotrophs half-saturation Cste for Iron |
|---|
| 225 | hratio = 1.e+7 ! Fe to 3He ratio assumed for vent iron supply |
|---|
| 226 | fep_rats = 1. ! Fep/Fer ratio from sed sources |
|---|
| 227 | fep_rath = 1. ! Fep/Fer ratio from sed hydro sources |
|---|
| 228 | / |
|---|
| 229 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 230 | &nampislig ! Namelist parameters for ligands, nampislig |
|---|
| 231 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 232 | rfep = 0.001 ! Dissolution rate of FeP |
|---|
| 233 | rlgw = 1. ! Lifetime (years) of weak ligands |
|---|
| 234 | rlig = 1.E-4 ! Remin ligand production per unit C |
|---|
| 235 | prlgw = 1.E-4 ! Photolysis of weak ligand |
|---|
| 236 | rlgs = 1000. ! Lifetime (years) of strong ligands |
|---|
| 237 | / |
|---|
| 238 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 239 | &nampisice ! Prescribed sea ice tracers |
|---|
| 240 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 241 | ! constant ocean tracer concentrations are defined in trcice_pisces.F90 (Global, Arctic, Antarctic and Baltic) |
|---|
| 242 | ! trc_ice_ratio * betw 0 and 1: prescribed ice/ocean tracer concentration ratio |
|---|
| 243 | ! * -1 => the ice-ocean tracer concentration ratio follows the |
|---|
| 244 | ! ice-ocean salinity ratio |
|---|
| 245 | ! * -2 => tracer concentration in sea ice is prescribed and |
|---|
| 246 | ! trc_ice_prescr is used |
|---|
| 247 | ! trc_ice_prescr * prescribed tracer concentration. used only if |
|---|
| 248 | ! trc_ice_ratio = -2. equals -99 if not used. |
|---|
| 249 | ! cn_trc_o * 'GL' use global ocean values making the Baltic distinction only |
|---|
| 250 | ! 'AA' use specific Arctic/Antarctic/Baltic values |
|---|
| 251 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 252 | ! sn_tri_ ! trc_ice_ratio ! trc_ice_prescr ! cn_trc_o |
|---|
| 253 | sn_tri_dic = -1., -99., 'AA' |
|---|
| 254 | sn_tri_doc = 0., -99., 'AA' |
|---|
| 255 | sn_tri_tal = -1., -99., 'AA' |
|---|
| 256 | sn_tri_oxy = -1., -99., 'AA' |
|---|
| 257 | sn_tri_cal = 0., -99., 'AA' |
|---|
| 258 | sn_tri_po4 = -1., -99., 'AA' |
|---|
| 259 | sn_tri_poc = 0., -99., 'AA' |
|---|
| 260 | sn_tri_goc = 0., -99., 'AA' |
|---|
| 261 | sn_tri_bfe = 0., -99., 'AA' |
|---|
| 262 | sn_tri_num = 0., -99., 'AA' |
|---|
| 263 | sn_tri_sil = -1., -99., 'AA' |
|---|
| 264 | sn_tri_dsi = 0., -99., 'AA' |
|---|
| 265 | sn_tri_gsi = 0., -99., 'AA' |
|---|
| 266 | sn_tri_phy = 0., -99., 'AA' |
|---|
| 267 | sn_tri_dia = 0., -99., 'AA' |
|---|
| 268 | sn_tri_zoo = 0., -99., 'AA' |
|---|
| 269 | sn_tri_mes = 0., -99., 'AA' |
|---|
| 270 | sn_tri_fer = -2., 15E-9, 'AA' |
|---|
| 271 | sn_tri_sfe = 0., -99., 'AA' |
|---|
| 272 | sn_tri_dfe = 0., -99., 'AA' |
|---|
| 273 | sn_tri_nfe = 0., -99., 'AA' |
|---|
| 274 | sn_tri_nch = 0., -99., 'AA' |
|---|
| 275 | sn_tri_dch = 0., -99., 'AA' |
|---|
| 276 | sn_tri_no3 = -1., -99., 'AA' |
|---|
| 277 | sn_tri_nh4 = 1., -99., 'AA' |
|---|
| 278 | / |
|---|
| 279 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 280 | &nampiskrp ! Kriest parameterization : parameters "key_kriest" |
|---|
| 281 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 282 | xkr_eta = 1.17 ! Sinking exponent |
|---|
| 283 | xkr_zeta = 2.28 ! N content exponent |
|---|
| 284 | xkr_ncontent = 5.7E-6 ! N content factor |
|---|
| 285 | xkr_mass_min = 0.0002 ! Minimum mass for Aggregates |
|---|
| 286 | xkr_mass_max = 1. ! Maximum mass for Aggregates |
|---|
| 287 | / |
|---|
| 288 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 289 | &nampiskrs ! Kriest parameterization : size classes "key_kriest" |
|---|
| 290 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 291 | xkr_sfact = 942. ! Sinking factor |
|---|
| 292 | xkr_stick = 0.5 ! Stickiness |
|---|
| 293 | xkr_nnano = 2.337 ! Nbr of cell in nano size class |
|---|
| 294 | xkr_ndiat = 3.718 ! Nbr of cell in diatoms size class |
|---|
| 295 | xkr_nmeso = 7.147 ! Nbr of cell in mesozoo size class |
|---|
| 296 | xkr_naggr = 9.877 ! Nbr of cell in aggregates size class |
|---|
| 297 | / |
|---|
| 298 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 299 | &nampisdia ! additional 2D/3D tracers diagnostics |
|---|
| 300 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 301 | ! ! name ! title of the field ! units ! |
|---|
| 302 | ! ! ! ! ! |
|---|
| 303 | pisdia2d(1) = 'Cflx ' , 'DIC flux ', 'molC/m2/s ' |
|---|
| 304 | pisdia2d(2) = 'Oflx ' , 'Oxygen flux ', 'molC/m2/s ' |
|---|
| 305 | pisdia2d(3) = 'Kg ' , 'Gas transfer ', 'mol/m2/s/uatm' |
|---|
| 306 | pisdia2d(4) = 'Delc ' , 'Delta CO2 ', 'uatm ' |
|---|
| 307 | pisdia2d(5) = 'PMO ' , 'POC export ', 'molC/m2/s ' |
|---|
| 308 | pisdia2d(6) = 'PMO2 ' , 'GOC export ', 'molC/m2/s ' |
|---|
| 309 | pisdia2d(7) = 'ExpFe1 ' , 'Nano iron export ', 'molFe/m2/s ' |
|---|
| 310 | pisdia2d(8) = 'ExpFe2 ' , 'Diatoms iron export ', 'molFe/m2/s ' |
|---|
| 311 | pisdia2d(9) = 'ExpSi ' , 'Silicate export ', 'molSi/m2/s ' |
|---|
| 312 | pisdia2d(10) = 'ExpCaCO3 ' , 'Calcite export ', 'molC/m2/s ' |
|---|
| 313 | pisdia2d(11) = 'heup ' , 'euphotic layer depth ', 'm ' |
|---|
| 314 | pisdia2d(12) = 'Fedep ' , 'Iron dep ', 'molFe/m2/s ' |
|---|
| 315 | pisdia2d(13) = 'Nfix ' , 'Nitrogen Fixation ', 'molN/m2/s ' |
|---|
| 316 | pisdia3d(1) = 'PH ' , 'PH ', '- ' |
|---|
| 317 | pisdia3d(2) = 'CO3 ' , 'Bicarbonates ', 'mol/l ' |
|---|
| 318 | pisdia3d(3) = 'CO3sat ' , 'CO3 saturation ', 'mol/l ' |
|---|
| 319 | pisdia3d(4) = 'PAR ' , 'light penetration ', 'W/m2 ' |
|---|
| 320 | pisdia3d(5) = 'PPPHY ' , 'Primary production of nanophyto ', 'molC/m3/s ' |
|---|
| 321 | pisdia3d(6) = 'PPPHY2 ' , 'Primary production of diatoms ', 'molC/m3/s ' |
|---|
| 322 | pisdia3d(7) = 'PPNEWN ' , 'New Primary production of nano ', 'molC/m3/s ' |
|---|
| 323 | pisdia3d(8) = 'PPNEWD ' , 'New Primary production of diat ', 'molC/m3/s ' |
|---|
| 324 | pisdia3d(9) = 'PBSi ' , 'Primary production of Si diatoms ', 'molSi/m3/s ' |
|---|
| 325 | pisdia3d(10) = 'PFeN ' , 'Primary production of nano iron ', 'molFe/m3/s ' |
|---|
| 326 | pisdia3d(11) = 'PFeD ' , 'Primary production of diatoms iron', 'molFe/m3/s ' |
|---|
| 327 | / |
|---|
| 328 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 329 | &nampisdmp ! Damping |
|---|
| 330 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 331 | ln_pisdmp = .true. ! Relaxation fo some tracers to a mean value |
|---|
| 332 | nn_pisdmp = 5475 ! Frequency of Relaxation |
|---|
| 333 | / |
|---|
| 334 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 335 | &nampismass ! Mass conservation |
|---|
| 336 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 337 | ln_check_mass = .false. ! Check mass conservation |
|---|
| 338 | / |
|---|
| 339 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 340 | !! PISCES reduced (key_pisces_reduced, ex LOBSTER) : namelists |
|---|
| 341 | !! 1 - biological parameters for phytoplankton (namlobphy) |
|---|
| 342 | !! 2 - biological parameters for nutrients (namlobnut) |
|---|
| 343 | !! 3 - biological parameters for zooplankton (namlobzoo) |
|---|
| 344 | !! 4 - biological parameters for detritus (namlobdet) |
|---|
| 345 | !! 5 - biological parameters for DOM (namlobdom) |
|---|
| 346 | !! 6 - parameters from aphotic layers to sediment (namlobsed) |
|---|
| 347 | !! 7 - general coefficients (namlobrat) |
|---|
| 348 | !! 8 - optical parameters (namlobopt) |
|---|
| 349 | |
|---|
| 350 | !! 10 - biological diagnostics trends (namlobdbi) |
|---|
| 351 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 352 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 353 | &namlobphy ! biological parameters for phytoplankton |
|---|
| 354 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 355 | tmumax = 1.21e-5 ! maximal phytoplankton growth rate [s-1] |
|---|
| 356 | rgamma = 0.05 ! phytoplankton exudation fraction [%] |
|---|
| 357 | fphylab = 0.75 ! NH4 fraction of phytoplankton exsudation |
|---|
| 358 | tmminp = 5.8e-7 ! minimal phytoplancton mortality rate [0.05/86400 s-1=20 days] |
|---|
| 359 | aki = 33. ! light photosynthesis half saturation constant[W/m2] |
|---|
| 360 | / |
|---|
| 361 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 362 | &namlobnut ! biological parameters for nutrients |
|---|
| 363 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 364 | akno3 = 0.7 ! nitrate limitation half-saturation value [mmol/m3] |
|---|
| 365 | aknh4 = 0.001 ! ammonium limitation half-saturation value [mmol/m3] |
|---|
| 366 | taunn = 5.80e-7 ! nitrification rate [s-1] |
|---|
| 367 | psinut = 3. ! inhibition of nitrate uptake by ammonium |
|---|
| 368 | / |
|---|
| 369 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 370 | &namlobzoo ! biological parameters for zooplankton |
|---|
| 371 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 372 | rppz = 0.8 ! zooplankton nominal preference for phytoplancton food [%] |
|---|
| 373 | taus = 9.26E-6 ! specific zooplankton maximal grazing rate [s-1] |
|---|
| 374 | ! ! 0.75/86400 s-1=8.680555E-6 1/86400 = 1.15e-5 |
|---|
| 375 | aks = 1. ! half-saturation constant for total zooplankton grazing [mmolN.m-3] |
|---|
| 376 | rpnaz = 0.3 ! non-assimilated phytoplankton by zooplancton [%] |
|---|
| 377 | rdnaz = 0.3 ! non-assimilated detritus by zooplankton [%] |
|---|
| 378 | tauzn = 8.1e-7 ! zooplancton specific excretion rate [0.1/86400 s-1=10 days] |
|---|
| 379 | fzoolab = 0.5 ! NH4 fraction of zooplankton excretion |
|---|
| 380 | fdbod = 0.5 ! zooplankton mortality fraction that goes to detritus |
|---|
| 381 | tmminz = 2.31e-6 ! minimal zooplankton mortality rate [(mmolN/m3)-1 d-1] |
|---|
| 382 | / |
|---|
| 383 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 384 | &namlobdet ! biological parameters for detritus |
|---|
| 385 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 386 | taudn = 5.80e-7 ! detritus breakdown rate [0.1/86400 s-1=10 days] |
|---|
| 387 | fdetlab = 0. ! NH4 fraction of detritus dissolution |
|---|
| 388 | / |
|---|
| 389 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 390 | &namlobdom ! biological parameters for DOM |
|---|
| 391 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 392 | taudomn = 6.43e-8 ! DOM breakdown rate [s-1] |
|---|
| 393 | ! ! slow remineralization rate of semi-labile dom to nh4 (1 month) |
|---|
| 394 | / |
|---|
| 395 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 396 | &namlobsed ! parameters from aphotic layers to sediment |
|---|
| 397 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 398 | sedlam = 3.86e-7 ! time coefficient of POC remineralization in sediments [s-1] |
|---|
| 399 | sedlostpoc = 0. ! mass of POC lost in sediments |
|---|
| 400 | vsed = 3.47e-5 ! detritus sedimentation speed [m/s] |
|---|
| 401 | xhr = -0.858 ! coeff for martin''s remineralisation profile |
|---|
| 402 | / |
|---|
| 403 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 404 | &namlobrat ! general coefficients |
|---|
| 405 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 406 | rcchl = 60. ! Carbone/Chlorophyl ratio [mgC.mgChla-1] |
|---|
| 407 | redf = 6.56 ! redfield ratio (C:N) for phyto |
|---|
| 408 | reddom = 6.56 ! redfield ratio (C:N) for DOM |
|---|
| 409 | / |
|---|
| 410 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 411 | &namlobopt ! optical parameters |
|---|
| 412 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 413 | xkg0 = 0.0232 ! green absorption coefficient of water |
|---|
| 414 | xkr0 = 0.225 ! red absorption coefficent of water |
|---|
| 415 | xkgp = 0.074 ! green absorption coefficient of chl |
|---|
| 416 | xkrp = 0.037 ! red absorption coefficient of chl |
|---|
| 417 | xlg = 0.674 ! green chl exposant for absorption |
|---|
| 418 | xlr = 0.629 ! red chl exposant for absorption |
|---|
| 419 | rpig = 0.7 ! chla/chla+pheo ratio |
|---|
| 420 | / |
|---|
| 421 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 422 | &nampisdbi ! biological diagnostics trends |
|---|
| 423 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 424 | ! ! 2D bio diagnostics units : mmole/m2/s ("key_trdmld_trc") |
|---|
| 425 | ! ! name ! title of the field ! units ! |
|---|
| 426 | pisdiabio(1) = 'NO3PHY' , 'Flux from NO3 to PHY ', 'mmole/m3/s' |
|---|
| 427 | pisdiabio(2) = 'NH4PHY' , 'Flux from NH4 to PHY ', 'mmole/m3/s' |
|---|
| 428 | pisdiabio(3) = 'PHYNH4' , 'Flux from PHY to NH4 ', 'mmole/m3/s' |
|---|
| 429 | pisdiabio(4) = 'PHYDOM' , 'Flux from PHY to DOM ', 'mmole/m3/s' |
|---|
| 430 | pisdiabio(5) = 'PHYZOO' , 'Flux from PHY to ZOO ', 'mmole/m3/s' |
|---|
| 431 | pisdiabio(6) = 'PHYDET' , 'Flux from PHY to DET ', 'mmole/m3/s' |
|---|
| 432 | pisdiabio(7) = 'DETZOO' , 'Flux from DET to ZOO ', 'mmole/m3/s' |
|---|
| 433 | pisdiabio(8) = 'DETSED' , 'Flux from DET to SED ', 'mmole/m3/s' |
|---|
| 434 | pisdiabio(9) = 'ZOODET' , 'Flux from ZOO to DET ', 'mmole/m3/s' |
|---|
| 435 | pisdiabio(10) = 'ZOOBOD' , 'Zooplankton closure ', 'mmole/m3/s' |
|---|
| 436 | pisdiabio(11) = 'ZOONH4' , 'Flux from ZOO to NH4 ', 'mmole/m3/s' |
|---|
| 437 | pisdiabio(12) = 'ZOODOM' , 'Flux from ZOO to DOM ', 'mmole/m3/s' |
|---|
| 438 | pisdiabio(13) = 'NH4NO3' , 'Flux from NH4 to NO3 ', 'mmole/m3/s' |
|---|
| 439 | pisdiabio(14) = 'DOMNH4' , 'Flux from DOM to NH4 ', 'mmole/m3/s' |
|---|
| 440 | pisdiabio(15) = 'DETNH4' , 'Flux from DET to NH4 ', 'mmole/m3/s' |
|---|
| 441 | pisdiabio(16) = 'DETDOM' , 'Flux from DET to DOM ', 'mmole/m3/s' |
|---|
| 442 | pisdiabio(17) = 'SEDNO3' , 'NO3 remineralization from SED', 'mmole/m3/s' |
|---|
| 443 | / |
|---|
