/[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.3 by jahn, Thu Sep 22 16:01:30 2011 UTC revision 1.6 by jahn, Fri Jun 21 12:49:47 2013 UTC
# Line 49  C============== Local variables ======== Line 49  C============== Local variables ========
49         INTEGER np, nz, npnum         INTEGER np, nz, npnum
50    
51    
       _BEGIN_MASTER(myThid)  
52  c  c
53           tmptotP=0. _d 0           tmptotP=0. _d 0
54           tmptotN=0. _d 0           tmptotN=0. _d 0
# Line 82  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 95  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 106  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 130  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 !!!!!!!!!!!
# Line 163  c       print*,'tmptotA',tmptotA Line 174  c       print*,'tmptotA',tmptotA
174  c       print*,'tmptotO',tmptotO  c       print*,'tmptotO',tmptotO
175  #endif  #endif
176    
       _END_MASTER(myThid)  
   
177         IF ( myProcId.EQ.0 .AND. myThid.EQ.1 ) THEN         IF ( myProcId.EQ.0 .AND. myThid.EQ.1 ) THEN
178          write(DAR_cons_unit1,*), myIter, tmptotP, tmptotP/voltot          write(DAR_cons_unitP,*), myIter, tmptotP, tmptotP/voltot
179          write(DAR_cons_unit2,*), myIter, tmptotN, tmptotN/voltot          write(DAR_cons_unitN,*), myIter, tmptotN, tmptotN/voltot
180          write(DAR_cons_unit3,*), myIter, tmptotFe, tmptotFe/voltot          write(DAR_cons_unitF,*), myIter, tmptotFe, tmptotFe/voltot
181          write(DAR_cons_unit4,*), myIter, tmptotSi, tmptotSi/voltot          write(DAR_cons_unitS,*), myIter, tmptotSi, tmptotSi/voltot
182  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
183          write(DAR_cons_unit5,*), myIter, tmptotC, tmptotC/voltot          write(DAR_cons_unitC,*), myIter, tmptotC, tmptotC/voltot
184          write(DAR_cons_unit6,*), myIter, tmptotA, tmptotA/voltot          write(DAR_cons_unitA,*), myIter, tmptotA, tmptotA/voltot
185          write(DAR_cons_unit7,*), myIter, tmptotO, tmptotO/voltot          write(DAR_cons_unitO,*), myIter, tmptotO, tmptotO/voltot
186  #endif  #endif
187        ENDIF        ENDIF
188    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22