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

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