/[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.9 by jmc, Tue Apr 28 23:27:24 2009 UTC revision 1.12 by jmc, Tue Nov 4 17:15:52 2014 UTC
# Line 70  C   KPO4        :: half saturation phosp Line 70  C   KPO4        :: half saturation phosp
70  C   KFE         :: half saturation fe constant (mol/m3)  C   KFE         :: half saturation fe constant (mol/m3)
71  CC Iron chemisty values  CC Iron chemisty values
72  C   alpfe       :: solubility of aeolian fe  C   alpfe       :: solubility of aeolian fe
73    C   fesedflux_pcm :: ratio of sediment iron to sinking organic matter
74    C   FeIntSec    :: y-axis crossing for Fe_flux = fesedflux_pcm*pflux + FeIntSec
75  C   freefemax   :: max solubility of free iron (mol/m3)  C   freefemax   :: max solubility of free iron (mol/m3)
76  CC Control variables  CC Control variables
77  C   KScav       :: iron scavenging rate QQ  C   KScav       :: iron scavenging rate QQ
# Line 83  C                  read in rainRatioUnif Line 85  C                  read in rainRatioUnif
85        NAMELIST /BIOTIC_PARMS/        NAMELIST /BIOTIC_PARMS/
86       & DOPfraction, KDOPRemin, KRemin, zcrit,       & DOPfraction, KDOPRemin, KRemin, zcrit,
87       & O2crit, R_OP, R_CP, R_NP, R_FeP, zca,       & O2crit, R_OP, R_CP, R_NP, R_FeP, zca,
88       & parfrac, k0, lit0, KPO4, KFE,       & parfrac, k0, lit0, KPO4, KFE, kchl,
89       & alpfe, freefemax,       & alpfe, fesedflux_pcm, FeIntSec, freefemax,
90       & KScav, ligand_stab, ligand_tot,       & KScav, ligand_stab, ligand_tot,
91       & alphaUniform, rainRatioUniform       & alphaUniform, rainRatioUniform
92  #endif  #endif
93    
94        NAMELIST /DIC_FORCING/        NAMELIST /DIC_FORCING/
95       &          DIC_windFile, DIC_atmospFile, DIC_iceFile,       &          DIC_windFile, DIC_atmospFile, DIC_iceFile,
96       &          DIC_ironFile, DIC_silicaFile,       &          DIC_ironFile, DIC_silicaFile, DIC_parFile,
97         &          DIC_chlaFile,
98       &          DIC_forcingPeriod, DIC_forcingCycle,       &          DIC_forcingPeriod, DIC_forcingCycle,
99       &          dic_int1, dic_int2, dic_int3, dic_int4, dic_pCO2       &          dic_int1, dic_int2, dic_int3, dic_int4, dic_pCO2
100    
# Line 114  C---+----1----+----2----+----3----+----4 Line 117  C---+----1----+----2----+----3----+----4
117         zca         = 3500. _d 0         zca         = 3500. _d 0
118         parfrac     = 0.4 _d 0         parfrac     = 0.4 _d 0
119         k0          = 0.02 _d 0         k0          = 0.02 _d 0
120           kchl        = 0.02 _d 0
121         lit0        = 30. _d 0         lit0        = 30. _d 0
122         KPO4        = 5. _d -4         KPO4        = 5. _d -4
123         KFE         = 1.2 _d -7         KFE         = 1.2 _d -7
124         alpfe       = 0.01 _d 0         alpfe       = 0.01 _d 0
125           fesedflux_pcm = 6.8 _d -4 * 106. _d 0
126           FeIntSec    = 0.5 _d -6 / 86400. _d 0
127         freefemax   = 3. _d -7         freefemax   = 3. _d -7
128         KScav       = 0.19 _d 0/(360. _d 0*86400. _d 0)         KScav       = 0.19 _d 0/(360. _d 0*86400. _d 0)
129         ligand_stab = 1. _d 8         ligand_stab = 1. _d 8
# Line 130  C---+----1----+----2----+----3----+----4 Line 136  C---+----1----+----2----+----3----+----4
136         DIC_iceFile   = ' '         DIC_iceFile   = ' '
137         DIC_ironFile  = ' '         DIC_ironFile  = ' '
138         DIC_silicaFile= ' '         DIC_silicaFile= ' '
139           DIC_parFile   = ' '
140           DIC_chlaFile  = ' '
141         dic_int1    = 0         dic_int1    = 0
142         dic_int2    = 0         dic_int2    = 0
143         dic_int3    = 0         dic_int3    = 0
144         dic_int4    = 0         dic_int4    = 0
145         dic_pCO2    = 0. _d 0         dic_pCO2    = 278. _d -6
146  c default periodic forcing to same as for physics  c default periodic forcing to same as for physics
147         DIC_forcingPeriod = externForcingPeriod         DIC_forcingPeriod = externForcingPeriod
148         DIC_forcingCycle  = externForcingCycle         DIC_forcingCycle  = externForcingCycle
# Line 226  C- namelist BIOTIC_PARMS Line 234  C- namelist BIOTIC_PARMS
234         CALL WRITE_0D_RL( parfrac, INDEX_NONE,'parfrac =',         CALL WRITE_0D_RL( parfrac, INDEX_NONE,'parfrac =',
235       &  ' /* Fraction of Qsw that is PAR */')       &  ' /* Fraction of Qsw that is PAR */')
236         CALL WRITE_0D_RL( k0, INDEX_NONE,'k0 =',         CALL WRITE_0D_RL( k0, INDEX_NONE,'k0 =',
237       &  ' /* Light attentuation coefficient (1/m) */')       &  ' /* Light attentuation coefficient, water (1/m) */')
238           CALL WRITE_0D_RL( kchl, INDEX_NONE,'kchl =',
239         &  ' /* Light attentuation coefficient, chlorophyll (m2/mg) */')
240         CALL WRITE_0D_RL( lit0, INDEX_NONE,'lit0 =',         CALL WRITE_0D_RL( lit0, INDEX_NONE,'lit0 =',
241       &  ' /* Half saturation light constant (W/m2) */')       &  ' /* Half saturation light constant (W/m2) */')
242         CALL WRITE_0D_RL( KPO4, INDEX_NONE,'KPO4 =',         CALL WRITE_0D_RL( KPO4, INDEX_NONE,'KPO4 =',
# Line 235  C- namelist BIOTIC_PARMS Line 245  C- namelist BIOTIC_PARMS
245       &  ' /* Half saturation fe constant (mol/m3) */')       &  ' /* Half saturation fe constant (mol/m3) */')
246         CALL WRITE_0D_RL( alpfe, INDEX_NONE,'alpfe =',         CALL WRITE_0D_RL( alpfe, INDEX_NONE,'alpfe =',
247       &  ' /* Solubility of aeolian fe */')       &  ' /* Solubility of aeolian fe */')
248           CALL WRITE_0D_RL( fesedflux_pcm, INDEX_NONE,'fesedflux_pcm =',
249         &  ' /* Sediment Fe flux = fesedflux_pcm*pflux+FeIntSec */')
250           CALL WRITE_0D_RL( FeIntSec, INDEX_NONE,'FeIntSec =',
251         &  ' /* Sediment Fe flux = fesedflux_pcm * pflux + FeIntSec */')
252         CALL WRITE_0D_RL( freefemax, INDEX_NONE,'freefemax =',         CALL WRITE_0D_RL( freefemax, INDEX_NONE,'freefemax =',
253       &  ' /* Max solubility of free iron (mol/m3) */')       &  ' /* Max solubility of free iron (mol/m3) */')
254         CALL WRITE_0D_RL( KScav, INDEX_NONE,'KScav =',         CALL WRITE_0D_RL( KScav, INDEX_NONE,'KScav =',
# Line 263  C- namelist DIC_FORCING Line 277  C- namelist DIC_FORCING
277       & '  /* File name of aeolian iron flux */')       & '  /* File name of aeolian iron flux */')
278         CALL WRITE_0D_C( DIC_silicaFile, -1,INDEX_NONE,'DIC_silicaFile=',         CALL WRITE_0D_C( DIC_silicaFile, -1,INDEX_NONE,'DIC_silicaFile=',
279       & '  /* File name of surface silica */')       & '  /* File name of surface silica */')
280           CALL WRITE_0D_C( DIC_parFile, -1,INDEX_NONE,'DIC_parFile=',
281         & '  /* File name of photosynthetically available radiation */')
282           CALL WRITE_0D_C( DIC_chlaFile, -1,INDEX_NONE,'DIC_chlaFile=',
283         & '  /* File name of chlorophyll climatology */')
284         CALL WRITE_0D_RL( DIC_forcingPeriod,         CALL WRITE_0D_RL( DIC_forcingPeriod,
285       &   INDEX_NONE,'DIC_forcingPeriod =',       &   INDEX_NONE,'DIC_forcingPeriod =',
286       &  ' /* Periodic forcing parameter specific for DIC (s) */')       &  ' /* Periodic forcing parameter specific for DIC (s) */')
# Line 282  C- namelist DIC_FORCING Line 300  C- namelist DIC_FORCING
300    
301  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
302    
303          IF ( dic_int1.EQ.0 .AND. dic_pCO2.NE.278. _d -6 ) THEN
304            WRITE(msgBuf,'(A)')
305         &    'DIC_READPARMS: cannot change default dic_pCO2 if dic_int1=0'
306            CALL PRINT_ERROR( msgBuf, myThid )
307            STOP 'ABNORMAL END: S/R DIC_READPARMS: dic_pCO2 error'
308          ENDIF
309    
310        _END_MASTER(myThid)        _END_MASTER(myThid)
311    
312  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.9  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22