/[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.4 by jahn, Wed Nov 9 23:34:17 2011 UTC revision 1.7 by stephd, Thu Jul 26 18:01:22 2012 UTC
# Line 106  c depth for denitrification to start Line 106  c depth for denitrification to start
106    
107  c critical oxygen for O2/NO3 remineralization  c critical oxygen for O2/NO3 remineralization
108        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)
109  c ratio of no3 to p in denitrification process  c ratio of n to p in denitrification process
110        denit_np = 120.0 _d 0        denit_np = 120.0 _d 0
111    c ratio no3 used relative to all n in denitrification process
112          denit_no3 = 104.0 _d 0
113    c critical nitrate below which no denit (or remin) happens
114          no3crit = 1. _d -2
115  c  c
116  c oxidation rates for ammonium and nitrite  c oxidation rates for ammonium and nitrite
117  c i.e. Knita ...    NH4 -> NO2  c i.e. Knita ...    NH4 -> NO2
# Line 307  c sinking rate for particulate matter (m Line 310  c sinking rate for particulate matter (m
310           Pa2Atm      = 1.01325 _d 5           Pa2Atm      = 1.01325 _d 5
311  #endif  #endif
312    
313    #ifdef ALLOW_CDOM
314             fraccdom=2. _d 0 / 100. _d 0
315             cdomdegrd= 1. _d 0 / (200 _d 0 *pday)
316             cdombleach = 1. _d 0 / (15 _d 0 *pday)
317             PARcdom  = 20. _d 0
318             rnp_cdom = 16. _d 0
319             rfep_cdom = 1. _d -3
320             rcp_cdom  = 120. _d 0
321             cdomcoeff = .1 _d -1  / 1.d -4
322    #endif
323    
324  C make sure we have reserved enough space in Ptracers  C make sure we have reserved enough space in Ptracers
325        IF ( nCompZooMax .LT. 4 ) THEN        IF ( nCompZooMax .LT. 4 ) THEN
326          WRITE(msgBuf,'(A,A,I3)')          WRITE(msgBuf,'(A,A,I3)')
# Line 317  C make sure we have reserved enough spac Line 331  C make sure we have reserved enough spac
331        ENDIF        ENDIF
332        DO nz = 1,nzmax        DO nz = 1,nzmax
333          iZooP (nz)  = iZoo +                   (nz-1)*strideTypeZoo          iZooP (nz)  = iZoo +                   (nz-1)*strideTypeZoo
334    #ifdef ONLY_P_CYCLE
335            iZooN (nz)  = nptot
336            iZooFe(nz)  = nptot
337            iZooSi(nz)  = nptot
338    #else
339          iZooN (nz)  = iZoo + 1*strideCompZoo + (nz-1)*strideTypeZoo          iZooN (nz)  = iZoo + 1*strideCompZoo + (nz-1)*strideTypeZoo
340          iZooFe(nz)  = iZoo + 2*strideCompZoo + (nz-1)*strideTypeZoo          iZooFe(nz)  = iZoo + 2*strideCompZoo + (nz-1)*strideTypeZoo
341          iZooSi(nz)  = iZoo + 3*strideCompZoo + (nz-1)*strideTypeZoo          iZooSi(nz)  = iZoo + 3*strideCompZoo + (nz-1)*strideTypeZoo
342    #endif
343        ENDDO        ENDDO
344  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
345         DO nz = 1,nzmax         DO nz = 1,nzmax
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22