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

Contents of /MITgcm/pkg/gchem/gchem_add_tendency.F

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

Revision 1.2 - (show annotations) (download)
Mon Nov 29 03:34:55 2004 UTC (19 years, 7 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint58e_post, checkpoint57v_post, checkpoint58u_post, checkpoint58w_post, checkpoint57m_post, checkpoint57s_post, checkpoint57k_post, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint58r_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint58n_post, checkpoint58x_post, checkpoint57g_pre, checkpoint58t_post, checkpoint58h_post, checkpoint56c_post, checkpoint57y_pre, checkpoint57f_pre, checkpoint57a_post, checkpoint58q_post, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59i, checkpoint59h, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint57a_pre, checkpoint57, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint57x_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, checkpoint58o_post, checkpoint57z_post, checkpoint57c_post, checkpoint58y_post, checkpoint58k_post, checkpoint58v_post, checkpoint58s_post, checkpoint58p_post, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint57h_post
Changes since 1.1: +1 -22 lines 
o remove obsolete files
o remove unnecessary comments
1 C $Header: /u/gcmpack/MITgcm/pkg/gchem/gchem_add_tendency.F,v 1.1 2004/11/28 23:48:31 mlosch Exp $
2 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 C In the case of GCHEM_SEPARATE_FORCING undefined, GCHEM_ADD_TENDENCY
15 C will update passive tracer tendencies gPtr with gchemTendency,
16 C the latter having been computed in GCHEM_CALC_TENDENDY, so that
17 C they will be incorporated into regular timestepping in
18 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 #include "PTRACERS.h"
30 #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