/[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.6 - (show annotations) (download)
Tue Nov 23 20:20:46 2004 UTC (19 years, 6 months ago) by mlosch
Branch: MAIN
Changes since 1.5: +3 -1 lines 
  - renamed PTRACERS_SEPARATE_FORCING to GCHEM_SEPARATE_FORCING;
    this flag is now only visible within GCHEM-pkg routines.
  - moved gchem_forcing_int temporarily into ptracers_forcing, this will
    be replaced by a 3D array that is computed before thermodynamics; for
    this, the tracer ID iTracer had to be included into the parameter list
    of ptracers_forcing
  - cleaning up, in particular CPP-flags and call to ptracer_statvars
1 C $Header: /u/gcmpack/MITgcm/pkg/gchem/gchem_forcing_int.F,v 1.5 2004/11/19 21:28:47 mlosch Exp $
2 C $Name: $
3
4 #include "GCHEM_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: GCHEM_FORCING_INT
8
9 C !INTERFACE: ==========================================================
10 SUBROUTINE GCHEM_FORCING_INT(bi,bj,iMin,iMax,jMin,jMax,k,
11 & iTracer, myTime,myIter, myThid )
12
13 C !DESCRIPTION:
14 C calls subroutine that will update tracer tendencies
15 C this will be incorporated into regular timestepping in
16 C ptracers_intergrate.F
17
18 C !USES: ===============================================================
19 IMPLICIT NONE
20 #include "SIZE.h"
21 #include "GRID.h"
22 #include "DYNVARS.h"
23 #include "EEPARAMS.h"
24 #include "PARAMS.h"
25 #include "PTRACERS_SIZE.h"
26 #include "PTRACERS.h"
27
28 C !INPUT PARAMETERS: ===================================================
29 C myThid :: thread number
30 C myIter :: current timestep
31 C myTime :: current time
32 C iTracer :: ptracer number
33 C bi,bj :: tile indices
34 C k :: vertical level
35 INTEGER myThid, myIter
36 _RL myTime
37 INTEGER iTracer
38 INTEGER bi,bj,imin,imax,jmin,jmax,k
39
40 C !LOCAL VARIABLES: ====================================================
41 C i,j :: loop indices
42 INTEGER i,j
43 INTEGER niter
44 CEOP
45
46 #ifdef ALLOW_GCHEM
47 # ifndef GCHEM_SEPARATE_FORCING
48 c
49
50 ccccccccccccccccccccccccc
51 c chemical forcing c
52 ccccccccccccccccccccccccc
53 ccccccccccccccccccccccccccc CFC cccccccccccccccccccccccccccccccc
54 c
55 #ifdef ALLOW_CFC
56 if (iTracer.eq.1.and.k.eq.1)
57 & call CFC11_Forcing( Ptracer(1-Olx,1-Oly,1,bi,bj,1),
58 & gPtr(1-Olx,1-Oly,1,bi,bj,1),
59 & bi,bj,imin,imax,jmin,jmax,
60 & myIter,myTime,myThid)
61 if (iTracer.eq.2.and.k.eq.1)
62 & call CFC12_Forcing( Ptracer(1-Olx,1-Oly,1,bi,bj,2),
63 & gPtr(1-Olx,1-Oly,1,bi,bj,2),
64 & bi,bj,imin,imax,jmin,jmax,
65 & myIter,myTime,myThid)
66 #endif /* ALLOW_CFC */
67 cccccccccccccccccccccccccc END CFC cccccccccccccccccccccccccccccccccc
68 c
69 # endif /* GCHEM_SEPARATE_FORCING */
70 #endif /* ALLOW_GCHEM */
71
72 RETURN
73 END
  ViewVC Help
Powered by ViewVC 1.1.22