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 @ 433

Last change on this file since 433 was 341, checked in by opalod, 19 years ago

nemo_v1_update_028 : CT : add missing headers
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 10.3 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/natgas/ gasfac, igaswind, icice
37      NAMELIST/natext/ atcco2
38      NAMELIST/natbio/caco3r,                                          &
39         &          dispo0,conc0,oxymin,grosip, sedlam,              &
40         &          sedlostpoc,sedlostcal,sedlostsil,nrdttrc,        &
41         &          pislope, excret,wsbio,wchl,wchld,resrat,mprat,mzrat,   &
42         &          grazrat,xprefc,xprefp,unass,xkgraz,xkmort,xksi1, &
43         &          xksi2,xremip,xremik,xsirem,xkdoc1,xkdoc2, &
44         &          excret2,resrat2,mprat2,mpratm,mzrat2,grazrat2,   &
45         &          xprefz,xprefpoc,unass2,xkgraz2,xlam1,      &
46         &          ferat3,conc1,conc2,conc3,concnnh4,concdnh4,      &
47         &          nitrif,epsher,epsher2,pislope2,wsbio2,sigma1,    &
48         &          sigma2, zprefc, zprefp, zprefd,fecnm,fecdm,      &
49         &          chlcnm,chlcdm, sedfeinput
50      NAMELIST/natsms/bdustfer, briver, bndepo, bsedinput
51
52      ! initialize the number of LOGICAL UNIT used
53      ! ------------------------------------------
54
55      IF(lwp) THEN
56         WRITE(numout,*) ' '
57         WRITE(numout,*) ' ROUTINE trclec'
58         WRITE(numout,*) ' **************'
59         WRITE(numout,*) ' '
60         WRITE(numout,*) ' namelist for PISCES model'
61         WRITE(numout,*) ' ***********************'
62         WRITE(numout,*) ' '
63      ENDIF
64
65      numnat=80
66      clname ='namelist.trc.sms'
67      OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old')
68
69
70      ! 1 Namelist natgas :
71      ! -------------------
72
73      READ(numnat,natgas)
74
75      IF(lwp) THEN
76         WRITE(numout,*) ' '
77         WRITE(numout,*) 'natgas'
78         write(numout,*) 'gasfac = ',gasfac
79         WRITE(numout,*) ' '
80         write(numout,*) 'igaswind = ',igaswind
81         WRITE(numout,*) ' '
82         write(numout,*) 'icice = ',icice
83         WRITE(numout,*) ' '
84      ENDIF
85
86      ! 2 Namelist natext :
87      ! -------------------
88      READ(numnat,natext)
89
90      IF(lwp) THEN
91         WRITE(numout,*) ' '
92         WRITE(numout,*) 'natext'
93         WRITE(numout,*) ' '
94         WRITE(numout,*) 'atmospheric pCO2= ',atcco2
95         WRITE(numout,*) ' '
96      ENDIF
97
98
99      READ(numnat,natbio)
100      IF(lwp) THEN
101         WRITE(numout,*) 'natbio'
102         WRITE(numout,*) ' '
103         WRITE(numout,*)      &
104            &   ' mean rainratio                             =', caco3r
105         WRITE(numout,*)      &
106            &   ' mean Si/C ratio                            =', grosip
107         WRITE(numout,*)      &
108            &   ' Calcite dissolution half saturation        =', dispo0
109         WRITE(numout,*)      &
110            &   ' Phosphate half saturation                  =', conc0
111         WRITE(numout,*)      &
112            &   ' Sediment bioturbation factor               =', sedlam
113         WRITE(numout,*)      &
114            &   ' Sediment burying ratio for POC         =', sedlostpoc
115         WRITE(numout,*)      &
116            &   ' Sediment burying ratio for CACO3       =', sedlostcal
117         WRITE(numout,*)      &
118            &   ' Sediment burying ratio for SI          =', sedlostsil
119         WRITE(numout,*)      &
120            &   ' frequence pour la biologie                 =', nrdttrc
121         WRITE(numout,*)      &
122            &   ' P-I slope                                  =', pislope
123         WRITE(numout,*)      &
124            &   ' excretion ratio of phytoplankton           =', excret
125         WRITE(numout,*)      &
126            &   ' POC sinking speed                          =', wsbio
127         WRITE(numout,*)      &
128            &   ' exsudation rate of zooplankton             =', resrat
129         WRITE(numout,*)      &
130            &   ' phytoplankton mortality rate               =', mprat
131         WRITE(numout,*)      &
132            &   ' zooplankton mortality rate                 =', mzrat
133         WRITE(numout,*)      &
134            &   ' zoo preference for phyto                   =', xprefc
135         WRITE(numout,*)      &
136            &   ' zoo preference for POC                     =', xprefp
137         WRITE(numout,*)      &
138            &   ' maximal zoo grazing rate                   =', grazrat
139         WRITE(numout,*)      &
140            &   ' non assimilated fraction of phyto by zoo   =', unass
141         WRITE(numout,*)      &
142            &   ' half sturation constant for grazing        =', xkgraz
143         WRITE(numout,*)      &
144            &   ' half saturation constant for mortality     =', xkmort
145         WRITE(numout,*)      &
146            &   ' half saturation constant for Si uptake     =', xksi1
147         WRITE(numout,*)      &
148            &   ' half saturation constant for Si/C          =', xksi2
149         WRITE(numout,*)      &
150            &   ' remineralisation rate of POC               =', xremip
151         WRITE(numout,*)      &
152            &   ' remineralization rate of DOC               =', xremik
153         WRITE(numout,*)      &
154            &   ' remineralization rate of Si                =', xsirem
155         WRITE(numout,*)      &
156            &   ' 1st half-sat. of DOC remineralization      =', xkdoc1
157         WRITE(numout,*)      &
158            &   ' 2nd half-sat. of DOC remineralization      =', xkdoc2
159         WRITE(numout,*)      &
160            &   ' excretion ratio of diatoms                 =', excret2
161         WRITE(numout,*)      &
162            &   ' exsudation rate of mesozooplankton         =', resrat2
163         WRITE(numout,*)      &
164            &   ' Diatoms mortality rate                     =', mprat2
165         WRITE(numout,*)      &
166            &   ' Phytoplankton minimum mortality rate       =', mpratm
167         WRITE(numout,*)      &
168            &   ' mesozooplankton mortality rate             =', mzrat2
169         WRITE(numout,*)      &
170            &   ' zoo preference for zoo                     =', xprefz
171         WRITE(numout,*)      &
172            &   ' zoo preference for poc                   =', xprefpoc
173         WRITE(numout,*)      &
174            &   ' maximal mesozoo grazing rate               =', grazrat2
175         WRITE(numout,*)      &
176            &   ' non assimilated fraction of P by mesozoo   =', unass2
177         WRITE(numout,*)      &
178            &   ' Efficicency of Mesozoo growth              =', epsher2
179         WRITE(numout,*)      &
180            &   ' Efficiency of microzoo growth              =', epsher
181         WRITE(numout,*)      &
182            &   ' half sturation constant for grazing 2      =', xkgraz2
183         WRITE(numout,*)      &
184            &   ' Maximum aggregation rate for diatoms       =', wchld
185         WRITE(numout,*)      &
186            &   ' scavenging rate of Iron                    =', xlam1
187         WRITE(numout,*)      &
188            &   ' Fe/C in zooplankton                        =', ferat3
189         WRITE(numout,*)      &
190            &   ' Phosphate half saturation for diatoms      =', conc1
191         WRITE(numout,*)      &
192            &   ' Iron half saturation for phyto             =', conc2
193         WRITE(numout,*)      &
194            &   ' Iron half saturation for diatoms           =', conc3
195         WRITE(numout,*)      &
196            &   ' NH4 half saturation for phyto              =', concnnh4
197         WRITE(numout,*)      &
198            &   ' NH4 half saturation for diatoms            =', concdnh4
199         WRITE(numout,*)      &
200            &   ' NH4 nitrification rate                     =', nitrif
201         WRITE(numout,*)      &
202            &   ' P-I slope  for diatoms                     =', pislope2
203         WRITE(numout,*)      &
204            &   ' Big particles sinking speed                =', wsbio2
205         WRITE(numout,*)      &
206            &   ' Fraction of microzoo excretion as DOM      =', sigma1
207         WRITE(numout,*)      &
208            &   ' Fraction of mesozoo excretion as DOM       =', sigma2
209         WRITE(numout,*)      &
210            &   ' Microzoo preference for POM                =', zprefc
211         WRITE(numout,*)      &
212            &   ' Microzoo preference for Nanophyto          =', zprefp
213         WRITE(numout,*)      &
214            &   ' Microzoo preference for Diatoms          =', zprefd
215         WRITE(numout,*)      &
216            &   ' Minimum Chl/C in nanophytoplankton         =', chlcnm
217         WRITE(numout,*)      &
218            &   ' Minimum Chl/C in diatoms                   =', chlcdm
219         WRITE(numout,*)      &
220            &   ' Maximum Fe/C in nanophytoplankton          =', fecnm
221         WRITE(numout,*)      &
222            &   ' Minimum Fe/C in diatoms                    =', fecdm
223         WRITE(numout,*)      &
224            &   ' Coastal release of Iron                 =', sedfeinput
225      ENDIF
226
227      READ(numnat,natsms)
228      IF(lwp) THEN
229         WRITE(numout,*) ' '
230         WRITE(numout,*) 'natsms'
231         WRITE(numout,*) ' '
232         WRITE(numout,*) 'Dust input from the atmosphere : ', bdustfer
233         WRITE(numout,*) ' '
234         WRITE(numout,*) 'River input of nutrients : ', briver
235         WRITE(numout,*) ' '
236         WRITE(numout,*) 'Atmospheric deposition of N : ', bndepo
237         WRITE(numout,*) ' '
238         WRITE(numout,*) 'Fe input from sediments : ', bsedinput
239         WRITE(numout,*) ' '
240      ENDIF
241
242   END SUBROUTINE trc_lsm
Note: See TracBrowser for help on using the repository browser.