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