/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_init_fixed.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/monod/monod_init_fixed.F

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

revision 1.3 by stephd, Fri Oct 21 20:56:53 2011 UTC revision 1.11 by jahn, Thu May 1 16:19:32 2014 UTC
# Line 39  C     /--------------------------------- Line 39  C     /---------------------------------
39  C     |  initialise common block biochemical parameters               |  C     |  initialise common block biochemical parameters               |
40  C     \--------------------------------------------------------------/  C     \--------------------------------------------------------------/
41    
42          _BEGIN_MASTER(myThid)
43    
44        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
45       &'// ======================================================='       &'// ======================================================='
46        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
# Line 106  c depth for denitrification to start Line 108  c depth for denitrification to start
108    
109  c critical oxygen for O2/NO3 remineralization  c critical oxygen for O2/NO3 remineralization
110        O2crit = 6.0 _d 0    !(Lipschultz et al 1990, DSR 37, 1513-1541)        O2crit = 6.0 _d 0    !(Lipschultz et al 1990, DSR 37, 1513-1541)
111  c ratio of no3 to p in denitrification process  c ratio of n to p in denitrification process
112        denit_np = 120.0 _d 0        denit_np = 120.0 _d 0
113    c ratio no3 used relative to all n in denitrification process
114          denit_no3 = 104.0 _d 0
115    c critical nitrate below which no denit (or remin) happens
116          no3crit = 1. _d -2
117  c  c
118  c oxidation rates for ammonium and nitrite  c oxidation rates for ammonium and nitrite
119  c i.e. Knita ...    NH4 -> NO2  c i.e. Knita ...    NH4 -> NO2
# Line 218  c ammonia and nitrite inhibition Line 223  c ammonia and nitrite inhibition
223  c set si half sat  c set si half sat
224         val_ksatsi=1. _d 0         val_ksatsi=1. _d 0
225  c set nutrient ratios for phyto  c set nutrient ratios for phyto
226    #ifdef ONLY_P_CYCLE
227           val_R_SiP_diatom=0.0 _d 0
228           val_R_NP=0.0 _d 0
229           val_RFeP=0.0 _d -3
230           val_R_NP_diaz=0.0 _d  0
231           val_RFeP_diaz=0.0 _d 0 * val_RFeP
232           val_R_PC=0.0 _d 0
233           val_R_PICPOC=0.0 _d 0
234    #else
235         val_R_SiP_diatom=16.0 _d 0  ! 32 for Fanny's runs         val_R_SiP_diatom=16.0 _d 0  ! 32 for Fanny's runs
236         val_R_NP=16.0 _d 0         val_R_NP=16.0 _d 0
237         val_RFeP=1.0 _d -3         val_RFeP=1.0 _d -3
# Line 225  c set nutrient ratios for phyto Line 239  c set nutrient ratios for phyto
239         val_RFeP_diaz=30.0 _d 0 * val_RFeP         val_RFeP_diaz=30.0 _d 0 * val_RFeP
240         val_R_PC=120.0 _d 0         val_R_PC=120.0 _d 0
241         val_R_PICPOC=0.8 _d 0         val_R_PICPOC=0.8 _d 0
242    #endif
243    #ifdef FIX_ZOO_QUOTAS
244           val_R_NP_zoo=val_R_NP
245           val_RFeP_zoo=val_RFeP
246           val_R_SiP_zoo=0.0 _d 0
247           val_R_PC_zoo=val_R_PC
248    #endif
249    c
250  #ifdef OLD_GRAZE  #ifdef OLD_GRAZE
251  c grazing hlaf saturation  c grazing hlaf saturation
252         kgrazesat = 0.1 _d 0         kgrazesat = 0.1 _d 0
# Line 307  c sinking rate for particulate matter (m Line 329  c sinking rate for particulate matter (m
329           Pa2Atm      = 1.01325 _d 5           Pa2Atm      = 1.01325 _d 5
330  #endif  #endif
331    
332    #ifdef ALLOW_CDOM
333             fraccdom=2. _d 0 / 100. _d 0
334             cdomdegrd= 1. _d 0 / (200 _d 0 *pday)
335             cdombleach = 1. _d 0 / (15 _d 0 *pday)
336             PARcdom  = 20. _d 0
337             rnp_cdom = 16. _d 0
338             rfep_cdom = 1. _d -3
339             rcp_cdom  = 120. _d 0
340             cdomcoeff = .1 _d -1  / 1.d -4
341    #endif
342    
343  C make sure we have reserved enough space in Ptracers  C make sure we have reserved enough space in Ptracers
344        IF ( nCompZooMax .LT. 4 ) THEN        IF ( nCompZooMax .LT. 4 ) THEN
345    #ifndef ONLY_P_CYCLE
346          WRITE(msgBuf,'(A,A,I3)')          WRITE(msgBuf,'(A,A,I3)')
347       &    'MONOD_INIT_FIXED: ERROR: 4 zooplankton components, but ',       &    'MONOD_INIT_FIXED: ERROR: 4 zooplankton components, but ',
348       &    'nCompZooMax = ', nCompZooMax       &    'nCompZooMax = ', nCompZooMax
349          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , 1)
350          STOP 'ABNORMAL END: S/R MONOD_INIT_FIXED'          STOP 'ABNORMAL END: S/R MONOD_INIT_FIXED'
351    #endif
352        ENDIF        ENDIF
353        DO nz = 1,nzmax        DO nz = 1,nzmax
354          iZooP (nz)  = iZoo +                   (nz-1)*strideTypeZoo          iZooP (nz)  = iZoo +                   (nz-1)*strideTypeZoo
355    #ifdef ONLY_P_CYCLE
356            iZooN (nz)  = nptot
357            iZooFe(nz)  = nptot
358            iZooSi(nz)  = nptot
359    #else
360          iZooN (nz)  = iZoo + 1*strideCompZoo + (nz-1)*strideTypeZoo          iZooN (nz)  = iZoo + 1*strideCompZoo + (nz-1)*strideTypeZoo
361          iZooFe(nz)  = iZoo + 2*strideCompZoo + (nz-1)*strideTypeZoo          iZooFe(nz)  = iZoo + 2*strideCompZoo + (nz-1)*strideTypeZoo
362          iZooSi(nz)  = iZoo + 3*strideCompZoo + (nz-1)*strideTypeZoo          iZooSi(nz)  = iZoo + 3*strideCompZoo + (nz-1)*strideTypeZoo
363    #endif
364        ENDDO        ENDDO
365  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
366         DO nz = 1,nzmax         DO nz = 1,nzmax
# Line 338  c diver3 - fraction of biomass to count Line 379  c diver3 - fraction of biomass to count
379          diver_thresh3=.999 _d 0          diver_thresh3=.999 _d 0
380  c diver4 - fraction of maximum species  c diver4 - fraction of maximum species
381          diver_thresh4=1 _d -5          diver_thresh4=1 _d -5
382    c threshold on total biomass for contributing to Shannon and Simpson ind
383    c (these become large at very small biomass)
384            shannon_thresh = 1 _d -8   ! mmol P m-3
385  #endif  #endif
386    
387  c set up diagnostics  c set up diagnostics
# Line 361  COJ set up diagnostics Line 405  COJ set up diagnostics
405  #endif /* ALLOW_DIAGNOSTICS */  #endif /* ALLOW_DIAGNOSTICS */
406  COJ  COJ
407    
408          _END_MASTER(myThid)
409    
410          RETURN
411          RETURN        END
         END  
412  C============================================================================  C============================================================================
413  #endif  #endif
414  #endif  #endif
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22