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

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

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


Revision 1.3 - (show annotations) (download)
Tue Jul 13 18:04:42 2004 UTC (19 years, 11 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint54f_post, checkpoint55j_post, checkpoint56, checkpoint55g_post, checkpoint55c_post, checkpoint55f_post, checkpoint55, checkpoint55h_post, checkpoint55e_post, checkpoint54c_post, checkpoint54d_post, checkpoint54e_post, checkpoint55b_post, checkpoint55i_post, checkpoint55d_pre, checkpoint55a_post, checkpoint55d_post
Changes since 1.2: +4 -0 lines
max number of passive tracers is now defined in PTRACERS_SIZE.h

1 C $Header: $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5 #include "PTRACERS_OPTIONS.h"
6 #include "GCHEM_OPTIONS.h"
7
8 CBOP
9 C !ROUTINE: GCHEM_FORCING_INT
10
11 C !INTERFACE: ==========================================================
12 SUBROUTINE GCHEM_FORCING_INT(bi,bj,iMin,iMax,jMin,jMax,k,
13 & iTracer, myTime,myIter, myThid )
14
15 C !DESCRIPTION:
16 C calls subroutine that will update tracer tendencies
17 C this will be incorporated into regular timestepping in
18 C ptracers_intergrate.F
19
20 C !USES: ===============================================================
21 IMPLICIT NONE
22 #include "SIZE.h"
23 #include "GRID.h"
24 #include "DYNVARS.h"
25 #include "EEPARAMS.h"
26 #include "PARAMS.h"
27 #ifdef ALLOW_PTRACERS
28 #include "PTRACERS_SIZE.h"
29 #include "PTRACERS.h"
30 #endif
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_PTRACERS
51 #ifdef ALLOW_GCHEM
52 c
53
54 ccccccccccccccccccccccccc
55 c chemical forcing c
56 ccccccccccccccccccccccccc
57 ccccccccccccccccccccccccccc DIC cccccccccccccccccccccccccccccccc
58 c
59 #ifdef ALLOW_CFC
60 if (iTracer.eq.1.and.k.eq.1)
61 & call CFC11_Forcing( Ptracer(1-Olx,1-Oly,1,bi,bj,1),
62 & gPtr(1-Olx,1-Oly,1,bi,bj,1),
63 & bi,bj,imin,imax,jmin,jmax,
64 & myIter,myTime,myThid)
65 if (iTracer.eq.2.and.k.eq.1)
66 & call CFC12_Forcing( Ptracer(1-Olx,1-Oly,1,bi,bj,2),
67 & gPtr(1-Olx,1-Oly,1,bi,bj,2),
68 & bi,bj,imin,imax,jmin,jmax,
69 & myIter,myTime,myThid)
70 #endif
71 cccccccccccccccccccccccccc END DIC cccccccccccccccccccccccccccccccccc
72 c
73 #endif
74 #endif /* ALLOW_PTRACERS */
75
76 RETURN
77 END

  ViewVC Help
Powered by ViewVC 1.1.22