/[MITgcm]/MITgcm/pkg/gchem/gchem_forcing_sep.F
ViewVC logotype

Diff of /MITgcm/pkg/gchem/gchem_forcing_sep.F

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

revision 1.22 by dfer, Fri Apr 4 21:43:18 2008 UTC revision 1.26 by jahn, Tue Jun 30 16:43:00 2009 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "GCHEM_OPTIONS.h"  #include "GCHEM_OPTIONS.h"
5    #ifdef ALLOW_DIC
6    #include "DIC_OPTIONS.h"
7    #endif
8    #ifdef ALLOW_DARWIN
9    #include "DARWIN_OPTIONS.h"
10    #endif
11    
12    
13  CBOP  CBOP
14  C !ROUTINE: GCHEM_FORCING_SEP  C !ROUTINE: GCHEM_FORCING_SEP
# Line 23  C !USES: =============================== Line 30  C !USES: ===============================
30  #include "EEPARAMS.h"  #include "EEPARAMS.h"
31  #include "PARAMS.h"  #include "PARAMS.h"
32  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
 c#include "PTRACERS_PARAMS.h"  
33  #include "PTRACERS_FIELDS.h"  #include "PTRACERS_FIELDS.h"
34    #include "GCHEM.h"
35  #ifdef ALLOW_DIC  #ifdef ALLOW_DIC
36  #include "DIC_VARS.h"  #include "DIC_VARS.h"
37  #endif /* ALLOW_DIC */  #endif /* ALLOW_DIC */
# Line 55  c find global surface averages Line 62  c find global surface averages
62         gsm_s = 0. _d 0         gsm_s = 0. _d 0
63         gsm_dic = 0. _d 0         gsm_dic = 0. _d 0
64         gsm_alk = 0. _d 0         gsm_alk = 0. _d 0
65         call tracer_meanarea(myThid,salt, 1,gsm_s)         call tracer_meanarea(salt, 1,gsm_s,myThid)
66         call tracer_meanarea(myThid,         call tracer_meanarea(
67       &               ptracer(1-Olx,1-Oly,1,1,1,1), 1, gsm_dic)       &             ptracer(1-Olx,1-Oly,1,1,1,1), 1, gsm_dic, myThid )
68         print*,'mean surface dic', gsm_dic,gsm_s         print*,'mean surface dic', gsm_dic,gsm_s
69         call tracer_meanarea(myThid,         call tracer_meanarea(
70       &               ptracer(1-Olx,1-Oly,1,1,1,2), 1, gsm_alk)       &             ptracer(1-Olx,1-Oly,1,1,1,2), 1, gsm_alk, myThid )
71  #endif /* ALLOW_OLD_VIRTUALFLUX */  #endif /* ALLOW_OLD_VIRTUALFLUX */
72  #endif  #endif
73  ccccccccccccccccccccccccccccccccccccccccccc  ccccccccccccccccccccccccccccccccccccccccccc
# Line 80  C$taf loop = parallel Line 87  C$taf loop = parallel
87          iMax=sNx          iMax=sNx
88  c  c
89  ccccccccccccccccccccccccccc DIC cccccccccccccccccccccccccccccccc  ccccccccccccccccccccccccccc DIC cccccccccccccccccccccccccccccccc
90  c  
91  #ifdef ALLOW_DIC  #ifdef ALLOW_DIC
92  #ifdef ALLOW_FE  #ifdef ALLOW_FE
93            call DIC_BIOTIC_Forcing( Ptracer(1-Olx,1-Oly,1,bi,bj,1),            call DIC_BIOTIC_Forcing( Ptracer(1-Olx,1-Oly,1,bi,bj,1),
# Line 112  c Line 119  c
119    
120  #endif  #endif
121  cccccccccccccccccccccccccc END DIC cccccccccccccccccccccccccccccccccc  cccccccccccccccccccccccccc END DIC cccccccccccccccccccccccccccccccccc
122    #ifdef ALLOW_DARWIN
123            IF ( useDARWIN ) THEN
124    #ifdef NUT_SUPPLY
125    c articficial supply of nutrients
126              call DARWIN_Nut_Supply( Ptracer(1-Olx,1-Oly,1,bi,bj,1),
127         &                          bi,bj,imin,imax,jmin,jmax,
128         &                          myIter,myTime,myThid)
129              call DARWIN_Nut_Supply( Ptracer(1-Olx,1-Oly,1,bi,bj,2),
130         &                          bi,bj,imin,imax,jmin,jmax,
131         &                          myIter,myTime,myThid)
132              call DARWIN_Nut_Supply( Ptracer(1-Olx,1-Oly,1,bi,bj,3),
133         &                          bi,bj,imin,imax,jmin,jmax,
134         &                          myIter,myTime,myThid)
135              call DARWIN_Nut_Supply( Ptracer(1-Olx,1-Oly,1,bi,bj,4),
136         &                          bi,bj,imin,imax,jmin,jmax,
137         &                          myIter,myTime,myThid)
138    #endif
139    ccccccccccccccc
140              call DARWIN_cons(myIter,myTime,myThid)
141    C darwin_forcing operates on bi,bj part only, but needs to get full
142    C array because of last (iPtr) index
143              call DARWIN_Forcing(  Ptracer(1-Olx,1-Oly,1,1,1,1),
144         &                          bi,bj,imin,imax,jmin,jmax,
145         &                          myIter,myTime,myThid)
146              call DARWIN_cons(myIter,myTime,myThid)
147    #ifdef ALLOW_CARBON
148           CALL DIC_ATMOS( 1, myTime, myIter, myThid )
149    #endif
150            ENDIF
151    #endif
152    
 c  
153          ENDDO          ENDDO
154         ENDDO         ENDDO
155    
156  #ifdef ALLOW_DIC  #ifdef ALLOW_DIC
157         CALL DIC_ATMOS( 1, myTime, myIter, myThid )         CALL DIC_ATMOS( 1, myTime, myIter, myThid )
158           CALL DIC_STORE_FLUXCO2( myTime, myIter, myThid )
159  #endif  #endif
160    
161  #endif /* GCHEM_SEPARATE_FORCING */  #endif /* GCHEM_SEPARATE_FORCING */
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22