/[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.3 by benw, Tue May 19 14:32:43 2015 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 43  c Line 47  c
47         _RL felimit(npmax)         _RL felimit(npmax)
48         _RL alpha_fe         _RL alpha_fe
49         _RL up_inorg(iimax,npmax)         _RL up_inorg(iimax,npmax)
        _RL N_uptake(npmax)  
50         _RL photo_Tempfunction         _RL photo_Tempfunction
51  c  c
52         _RL dchldt(npmax)         _RL dchldt(npmax)
# Line 51  c Line 54  c
54         _RL E0         _RL E0
55         _RL C_biomass         _RL C_biomass
56         _RL chl         _RL chl
        _RL VCref  
57  c  c
58         _RL Chl2C         _RL Chl2C
59         _RL PCmax         _RL PCmax
60         _RL PCPhot(npmax)         _RL PCPhot(npmax)
61         _RL rhochl         _RL rhochl
62         _RL chlsynth(npmax)         _RL chlsynth(npmax)
        _RL VCmax  
63         _RL VCN         _RL VCN
64           _RL PP
65    #ifdef QUOTA_DIAG_LIMIT
66           _RL Ilim(npmax)
67    #endif
68  c  c
69  c Geider, MacIntyre and Kana (1998) photosynthesis model  c Geider, MacIntyre and Kana (1998) photosynthesis model
70  c adapted for multiple nutrients following Moore et al (2002)  c adapted for multiple nutrients following Moore et al (2002)
71  c  c
72         E0 = PARlocal ! muEin m^-2 s^-1         E0 = PARlocal ! muEin m^-2 s^-1
73           PP = 0.0 _d 0
74         do jp=1,npmax         do jp=1,npmax
75  #ifdef FQUOTA  #ifdef FQUOTA
76           alpha_fe = alphachl(jp) * felimit(jp)           alpha_fe = alphachl(jp) * felimit(jp)
# Line 96  c Line 102  c
102  c  c
103               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
104                 PCPhot(jp) = PCmax                      ! s^-1                 PCPhot(jp) = PCmax                      ! s^-1
105       &          *(1-exp(-alpha_fe*Chl2C*E0/PCmax))       &          *(1. _d 0 - exp(-alpha_fe*Chl2C*E0/PCmax))
106  c  c
107                 rhochl = Chl2Nmax * PCPhot(jp)          ! mg chl (mmol N)^-1                 rhochl = Chl2Nmax * PCPhot(jp)          ! mg chl (mmol N)^-1
108       &                /(alphachl(jp)*Chl2C*E0)       &                /(alphachl(jp)*Chl2C*E0)
# Line 112  c Line 118  c
118             endif             endif
119  c-----------------------------------------------------------------  c-----------------------------------------------------------------
120             ! for passing back to quota_plankton             ! for passing back to quota_plankton
121               PP = PP + PCPhot(jp) * C_biomass
122             up_inorg(iDIC,jp) = PCPhot(jp) - biosynth*VCN  ! s^-1             up_inorg(iDIC,jp) = PCPhot(jp) - biosynth*VCN  ! s^-1
123             dchldt(jp)        = chlsynth(jp)  ! mg chl m^-3 s^-1             dchldt(jp)        = chlsynth(jp)  ! mg chl m^-3 s^-1
124           else ! else if jp is a heterotroph           else ! else if jp is a heterotroph
# Line 119  c--------------------------------------- Line 126  c---------------------------------------
126             dchldt(jp)        = 0. _d 0       ! mg chl m^-3 s^-1             dchldt(jp)        = 0. _d 0       ! mg chl m^-3 s^-1
127  c-----------------------------------------------------------------  c-----------------------------------------------------------------
128           endif           endif
129    #ifdef QUOTA_DIAG_LIMIT
130             Ilim(jp) = 1. _d 0 - exp(-alphachl(jp)*1.59*E0/vmaxi(iDIC,jp))
131    #endif
132    !      print*,"vmax",PP,PCPhot(jp),biomass(iCarb,jp),
133    !     &       E0,vmaxi(iDIC,jp),chl,C_biomass
134         enddo         enddo
135    
136         RETURN         RETURN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22