Changeset 336 for trunk/NEMO/TOP_SRC/SMS/trclsm.pisces.h90
- Timestamp:
- 2005-11-14T13:13:26+01:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/SMS/trclsm.pisces.h90
r274 r336 1 !!2 !! trclsm.pisces.h3 !! ****************4 !!5 !! PURPOSE :6 !! ---------7 !! READs and PRINT options for PISCES namelist8 !!9 !! MODIFICATIONS:10 !! --------------11 !! original : 99-10 (M.A. Foujols, M. Levy) passive tracer12 !! addition : 00-01 (L. Bopp) hamocc3,p3zd13 !!14 1 !!---------------------------------------------------------------------- 2 !! *** trclsm.pisces.h90 *** 15 3 !!---------------------------------------------------------------------- 16 !! local declarations 17 !! ================== 18 19 #if defined key_passivetrc && defined key_trc_pisces 4 CONTAINS 5 6 SUBROUTINE trc_lsm 7 !!---------------------------------------------------------------------- 8 !! 9 !! trclsm.pisces.h 10 !! **************** 11 !! 12 !! PURPOSE : 13 !! --------- 14 !! READs and PRINT options for PISCES namelist 15 !! 16 !! MODIFICATIONS: 17 !! -------------- 18 !! original : 99-10 (M.A. Foujols, M. Levy) passive tracer 19 !! addition : 00-01 (L. Bopp) hamocc3,p3zd 20 !! 21 !!---------------------------------------------------------------------- 22 !!---------------------------------------------------------------------- 23 !! local declarations 24 !! ================== 20 25 CHARACTER (len=32) clname 21 26 22 !!--------------------------------------------------------------------- 23 !! TOP 1.0, LOCEAN-IPSL (2005) 24 !! $Header$ 25 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 26 !!--------------------------------------------------------------------- 27 28 ! 0. initializations 29 ! ------------------ 30 ! 31 NAMELIST/natgas/ gasfac, igaswind, icice 32 NAMELIST/natext/ atcco2 33 NAMELIST/natbio/caco3r, & 34 & dispo0,conc0,oxymin,grosip, sedlam, & 35 & sedlostpoc,sedlostcal,sedlostsil,nrdttrc, & 36 & pislope, excret,wsbio,wchl,resrat,mprat,mzrat, & 37 & grazrat,xprefc,xprefp,unass,xkgraz,xkmort,xksi1, & 38 & xksi2,sicmax,xremip,xremik,xsirem,xkdoc1,xkdoc2, & 39 & excret2,resrat2,mprat2,mpratm,mzrat2,grazrat2, & 40 & xprefz,xprefpoc,unass2,wchl2,xkgraz2,xlam1, & 41 & ferat3,conc1,conc2,conc3,concnnh4,concdnh4, & 42 & nitrif,epsher,epsher2,pislope2,wsbio2,sigma1, & 43 & sigma2, zprefc, zprefp, zprefd,sedfeinput 44 NAMELIST/natsms/bdustfer, briver, bndepo, bsedinput 45 46 ! initialize the number of LOGICAL UNIT used 47 ! ------------------------------------------ 48 49 IF(lwp) THEN 50 WRITE(numout,*) ' ' 51 WRITE(numout,*) ' ROUTINE trclec' 52 WRITE(numout,*) ' **************' 53 WRITE(numout,*) ' ' 54 WRITE(numout,*) ' namelist for PISCES model' 55 WRITE(numout,*) ' ***********************' 56 WRITE(numout,*) ' ' 27 !!--------------------------------------------------------------------- 28 !! OPA8, LODYC (15/11/96) 29 !!--------------------------------------------------------------------- 30 31 ! 0. initializations 32 ! ------------------ 33 ! 34 NAMELIST/natgas/ gasfac, igaswind, icice 35 NAMELIST/natext/ atcco2 36 NAMELIST/natbio/caco3r, & 37 & dispo0,conc0,oxymin,grosip, sedlam, & 38 & sedlostpoc,sedlostcal,sedlostsil,nrdttrc, & 39 & pislope, excret,wsbio,wchl,wchld,resrat,mprat,mzrat, & 40 & grazrat,xprefc,xprefp,unass,xkgraz,xkmort,xksi1, & 41 & xksi2,xremip,xremik,xsirem,xkdoc1,xkdoc2, & 42 & excret2,resrat2,mprat2,mpratm,mzrat2,grazrat2, & 43 & xprefz,xprefpoc,unass2,xkgraz2,xlam1, & 44 & ferat3,conc1,conc2,conc3,concnnh4,concdnh4, & 45 & nitrif,epsher,epsher2,pislope2,wsbio2,sigma1, & 46 & sigma2, zprefc, zprefp, zprefd,fecnm,fecdm, & 47 & chlcnm,chlcdm, sedfeinput 48 NAMELIST/natsms/bdustfer, briver, bndepo, bsedinput 49 50 ! initialize the number of LOGICAL UNIT used 51 ! ------------------------------------------ 52 53 IF(lwp) THEN 54 WRITE(numout,*) ' ' 55 WRITE(numout,*) ' ROUTINE trclec' 56 WRITE(numout,*) ' **************' 57 WRITE(numout,*) ' ' 58 WRITE(numout,*) ' namelist for PISCES model' 59 WRITE(numout,*) ' ***********************' 60 WRITE(numout,*) ' ' 57 61 ENDIF 58 62 … … 62 66 63 67 64 ! 1 Namelist natgas :65 ! -------------------68 ! 1 Namelist natgas : 69 ! ------------------- 66 70 67 71 READ(numnat,natgas) 68 72 69 73 IF(lwp) THEN 70 71 72 73 74 75 76 77 78 ENDIF 79 80 ! 2 Namelist natext :81 ! -------------------74 WRITE(numout,*) ' ' 75 WRITE(numout,*) 'natgas' 76 write(numout,*) 'gasfac = ',gasfac 77 WRITE(numout,*) ' ' 78 write(numout,*) 'igaswind = ',igaswind 79 WRITE(numout,*) ' ' 80 write(numout,*) 'icice = ',icice 81 WRITE(numout,*) ' ' 82 ENDIF 83 84 ! 2 Namelist natext : 85 ! ------------------- 82 86 READ(numnat,natext) 83 87 84 88 IF(lwp) THEN 85 86 87 88 89 89 WRITE(numout,*) ' ' 90 WRITE(numout,*) 'natext' 91 WRITE(numout,*) ' ' 92 WRITE(numout,*) 'atmospheric pCO2= ',atcco2 93 WRITE(numout,*) ' ' 90 94 ENDIF 91 95 … … 93 97 READ(numnat,natbio) 94 98 IF(lwp) THEN 95 WRITE(numout,*) 'natbio' 96 WRITE(numout,*) ' ' 97 WRITE(numout,*) & 98 & ' mean rainratio =', caco3r 99 WRITE(numout,*) & 100 & ' mean Si/C ratio =', grosip 101 WRITE(numout,*) & 102 & ' Calcite dissolution half saturation =', dispo0 103 WRITE(numout,*) & 104 & ' Phosphate half saturation =', conc0 105 WRITE(numout,*) & 106 & ' Sediment bioturbation factor =', sedlam 107 WRITE(numout,*) & 108 & ' Sediment burying ratio for POC =', sedlostpoc 109 WRITE(numout,*) & 110 & ' Sediment burying ratio for CACO3 =', sedlostcal 111 WRITE(numout,*) & 112 & ' Sediment burying ratio for SI =', sedlostsil 113 WRITE(numout,*) & 114 & ' frequence pour la biologie =', nrdttrc 115 WRITE(numout,*) & 116 & ' P-I slope =', pislope 117 WRITE(numout,*) & 118 & ' excretion ratio of phytoplankton =', excret 119 WRITE(numout,*) & 120 & ' POC sinking speed =', wsbio 121 WRITE(numout,*) & 122 & ' exsudation rate of zooplankton =', resrat 123 WRITE(numout,*) & 124 & ' phytoplankton mortality rate =', mprat 125 WRITE(numout,*) & 126 & ' zooplankton mortality rate =', mzrat 127 WRITE(numout,*) & 128 & ' zoo preference for phyto =', xprefc 129 WRITE(numout,*) & 130 & ' zoo preference for POC =', xprefp 131 WRITE(numout,*) & 132 & ' maximal zoo grazing rate =', grazrat 133 WRITE(numout,*) & 134 & ' non assimilated fraction of phyto by zoo =', unass 135 WRITE(numout,*) & 136 & ' half sturation constant for grazing =', xkgraz 137 WRITE(numout,*) & 138 & ' half saturation constant for mortality =', xkmort 139 WRITE(numout,*) & 140 & ' half saturation constant for Si uptake =', xksi1 141 WRITE(numout,*) & 142 & ' half saturation constant for Si/C =', xksi2 143 WRITE(numout,*) & 144 & ' maximum Si/C =', sicmax 145 WRITE(numout,*) & 146 & ' remineralisation rate of POC =', xremip 147 WRITE(numout,*) & 148 & ' remineralization rate of DOC =', xremik 149 WRITE(numout,*) & 150 & ' remineralization rate of Si =', xsirem 151 WRITE(numout,*) & 152 & ' 1st half-sat. of DOC remineralization =', xkdoc1 153 WRITE(numout,*) & 154 & ' 2nd half-sat. of DOC remineralization =', xkdoc2 155 WRITE(numout,*) & 156 & ' excretion ratio of diatoms =', excret2 157 WRITE(numout,*) & 158 & ' exsudation rate of mesozooplankton =', resrat2 159 WRITE(numout,*) & 160 & ' Diatoms mortality rate =', mprat2 161 WRITE(numout,*) & 162 & ' Phytoplankton minimum mortality rate =', mpratm 163 WRITE(numout,*) & 164 & ' mesozooplankton mortality rate =', mzrat2 165 WRITE(numout,*) & 166 & ' zoo preference for zoo =', xprefz 167 WRITE(numout,*) & 168 & ' zoo preference for poc =', xprefpoc 169 WRITE(numout,*) & 170 & ' maximal mesozoo grazing rate =', grazrat2 171 WRITE(numout,*) & 172 & ' non assimilated fraction of P by mesozoo =', unass2 173 WRITE(numout,*) & 174 & ' Efficicency of Mesozoo growth =', epsher2 175 WRITE(numout,*) & 176 & ' Efficiency of microzoo growth =', epsher 177 WRITE(numout,*) & 178 & ' half sturation constant for grazing 2 =', xkgraz2 179 WRITE(numout,*) & 180 & ' Aggregation rate for diatoms =', wchl2 181 WRITE(numout,*) & 182 & ' scavenging rate of Iron =', xlam1 183 WRITE(numout,*) & 184 & ' Fe/C in zooplankton =', ferat3 185 WRITE(numout,*) & 186 & ' Phosphate half saturation for diatoms =', conc1 187 WRITE(numout,*) & 188 & ' Iron half saturation for phyto =', conc2 189 WRITE(numout,*) & 190 & ' Iron half saturation for diatoms =', conc3 191 WRITE(numout,*) & 192 & ' NH4 half saturation for phyto =', concnnh4 193 WRITE(numout,*) & 194 & ' NH4 half saturation for diatoms =', concdnh4 195 WRITE(numout,*) & 196 & ' NH4 nitrification rate =', nitrif 197 WRITE(numout,*) & 198 & ' P-I slope for diatoms =', pislope2 199 WRITE(numout,*) & 200 & ' Big particles sinking speed =', wsbio2 201 WRITE(numout,*) & 202 & ' Fraction of microzoo excretion as DOM =', sigma1 203 WRITE(numout,*) & 204 & ' Fraction of mesozoo excretion as DOM =', sigma2 205 WRITE(numout,*) & 206 & ' Microzoo preference for POM =', zprefc 207 WRITE(numout,*) & 208 & ' Microzoo preference for Nanophyto =', zprefp 209 WRITE(numout,*) & 210 & ' Microzoo preference for Diatoms =', zprefd 211 WRITE(numout,*) & 212 & ' Coastal release of Iron =', sedfeinput 99 WRITE(numout,*) 'natbio' 100 WRITE(numout,*) ' ' 101 WRITE(numout,*) & 102 & ' mean rainratio =', caco3r 103 WRITE(numout,*) & 104 & ' mean Si/C ratio =', grosip 105 WRITE(numout,*) & 106 & ' Calcite dissolution half saturation =', dispo0 107 WRITE(numout,*) & 108 & ' Phosphate half saturation =', conc0 109 WRITE(numout,*) & 110 & ' Sediment bioturbation factor =', sedlam 111 WRITE(numout,*) & 112 & ' Sediment burying ratio for POC =', sedlostpoc 113 WRITE(numout,*) & 114 & ' Sediment burying ratio for CACO3 =', sedlostcal 115 WRITE(numout,*) & 116 & ' Sediment burying ratio for SI =', sedlostsil 117 WRITE(numout,*) & 118 & ' frequence pour la biologie =', nrdttrc 119 WRITE(numout,*) & 120 & ' P-I slope =', pislope 121 WRITE(numout,*) & 122 & ' excretion ratio of phytoplankton =', excret 123 WRITE(numout,*) & 124 & ' POC sinking speed =', wsbio 125 WRITE(numout,*) & 126 & ' exsudation rate of zooplankton =', resrat 127 WRITE(numout,*) & 128 & ' phytoplankton mortality rate =', mprat 129 WRITE(numout,*) & 130 & ' zooplankton mortality rate =', mzrat 131 WRITE(numout,*) & 132 & ' zoo preference for phyto =', xprefc 133 WRITE(numout,*) & 134 & ' zoo preference for POC =', xprefp 135 WRITE(numout,*) & 136 & ' maximal zoo grazing rate =', grazrat 137 WRITE(numout,*) & 138 & ' non assimilated fraction of phyto by zoo =', unass 139 WRITE(numout,*) & 140 & ' half sturation constant for grazing =', xkgraz 141 WRITE(numout,*) & 142 & ' half saturation constant for mortality =', xkmort 143 WRITE(numout,*) & 144 & ' half saturation constant for Si uptake =', xksi1 145 WRITE(numout,*) & 146 & ' half saturation constant for Si/C =', xksi2 147 WRITE(numout,*) & 148 & ' remineralisation rate of POC =', xremip 149 WRITE(numout,*) & 150 & ' remineralization rate of DOC =', xremik 151 WRITE(numout,*) & 152 & ' remineralization rate of Si =', xsirem 153 WRITE(numout,*) & 154 & ' 1st half-sat. of DOC remineralization =', xkdoc1 155 WRITE(numout,*) & 156 & ' 2nd half-sat. of DOC remineralization =', xkdoc2 157 WRITE(numout,*) & 158 & ' excretion ratio of diatoms =', excret2 159 WRITE(numout,*) & 160 & ' exsudation rate of mesozooplankton =', resrat2 161 WRITE(numout,*) & 162 & ' Diatoms mortality rate =', mprat2 163 WRITE(numout,*) & 164 & ' Phytoplankton minimum mortality rate =', mpratm 165 WRITE(numout,*) & 166 & ' mesozooplankton mortality rate =', mzrat2 167 WRITE(numout,*) & 168 & ' zoo preference for zoo =', xprefz 169 WRITE(numout,*) & 170 & ' zoo preference for poc =', xprefpoc 171 WRITE(numout,*) & 172 & ' maximal mesozoo grazing rate =', grazrat2 173 WRITE(numout,*) & 174 & ' non assimilated fraction of P by mesozoo =', unass2 175 WRITE(numout,*) & 176 & ' Efficicency of Mesozoo growth =', epsher2 177 WRITE(numout,*) & 178 & ' Efficiency of microzoo growth =', epsher 179 WRITE(numout,*) & 180 & ' half sturation constant for grazing 2 =', xkgraz2 181 WRITE(numout,*) & 182 & ' Maximum aggregation rate for diatoms =', wchld 183 WRITE(numout,*) & 184 & ' scavenging rate of Iron =', xlam1 185 WRITE(numout,*) & 186 & ' Fe/C in zooplankton =', ferat3 187 WRITE(numout,*) & 188 & ' Phosphate half saturation for diatoms =', conc1 189 WRITE(numout,*) & 190 & ' Iron half saturation for phyto =', conc2 191 WRITE(numout,*) & 192 & ' Iron half saturation for diatoms =', conc3 193 WRITE(numout,*) & 194 & ' NH4 half saturation for phyto =', concnnh4 195 WRITE(numout,*) & 196 & ' NH4 half saturation for diatoms =', concdnh4 197 WRITE(numout,*) & 198 & ' NH4 nitrification rate =', nitrif 199 WRITE(numout,*) & 200 & ' P-I slope for diatoms =', pislope2 201 WRITE(numout,*) & 202 & ' Big particles sinking speed =', wsbio2 203 WRITE(numout,*) & 204 & ' Fraction of microzoo excretion as DOM =', sigma1 205 WRITE(numout,*) & 206 & ' Fraction of mesozoo excretion as DOM =', sigma2 207 WRITE(numout,*) & 208 & ' Microzoo preference for POM =', zprefc 209 WRITE(numout,*) & 210 & ' Microzoo preference for Nanophyto =', zprefp 211 WRITE(numout,*) & 212 & ' Microzoo preference for Diatoms =', zprefd 213 WRITE(numout,*) & 214 & ' Minimum Chl/C in nanophytoplankton =', chlcnm 215 WRITE(numout,*) & 216 & ' Minimum Chl/C in diatoms =', chlcdm 217 WRITE(numout,*) & 218 & ' Maximum Fe/C in nanophytoplankton =', fecnm 219 WRITE(numout,*) & 220 & ' Minimum Fe/C in diatoms =', fecdm 221 WRITE(numout,*) & 222 & ' Coastal release of Iron =', sedfeinput 213 223 ENDIF 214 224 215 225 READ(numnat,natsms) 216 226 IF(lwp) THEN 217 WRITE(numout,*) ' ' 218 WRITE(numout,*) 'natsms' 219 WRITE(numout,*) ' ' 220 WRITE(numout,*) 'Dust input from the atmosphere : ', bdustfer 221 WRITE(numout,*) ' ' 222 WRITE(numout,*) 'River input of nutrients : ', briver 223 WRITE(numout,*) ' ' 224 WRITE(numout,*) 'Atmospheric deposition of N : ', bndepo 225 WRITE(numout,*) ' ' 226 WRITE(numout,*) 'Fe input from sediments : ', bsedinput 227 WRITE(numout,*) ' ' 228 ENDIF 229 230 231 #else 232 233 ! no passive tracers 234 235 #endif 227 WRITE(numout,*) ' ' 228 WRITE(numout,*) 'natsms' 229 WRITE(numout,*) ' ' 230 WRITE(numout,*) 'Dust input from the atmosphere : ', bdustfer 231 WRITE(numout,*) ' ' 232 WRITE(numout,*) 'River input of nutrients : ', briver 233 WRITE(numout,*) ' ' 234 WRITE(numout,*) 'Atmospheric deposition of N : ', bndepo 235 WRITE(numout,*) ' ' 236 WRITE(numout,*) 'Fe input from sediments : ', bsedinput 237 WRITE(numout,*) ' ' 238 ENDIF 239 240 END SUBROUTINE trc_lsm
Note: See TracChangeset
for help on using the changeset viewer.