/[MITgcm]/MITgcm/pkg/dic/dic_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/dic_readparms.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by jmc, Wed Apr 9 00:44:03 2008 UTC revision 1.4 by dfer, Wed Apr 9 16:22:54 2008 UTC
# Line 57  C             1 Atm = 1.01325e5 Pa = 101 Line 57  C             1 Atm = 1.01325e5 Pa = 101
57    
58  C-- Biotic dic parameters:  C-- Biotic dic parameters:
59  C   DOPfraction :: fraction of new production going to DOP  C   DOPfraction :: fraction of new production going to DOP
60  C   KDOPRemin   :: DOP remineralization rate (s) = 1/(6 month)  C   KDOPRemin   :: DOP remineralization rate (1/s) = 1/(6 month)
61  C   KRemin      :: remin power law coeff  C   KRemin      :: remin power law coeff
62  C   zcrit       :: Minimum Depth (m) over which biological activity  C   zcrit       :: Minimum Depth (m) over which biological activity
63  C                  is computed --> determines nlev as the indice of the  C                  is computed --> determines nlev as the indice of the
# Line 186  C taken into account Line 186  C taken into account
186  #endif  #endif
187  #endif /* DIC_BIOTIC */  #endif /* DIC_BIOTIC */
188    
189    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
190    C--   Print out parameter values :
191    
192          iUnit = standardMessageUnit
193          WRITE(msgBuf,'(A)') ' '
194          CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
195          WRITE(msgBuf,'(A)') '// ==================================='
196          CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
197          WRITE(msgBuf,'(A)') '// Dic package parameters :'
198          CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
199          WRITE(msgBuf,'(A)') '// ==================================='
200          CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
201    
202    C- namelist ABIOTIC_PARMS
203           CALL WRITE_0D_R8( permil, INDEX_NONE,'permil =',
204         &  ' /* Ref. density to convert mol/m3 to mol/kg */')
205           CALL WRITE_0D_R8( Pa2Atm, INDEX_NONE,'Pa2Atm =',
206         &  ' /* Atmosph. pressure conversion coeff (to Atm) */')
207    
208    C- namelist BIOTIC_PARMS
209           CALL WRITE_0D_R8( DOPfraction, INDEX_NONE,'DOPfraction =',
210         &  ' /* Fraction of new production going to DOP */')
211           CALL WRITE_0D_R8( KDOPRemin, INDEX_NONE,'KDOPRemin =',
212         &  ' /* DOP remineralization rate (1/s) */')
213           CALL WRITE_0D_R8( KRemin, INDEX_NONE,'KRemin =',
214         &  ' /* Remin power law coeff. */')
215           CALL WRITE_0D_R8( zcrit, INDEX_NONE,'zcrit =',
216         &  ' /* Minimum depth for biological activity (m) */')
217           CALL WRITE_0D_R8( O2crit, INDEX_NONE,'O2crit =',
218         &  ' /* Critical oxygen level (mol/m3) */')
219           CALL WRITE_0D_R8( R_OP, INDEX_NONE,'R_OP =',
220         &  ' /* Stochiometric ratio R_OP */')
221           CALL WRITE_0D_R8( R_CP, INDEX_NONE,'R_CP =',
222         &  ' /* Stochiometric ratio R_CP */')
223           CALL WRITE_0D_R8( R_NP, INDEX_NONE,'R_NP =',
224         &  ' /* Stochiometric ratio R_NP */')
225           CALL WRITE_0D_R8( R_FeP, INDEX_NONE,'R_FeP =',
226         &  ' /* Stochiometric ratio R_FeP */')
227           CALL WRITE_0D_R8( zca, INDEX_NONE,'zca =',
228         &  ' /* Scale depth for CaCO3 remineralization (m) */')
229           CALL WRITE_0D_R8( parfrac, INDEX_NONE,'parfrac =',
230         &  ' /* Fraction of Qsw that is PAR */')
231           CALL WRITE_0D_R8( k0, INDEX_NONE,'k0 =',
232         &  ' /* Light attentuation coefficient (1/m) */')
233           CALL WRITE_0D_R8( lit0, INDEX_NONE,'lit0 =',
234         &  ' /* Half saturation light constant (W/m2) */')
235           CALL WRITE_0D_R8( KPO4, INDEX_NONE,'KPO4 =',
236         &  ' /* Half saturation phosphate constant (mol/m3) */')
237           CALL WRITE_0D_R8( KFE, INDEX_NONE,'KFE =',
238         &  ' /* Half saturation fe constant (mol/m3) */')
239           CALL WRITE_0D_R8( alpfe, INDEX_NONE,'alpfe =',
240         &  ' /* Solubility of aeolian fe */')
241           CALL WRITE_0D_R8( freefemax, INDEX_NONE,'freefemax =',
242         &  ' /* Max solubility of free iron (mol/m3) */')
243           CALL WRITE_0D_R8( KScav, INDEX_NONE,'KScav =',
244         &  ' /* Iron scavenging rate */')
245           CALL WRITE_0D_R8( ligand_stab, INDEX_NONE,'ligand_stab =',
246         &  ' /* Ligand-free iron stability constant (m3/mol) */')
247           CALL WRITE_0D_R8( ligand_tot, INDEX_NONE,'ligand_tot =',
248         &  ' /* Total free ligand  (mol/m3) */')
249           CALL WRITE_0D_R8( alphaloc, INDEX_NONE,'alphaloc =',
250         &  ' /* Timescape for biological activity */')
251           CALL WRITE_0D_R8( rain_ratioloc, INDEX_NONE,'rain_ratioloc =',
252         &  ' /* Inorganic/organic carbon rain ratio */')
253           CALL WRITE_0D_C( WindFileDic, 0, INDEX_NONE, 'WindFileDic =',
254         & '  /* File name of wind speeds */')
255           CALL WRITE_0D_C( AtmospFileDic, 0, INDEX_NONE, 'AtmospFileDic =',
256         & '  /* File name of atmospheric pressure*/')
257           CALL WRITE_0D_C( IceFileDic, 0, INDEX_NONE, 'IceFileDic =',
258         & '  /* File name of seaice fraction */')
259           CALL WRITE_0D_C( IronFileDic, 0, INDEX_NONE, 'IronFileDic =',
260         & '  /* File name of aeolian iron flux */')
261           CALL WRITE_0D_C( SilicaFileDic, 0, INDEX_NONE, 'SilicaFileDic =',
262         & '  /* File name of surface silica */')
263           CALL WRITE_0D_R8( dic_ForcingPeriod,
264         &   INDEX_NONE,'dic_ForcingPeriod =',
265         &  ' /* Periodic forcing parameter specific for dic (s) */')
266           CALL WRITE_0D_R8( dic_ForcingCycle,
267         &   INDEX_NONE,'dic_ForcingCycle =',
268         &  ' /* Periodic forcing parameter specific for dic (s) */')
269           CALL WRITE_0D_I( dic_int1, INDEX_NONE, 'dic_int1 =',
270         &  '  /*  */')
271           CALL WRITE_0D_I( dic_int2, INDEX_NONE, 'dic_int2 =',
272         &  '  /*  */')
273           CALL WRITE_0D_I( dic_int3, INDEX_NONE, 'dic_int3 =',
274         &  '  /*  */')
275           CALL WRITE_0D_I( dic_int4, INDEX_NONE, 'dic_int4 =',
276         &  '  /*  */')
277           CALL WRITE_0D_R8( dic_pCO2, INDEX_NONE,'dic_pCO2 =',
278         &  ' /* Atmospheric pCO2 to be read in data.dic */')
279           CALL WRITE_0D_L( QSW_underice, INDEX_NONE, 'QSW_underice  =',
280         &  '  /* Flag for Qsw under Sea-Ice (i.e. SI fract included) */')
281    
282    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
283    
284        _END_MASTER(myThid)        _END_MASTER(myThid)
285    
286  C--   Everyone else must wait for the parameters to be loaded  C--   Everyone else must wait for the parameters to be loaded
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22