/[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.11 by stephd, Fri Jun 8 20:31:10 2012 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
# Line 92  C                  read in rainRatioUnif Line 94  C                  read in rainRatioUnif
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_parFile,       &          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 131  C---+----1----+----2----+----3----+----4 Line 137  C---+----1----+----2----+----3----+----4
137         DIC_ironFile  = ' '         DIC_ironFile  = ' '
138         DIC_silicaFile= ' '         DIC_silicaFile= ' '
139         DIC_parFile   = ' '         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
# Line 227  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 236  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 266  C- namelist DIC_FORCING Line 279  C- namelist DIC_FORCING
279       & '  /* File name of surface silica */')       & '  /* File name of surface silica */')
280         CALL WRITE_0D_C( DIC_parFile, -1,INDEX_NONE,'DIC_parFile=',         CALL WRITE_0D_C( DIC_parFile, -1,INDEX_NONE,'DIC_parFile=',
281       & '  /* File name of photosynthetically available radiation */')       & '  /* 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) */')
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22