/[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.2 by benw, Mon Jul 2 09:44:24 2012 UTC revision 1.3 by benw, Tue May 19 14:32:43 2015 UTC
# Line 47  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 55  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         _RL PP
65  #ifdef QUOTA_DIAG_LIMIT  #ifdef QUOTA_DIAG_LIMIT
# Line 105  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. _d 0-exp(-alpha_fe*Chl2C*E0/PCmax))       &          *(1. _d 0 - exp(-alpha_fe*Chl2C*E0/PCmax))
 c  
 #ifdef QUOTA_DIAG_LIMIT  
                Ilim(jp)=1. _d 0  
      &                 -exp(-alphachl(jp)*Chl2C*E0/vmaxi(iDIC,jp))  
 #endif  
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 134  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.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22