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

source: trunk/NEMO/TOP_SRC/SMS/trclsm.hamocc3.h90 @ 260

Last change on this file since 260 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: 6.9 KB
Line 
1!!! $Header$
2!!!
3!!!                       trclsm.hamocc3.h
4!!!                       ****************
5!!!
6!!!  PURPOSE :
7!!!  ---------
8!!!     READs and PRINT options for HAMOCC3 or P3ZD NAMELIST
9!!!
10!!   METHOD :                   : no
11!!   -------
12!!
13!!   MODIFICATIONS:
14!!   --------------
15!!      original  : 99-10 (M.A. Foujols, M. Levy) passive tracer
16!!      addition  : 00-01 (L. Bopp) hamocc3,p3zd
17!!     
18!!----------------------------------------------------------------------
19!! local declarations
20!! ==================
21
22#if defined key_passivetrc && defined key_trc_hamocc3
23      INTEGER :: ji
24      CHARACTER (len=32) :: clname
25
26!!!---------------------------------------------------------------------
27!!!  OPA8, LODYC (15/11/96)
28!!!---------------------------------------------------------------------
29
30! 0. initializations
31! ------------------
32
33       NAMELIST/natgas/ gasfac, igaswind, icice
34       NAMELIST/natext/ ilecvit,atcco2
35#    if defined key_trc_p3zd
36       NAMELIST/natbio/silic0,pendec,pendes,fluexp,caco3r,dispo0,conc0,   &
37       &               oxymin,grosip,jpkb,sedlam,sedlostpoc,sedlostsil,   &
38       &               sedlostcal,pislope,betslope,ekc,ekw1,ekw2,excret,  &
39       &               epsbio,wsbio,wchl,resrat,mprat,mzrat,grazrat,      &
40       &               xprefc,xprefp,unass,xkgraz,xkmort,xksi1,xksi2,sicmax, &
41       &               xremip,xremik,smax,xkdoc1,xkdoc2
42#else
43       NAMELIST/natbio/silic0,pendec,pendes,fluexp,caco3r,dispo0,conc0,   &
44       &               oxymin,grosip,jpkb,sedlam,sedlostpoc,sedlostsil,   &
45       &               sedlostcal
46#endif
47
48
49#    if defined key_trc_p3zd
50      IF(lwp) THEN
51          WRITE(numout,*) ' '
52          WRITE(numout,*) ' ROUTINE trclec'
53          WRITE(numout,*) ' **************'
54          WRITE(numout,*) ' '
55          WRITE(numout,*) ' NAMELIST for P3ZD model'
56          WRITE(numout,*) ' ***********************'
57          WRITE(numout,*) ' '
58      ENDIF
59#    else
60      IF(lwp) THEN
61          WRITE(numout,*) ' '
62          WRITE(numout,*) ' ROUTINE trclec'
63          WRITE(numout,*) ' **************'
64          WRITE(numout,*) ' '
65          WRITE(numout,*) ' NAMELIST for HAMOCC3 model'
66          WRITE(numout,*) ' ***********************'
67          WRITE(numout,*) ' '
68      ENDIF
69#    endif
70
71      numnat=80
72      clname ='namelist.trc.sms'
73      OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old')
74
75! 1 Namelist natgas :
76! -------------------
77
78      READ(numnat,natgas)
79
80      IF(lwp) THEN
81          WRITE(numout,*) ' '
82          WRITE(numout,*) 'natgas'
83          write(numout,*) 'gasfac = ',gasfac
84          WRITE(numout,*) ' '
85          write(numout,*) 'igaswind = ',igaswind
86          WRITE(numout,*) ' '
87          write(numout,*) 'icice = ',icice
88          WRITE(numout,*) ' '
89
90      ENDIF
91
92! 2 Namelist natext :
93! -------------------
94
95      READ(numnat,natext)
96
97      IF(lwp) THEN
98          WRITE(numout,*) ' '
99          WRITE(numout,*) 'natext'
100          write(numout,*)'ilecvit = ',ilecvit
101          WRITE(numout,*) ' '
102          WRITE(numout,*) 'atmospheric pCO2= ',atcco2
103          WRITE(numout,*) ' '
104      ENDIF
105
106      READ(numnat,natbio)
107      IF(lwp) THEN
108          WRITE(numout,*) 'natbio'
109          WRITE(numout,*) ' '
110          WRITE(numout,*)    &
111          &   ' Silicon half saturation constant           =', silic0
112          WRITE(numout,*)    &
113          &   ' export length for calcite                  =', pendec
114          WRITE(numout,*)    &
115          &   ' export length for biogenic silicon         =', pendes
116          WRITE(numout,*)    &
117          &   ' power constant for POC export              =', fluexp
118          WRITE(numout,*)    &
119          &   ' mean rainratio                             =', caco3r
120          WRITE(numout,*)    &
121          &   ' Calcite dissolution half saturation        =', dispo0
122          WRITE(numout,*)    &
123          &   ' Phosphate half saturation                  =', conc0
124          WRITE(numout,*)    &
125          &   ' Productive zone vertical indice            =', jpkb
126          WRITE(numout,*)    &
127          &   ' Sediment bioturbation constant             =', sedlam
128          WRITE(numout,*)    &
129          &   ' Sediment geol loss for POC            =', sedlostpoc
130          WRITE(numout,*)    &
131          &   ' Sediment geol loss for SIL             =', sedlostsil
132          WRITE(numout,*)    &
133          &   ' Sediment geol loss for CAL             =', sedlostcal
134#    if defined key_trc_p3zd
135          WRITE(numout,*)    &
136          &   ' P-I slope                                  =', pislope
137          WRITE(numout,*)    &
138          &   ' Light inhibition factor                    =', betslope
139          WRITE(numout,*)    &
140          &   ' pigment light absorption coeff.            =', ekc
141          WRITE(numout,*)    &
142          &   ' green light absorption coeff. of water     =', ekw1
143          WRITE(numout,*)    &
144          &   ' red light absorption coeff. of water       =', ekw2
145          WRITE(numout,*)    &
146          &   ' excretion ratio of phytoplankton           =', excret
147          WRITE(numout,*)    &
148          &   ' fraction of dead zooplankton --> POC       =', epsbio
149          WRITE(numout,*)    &
150          &   ' POC sinking speed                          =', wsbio
151          WRITE(numout,*)    &
152          &   ' exsudation rate of zooplankton             =', resrat
153          WRITE(numout,*)    &
154          &   ' phytoplankton mortality rate               =', mprat
155          WRITE(numout,*)    &
156          &   ' zooplankton mortality rate                 =', mzrat
157          WRITE(numout,*)    &
158          &   ' zoo preference for phyto                   =', xprefc
159          WRITE(numout,*)    &
160          &   ' zoo preference for POC                     =', xprefp
161          WRITE(numout,*)    &
162          &   ' maximal zoo grazing rate                   =', grazrat
163          WRITE(numout,*)    &
164          &   ' non assimilated fraction of phyto by zoo   =', unass
165          WRITE(numout,*)    &
166          &   ' half saturation constant for grazing        =', xkgraz
167          WRITE(numout,*)    &
168          &   ' half saturation constant for mortality     =', xkmort
169          WRITE(numout,*)    &
170          &   ' half saturation constant for Si uptake     =', xksi1
171          WRITE(numout,*)    &
172          &   ' half saturation constant for Si/C          =', xksi2
173          WRITE(numout,*)    &
174          &   ' maximum Si/C                               =', sicmax
175          WRITE(numout,*)    &
176          &   ' remineralisation rate of POC               =', xremip
177          WRITE(numout,*)    &
178          &   ' remineralization rate of DOC               =', xremik
179          WRITE(numout,*)    &
180          &   ' maximum sinking rate of phyto              =', smax
181          WRITE(numout,*)    &
182          &   ' 1st half-sat. of DOC remineralization      =', xkdoc1
183          WRITE(numout,*)    &
184          &   ' 2nd half-sat. of DOC remineralization      =', xkdoc2
185          WRITE(numout,*) ' '
186#    endif
187      ENDIF
188
189#else
190
191! no passive tracers
192
193#endif
Note: See TracBrowser for help on using the repository browser.