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

source: tags/nemo_v1_04/NEMO/TOP_SRC/SMS/trclsm.npzd.h90 @ 280

Last change on this file since 280 was 274, checked in by opalod, 19 years ago

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