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.npzd.h90 in trunk/NEMO/TOP_SRC/SMS – NEMO

source: trunk/NEMO/TOP_SRC/SMS/trclsm.npzd.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:keywords set to Author Date Id Revision
File size: 8.1 KB
Line 
1! $Id$
2!!!
3!!!                       trclsm.npzd.h
4!!!                       *************
5!!!
6!!!  PURPOSE :
7!!!  ---------
8!!!     READs and PRINT options for NPZD namelist
9!!!
10!!   MODIFICATIONS:
11!!   --------------
12!!      original  : 99-10 (M.A. Foujols, M. Levy) passive tracer
13!!      additions : 01-01 (E. Kestenare) add sediments
14!!----------------------------------------------------------------------
15!! local declarations
16!! ==================
17
18#if defined key_passivetrc && defined key_trc_npzd
19      INTEGER :: ji
20      CHARACTER (len=32) :: clname
21
22
23!!!---------------------------------------------------------------------
24!!!  OPA8, LODYC (15/11/96)
25!!!---------------------------------------------------------------------
26
27! 0. initializations
28! ------------------
29
30      NAMELIST/natbio/apmin,azmin,anmin,admin,      &
31      &   redf,slopet,toptp,aknut,rcchl,            &
32      &   rgamma,toptgz,tmaxgz,rgz,                 &
33      &   rppz,taus,aks,filmax,rpnaz,rdnaz,eggzoo,tauzn,   &
34      &   tmmaxp,tmminp,tmmaxz,tmminz,anumin,afdmin,taudn, &
35      &   vsed,tmumax,aki,tmaxr,tminr,sedlam,sedlostpoc
36      NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig
37#if defined key_trc_diabio
38      NAMELIST/natdbi/ctrbio,ctrbil,ctrbiu,nwritebio
39#endif
40
41
42      IF(lwp) THEN
43          WRITE(numout,*) ' '
44          WRITE(numout,*) ' ROUTINE trclec'
45          WRITE(numout,*) ' **************'
46          WRITE(numout,*) ' '
47          WRITE(numout,*) ' namelist for NPZD model'
48          WRITE(numout,*) ' ***********************'
49          WRITE(numout,*) ' '
50      ENDIF
51
52      numnat=80
53      clname ='namelist.trc.sms'
54      OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old')
55
56! 1.4 namelist natbio : biological parameters
57! -------------------------------------------
58
59      apmin = 0
60      azmin = 0
61      anmin = 0
62      admin = 0
63      redf  = 0
64      slopet= 0
65      toptp = 0
66      aknut = 0
67      rcchl = 0
68      rgamma= 0
69      toptgz= 0
70      tmaxgz= 0
71      rgz   = 0
72      rppz  = 0
73      taus  = 0
74      aks   = 0
75      filmax= 0
76      rpnaz = 0
77      rdnaz = 0
78      eggzoo= 0
79      tauzn = 0
80      tmmaxp= 0
81      tmminp= 0
82      tmmaxz= 0
83      tmminz= 0
84      anumin= 0
85      afdmin= 0
86      taudn = 0
87      vsed  = 0
88      tmumax= 0
89      aki   = 0
90      tmaxr   = 1./(     4.*rday)*0.
91      tminr   = 1./(24.*30.*rday)*0.
92      sedlam=0
93      sedlostpoc=0
94
95      READ(numnat,natbio)
96
97      IF(lwp) THEN
98          WRITE(numout,*) 'natbio'
99          WRITE(numout,*) ' '
100          WRITE(numout,*)   &
101          &   ' minimum phytoplancton concentration  apmin =', apmin
102          WRITE(numout,*)   &
103          &   ' minimum zooplancton   concentration  azmin =', azmin
104          WRITE(numout,*)   &
105          &   ' minimum nutrients     concentration  anmin =', anmin
106          WRITE(numout,*)   &
107          &   ' minimum detritus      concentration  admin =', admin
108          WRITE(numout,*)   &
109          &   ' redfield ratio  c:n                   redf =', redf
110          WRITE(numout,*)   &
111          &   ' van t hoff coefficient              slopet =', slopet
112          WRITE(numout,*)   &
113          &   ' optimal photosynthesis temperature   toptp =', toptp
114          WRITE(numout,*)   &
115          &   ' half-saturation nutrient             aknut =', aknut
116          WRITE(numout,*)   &
117          &   ' carbone/chlorophyl ratio             rcchl =', rcchl
118          WRITE(numout,*)   &
119          &   ' phytoplankton exudation fraction    rgamma =', rgamma
120          WRITE(numout,*)   &
121          &   ' optimal temperature for zoo growth  toptgz =', toptgz
122          WRITE(numout,*)   &
123          &   ' maximal temperature for zoo growth  tmaxgz =', tmaxgz
124          WRITE(numout,*)   &
125          &   ' widtht of zoo temperature FUNCTION     rgz =', rgz
126          WRITE(numout,*)   &
127          &   ' zoo preference for phyto              rppz =', rppz
128          WRITE(numout,*)   &
129          &   ' maximal zoo grazing rate              taus =',86400*taus
130          WRITE(numout,*)   &
131          &   ' half saturation constant for zoo food  aks =', aks
132          WRITE(numout,*)   &
133          &   ' maximal mass clearance rate for zoo filmax =', filmax
134          WRITE(numout,*)   &
135          &   ' non-assimilated phyto by zoo         rpnaz =', rpnaz
136          WRITE(numout,*)   &
137          &   ' non-assimilated detritus by zoo      rdnaz =', rdnaz
138          WRITE(numout,*)   &
139          &   ' minimum  for zoo concentration      eggzoo =', eggzoo
140          WRITE(numout,*)   &
141          &   ' zoo specific excretion rate          tauzn =',86400   &
142          &   *tauzn
143          WRITE(numout,*)   &
144          &   ' maximal phyto mortality rate        tmmaxp =',86400   &
145          &   *tmmaxp
146          WRITE(numout,*)   &
147          &   ' minimal phyto mortality rate        tmminp =',86400   &
148          &   *tmminp
149          WRITE(numout,*)   &
150          &   ' maximal zoo mortality rate          tmmaxz =',86400   &
151          &   *tmmaxz
152          WRITE(numout,*)   &
153          &   ' minimal zoo mortality rate          tmminz =',86400   &
154          &   *tmminz
155          WRITE(numout,*)   &
156          &   ' nutrient threshold for phyto mort   anumin =', anumin
157          WRITE(numout,*)   &
158          &   ' food threshold for zoo mort         afdmin =', afdmin
159          WRITE(numout,*)   &
160          &   ' detrital breakdown rate              taudn =',86400   &
161          &   *taudn
162          WRITE(numout,*)   &
163          &   ' detritus sedimentation speed          vsed =',86400*vsed
164          WRITE(numout,*)   &
165          &   ' phyto max growth rate               tmumax =',86400   &
166          &   *tmumax
167          WRITE(numout,*)   &
168          &   ' light hlaf saturation constant         aki =', aki
169          WRITE(numout,*)   &
170          &   ' maximum damping for d z or p         tmaxr =', tmaxr
171          WRITE(numout,*)   &
172          &   ' minimum damping for d z or p         tminr =', tminr
173          WRITE(numout,*)   &
174          &   ' time coeff of POC in sediments      sedlam =', sedlam
175          WRITE(numout,*)   &
176          &   ' Sediment geol loss for POC  sedlostpoc =', sedlostpoc
177      ENDIF
178
179! 1.5 namelist natopt : parameters for optic
180! ------------------------------------------
181
182      xkg0  = 0.
183      xkr0  = 0.
184      xkgp  = 0.
185      xkrp  = 0.
186      xlg   = 0.
187      xlr   = 0.
188      rpig  = 0.
189
190      READ(numnat,natopt)
191
192      IF(lwp) THEN
193          WRITE(numout,*) 'natopt'
194          WRITE(numout,*) ' '
195          WRITE(numout,*) ' green   water absorption coeff  xkg0  = '  &
196          &   ,xkg0
197          WRITE(numout,*) ' red water absorption coeff      xkr0  = '  &
198          &   ,xkr0
199          WRITE(numout,*) ' pigment red absorption coeff    xkrp  = '  &
200          &   ,xkrp
201          WRITE(numout,*) ' pigment green absorption coeff  xkgp  = '  &
202          &   ,xkgp
203          WRITE(numout,*) ' green chl exposant              xlg   = '  &
204          &   ,xlg
205          WRITE(numout,*) ' red   chl exposant              xlr   = '  &
206          &   ,xlr
207          WRITE(numout,*) ' chla/chla+phea ratio            rpig  = '  &
208          &   ,rpig
209          WRITE(numout,*) ' '
210      ENDIF
211
212#if defined key_trc_diabio
213
214! NAMELIST : natdbi
215
216! default name for biological trends : short and long name, units
217
218      DO ji=1,jpdiabio
219        IF (ji < 10) THEN
220            WRITE (ctrbio(ji),'("BIO_",I1)') ji
221        ELSE IF (ji < 100) THEN
222            WRITE (ctrbio(ji),'("BIO_",I2)') ji
223        ELSE
224            WRITE (ctrbio(ji),'("BIO_",I3)') ji
225        ENDIF
226        WRITE (ctrbil(ji),'("BIOLOGICAL TREND NUMBER ",I2)') ji
227        ctrbiu(ji)='mmoleN/m3/s'
228      END DO
229
230      nwritebio = 10
231
232      READ(numnat,natdbi)
233
234      IF(lwp) THEN
235          WRITE(numout,*) 'natdbi'
236          WRITE(numout,*) ' '
237          WRITE(numout,*)        &
238          &   ' frequency of outputs for biological outputs = '  &
239          &   ,nwritebio
240          WRITE(numout,*) ' '
241          DO ji=1,jpdiabio
242            WRITE(numout,*)    &
243            &   'name of biological trend number :',ji,' : ',ctrbio(ji) 
244            WRITE(numout,*) ctrbil(ji) 
245            WRITE(numout,*) ' in unit = ',ctrbiu(ji)
246          END DO
247      END IF
248
249#endif
250
251#else
252
253! no passive tracers
254
255#endif
Note: See TracBrowser for help on using the repository browser.