/[MITgcm]/MITgcm_contrib/darwin2/pkg/quota/quota_plankton.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/quota/quota_plankton.F

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

revision 1.3 by benw, Tue May 19 14:32:43 2015 UTC revision 1.4 by benw, Tue May 19 15:23:46 2015 UTC
# Line 288  c             apply ammonium inhibition Line 288  c             apply ammonium inhibition
288                up_inorg(ii,jp) = 0. _d 0                up_inorg(ii,jp) = 0. _d 0
289              endif              endif
290            enddo ! ii            enddo ! ii
 #ifdef FQUOTA  
           up_inorg(iNO3,jp) = up_inorg(iNO3,jp) * felimit(jp)  
 #endif  
291          else ! if autotrophy(jp).eq.0          else ! if autotrophy(jp).eq.0
292            do ii=1,iimax            do ii=1,iimax
293              up_inorg(ii,jp) = 0. _d 0              up_inorg(ii,jp) = 0. _d 0
# Line 349  c PRE-ASSIMILATION grazing of type jpred Line 346  c PRE-ASSIMILATION grazing of type jpred
346              if (food1.gt.0. _d 0) then              if (food1.gt.0. _d 0) then
347                grazing(iCarb,jpred,jprey)                 ! d^-1                grazing(iCarb,jpred,jprey)                 ! d^-1
348       &          = tmp1 ! grazing effort       &          = tmp1 ! grazing effort
349  !#ifdef ONEGRAZER  #ifdef ONEGRAZER
350  !     &          * graz(jprey) ! prey dependent maximum rate       &          * graz(jprey) ! prey dependent maximum rate
351  !#else  #else
352       &          * graz(jpred) ! predator dependent maximum rate       &          * graz(jpred) ! predator dependent maximum rate
353  !#endif  #endif
354  #ifdef SWITCHING  #ifdef SWITCHING
355       &          *(graz_pref(jpred,jprey)*biomass(iCarb,jprey))**ns/food2       &          *(graz_pref(jpred,jprey)*biomass(iCarb,jprey))**ns/food2
356  #else  #else
# Line 478  c           Exudation of elemental reser Line 475  c           Exudation of elemental reser
475              endif              endif
476            endif            endif
477  !  !
478  ! calculate temperature adjusted mortality rates  ! calculate (temperature adjusted?) mortality rates
479            m_ref(jp) = kmort(jp) !* activ_Tempfunction            m_ref(jp) = kmort(jp) !* activ_Tempfunction
 #ifdef ALLOWPFT  
           ! Z mortality is t dependent - a la Moore 2002  
           if (pft(jp).eq.6) then  
             m_ref(jp) = kmort(jp) * activ_Tempfunction  
           endif  
 #endif  
480  !  !
481  ! Loss and sinking terms - include silicate and chlorophyll  ! Loss and sinking terms - include silicate and chlorophyll
482              dbiomassdt(io,jp) = dbiomassdt(io,jp)              dbiomassdt(io,jp) = dbiomassdt(io,jp)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22