/[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.1 by jahn, Wed Apr 13 18:56:25 2011 UTC revision 1.4 by stephd, Fri Oct 21 20:56:53 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       &           ( tmpz/  #endif
974       &             (tmpz+kgrazesat) )       &           ( tmpz**hollexp/
975         &             (tmpz**hollexp+kgrazesat**hollexp) )
976             enddo             enddo
977            enddo            enddo
978            if (debug.eq.2) print*,'allphyto',allphyto            if (debug.eq.2) print*,'allphyto',allphyto
# Line 1107  cccccccccccccccccccccccccccccccccccccccc Line 1123  cccccccccccccccccccccccccccccccccccccccc
1123  c accumulate particulate and dissolved detritus  c accumulate particulate and dissolved detritus
1124              do nz=1, nzmax              do nz=1, nzmax
1125                 totzoo_pop=totzoo_pop+                 totzoo_pop=totzoo_pop+
1126       &                     ExportFracZ(nz)*mortzoo(nz)*zooP(nz)       &                     ExportFracZ(nz)*( mortzoo(nz)*zooP(nz)
1127         &                     +  mortzoo2(nz)*zooP(nz)**2 )
1128                 totzoo_dop=totzoo_dop+                 totzoo_dop=totzoo_dop+
1129       &                    (1. _d 0-ExportFracZ(nz))*mortzoo(nz)*zooP(nz)       &                    (1. _d 0-ExportFracZ(nz))*(
1130                 totzoo_pon=totzoo_pon+       &                     mortzoo(nz)*zooP(nz)+
1131       &                     ExportFracZ(nz)*mortzoo(nz)*zooN(nz)       &                     mortzoo2(nz)*zooP(nz)**2 )
1132                   totzoo_pon=totzoo_pon+
1133         &                     ExportFracZ(nz)*( mortzoo(nz)*zooN(nz)
1134         &                     +  mortzoo2(nz)*zooN(nz)**2 )
1135                 totzoo_don=totzoo_don+                 totzoo_don=totzoo_don+
1136       &                    (1. _d 0-ExportFracZ(nz))*mortzoo(nz)*zooN(nz)       &                    (1. _d 0-ExportFracZ(nz))*(
1137         &                     mortzoo(nz)*zooN(nz)+
1138         &                     mortzoo2(nz)*zooN(nz)**2 )
1139                 totzoo_pofe=totzoo_pofe+                 totzoo_pofe=totzoo_pofe+
1140       &                     ExportFracZ(nz)*mortzoo(nz)*zooFe(nz)       &                     ExportFracZ(nz)*( mortzoo(nz)*zooFe(nz)
1141         &                     +  mortzoo2(nz)*zooFe(nz)**2 )
1142                 totzoo_dofe=totzoo_dofe+                 totzoo_dofe=totzoo_dofe+
1143       &                   (1. _d 0-ExportFracZ(nz))*mortzoo(nz)*zooFe(nz)       &                   (1. _d 0-ExportFracZ(nz))*(
1144                 totzoo_posi=totzoo_posi+       &                    mortzoo(nz)*zooFe(nz) +
1145       &                     mortzoo(nz)*zooSi(nz)       &                    mortzoo2(nz)*zooFe(nz)**2 )
1146                   totzoo_posi=totzoo_posi+
1147         &                     ( mortzoo(nz)*zooSi(nz)+
1148         &                       mortzoo2(nz)*zooSi(nz)**2 )
1149  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
1150                 totzoo_poc=totzoo_poc+                 totzoo_poc=totzoo_poc+
1151       &                    ExportFracZ(nz)*mortzoo(nz)*zooClocal(nz)       &                    ExportFracZ(nz)*( mortzoo(nz)*zooClocal(nz)
1152         &                      +  mortzoo2(nz)*zooClocal(nz)**2 )
1153                 totzoo_doc=totzoo_doc+                 totzoo_doc=totzoo_doc+
1154       &           (1. _d 0-ExportFracZ(nz))*mortzoo(nz)*zooClocal(nz)       &           (1. _d 0-ExportFracZ(nz))*( mortzoo(nz)*zooClocal(nz)
1155         &                      +  mortzoo2(nz)*zooClocal(nz)**2 )
1156  #endif  #endif
1157              enddo              enddo
1158    
# Line 1420  c zooplankton mortality Line 1448  c zooplankton mortality
1448  c zoo in P currency  c zoo in P currency
1449                dzooPdt(nz)  = dzooPdt(nz)                dzooPdt(nz)  = dzooPdt(nz)
1450       &                           - mortzoo(nz)*zooP(nz)       &                           - mortzoo(nz)*zooP(nz)
1451         &                           - mortzoo2(nz)*zooP(nz)**2
1452  c zooplankton in other currencies  c zooplankton in other currencies
1453  C zooplankton stoichiometry varies according to food source  C zooplankton stoichiometry varies according to food source
1454                dzooNdt(nz)  = dzooNdt(nz)                dzooNdt(nz)  = dzooNdt(nz)
1455       &                           - mortzoo(nz)*zooN(nz)       &                           - mortzoo(nz)*zooN(nz)
1456         &                           - mortzoo2(nz)*zooN(nz)**2
1457                dzooFedt(nz) = dzooFedt(nz)                dzooFedt(nz) = dzooFedt(nz)
1458       &                           - mortzoo(nz)*zooFe(nz)       &                           - mortzoo(nz)*zooFe(nz)
1459         &                           - mortzoo2(nz)*zooFe(nz)**2
1460                dzooSidt(nz) = dzooSidt(nz)                dzooSidt(nz) = dzooSidt(nz)
1461       &                           - mortzoo(nz)*zooSi(nz)       &                           - mortzoo(nz)*zooSi(nz)
1462         &                           - mortzoo2(nz)*zooSi(nz)**2
1463              enddo              enddo
1464    
1465    
# Line 1486  c loss of O2 by remineralization Line 1518  c loss of O2 by remineralization
1518              do nz=1,nzmax              do nz=1,nzmax
1519                dzooCdt(nz)  = grazingC(nz)*zooClocal(nz)                dzooCdt(nz)  = grazingC(nz)*zooClocal(nz)
1520       &                           - mortzoo(nz)*zooClocal(nz)       &                           - mortzoo(nz)*zooClocal(nz)
1521         &                           - mortzoo2(nz)*zooClocal(nz)**2
1522              enddo              enddo
1523  #else  #else
1524              do nz=1,nzmax              do nz=1,nzmax
1525                dzooCdt(nz)  = sumgrazzooc(nz)                dzooCdt(nz)  = sumgrazzooc(nz)
1526       &                           - mortzoo(nz)*zooClocal(nz)       &                           - mortzoo(nz)*zooClocal(nz)
1527         &                           - mortzoo2(nz)*zooClocal(nz)**2
1528              enddo              enddo
1529  #endif  #endif
1530    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22