/[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.3 by stephd, Wed Oct 5 20:43:39 2011 UTC revision 1.4 by stephd, Fri Oct 21 20:56:53 2011 UTC
# Line 1123  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 1436  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 1502  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.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22