/[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.13 by stephd, Thu Mar 28 18:14:17 2013 UTC revision 1.15 by jahn, Thu May 1 16:19:32 2014 UTC
# Line 1080  c sum over phy for impact on zoo, and al Line 1080  c sum over phy for impact on zoo, and al
1080             sumgrazlossPIC(nz)=0. _d 0             sumgrazlossPIC(nz)=0. _d 0
1081  #endif  #endif
1082             do np=1,npmax             do np=1,npmax
1083    #ifdef FIX_ZOO_QUOTAS
1084                  sumgrazzoo(nz)=sumgrazzoo(nz)+
1085         &             asseff(np,nz)*grazphy(np,nz)*zooP(nz)
1086                  sumgrazloss(nz)=sumgrazloss(nz)+
1087         &             (1. _d 0-asseff(np,nz))*grazphy(np,nz)*zooP(nz)
1088                  sumgrazzooN(nz)=sumgrazzooN(nz)+
1089         &             asseff(np,nz)*grazphy(np,nz)*zooP(nz)*R_NP_zoo(nz)
1090                  sumgrazlossN(nz)=sumgrazlossN(nz)+
1091         &             (R_NP(np)-asseff(np,nz)*R_NP_zoo(nz))*
1092         &             grazphy(np,nz)*zooP(nz)
1093                  sumgrazzooFe(nz)=sumgrazzooFe(nz)+
1094         &             asseff(np,nz)*grazphy(np,nz)*
1095         &                                      zooP(nz)*R_FeP_zoo(nz)
1096                  sumgrazlossFe(nz)=sumgrazlossFe(nz)+
1097         &             (R_FeP(np)-asseff(np,nz)*R_FeP_zoo(nz))*
1098         &             grazphy(np,nz)*zooP(nz)
1099                  sumgrazzooSi(nz)=sumgrazzooSi(nz)+
1100         &             asseff(np,nz)*grazphy(np,nz)*
1101         &                                     zooP(nz)*R_SiP_zoo(nz)
1102                  sumgrazlossSi(nz)=sumgrazlossSi(nz)+
1103         &             (R_SiP(np)-asseff(np,nz)*R_SiP_zoo(nz))*
1104         &             grazphy(np,nz)*zooP(nz)
1105    #ifdef ALLOW_CARBON
1106                  sumgrazzooC(nz)=sumgrazzooC(nz)+
1107         &             asseff(np,nz)*grazphy(np,nz)*zooP(nz)*R_PC_zoo(nz)
1108                  sumgrazlossC(nz)=sumgrazlossC(nz)+
1109         &             (R_PC(np)-asseff(np,nz)*R_PC_zoo(nz))*
1110         &             grazphy(np,nz)*zooP(nz)
1111                  sumgrazlossPIC(nz)=sumgrazlossPIC(nz)+
1112         &             (1. _d 0)*grazphy(np,nz)*
1113         &                           zooP(nz)*R_PC(np)*R_PICPOC(np)
1114    #endif
1115    #else
1116                sumgrazzoo(nz)=sumgrazzoo(nz)+                sumgrazzoo(nz)=sumgrazzoo(nz)+
1117       &             asseff(np,nz)*grazphy(np,nz)*zooP(nz)       &             asseff(np,nz)*grazphy(np,nz)*zooP(nz)
1118                sumgrazloss(nz)=sumgrazloss(nz)+                sumgrazloss(nz)=sumgrazloss(nz)+
# Line 1111  c sum over phy for impact on zoo, and al Line 1144  c sum over phy for impact on zoo, and al
1144       &             (1. _d 0)*grazphy(np,nz)*       &             (1. _d 0)*grazphy(np,nz)*
1145       &                           zooP(nz)*R_PC(np)*R_PICPOC(np)       &                           zooP(nz)*R_PC(np)*R_PICPOC(np)
1146  #endif  #endif
1147    #endif
1148             enddo             enddo
1149            enddo            enddo
1150            if (debug.eq.2) print*,'sumgrazzoo',sumgrazzoo            if (debug.eq.2) print*,'sumgrazzoo',sumgrazzoo
# Line 1797  c loss of O2 by remineralization Line 1831  c loss of O2 by remineralization
1831              enddo              enddo
1832  #endif  #endif
1833    
1834    #endif /* ALLOW_CARBON */
1835    
1836  #ifdef ALLOW_CDOM  #ifdef ALLOW_CDOM
1837              dcdomdt = totphy_cdomp + totzoo_cdomp +  preminP_cdom              dcdomdt = totphy_cdomp + totzoo_cdomp +  preminP_cdom
1838       &                  -cdomp_degrd       &                  -cdomp_degrd
1839  #endif  #endif
1840    
 #endif  
   
1841              if (debug.eq.10) print*,'dDOPdt', dDOPdt              if (debug.eq.10) print*,'dDOPdt', dDOPdt
1842              if (debug.eq.10) print*,'dpopdt',dpopdt              if (debug.eq.10) print*,'dpopdt',dpopdt
1843              if (debug.eq.10) print*,'dDONdt',dDONdt              if (debug.eq.10) print*,'dDONdt',dDONdt
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22