/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_check_chl.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/monod/monod_check_chl.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 jahn, Fri Jan 25 21:59:38 2013 UTC
# Line 88  c if spectral use wavebands to give a si Line 88  c if spectral use wavebands to give a si
88    
89               DO k=1,nR               DO k=1,nR
90                if (HFacC(i,j,k,bi,bj).gt.0. _d 0) then                if (HFacC(i,j,k,bi,bj).gt.0. _d 0) then
91                   atten= atten + (k0 )*5. _d -1*drF(k)
92                   if (k.gt.1)then
93                    atten = atten + (k0)*5. _d -1*drF(k-1)
94                   endif
95                   PARlocal = lite*exp(-atten)
96                 DO np=1,npmax                 DO np=1,npmax
97                  chltmp=Ptracer(i,j,k,bi,bj,iChl+np-1)                  chltmp=Ptracer(i,j,k,bi,bj,iChl+np-1)
98                  phytmp=Ptracer(i,j,k,bi,bj,iPhy+np-1)                  phytmp=Ptracer(i,j,k,bi,bj,iPhy+np-1)
# Line 95  c if initialized to zero, set to balance Line 100  c if initialized to zero, set to balance
100  c with no nutrient or temperature limitation  c with no nutrient or temperature limitation
101  c (Eq. A15 Geider et al 1997)  c (Eq. A15 Geider et al 1997)
102                  if (chltmp.eq.0. _d 0) then                  if (chltmp.eq.0. _d 0) then
                      atten= atten + (k0 )*5. _d -1*drF(k)  
                      if (k.gt.1)then  
                       atten = atten + (k0)*5. _d -1*drF(k-1)  
                      endif  
                      PARlocal = lite*exp(-atten)  
103                       pcm=pcmax(np)  !*limit(np)*phytoTempFunction(np)                       pcm=pcmax(np)  !*limit(np)*phytoTempFunction(np)
104                       chl2ctmp=chl2cmax(np)/                       chl2ctmp=chl2cmax(np)/
105       &                       (1+(chl2cmax(np)*alphachl(np)*PARlocal)/       &                       (1+(chl2cmax(np)*alphachl(np)*PARlocal)/
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22