/[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.11 by stephd, Thu Aug 18 19:29:51 2005 UTC revision 1.21 by dfer, Fri Jan 18 01:06:29 2008 UTC
# Line 11  C !INTERFACE: ========================== Line 11  C !INTERFACE: ==========================
11  C !DESCRIPTION:  C !DESCRIPTION:
12  C     calls subroutine that will update passive tracers values  C     calls subroutine that will update passive tracers values
13  C     with a separate timestep. Since GCHEM_FORCING_SEP is now  C     with a separate timestep. Since GCHEM_FORCING_SEP is now
14  C     called before DO_FIELDS_BLOCKING_EXCHANGES, the passive  C     called before DO_FIELDS_BLOCKING_EXCHANGES, the passive
15  C     tracer values in the halo regions are not up to date and  C     tracer values in the halo regions are not up to date and
16  C     must not be used.  C     must not be used.
17    
18  C !USES: ===============================================================  C !USES: ===============================================================
# Line 23  C !USES: =============================== Line 23  C !USES: ===============================
23  #include "EEPARAMS.h"  #include "EEPARAMS.h"
24  #include "PARAMS.h"  #include "PARAMS.h"
25  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
26  #include "PTRACERS.h"  c#include "PTRACERS_PARAMS.h"
27    #include "PTRACERS_FIELDS.h"
28  #ifdef ALLOW_DIC  #ifdef ALLOW_DIC
29  #include "DIC_BIOTIC.h"  #include "DIC_BIOTIC.h"
30  #include "DIC_ABIOTIC.h"  #include "DIC_ABIOTIC.h"
# Line 46  CEOP Line 47  CEOP
47  #ifdef ALLOW_GCHEM  #ifdef ALLOW_GCHEM
48  #ifdef GCHEM_SEPARATE_FORCING  #ifdef GCHEM_SEPARATE_FORCING
49    
        _BARRIER  
        _BEGIN_MASTER( mythid )  
   
50  ccccccccccccccccccccccccc  ccccccccccccccccccccccccc
51  c global calculations   c  c global calculations   c
52  ccccccccccccccccccccccccc  ccccccccccccccccccccccccc
53  #ifdef ALLOW_DIC  #ifdef ALLOW_DIC
54    #ifdef ALLOW_OLD_VIRTUALFLUX
55  c find global surface averages  c find global surface averages
56         gsm_s = 0.         gsm_s = 0. _d 0
57         gsm_dic = 0.         gsm_dic = 0. _d 0
58         gsm_alk = 0.         gsm_alk = 0. _d 0
59         call tracer_meanarea(myThid,salt, 1,gsm_s)         call tracer_meanarea(myThid,salt, 1,gsm_s)
60         call tracer_meanarea(myThid,         call tracer_meanarea(myThid,
61       &               ptracer(1-Olx,1-Oly,1,1,1,1), 1, gsm_dic)       &               ptracer(1-Olx,1-Oly,1,1,1,1), 1, gsm_dic)
62         print*,'mean surface dic', gsm_dic,gsm_s         print*,'mean surface dic', gsm_dic,gsm_s
63         call tracer_meanarea(myThid,         call tracer_meanarea(myThid,
64       &               ptracer(1-Olx,1-Oly,1,1,1,2), 1, gsm_alk)       &               ptracer(1-Olx,1-Oly,1,1,1,2), 1, gsm_alk)
65    #endif /* ALLOW_OLD_VIRTUALFLUX */
66  #endif  #endif
67  ccccccccccccccccccccccccccccccccccccccccccc  ccccccccccccccccccccccccccccccccccccccccccc
68    
       _END_MASTER( mythid )  
       _BARRIER  
   
69    
70  ccccccccccccccccccccccccc  ccccccccccccccccccccccccc
71  c chemical forcing      c  c chemical forcing      c
# Line 96  c Line 93  c
93       &                          bi,bj,imin,imax,jmin,jmax,       &                          bi,bj,imin,imax,jmin,jmax,
94       &                          myIter,myTime,myThid)       &                          myIter,myTime,myThid)
95  #else  #else
96    #ifdef ALLOW_O2
97            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),
98       &                          Ptracer(1-Olx,1-Oly,1,bi,bj,2),       &                          Ptracer(1-Olx,1-Oly,1,bi,bj,2),
99       &                          Ptracer(1-Olx,1-Oly,1,bi,bj,3),       &                          Ptracer(1-Olx,1-Oly,1,bi,bj,3),
# Line 103  c Line 101  c
101       &                          Ptracer(1-Olx,1-Oly,1,bi,bj,5),       &                          Ptracer(1-Olx,1-Oly,1,bi,bj,5),
102       &                          bi,bj,imin,imax,jmin,jmax,       &                          bi,bj,imin,imax,jmin,jmax,
103       &                          myIter,myTime,myThid)       &                          myIter,myTime,myThid)
104    #else
105              call DIC_BIOTIC_Forcing( Ptracer(1-Olx,1-Oly,1,bi,bj,1),
106         &                          Ptracer(1-Olx,1-Oly,1,bi,bj,2),
107         &                          Ptracer(1-Olx,1-Oly,1,bi,bj,3),
108         &                          Ptracer(1-Olx,1-Oly,1,bi,bj,4),
109         &                          bi,bj,imin,imax,jmin,jmax,
110         &                          myIter,myTime,myThid)
111  #endif  #endif
112  #endif  #endif
113    
114    #endif
115  cccccccccccccccccccccccccc END DIC cccccccccccccccccccccccccccccccccc  cccccccccccccccccccccccccc END DIC cccccccccccccccccccccccccccccccccc
116    
117  c  c
118          ENDDO          ENDDO
119         ENDDO         ENDDO
120    
121    #ifdef ALLOW_DIC
122           CALL DIC_ATMOS( 1, myTime, myIter, myThid )
123    #endif
124    
125  #endif /* GCHEM_SEPARATE_FORCING */  #endif /* GCHEM_SEPARATE_FORCING */
126  #endif /* ALLOW_GCHEM */  #endif /* ALLOW_GCHEM */
127    
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22