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

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

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

revision 1.4 by jahn, Wed Dec 7 20:04:00 2011 UTC revision 1.5 by stephd, Tue Jul 24 18:27:28 2012 UTC
# Line 81  c Phosphorus budget Line 81  c Phosphorus budget
81                npnum=iPhy+np-1                npnum=iPhy+np-1
82                tmptotP=tmptotP+Ptracer(i,j,k,bi,bj,npnum)*vol  !Phy                tmptotP=tmptotP+Ptracer(i,j,k,bi,bj,npnum)*vol  !Phy
83             enddo             enddo
84    #ifdef ALLOW_CDOM
85               tmptotP=tmptotP+Ptracer(i,j,k,bi,bj,iCDOM)*vol  !CDOM
86    #endif
87  c Nitrogen budget  c Nitrogen budget
88             tmptotN=tmptotN+Ptracer(i,j,k,bi,bj,iNO3)*vol    !NO3             tmptotN=tmptotN+Ptracer(i,j,k,bi,bj,iNO3)*vol    !NO3
89             tmptotN=tmptotN+Ptracer(i,j,k,bi,bj,iNH4)*vol   !NH4             tmptotN=tmptotN+Ptracer(i,j,k,bi,bj,iNH4)*vol   !NH4
# Line 94  c Nitrogen budget Line 97  c Nitrogen budget
97                npnum=iPhy+np-1                npnum=iPhy+np-1
98                tmptotN=tmptotN+Ptracer(i,j,k,bi,bj,npnum)*R_NP(np)*vol  !Phy                tmptotN=tmptotN+Ptracer(i,j,k,bi,bj,npnum)*R_NP(np)*vol  !Phy
99             enddo             enddo
100    #ifdef ALLOW_CDOM
101               tmptotN=tmptotN+Ptracer(i,j,k,bi,bj,iCDOM)*rnp_cdom*vol  !CDOM
102    #endif
103  c Iron budget  c Iron budget
104             tmptotFe=tmptotFe+Ptracer(i,j,k,bi,bj,iFeT)*vol    !Fet             tmptotFe=tmptotFe+Ptracer(i,j,k,bi,bj,iFeT)*vol    !Fet
105             tmptotFe=tmptotFe+Ptracer(i,j,k,bi,bj,iDOFe)*vol    !DOFe             tmptotFe=tmptotFe+Ptracer(i,j,k,bi,bj,iDOFe)*vol    !DOFe
# Line 105  c Iron budget Line 111  c Iron budget
111                npnum=iPhy+np-1                npnum=iPhy+np-1
112                tmptotFe=tmptotFe+Ptracer(i,j,k,bi,bj,npnum)*R_FeP(np)*vol  !Phy                tmptotFe=tmptotFe+Ptracer(i,j,k,bi,bj,npnum)*R_FeP(np)*vol  !Phy
113             enddo             enddo
114    #ifdef ALLOW_CDOM
115               tmptotFe=tmptotFe+Ptracer(i,j,k,bi,bj,iCDOM)*rfep_cdom*vol  !CDOM
116    #endif
117  c Silica budget  c Silica budget
118             tmptotSi=tmptotSi+Ptracer(i,j,k,bi,bj,iSi)*vol    !Si             tmptotSi=tmptotSi+Ptracer(i,j,k,bi,bj,iSi)*vol    !Si
119             tmptotSi=tmptotSi+Ptracer(i,j,k,bi,bj,iPOSi)*vol   !POSi             tmptotSi=tmptotSi+Ptracer(i,j,k,bi,bj,iPOSi)*vol   !POSi
# Line 129  c carbon budget Line 138  c carbon budget
138                npnum=iPhy+np-1                npnum=iPhy+np-1
139                tmptotC=tmptotC+Ptracer(i,j,k,bi,bj,npnum)*R_PC(np)*vol  !Phy                tmptotC=tmptotC+Ptracer(i,j,k,bi,bj,npnum)*R_PC(np)*vol  !Phy
140             enddo             enddo
141    #ifdef ALLOW_CDOM
142               tmptotC=tmptotC+Ptracer(i,j,k,bi,bj,iCDOM)*rcp_cdom*vol  !CDOM
143    #endif
144  c alkalinity budget !!!!! NOT FULL BUDGET!!!!!  c alkalinity budget !!!!! NOT FULL BUDGET!!!!!
145             tmptotA=tmptotA+Ptracer(i,j,k,bi,bj,iALK)*vol    !ALK             tmptotA=tmptotA+Ptracer(i,j,k,bi,bj,iALK)*vol    !ALK
146  c oxygem budget !!!!!!!! NOT FULL BUDGET !!!!!!!!!!!  c oxygem budget !!!!!!!! NOT FULL BUDGET !!!!!!!!!!!
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22