New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
trclsm.pisces.h90 in tags/nemo_v2_2/NEMO/TOP_SRC/SMS – NEMO

source: tags/nemo_v2_2/NEMO/TOP_SRC/SMS/trclsm.pisces.h90 @ 630

Last change on this file since 630 was 617, checked in by opalod, 17 years ago

* empty log message *
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 12.6 KB
Line 
1!!----------------------------------------------------------------------
2!!                    ***  trclsm.pisces.h90 ***
3!!----------------------------------------------------------------------
4CONTAINS
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      !! ==================
25      CHARACTER (len=32) clname
26
27      !!---------------------------------------------------------------------
28      !!  TOP 1.0 , LOCEAN-IPSL (2005)
29   !! $Header$
30   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
31      !!---------------------------------------------------------------------
32
33      ! 0. initializations
34      ! ------------------
35      !
36      NAMELIST/natext/ atcco2
37      NAMELIST/natbio/caco3r, kdca, nca, part,                       &
38         &          dispo0,conc0,oxymin,grosip, sedlam,              &
39         &          sedlostpoc,sedlostcal,sedlostsil,nrdttrc,        &
40         &          pislope, excret,wsbio,wchl,wchld,resrat,mprat,mzrat,   &
41         &          grazrat,xprefc,xprefp,unass,xkgraz,xkmort,xksi1, &
42         &          xksi2,xremip,xremik,xsirem,xkdoc1,xkdoc2, &
43         &          excret2,resrat2,mprat2,mpratm,mzrat2,grazrat2,   &
44         &          xprefz,xprefpoc,unass2,xkgraz2,xlam1,      &
45         &          ferat3,conc1,conc2,conc3,concnnh4,concdnh4,      &
46         &          nitrif,epsher,epsher2,pislope2,wsbio2,sigma1,    &
47         &          sigma2, zprefc, zprefp, zprefd,fecnm,fecdm,      &
48         &          chlcnm,chlcdm, sedfeinput
49      NAMELIST/natsms/bdustfer, briver, bndepo, bsedinput
50#if defined key_trc_kriest
51      NAMELIST/natkriest/xkr_eta,xkr_zeta,xkr_sfact,xkr_mass_min,xkr_mass_max,&
52         &               xkr_dnano,xkr_ddiat,xkr_dmeso,xkr_daggr,xkr_stick
53#endif
54
55      ! initialize the number of LOGICAL UNIT used
56      ! ------------------------------------------
57
58      IF(lwp) THEN
59         WRITE(numout,*) ' '
60         WRITE(numout,*) ' ROUTINE trclec'
61         WRITE(numout,*) ' **************'
62         WRITE(numout,*) ' '
63         WRITE(numout,*) ' namelist for PISCES model'
64         WRITE(numout,*) ' ***********************'
65         WRITE(numout,*) ' '
66      ENDIF
67
68      numnat=80
69      clname ='namelist.trc.sms'
70      OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old')
71
72
73      ! 1 Namelist natext :
74      ! -------------------
75      READ(numnat,natext)
76
77      IF(lwp) THEN
78         WRITE(numout,*) ' '
79         WRITE(numout,*) 'natext'
80         WRITE(numout,*) ' '
81         WRITE(numout,*) 'atmospheric pCO2= ',atcco2
82         WRITE(numout,*) ' '
83      ENDIF
84
85
86      READ(numnat,natbio)
87      IF(lwp) THEN
88         WRITE(numout,*) 'natbio'
89         WRITE(numout,*) ' '
90         WRITE(numout,*)      &
91            &   ' mean rainratio                             =', caco3r
92         WRITE(numout,*)      &
93            &   ' diss. rate constant calcite (per month)    =', kdca
94         WRITE(numout,*)      &
95            &   ' order of reaction for calcite dissolution  =', nca
96         WRITE(numout,*)      &
97            &   ' part of calcite not dissolved in guts      =', part
98         WRITE(numout,*)      &
99            &   ' mean Si/C ratio                            =', grosip
100         WRITE(numout,*)      &
101            &   ' Calcite dissolution half saturation        =', dispo0
102         WRITE(numout,*)      &
103            &   ' Phosphate half saturation                  =', conc0
104         WRITE(numout,*)      &
105            &   ' Sediment bioturbation factor               =', sedlam
106         WRITE(numout,*)      &
107            &   ' Sediment burying ratio for POC             =', sedlostpoc
108         WRITE(numout,*)      &
109            &   ' Sediment burying ratio for CACO3           =', sedlostcal
110         WRITE(numout,*)      &
111            &   ' Sediment burying ratio for SI              =', sedlostsil
112         WRITE(numout,*)      &
113            &   ' frequence pour la biologie                 =', nrdttrc
114         WRITE(numout,*)      &
115            &   ' P-I slope                                  =', pislope
116         WRITE(numout,*)      &
117            &   ' excretion ratio of phytoplankton           =', excret
118         WRITE(numout,*)      &
119            &   ' POC sinking speed                          =', wsbio
120         WRITE(numout,*)      &
121            &   ' quadratic mortality of phytoplankton       =', wchl
122         WRITE(numout,*)      &
123            &   ' maximum quadratic mortality of diatoms     =', wchld
124         WRITE(numout,*)      &
125            &   ' exsudation rate of zooplankton             =', resrat
126         WRITE(numout,*)      &
127            &   ' phytoplankton mortality rate               =', mprat
128         WRITE(numout,*)      &
129            &   ' zooplankton mortality rate                 =', mzrat
130         WRITE(numout,*)      &
131            &   ' zoo preference for phyto                   =', xprefc
132         WRITE(numout,*)      &
133            &   ' zoo preference for POC                     =', xprefp
134         WRITE(numout,*)      &
135            &   ' maximal zoo grazing rate                   =', grazrat
136         WRITE(numout,*)      &
137            &   ' non assimilated fraction of phyto by zoo   =', unass
138         WRITE(numout,*)      &
139            &   ' half sturation constant for grazing        =', xkgraz
140         WRITE(numout,*)      &
141            &   ' half saturation constant for mortality     =', xkmort
142         WRITE(numout,*)      &
143            &   ' half saturation constant for Si uptake     =', xksi1
144         WRITE(numout,*)      &
145            &   ' half saturation constant for Si/C          =', xksi2
146         WRITE(numout,*)      &
147            &   ' remineralisation rate of POC               =', xremip
148         WRITE(numout,*)      &
149            &   ' remineralization rate of DOC               =', xremik
150         WRITE(numout,*)      &
151            &   ' remineralization rate of Si                =', xsirem
152         WRITE(numout,*)      &
153            &   ' 1st half-sat. of DOC remineralization      =', xkdoc1
154         WRITE(numout,*)      &
155            &   ' 2nd half-sat. of DOC remineralization      =', xkdoc2
156         WRITE(numout,*)      &
157            &   ' excretion ratio of diatoms                 =', excret2
158         WRITE(numout,*)      &
159            &   ' exsudation rate of mesozooplankton         =', resrat2
160         WRITE(numout,*)      &
161            &   ' Diatoms mortality rate                     =', mprat2
162         WRITE(numout,*)      &
163            &   ' Phytoplankton minimum mortality rate       =', mpratm
164         WRITE(numout,*)      &
165            &   ' mesozooplankton mortality rate             =', mzrat2
166         WRITE(numout,*)      &
167            &   ' zoo preference for zoo                     =', xprefz
168         WRITE(numout,*)      &
169            &   ' zoo preference for poc                   =', xprefpoc
170         WRITE(numout,*)      &
171            &   ' maximal mesozoo grazing rate               =', grazrat2
172         WRITE(numout,*)      &
173            &   ' non assimilated fraction of P by mesozoo   =', unass2
174         WRITE(numout,*)      &
175            &   ' Efficicency of Mesozoo growth              =', epsher2
176         WRITE(numout,*)      &
177            &   ' Efficiency of microzoo growth              =', epsher
178         WRITE(numout,*)      &
179            &   ' half sturation constant for grazing 2      =', xkgraz2
180         WRITE(numout,*)      &
181            &   ' Maximum aggregation rate for diatoms       =', wchld
182         WRITE(numout,*)      &
183            &   ' scavenging rate of Iron                    =', xlam1
184         WRITE(numout,*)      &
185            &   ' Fe/C in zooplankton                        =', ferat3
186         WRITE(numout,*)      &
187            &   ' Phosphate half saturation for diatoms      =', conc1
188         WRITE(numout,*)      &
189            &   ' Iron half saturation for phyto             =', conc2
190         WRITE(numout,*)      &
191            &   ' Iron half saturation for diatoms           =', conc3
192         WRITE(numout,*)      &
193            &   ' NH4 half saturation for phyto              =', concnnh4
194         WRITE(numout,*)      &
195            &   ' NH4 half saturation for diatoms            =', concdnh4
196         WRITE(numout,*)      &
197            &   ' NH4 nitrification rate                     =', nitrif
198         WRITE(numout,*)      &
199            &   ' P-I slope  for diatoms                     =', pislope2
200         WRITE(numout,*)      &
201            &   ' Big particles sinking speed                =', wsbio2
202         WRITE(numout,*)      &
203            &   ' Fraction of microzoo excretion as DOM      =', sigma1
204         WRITE(numout,*)      &
205            &   ' Fraction of mesozoo excretion as DOM       =', sigma2
206         WRITE(numout,*)      &
207            &   ' Microzoo preference for POM                =', zprefc
208         WRITE(numout,*)      &
209            &   ' Microzoo preference for Nanophyto          =', zprefp
210         WRITE(numout,*)      &
211            &   ' Microzoo preference for Diatoms          =', zprefd
212         WRITE(numout,*)      &
213            &   ' Minimum Chl/C in nanophytoplankton         =', chlcnm
214         WRITE(numout,*)      &
215            &   ' Minimum Chl/C in diatoms                   =', chlcdm
216         WRITE(numout,*)      &
217            &   ' Maximum Fe/C in nanophytoplankton          =', fecnm
218         WRITE(numout,*)      &
219            &   ' Minimum Fe/C in diatoms                    =', fecdm
220         WRITE(numout,*)      &
221            &   ' Coastal release of Iron                 =', sedfeinput
222      ENDIF
223
224      READ(numnat,natsms)
225      IF(lwp) THEN
226         WRITE(numout,*) ' '
227         WRITE(numout,*) 'natsms'
228         WRITE(numout,*) ' '
229         WRITE(numout,*) 'Dust input from the atmosphere : ', bdustfer
230         WRITE(numout,*) ' '
231         WRITE(numout,*) 'River input of nutrients : ', briver
232         WRITE(numout,*) ' '
233         WRITE(numout,*) 'Atmospheric deposition of N : ', bndepo
234         WRITE(numout,*) ' '
235         WRITE(numout,*) 'Fe input from sediments : ', bsedinput
236         WRITE(numout,*) ' '
237      ENDIF
238
239#if defined key_trc_kriest
240
241      READ(numnat,natkriest)
242      IF(lwp) THEN
243         WRITE(numout,*) ' '
244         WRITE(numout,*) 'natkriest'
245         WRITE(numout,*) ' '
246         WRITE(numout,*) 'Sinking  exponent                    xkr_eta       = ', xkr_eta
247         WRITE(numout,*) 'N content exponent                   xkr_zeta      = ', xkr_zeta
248         WRITE(numout,*) 'Sinking factor                       xkr_sfact     = ', xkr_sfact
249         WRITE(numout,*) 'Stickiness                           xkr_stick     = ', xkr_stick
250         WRITE(numout,*) 'Minimum mass for Aggregates          xkr_mass_min  = ', xkr_mass_min
251         WRITE(numout,*) 'Maximum mass for Aggregates          xkr_mass_max  = ', xkr_mass_max
252         WRITE(numout,*) 'Size of particles in nano pool       xkr_dnano     = ', xkr_dnano
253         WRITE(numout,*) 'Size of particles in diatoms pool    xkr_ddiat     = ', xkr_ddiat
254         WRITE(numout,*) 'Size of particles in mesozoo pool    xkr_dmeso     = ', xkr_dmeso
255         WRITE(numout,*) 'Size of particles in aggregates pool xkr_daggr     = ', xkr_daggr
256     ENDIF
257
258
259     ! Computation of some variables
260     xkr_massp     = 5.7E-6 * 7.6 * xkr_mass_min**xkr_zeta
261     ! max and min vertical particle speed
262     xkr_wsbio_min = xkr_sfact * xkr_mass_min**xkr_eta
263     xkr_wsbio_max = xkr_sfact * xkr_mass_max**xkr_eta
264     WRITE(numout,*) ' max and min vertical particle speed ',xkr_wsbio_min,xkr_wsbio_max
265     !
266     !    effect of the sizes of the different living pools on particle numbers
267     !    nano = 2um-20um -> mean size=6.32 um -> ws=2.596 -> xnum=xnnano=2.337
268     !    diat and microzoo = 10um-200um -> 44.7 -> 8.732 -> xnum=xndiat=3.718
269     !    mesozoo = 200um-2mm -> 632.45 -> 45.14 -> xnum=xnmeso=7.147
270     !    aggregates = 200um-10mm -> 1414 -> 74.34 -> xnum=xnaggr=9.877
271     !    doc aggregates = 1um
272     ! ----------------------------------------------------------
273
274     xkr_nnano = 1. / ( xkr_massp * xkr_dnano )
275     xkr_ndiat = 1. / ( xkr_massp * xkr_ddiat )
276     xkr_nmeso = 1. / ( xkr_massp * xkr_dmeso )
277     xkr_naggr = 1. / ( xkr_massp * xkr_daggr )
278   
279
280#endif
281
282   END SUBROUTINE trc_lsm
Note: See TracBrowser for help on using the repository browser.