/[MITgcm]/MITgcm/verification/aim.5l_cs/code/calc_gs.F
ViewVC logotype

Contents of /MITgcm/verification/aim.5l_cs/code/calc_gs.F

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


Revision 1.4 - (show annotations) (download)
Sat Jun 15 03:50:52 2002 UTC (21 years, 9 months ago) by jmc
Branch: MAIN
CVS Tags: HEAD
Changes since 1.3: +1 -1 lines
FILE REMOVED
T,S forcing outside Adams-Bashforh is now a standard option

1 C $Header: /u/gcmpack/MITgcm/verification/aim.5l_cs/code/calc_gs.F,v 1.3 2001/09/27 15:49:45 jmc Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: CALC_GS
8 C !INTERFACE:
9 SUBROUTINE CALC_GS(
10 I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
11 I xA,yA,uTrans,vTrans,rTrans,maskUp,
12 I KappaRS,
13 U fVerS,
14 I myTime,myIter,myThid )
15 C !DESCRIPTION: \bv
16 C *==========================================================*
17 C | SUBROUTINE CALC_GS
18 C | o Calculate the salt tendency terms.
19 C *==========================================================*
20 C | A procedure called EXTERNAL_FORCING_S is called from
21 C | here. These procedures can be used to add per problem
22 C | E-P flux source terms.
23 C | Note: Although it is slightly counter-intuitive the
24 C | EXTERNAL_FORCING routine is not the place to put
25 C | file I/O. Instead files that are required to
26 C | calculate the external source terms are generally
27 C | read during the model main loop. This makes the
28 C | logisitics of multi-processing simpler and also
29 C | makes the adjoint generation simpler. It also
30 C | allows for I/O to overlap computation where that
31 C | is supported by hardware.
32 C | Aside from the problem specific term the code here
33 C | forms the tendency terms due to advection and mixing
34 C | The baseline implementation here uses a centered
35 C | difference form for the advection term and a tensorial
36 C | divergence of a flux form for the diffusive term. The
37 C | diffusive term is formulated so that isopycnal mixing and
38 C | GM-style subgrid-scale terms can be incorporated b simply
39 C | setting the diffusion tensor terms appropriately.
40 C *==========================================================*
41 C \ev
42
43 C !USES:
44 IMPLICIT NONE
45 C == GLobal variables ==
46 #include "SIZE.h"
47 #include "DYNVARS.h"
48 #include "EEPARAMS.h"
49 #include "PARAMS.h"
50 #include "GAD.h"
51
52 C !INPUT/OUTPUT PARAMETERS:
53 C == Routine arguments ==
54 C fVerS :: Flux of salt (S) in the vertical
55 C direction at the upper(U) and lower(D) faces of a cell.
56 C maskUp :: Land mask used to denote base of the domain.
57 C xA :: Tracer cell face area normal to X
58 C yA :: Tracer cell face area normal to X
59 C uTrans :: Zonal volume transport through cell face
60 C vTrans :: Meridional volume transport through cell face
61 C rTrans :: Vertical volume transport through cell face
62 C bi, bj, iMin, iMax, jMin, jMax :: Range of points for which calculation
63 C results will be set.
64 C myThid :: Instance number for this innvocation of CALC_GT
65 _RL fVerS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
66 _RS xA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
67 _RS yA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
68 _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
69 _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
70 _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
71 _RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
72 _RL KappaRS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
73 INTEGER k,kUp,kDown,kM1
74 INTEGER bi,bj,iMin,iMax,jMin,jMax
75 _RL myTime
76 INTEGER myIter
77 INTEGER myThid
78
79 CEOP
80
81 #ifdef ALLOW_AUTODIFF_TAMC
82 C-- only the kUp part of fverS is set in this subroutine
83 C-- the kDown is still required
84 fVerS(1,1,kDown) = fVerS(1,1,kDown)
85 #endif
86
87 CALL GAD_CALC_RHS(
88 I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
89 I xA,yA,uTrans,vTrans,rTrans,maskUp,
90 I diffKhS, diffK4S, KappaRS, Salt,
91 I GAD_SALINITY, saltAdvScheme,
92 U fVerS, gS,
93 I myThid )
94
95 C-- External forcing term(s)
96 c CALL EXTERNAL_FORCING_S(
97 c I iMin,iMax,jMin,jMax,bi,bj,k,
98 c I myTime,myThid)
99
100 IF ( saltAdvScheme.EQ.ENUM_CENTERED_2ND
101 & .OR.saltAdvScheme.EQ.ENUM_UPWIND_3RD
102 & .OR.saltAdvScheme.EQ.ENUM_CENTERED_4TH ) THEN
103 CALL ADAMS_BASHFORTH2(
104 I bi, bj, K,
105 U gS, gSnm1,
106 I myIter, myThid )
107 ENDIF
108
109 C-- External forcing term(s)
110 CALL EXTERNAL_FORCING_S(
111 I iMin,iMax,jMin,jMax,bi,bj,k,
112 I myTime,myThid)
113
114 #ifdef NONLIN_FRSURF
115 IF (nonlinFreeSurf.GT.0) THEN
116 CALL FREESURF_RESCALE_G(
117 I bi, bj, K,
118 U gS,
119 I myThid )
120 ENDIF
121 #endif /* NONLIN_FRSURF */
122
123 RETURN
124 END

  ViewVC Help
Powered by ViewVC 1.1.22