/[MITgcm]/MITgcm/pkg/bling/bling_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/bling/bling_readparms.F

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

revision 1.9 by jmc, Thu Feb 23 20:48:13 2017 UTC revision 1.10 by mmazloff, Thu Mar 16 17:03:26 2017 UTC
# Line 143  C ====================================== Line 143  C ======================================
143  #endif  #endif
144    
145  C ==========================================================  C ==========================================================
146  C     permil   :: set carbon mol/m3 <---> mol/kg conversion factor  C   secperday        :: seconds in a day = 24*60*60
147  C                 default permil = 1024.5 kg/m3  C   permil           :: set carbon mol/m3 <---> mol/kg conversion factor
148  C     Pa2Atm   :: Conversion factor for atmospheric pressure pLoad  C                       default permil = 1024.5 kg/m3
149  C                 (when coupled to atmospheric model) into Atm.  C   Pa2Atm           :: Conversion factor for atmospheric pressure pLoad
150  C                 Default assumes pLoad in Pascal  C                       (when coupled to atmospheric model) into Atm.  
151  C                 1 Atm = 1.01325e5 Pa = 1013.25 mb  C                       Default assumes pLoad in Pascal
152  C  C                       1 Atm = 1.01325e5 Pa = 1013.25 mb
153  C     CtoP            :: Carbon to phosphorus ratio in organic matter  C   pivotal          :: Pivotal phytoplankton biomass
154  C     NtoP            :: Nitrogen to phosphorus ratio in organic matter  C   Pc_0             :: Maximum phytoplankton carbon-specific growth rate at 0C
155  C     O2toP           :: Oxygen to phosphorus for biological activity  C   Pc_0_diaz        :: Maximum diazotroph carbon-specific growth rate at 0C
156  C     CatoP           :: Calcium to phosphorus uptake by small phyto  C   lambda_0         :: Carbon-specific phytoplankton mortality rate
157  C     NUTfac          :: If using nitrogen as macro-nutrient instead of ,  C   chl_min          :: minimum chlorophyll concentration [ug kg-1]
158  C                        phosphate, NUTfac is stoichiometric ratio N:P,  C   CtoN             :: Carbon to nitrogen ratio in organic matter
159  C                        otherwise equal to 1.  C   NO3toN           ::
160  C     Pstar           :: Pivotal phytoplankton biomass  C   HtoC             ::
161  C     Pc_0            :: Maximum carbon-specific growth rate at 0C  C   O2toN            :: Oxygen to nitrogen for biological activity
162  C     lambda_0        :: Carbon-specific phytoplankton mortality rate  C   CatoN            :: Calcium to nitrogen uptake by small phyto
163  C     chl_min         :: minimum chlorophyll concentration  C   masstoN          ::
164  C     theta_Fe_max_hi :: Maximum Chl:c ratio, abundant iron  C   alpha_photo      :: [g C g Chl-1 m2 W-1 s-1]
165  C     theta_Fe_max_lo :: Maximum Chl:c ratio, extreme iron limitation  C   theta_Fe_max_hi  :: Maximum Chl:c ratio, abundant iron
166  C     gamma_irr_mem   :: Photoadaptation time scale  C   theta_Fe_max_lo  :: Maximum Chl:c ratio, extreme iron limitation
167  C     gamma_DOM       :: Decay timescale of DOM  C   gamma_irr_mem    :: Photoadaptation time scale
168  C     gamma_POM       :: Remineralization rate of sinking POM  C   gamma_DON        :: Decay timescale of DON
169  C     k_Fe            :: Dissolved Fe uptake half-saturation constant  C   gamma_DOP        :: Decay timescale of DOP
170  C     k_O2            :: Half-saturation constant for aerobic respiration  C   gamma_POM        ::
171  C     k_NUT           :: Macro-nutrient uptake half-saturation constant  C   k_Fe             :: Dissolved Fe uptake half-saturation constant
172  C     k_FetoN         :: Half-saturation cellular Fe:N  C   k_Fe_diaz        :: Dissolved Fe uptake half-saturation constant for diazotrophs
173  C     FetoP_max       :: Maximum Fe:P uptake ratio  C   k_O2             :: Half-saturation constant for aerobic respiration
174  C     FetoP_sed       :: Fe:P in sediments  C   k_NO3            :: Nitrate uptake half-saturation constant
175  C     KFeLeq_max      :: Maximum Fe-ligand stability constant  C   k_PO4            :: Phosphate uptake half-saturation constant
176  C     KFeLeq_min      :: Minimum Fe-ligand stability constant  C   k_PtoN           :: Half-saturation cellular P:N
177  C     Fe_min          :: Constant having to do with photodissociation  C   k_FetoN          :: Half-saturation cellular Fe:N
178  C     IFeL            :: Iron ligand stability constant  C   kFe_eq_lig_max   :: Maximum Fe-ligand stability constant
179  C     kFe_org         :: Organic-matter dependent scavenging rate  C   kFe_eq_lig_min   :: Minimum Fe-ligand stability constant
180  C     kFe_inorg       :: Inorganic scavenging rate  C   kFe_eq_lig_Femin :: Constant having to do with photodissociation
181  C     remin_min       :: Minimum anaerobic respiration rate  C   kFe_eq_lig_irr   :: Iron ligand stability constant  
182  C     O2_min          :: Minimum O2 concentration for aerobic respiration  C   kFe_org          :: Organic-matter dependent scavenging rate
183  C     Ligand          :: Ligand concentration  C   kFe_inorg        :: Inorganic scavenging rate
184  C     kappa_eppley    :: Temperature dependence of growth  C   PtoN_min         :: Minimum P:N uptake ratio
185  C     kappa_remin     :: Temperature dependence of remineralization  C   PtoN_max         :: Maximum P:N uptake ratio
186  C     ca_remin_depth  :: CaCO3 remineralization lengthscale  C   FetoN_min        :: Minimum Fe:N uptake ratio
187  C     phi_DOM         :: Fraction of non-sinking production to DOM  C   FetoN_max        :: Maximum Fe:N uptake ratio
188  C     phi_sm          :: Fraction of small phytoplankton biomass converted to detritus  C   FetoC_sed        :: Fe:P in sediments
189  C     phi_lg          :: Fraction of large phytoplankton biomass converted to detritus  C   remin_min        :: Minimum anaerobic respiration rate
190  C     wsink0z         :: Depth at which sinking rate starts increasing  C   oxic_min         :: Minimum O2 concentration for aerobic respiration
191  C     wsink0          :: Initial sinking rate  C   ligand           :: Ligand concentration
192  C     wsinkacc        :: Acceleration rate of sinking with depth  C   kappa_eppley     :: Temperature dependence of growth
193  C     parfrac         :: fraction of Qsw avail for photosynthesis  C   kappa_eppley_diaz:: Temperature dependence of growth for diazotrophs
194  C     alpfe           :: solubility of aeolian iron  C   kappa_remin      :: Temperature dependence of remineralization
195  C     k0              :: Light attentuation coefficient  C   ca_remin_depth   :: CaCO3 remineralization lengthscale
196  C     epsln           :: a very small number  C   phi_DOM          :: Fraction of non-sinking production to DOM
197    C   phi_sm           :: Fraction of small phytoplankton biomass converted to detritus
198  C     secperday       :: seconds in a day = 24*60*60  C   phi_lg           :: Fraction of large phytoplankton biomass converted to detritus
199    C   phi_dvm          ::
200    C   sigma_dvm        ::
201    C   wsink0z          :: Depth at which sinking rate starts increasing  
202    C   wsink0           :: Initial sinking rate
203    C   wsinkacc         :: Acceleration rate of sinking with depth
204    C   parfrac          :: fraction of Qsw avail for photosynthesis
205    C   alpfe            :: solubility of aeolian iron
206    C   k0               :: Light attentuation coefficient
207    C   epsln            :: a very small number
208        
209        _RL secperday        _RL secperday
210        integer k        integer k
211  #ifdef USE_EXF_INTERPOLATION  #ifdef USE_EXF_INTERPOLATION
# Line 226  C     Default values Line 236  C     Default values
236        Pc_0_diaz            = 0.01 _d -5        Pc_0_diaz            = 0.01 _d -5
237        lambda_0             = 0.19 _d 0 / secperday        lambda_0             = 0.19 _d 0 / secperday
238        chl_min              = 1. _d -5        chl_min              = 1. _d -5
 Cxxx chl_min in ug kg-1  
239        alpha_photo          = 0.7 _d -5 * 2.77 _d 18 / 6.022 _d 17        alpha_photo          = 0.7 _d -5 * 2.77 _d 18 / 6.022 _d 17
 Cxxx alpha_photo in g C g Chl-1 m2 W-1 s-1  
240        theta_Fe_max_hi      = 0.04 _d 0        theta_Fe_max_hi      = 0.04 _d 0
241        theta_Fe_max_lo      = 0.01 _d 0        theta_Fe_max_lo      = 0.01 _d 0
242        gamma_irr_mem        = 1. _d 0 / secperday        gamma_irr_mem        = 1. _d 0 / secperday
243        gamma_DON            = 0.25 _d 0 / (365.25 _d 0 * secperday)        gamma_DON            = 0.25 _d 0 / (365.25 _d 0 * secperday)
244        gamma_DOP            = 0.5 _d 0 / (365.25 _d 0 * secperday)        gamma_DOP            = 0.5 _d 0 / (365.25 _d 0 * secperday)
245        gamma_POM            = 0.12 _d 0 / secperday        gamma_POM            = 0.12 _d 0 / secperday
246  cav      k_Fe                 = 1.6 _d -10 / permil        k_Fe                 = 1.6 _d -10 / permil
       k_Fe                 = 8. _d -10 / permil  
247        k_Fe_diaz            = 7. _d -10 / permil        k_Fe_diaz            = 7. _d -10 / permil
248        k_O2                 = 20. _d -6 / permil        k_O2                 = 20. _d -6 / permil
249        k_NO3                = 2. _d -6 / permil        k_NO3                = 2. _d -6 / permil
250  cav      k_PO4                = 1. _d -8 / permil        k_PO4                = 1. _d -8 / permil
       k_PO4                = 1. _d -7 / permil  
251        k_PtoN               = 1.5 _d -6 / permil        k_PtoN               = 1.5 _d -6 / permil
252        k_FetoN              = 8. _d -10 / permil        k_FetoN              = 8. _d -10 / permil
253        PtoN_min             = 1. / 28.        PtoN_min             = 1. / 28.
# Line 375  C taken into account Line 381  C taken into account
381         IF ( useThSIce ) QSW_underice = .TRUE.         IF ( useThSIce ) QSW_underice = .TRUE.
382  #endif  #endif
383    
 #ifdef USE_EXFCO2  
       IF ( useEXF .AND. apco2file .NE. ' ' ) THEN  
        CALL EXF_GETFFIELD_START( useExfYearlyFields,  
      I                    'bling', 'apco2', apco2period,  
      I                    apco2startdate1, apco2startdate2,  
      U                    apco2StartTime, errCount,  
      I                    myThid )  
       ENDIF  
 #endif /* USE_EXFCO2 */  
384        IF ( errCount.GE.1 ) THEN        IF ( errCount.GE.1 ) THEN
385         WRITE(msgBuf,'(A,I3,A)')         WRITE(msgBuf,'(A,I3,A)')
386       &     'BLING_READPARMS: detected', errCount,' fatal error(s)'       &     'BLING_READPARMS: detected', errCount,' fatal error(s)'
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22