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