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 trunk/NEMO/TOP_SRC/SMS – NEMO

source: trunk/NEMO/TOP_SRC/SMS/trclsm.pisces.h90 @ 262

Last change on this file since 262 was 260, checked in by opalod, 19 years ago

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