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

Diff of /MITgcm_contrib/darwin2/pkg/monod/monod_plankton.F

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

revision 1.2 by stephd, Wed Oct 5 20:37:23 2011 UTC revision 1.3 by stephd, Wed Oct 5 20:43:39 2011 UTC
# Line 280  c  variables for zooplankton grazing rat Line 280  c  variables for zooplankton grazing rat
280  c  variables for zooplankton grazing rates  c  variables for zooplankton grazing rates
281           _RL zooTempFunction(nzmax)           _RL zooTempFunction(nzmax)
282           _RL allphyto(nzmax)           _RL allphyto(nzmax)
283             _RL denphyto(nzmax)
284           _RL grazphy(npmax,nzmax)           _RL grazphy(npmax,nzmax)
285           _RL sumgrazphy(npmax)           _RL sumgrazphy(npmax)
286           _RL sumgrazzoo(nzmax)           _RL sumgrazzoo(nzmax)
# Line 951  c sum all palatability*phyto and find ph Line 952  c sum all palatability*phyto and find ph
952              allphyto(nz)=allphyto(nz)+palat(np,nz)*phyto(np)              allphyto(nz)=allphyto(nz)+palat(np,nz)*phyto(np)
953             enddo             enddo
954             if (allphyto(nz).le.0. _d 0) allphyto(nz)=phygrazmin             if (allphyto(nz).le.0. _d 0) allphyto(nz)=phygrazmin
955    #ifdef SER_GRAZ
956               denphyto(nz)=0. _d 0
957               do np=1,npmax
958                denphyto(nz)=denphyto(nz)+
959         &            (palat(np,nz)*phyto(np)/allphyto(nz))*phyto(np)
960               enddo
961               if (denphyto(nz).le.0. _d 0) denphyto(nz)=phygrazmin
962    #endif
963             do np=1,npmax             do np=1,npmax
964              tmpz=max(0. _d 0,(allphyto(nz)-phygrazmin) )              tmpz=max(0. _d 0,(allphyto(nz)-phygrazmin) )
965              grazphy(np,nz)=grazemax(nz)*              grazphy(np,nz)=grazemax(nz)*
966    #ifdef SER_GRAZ
967    c as in Vallina et al, 2011
968         &          (((palat(np,nz)*phyto(np)/allphyto(nz))*phyto(np))/
969         &           denphyto(nz)) *
970    #else
971    c as in Dutkiewicz et al, GBC, 2009
972       &           (palat(np,nz)*phyto(np)/allphyto(nz))*       &           (palat(np,nz)*phyto(np)/allphyto(nz))*
973    #endif
974       &           ( tmpz**hollexp/       &           ( tmpz**hollexp/
975       &             (tmpz**hollexp+kgrazesat**hollexp) )       &             (tmpz**hollexp+kgrazesat**hollexp) )
976             enddo             enddo
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22