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

source: tags/nemo_v1_04/NEMO/TOP_SRC/SMS/trclsm.hamocc3.h90 @ 8023

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