/[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.2 by stephd, Tue Oct 7 16:21:12 2003 UTC revision 1.5 by jmc, Tue Jul 13 18:04:42 2004 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5  #include "PTRACERS_OPTIONS.h"  #include "PTRACERS_OPTIONS.h"
6  #include "GCHEM_OPTIONS.h"  #include "GCHEM_OPTIONS.h"
7    
8    CBOP
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE GCHEM_FORCING_SEP(myTime,myIter, myThid )        SUBROUTINE GCHEM_FORCING_SEP(myTime,myIter, myThid )
11    
12    C !DESCRIPTION:
13    C calls subroutine that will update tracers values.
14    C this will be in a secondary timestep after that used for
15    C advection and diffusion in ptracers_intergrate.F
16    
17    C !USES: ===============================================================
18        IMPLICIT NONE        IMPLICIT NONE
19  #include "SIZE.h"  #include "SIZE.h"
20  #include "GRID.h"  #include "GRID.h"
# Line 12  C !INTERFACE: ========================== Line 22  C !INTERFACE: ==========================
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
23  #include "PARAMS.h"  #include "PARAMS.h"
24  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
25    #include "PTRACERS_SIZE.h"
26  #include "PTRACERS.h"  #include "PTRACERS.h"
27  #endif  #endif
28  #ifdef DIC_BIOTIC  #ifdef DIC_BIOTIC
# Line 23  C !INPUT PARAMETERS: =================== Line 34  C !INPUT PARAMETERS: ===================
34  C  myThid               :: thread number  C  myThid               :: thread number
35        INTEGER myThid, myIter        INTEGER myThid, myIter
36        _RL myTime        _RL myTime
37    
38     !LOCAL VARIABLES: ====================================================
39    C  i,j                  :: loop indices
40    C  bi,bj                :: tile indices
41    C  k                    :: vertical level
42        INTEGER bi,bj,imin,imax,jmin,jmax        INTEGER bi,bj,imin,imax,jmin,jmax
43        INTEGER i,j        INTEGER i,j
44        INTEGER niter        INTEGER niter
45    CEOP
46    
47  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
48  #ifdef ALLOW_GCHEM  #ifdef ALLOW_GCHEM
49    
50           _BARRIER         _BARRIER
51           _BEGIN_MASTER( mythid )         _BEGIN_MASTER( mythid )
52    
53  ccccccccccccccccccccccccc  ccccccccccccccccccccccccc
54  c global calculations   c  c global calculations   c
55  ccccccccccccccccccccccccc  ccccccccccccccccccccccccc
56  #ifdef DIC_BIOTIC  #ifdef DIC_BIOTIC
57  c find global surface averages  c find global surface averages
58           gsm_s = 0.
59           gsm_dic = 0.
60           gsm_alk = 0.
61         call tracer_meanarea(myThid,salt, 1,gsm_s)         call tracer_meanarea(myThid,salt, 1,gsm_s)
62         call tracer_meanarea(myThid,         call tracer_meanarea(myThid,
63       &               ptracer(1-Olx,1-Oly,1,1,1,1), 1, gsm_dic)       &               ptracer(1-Olx,1-Oly,1,1,1,1), 1, gsm_dic)
64         print*,'mean surface dic', gsm_dic,gsm_s         print*,'mean surface dic', gsm_dic,gsm_s
65         call tracer_meanarea(myThid,         call tracer_meanarea(myThid,
66       &               ptracer(1-Olx,1-Oly,1,1,1,2), 1, gsm_alk)       &               ptracer(1-Olx,1-Oly,1,1,1,2), 1, gsm_alk)
   
67  #endif  #endif
68  ccccccccccccccccccccccccccccccccccccccccccc  ccccccccccccccccccccccccccccccccccccccccccc
69    
# Line 55  cccccccccccccccccccccccccccccccccccccccc Line 74  cccccccccccccccccccccccccccccccccccccccc
74  ccccccccccccccccccccccccc  ccccccccccccccccccccccccc
75  c chemical forcing      c  c chemical forcing      c
76  ccccccccccccccccccccccccc  ccccccccccccccccccccccccc
77  C$TAF loop = parallel  C$taf loop = parallel
78         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
79  C$TAF loop = parallel  C$taf loop = parallel
80          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
81    
82          jMin=1-OLy          jMin=1-OLy
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22