/[MITgcm]/MITgcm_contrib/darwin/pkg/gchem/gchem_add_tendency.F
ViewVC logotype

Annotation of /MITgcm_contrib/darwin/pkg/gchem/gchem_add_tendency.F

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


Revision 1.3 - (hide annotations) (download)
Mon Jun 2 16:18:06 2008 UTC (17 years, 2 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin_ckpt61i_20090224, ctrb_darwin_ckpt61r_20090626, ctrb_darwin_ckpt61m_20090430, ctrb_darwin_ckpt61d_20081212, ctrb_darwin_ckpt61t_20090820, ctrb_darwin_ckpt59p_20080604, ctrb_darwin_ckpt61v_20091024, ctrb_darwin_ckpt61w_20091024, ctrb_darwin_ckpt61h_20090224, ctrb_darwin_ckpt62a_20100117, ctrb_darwin_ckpt61s_20090630, ctrb_darwin_ckpt61_20080822, ctrb_darwin_ckpt61y_20091120, ctrb_darwin_ckpt61q_20090626, ctrb_darwin_ckpt61_20080624, ctrb_darwin_ckpt61o_20090527, ctrb_darwin_ckpt61d_20081013, ctrb_darwin_ckpt61k_20090312, ctrb_darwin_ckpt61i_20090312, ctrb_darwin_ckpt62_20091227, ctrb_darwin_ckpt59q_20080605, ctrb_darwin_ckpt61o_20090610, ctrb_darwin_ckpt61x_20091024, ctrb_darwin_ckpt59o_20080602, ctrb_darwin_ckpt61p_20090610, ctrb_darwin_ckpt61b_20080822, ctrb_darwin_ckpt61u_20090825, ctrb_darwin_ckpt61l_20090408, ctrb_darwin_ckpt62b_20100201, ctrb_darwin_ckpt61z_20091207, ctrb_darwin_ckpt62c_20100303, ctrb_darwin_ckpt59r_20080606, ctrb_darwin_ckpt61n_20090519, ctrb_darwin_ckpt60_20080619, ctrb_darwin_ckpt61n_20090513
Changes since 1.2: +1 -2 lines
cleanup unused includes

1 jahn 1.3 C $Header: /u/gcmpack/MITgcm_contrib/darwin/pkg/gchem/gchem_add_tendency.F,v 1.2 2008/02/26 17:13:11 jahn Exp $
2 jahn 1.1 C $Name: $
3    
4     #include "GCHEM_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: GCHEM_ADD_TENDENCY
8    
9     C !INTERFACE: ==========================================================
10     SUBROUTINE GCHEM_ADD_TENDENCY(bi,bj,iMin,iMax,jMin,jMax,k,
11     & iTracer, myTime,myIter, myThid )
12    
13     C !DESCRIPTION:
14 jahn 1.2 C In the case of GCHEM_SEPARATE_FORCING undefined, GCHEM_ADD_TENDENCY
15 jahn 1.1 C will update passive tracer tendencies gPtr with gchemTendency,
16     C the latter having been computed in GCHEM_CALC_TENDENDY, so that
17 jahn 1.2 C they will be incorporated into regular timestepping in
18 jahn 1.1 C PTRACERS_INTERGRATE.
19     C IF GCHEM_SEPARATE_FORCING is defined, this routine is empty.
20    
21     C !USES: ===============================================================
22     IMPLICIT NONE
23     #include "SIZE.h"
24     #include "GRID.h"
25     #include "DYNVARS.h"
26     #include "EEPARAMS.h"
27     #include "PARAMS.h"
28     #include "PTRACERS_SIZE.h"
29 jahn 1.2 #include "PTRACERS_FIELDS.h"
30 jahn 1.1 #include "GCHEM_FIELDS.h"
31    
32     C !INPUT PARAMETERS: ===================================================
33     C myThid :: thread number
34     C myIter :: current timestep
35     C myTime :: current time
36     C iTracer :: ptracer number
37     C bi,bj :: tile indices
38     C k :: vertical level
39     INTEGER myThid, myIter
40     _RL myTime
41     INTEGER iTracer
42     INTEGER bi,bj,imin,imax,jmin,jmax,k
43    
44     C !LOCAL VARIABLES: ====================================================
45     C i,j :: loop indices
46     INTEGER i,j
47     INTEGER niter
48     CEOP
49    
50     #ifdef ALLOW_GCHEM
51     # ifndef GCHEM_SEPARATE_FORCING
52     C apply the tendencies computed in GCHEM_CALC_TENDENCY to the passive
53     C tracers tendendy terms. The part of the code could also go into
54     C ptracers_forcing and replace the call to gchem_add_tendency there,
55     C but this way, no gchem-related header files and CPP-flags are required
56     C within ptracers_forcing.
57     DO j=jMin,jMax
58     DO i=iMin,iMax
59     gPtr(I,J,K,bi,bj,iTracer) = gPtr(I,J,K,bi,bj,iTracer)
60     & + gchemTendency(I,J,K,bi,bj,iTracer)
61     ENDDO
62     ENDDO
63     # endif /* GCHEM_SEPARATE_FORCING */
64     #endif /* ALLOW_GCHEM */
65    
66     RETURN
67     END

  ViewVC Help
Powered by ViewVC 1.1.22