| 1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 2 | !! PISCES : 1 - air-sea exchange (nampisext) |
|---|
| 3 | !! namelists 2 - biological parameters (nampisbio) |
|---|
| 4 | !! 3 - parameters for nutrient limitations (nampislim) |
|---|
| 5 | !! 4 - parameters for phytoplankton (nampisprod,nampismort) |
|---|
| 6 | !! 5 - parameters for zooplankton (nampismes,nampiszoo) |
|---|
| 7 | !! 6 - parameters for remineralization (nampisrem) |
|---|
| 8 | !! 7 - parameters for calcite chemistry (nampiscal) |
|---|
| 9 | !! 8 - parameters for inputs deposition (nampissed) |
|---|
| 10 | !! 9 - parameters for Kriest parameterization (nampiskrp, nampiskrs) |
|---|
| 11 | !! 10 - additional 2D/3D diagnostics (nampisdia) |
|---|
| 12 | !! 11 - Damping (nampisdmp) |
|---|
| 13 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 14 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 15 | &nampisext ! air-sea exchange |
|---|
| 16 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 17 | ln_co2int = .false. ! read atm pco2 from a file (T) or constant (F) |
|---|
| 18 | atcco2 = 287. ! Constant value atmospheric pCO2 - ln_co2int = F |
|---|
| 19 | clname = 'atcco2.txt' ! Name of atm pCO2 file - ln_co2int = T |
|---|
| 20 | nn_offset = 0 ! Offset model-data start year - ln_co2int = T |
|---|
| 21 | ! ! If your model year is iyy, nn_offset=(years(1)-iyy) |
|---|
| 22 | ! ! then the first atmospheric CO2 record read is at years(1) |
|---|
| 23 | / |
|---|
| 24 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 25 | &nampisatm ! Atmospheric prrssure |
|---|
| 26 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 27 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! |
|---|
| 28 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
|---|
| 29 | sn_patm = 'presatm' , -1 , 'patm' , .true. , .true. , 'yearly' , '' , '' |
|---|
| 30 | cn_dir = './' ! root directory for the location of the dynamical files |
|---|
| 31 | / |
|---|
| 32 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 33 | &nampisbio ! biological parameters |
|---|
| 34 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 35 | nrdttrc = 1 ! time step frequency for biology |
|---|
| 36 | wsbio = 2. ! POC sinking speed |
|---|
| 37 | xkmort = 1.E-7 ! half saturation constant for mortality |
|---|
| 38 | ferat3 = 10.E-6 ! Fe/C in zooplankton |
|---|
| 39 | wsbio2 = 30. ! Big particles sinking speed |
|---|
| 40 | / |
|---|
| 41 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 42 | &nampislim ! parameters for nutrient limitations |
|---|
| 43 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 44 | conc0 = 2.e-6 ! Phosphate half saturation |
|---|
| 45 | conc1 = 8E-6 ! Phosphate half saturation for diatoms |
|---|
| 46 | conc2 = 1E-9 ! Iron half saturation for phyto |
|---|
| 47 | conc2m = 3E-9 ! Max iron half saturation for phyto |
|---|
| 48 | conc3 = 2E-9 ! Iron half saturation for diatoms |
|---|
| 49 | conc3m = 8E-9 ! Maxi iron half saturation for diatoms |
|---|
| 50 | concnnh4 = 1.E-7 ! NH4 half saturation for phyto |
|---|
| 51 | concdnh4 = 4.E-7 ! NH4 half saturation for diatoms |
|---|
| 52 | xksi1 = 2.E-6 ! half saturation constant for Si uptake |
|---|
| 53 | xksi2 = 3.33E-6 ! half saturation constant for Si/C |
|---|
| 54 | xkdoc = 417.E-6 ! half-saturation constant of DOC remineralization |
|---|
| 55 | concfebac = 3.E-11 ! Half-saturation for Fe limitation of Bacteria |
|---|
| 56 | qnfelim = 7.E-6 ! Optimal quota of phyto |
|---|
| 57 | qdfelim = 7.E-6 ! Optimal quota of diatoms |
|---|
| 58 | caco3r = 0.16 ! mean rain ratio |
|---|
| 59 | / |
|---|
| 60 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 61 | &nampisprod ! parameters for phytoplankton growth |
|---|
| 62 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 63 | pislope = 2. ! P-I slope |
|---|
| 64 | pislope2 = 2. ! P-I slope for diatoms |
|---|
| 65 | excret = 0.05 ! excretion ratio of phytoplankton |
|---|
| 66 | excret2 = 0.05 ! excretion ratio of diatoms |
|---|
| 67 | bresp = 0.00333 ! Basal respiration rate |
|---|
| 68 | chlcnm = 0.033 ! Minimum Chl/C in nanophytoplankton |
|---|
| 69 | chlcdm = 0.04 ! Minimum Chl/C in diatoms |
|---|
| 70 | chlcmin = 0.0033 ! Maximum Chl/c in phytoplankton |
|---|
| 71 | fecnm = 40E-6 ! Maximum Fe/C in nanophytoplankton |
|---|
| 72 | fecdm = 40E-6 ! Minimum Fe/C in diatoms |
|---|
| 73 | grosip = 0.151 ! mean Si/C ratio |
|---|
| 74 | / |
|---|
| 75 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 76 | &nampismort ! parameters for phytoplankton sinks |
|---|
| 77 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 78 | wchl = 0.001 ! quadratic mortality of phytoplankton |
|---|
| 79 | wchld = 0.02 ! maximum quadratic mortality of diatoms |
|---|
| 80 | mprat = 0.01 ! phytoplankton mortality rate |
|---|
| 81 | mprat2 = 0.01 ! Diatoms mortality rate |
|---|
| 82 | mpratm = 0.01 ! Phytoplankton minimum mortality rate |
|---|
| 83 | / |
|---|
| 84 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 85 | &nampismes ! parameters for mesozooplankton |
|---|
| 86 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 87 | part2 = 0.75 ! part of calcite not dissolved in mesozoo guts |
|---|
| 88 | grazrat2 = 0.9 ! maximal mesozoo grazing rate |
|---|
| 89 | resrat2 = 0.005 ! exsudation rate of mesozooplankton |
|---|
| 90 | mzrat2 = 0.04 ! mesozooplankton mortality rate |
|---|
| 91 | xprefc = 1. ! zoo preference for phyto |
|---|
| 92 | xprefp = 0.3 ! zoo preference for POC |
|---|
| 93 | xprefz = 1. ! zoo preference for zoo |
|---|
| 94 | xprefpoc = 0.3 ! zoo preference for poc |
|---|
| 95 | xthresh2 = 2E-7 ! Food threshold for grazing |
|---|
| 96 | xkgraz2 = 20.E-6 ! half sturation constant for meso grazing |
|---|
| 97 | epsher2 = 0.3 ! Efficicency of Mesozoo growth |
|---|
| 98 | sigma2 = 0.6 ! Fraction of mesozoo excretion as DOM |
|---|
| 99 | unass2 = 0.3 ! non assimilated fraction of P by mesozoo |
|---|
| 100 | grazflux = 3.e3 ! flux-feeding rate |
|---|
| 101 | / |
|---|
| 102 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 103 | &nampiszoo ! parameters for microzooplankton |
|---|
| 104 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 105 | part = 0.5 ! part of calcite not dissolved in microzoo gutsa |
|---|
| 106 | grazrat = 3.0 ! maximal zoo grazing rate |
|---|
| 107 | resrat = 0.03 ! exsudation rate of zooplankton |
|---|
| 108 | mzrat = 0.0 ! zooplankton mortality rate |
|---|
| 109 | xpref2c = 0.2 ! Microzoo preference for POM |
|---|
| 110 | xpref2p = 1. ! Microzoo preference for Nanophyto |
|---|
| 111 | xpref2d = 0.6 ! Microzoo preference for Diatoms |
|---|
| 112 | xthresh = 2E-7 ! Food threshold for feeding |
|---|
| 113 | xkgraz = 20.E-6 ! half sturation constant for grazing |
|---|
| 114 | epsher = 0.3 ! Efficiency of microzoo growth |
|---|
| 115 | sigma1 = 0.6 ! Fraction of microzoo excretion as DOM |
|---|
| 116 | unass = 0.3 ! non assimilated fraction of phyto by zoo |
|---|
| 117 | / |
|---|
| 118 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 119 | &nampisrem ! parameters for remineralization |
|---|
| 120 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 121 | xremik = 0.25 ! remineralization rate of DOC |
|---|
| 122 | xremip = 0.025 ! remineralisation rate of POC |
|---|
| 123 | nitrif = 0.05 ! NH4 nitrification rate |
|---|
| 124 | xsirem = 0.003 ! remineralization rate of Si |
|---|
| 125 | xsiremlab = 0.025 ! fast remineralization rate of Si |
|---|
| 126 | xsilab = 0.31 ! Fraction of labile biogenic silica |
|---|
| 127 | xlam1 = 0.005 ! scavenging rate of Iron |
|---|
| 128 | oxymin = 1.E-6 ! Half-saturation constant for anoxia |
|---|
| 129 | ligand = 0.6E-9 ! Ligands concentration |
|---|
| 130 | / |
|---|
| 131 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 132 | &nampiscal ! parameters for Calcite chemistry |
|---|
| 133 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 134 | kdca = 6. ! calcite dissolution rate constant (1/time) |
|---|
| 135 | nca = 1. ! order of dissolution reaction (dimensionless) |
|---|
| 136 | / |
|---|
| 137 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 138 | &nampissed ! parameters for inputs deposition |
|---|
| 139 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 140 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! |
|---|
| 141 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
|---|
| 142 | sn_dust = 'dust.orca' , -1 , 'dust' , .true. , .true. , 'yearly' , '' , '' |
|---|
| 143 | sn_riverdic = 'river.orca' , -12 , 'riverdic' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 144 | sn_riverdoc = 'river.orca' , -12 , 'riverdoc' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 145 | sn_ndepo = 'ndeposition.orca', -12 , 'ndep' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 146 | sn_ironsed = 'bathy.orca' , -12 , 'bathy' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 147 | ! |
|---|
| 148 | cn_dir = './' ! root directory for the location of the dynamical files |
|---|
| 149 | ln_dust = .true. ! boolean for dust input from the atmosphere |
|---|
| 150 | ln_river = .false. ! boolean for river input of nutrients |
|---|
| 151 | ln_ndepo = .true. ! boolean for atmospheric deposition of N |
|---|
| 152 | ln_ironsed = .true. ! boolean for Fe input from sediments |
|---|
| 153 | sedfeinput = 1E-9 ! Coastal release of Iron |
|---|
| 154 | dustsolub = 0.02 ! Solubility of the dust |
|---|
| 155 | wdust = 2.0 ! Dust sinking speed |
|---|
| 156 | nitrfix = 1E-7 ! Nitrogen fixation rate |
|---|
| 157 | diazolight = 50. ! Diazotrophs sensitivity to light (W/m2) |
|---|
| 158 | concfediaz = 1.E-10 ! Diazotrophs half-saturation Cste for Iron |
|---|
| 159 | / |
|---|
| 160 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 161 | &nampiskrp ! Kriest parameterization : parameters "key_kriest" |
|---|
| 162 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 163 | xkr_eta = 1.17 ! Sinking exponent |
|---|
| 164 | xkr_zeta = 2.28 ! N content exponent |
|---|
| 165 | xkr_mass_min = 0.0002 ! Minimum mass for Aggregates |
|---|
| 166 | xkr_mass_max = 1. ! Maximum mass for Aggregates |
|---|
| 167 | / |
|---|
| 168 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 169 | &nampiskrs ! Kriest parameterization : size classes "key_kriest" |
|---|
| 170 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 171 | xkr_sfact = 942. ! Sinking factor |
|---|
| 172 | xkr_stick = 0.5 ! Stickiness |
|---|
| 173 | xkr_nnano = 2.337 ! Nbr of cell in nano size class |
|---|
| 174 | xkr_ndiat = 3.718 ! Nbr of cell in diatoms size class |
|---|
| 175 | xkr_nmeso = 7.147 ! Nbr of cell in mesozoo size class |
|---|
| 176 | xkr_naggr = 9.877 ! Nbr of cell in aggregates size class |
|---|
| 177 | / |
|---|
| 178 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 179 | &nampisdia ! additional 2D/3D tracers diagnostics |
|---|
| 180 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 181 | nn_writedia = 5475 ! time step frequency for tracers diagnostics |
|---|
| 182 | ! ! name ! title of the field ! units ! |
|---|
| 183 | ! ! ! ! ! |
|---|
| 184 | pisdia2d(1) = 'Cflx ' , 'DIC flux ', 'molC/m2/s ' |
|---|
| 185 | pisdia2d(2) = 'Oflx ' , 'Oxygen flux ', 'molC/m2/s ' |
|---|
| 186 | pisdia2d(3) = 'Kg ' , 'Gas transfer ', 'mol/m2/s/uatm' |
|---|
| 187 | pisdia2d(4) = 'Delc ' , 'Delta CO2 ', 'uatm ' |
|---|
| 188 | pisdia2d(5) = 'PMO ' , 'POC export ', 'molC/m2/s ' |
|---|
| 189 | pisdia2d(6) = 'PMO2 ' , 'GOC export ', 'molC/m2/s ' |
|---|
| 190 | pisdia2d(7) = 'ExpFe1 ' , 'Nano iron export ', 'molFe/m2/s ' |
|---|
| 191 | pisdia2d(8) = 'ExpFe2 ' , 'Diatoms iron export ', 'molFe/m2/s ' |
|---|
| 192 | pisdia2d(9) = 'ExpSi ' , 'Silicate export ', 'molSi/m2/s ' |
|---|
| 193 | pisdia2d(10) = 'ExpCaCO3 ' , 'Calcite export ', 'molC/m2/s ' |
|---|
| 194 | pisdia2d(11) = 'heup ' , 'euphotic layer depth ', 'm ' |
|---|
| 195 | pisdia2d(12) = 'Fedep ' , 'Iron dep ', 'molFe/m2/s ' |
|---|
| 196 | pisdia2d(13) = 'Nfix ' , 'Nitrogen Fixation ', 'molN/m2/s ' |
|---|
| 197 | pisdia3d(1) = 'PH ' , 'PH ', '- ' |
|---|
| 198 | pisdia3d(2) = 'CO3 ' , 'Bicarbonates ', 'mol/l ' |
|---|
| 199 | pisdia3d(3) = 'CO3sat ' , 'CO3 saturation ', 'mol/l ' |
|---|
| 200 | pisdia3d(4) = 'PAR ' , 'light penetration ', 'W/m2 ' |
|---|
| 201 | pisdia3d(5) = 'PPPHY ' , 'Primary production of nanophyto ', 'molC/m3/s ' |
|---|
| 202 | pisdia3d(6) = 'PPPHY2 ' , 'Primary production of diatoms ', 'molC/m3/s ' |
|---|
| 203 | pisdia3d(7) = 'PPNEWN ' , 'New Primary production of nano ', 'molC/m3/s ' |
|---|
| 204 | pisdia3d(8) = 'PPNEWD ' , 'New Primary production of diat ', 'molC/m3/s ' |
|---|
| 205 | pisdia3d(9) = 'PBSi ' , 'Primary production of Si diatoms ', 'molSi/m3/s ' |
|---|
| 206 | pisdia3d(10) = 'PFeN ' , 'Primary production of nano iron ', 'molFe/m3/s ' |
|---|
| 207 | pisdia3d(11) = 'PFeD ' , 'Primary production of diatoms iron', 'molFe/m3/s ' |
|---|
| 208 | / |
|---|
| 209 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 210 | &nampisdmp ! Damping |
|---|
| 211 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 212 | ln_pisdmp = .true. ! Relaxation fo some tracers to a mean value |
|---|
| 213 | / |
|---|
