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

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

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

revision 1.1 by jahn, Wed Apr 13 18:56:25 2011 UTC revision 1.2 by benw, Mon Jul 2 09:44:24 2012 UTC
# Line 20  c ====================================== Line 20  c ======================================
20  #ifdef FQUOTA  #ifdef FQUOTA
21       I             felimit,       I             felimit,
22  #endif  #endif
23    #ifdef QUOTA_DIAG_LIMIT
24         O             Ilim,
25    #endif
26       I             up_inorg,       I             up_inorg,
27         O             PP,
28       I             photo_Tempfunction,       I             photo_Tempfunction,
29       O             dchldt, ! chlorophyll synthesis rate       O             dchldt, ! chlorophyll synthesis rate
30       I             myThid)       I             myThid)
# Line 60  c Line 64  c
64         _RL chlsynth(npmax)         _RL chlsynth(npmax)
65         _RL VCmax         _RL VCmax
66         _RL VCN         _RL VCN
67           _RL PP
68    #ifdef QUOTA_DIAG_LIMIT
69           _RL Ilim(npmax)
70    #endif
71  c  c
72  c Geider, MacIntyre and Kana (1998) photosynthesis model  c Geider, MacIntyre and Kana (1998) photosynthesis model
73  c adapted for multiple nutrients following Moore et al (2002)  c adapted for multiple nutrients following Moore et al (2002)
74  c  c
75         E0 = PARlocal ! muEin m^-2 s^-1         E0 = PARlocal ! muEin m^-2 s^-1
76           PP = 0.0 _d 0
77         do jp=1,npmax         do jp=1,npmax
78  #ifdef FQUOTA  #ifdef FQUOTA
79           alpha_fe = alphachl(jp) * felimit(jp)           alpha_fe = alphachl(jp) * felimit(jp)
# Line 96  c Line 105  c
105  c  c
106               if (PCmax.gt.0. _d 0.and.alpha_fe.gt.0. _d 0) then               if (PCmax.gt.0. _d 0.and.alpha_fe.gt.0. _d 0) then
107                 PCPhot(jp) = PCmax                      ! s^-1                 PCPhot(jp) = PCmax                      ! s^-1
108       &          *(1-exp(-alpha_fe*Chl2C*E0/PCmax))       &          *(1. _d 0-exp(-alpha_fe*Chl2C*E0/PCmax))
109    c
110    #ifdef QUOTA_DIAG_LIMIT
111                   Ilim(jp)=1. _d 0
112         &                 -exp(-alphachl(jp)*Chl2C*E0/vmaxi(iDIC,jp))
113    #endif
114  c  c
115                 rhochl = Chl2Nmax * PCPhot(jp)          ! mg chl (mmol N)^-1                 rhochl = Chl2Nmax * PCPhot(jp)          ! mg chl (mmol N)^-1
116       &                /(alphachl(jp)*Chl2C*E0)       &                /(alphachl(jp)*Chl2C*E0)
# Line 112  c Line 126  c
126             endif             endif
127  c-----------------------------------------------------------------  c-----------------------------------------------------------------
128             ! for passing back to quota_plankton             ! for passing back to quota_plankton
129               PP = PP + PCPhot(jp) * C_biomass
130             up_inorg(iDIC,jp) = PCPhot(jp) - biosynth*VCN  ! s^-1             up_inorg(iDIC,jp) = PCPhot(jp) - biosynth*VCN  ! s^-1
131             dchldt(jp)        = chlsynth(jp)  ! mg chl m^-3 s^-1             dchldt(jp)        = chlsynth(jp)  ! mg chl m^-3 s^-1
132           else ! else if jp is a heterotroph           else ! else if jp is a heterotroph
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22