/[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.6 by stephd, Tue Jul 24 16:01:56 2012 UTC
# Line 307  c sinking rate for particulate matter (m Line 307  c sinking rate for particulate matter (m
307           Pa2Atm      = 1.01325 _d 5           Pa2Atm      = 1.01325 _d 5
308  #endif  #endif
309    
310    #ifdef ALLOW_CDOM
311             fraccdom=2. _d 0 / 100. _d 0
312             cdomdegrd= 1. _d 0 / (200 _d 0 *pday)
313             cdombleach = 1. _d 0 / (15 _d 0 *pday)
314             PARcdom  = 20. _d 0
315             rnp_cdom = 16. _d 0
316             rfep_cdom = 1. _d -3
317             rcp_cdom  = 120. _d 0
318             cdomcoeff = .1 _d -1  / 1.d -4
319    #endif
320    
321  C make sure we have reserved enough space in Ptracers  C make sure we have reserved enough space in Ptracers
322        IF ( nCompZooMax .LT. 4 ) THEN        IF ( nCompZooMax .LT. 4 ) THEN
323          WRITE(msgBuf,'(A,A,I3)')          WRITE(msgBuf,'(A,A,I3)')
# Line 317  C make sure we have reserved enough spac Line 328  C make sure we have reserved enough spac
328        ENDIF        ENDIF
329        DO nz = 1,nzmax        DO nz = 1,nzmax
330          iZooP (nz)  = iZoo +                   (nz-1)*strideTypeZoo          iZooP (nz)  = iZoo +                   (nz-1)*strideTypeZoo
331    #ifdef ONLY_P_CYCLE
332            iZooN (nz)  = nptot
333            iZooFe(nz)  = nptot
334            iZooSi(nz)  = nptot
335    #else
336          iZooN (nz)  = iZoo + 1*strideCompZoo + (nz-1)*strideTypeZoo          iZooN (nz)  = iZoo + 1*strideCompZoo + (nz-1)*strideTypeZoo
337          iZooFe(nz)  = iZoo + 2*strideCompZoo + (nz-1)*strideTypeZoo          iZooFe(nz)  = iZoo + 2*strideCompZoo + (nz-1)*strideTypeZoo
338          iZooSi(nz)  = iZoo + 3*strideCompZoo + (nz-1)*strideTypeZoo          iZooSi(nz)  = iZoo + 3*strideCompZoo + (nz-1)*strideTypeZoo
339    #endif
340        ENDDO        ENDDO
341  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
342         DO nz = 1,nzmax         DO nz = 1,nzmax
# Line 338  c diver3 - fraction of biomass to count Line 355  c diver3 - fraction of biomass to count
355          diver_thresh3=.999 _d 0          diver_thresh3=.999 _d 0
356  c diver4 - fraction of maximum species  c diver4 - fraction of maximum species
357          diver_thresh4=1 _d -5          diver_thresh4=1 _d -5
358    c threshold on total biomass for contributing to Shannon and Simpson ind
359    c (these become large at very small biomass)
360            shannon_thresh = 1 _d -8   ! mmol P m-3
361  #endif  #endif
362    
363  c set up diagnostics  c set up diagnostics
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22