/[MITgcm]/MITgcm/verification/aim.5l_LatLon/code/calc_gt.F
ViewVC logotype

Contents of /MITgcm/verification/aim.5l_LatLon/code/calc_gt.F

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


Revision 1.1 - (show annotations) (download)
Thu Sep 27 15:36:29 2001 UTC (20 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint44h_post, checkpoint44b_post, checkpoint44e_post, release1_beta1, checkpoint44h_pre, release1_b1, checkpoint45a_post, release1-branch_tutorials, checkpoint43a-release1mods, release1-branch_branchpoint, chkpt44d_post, checkpoint42, release1_p12, release1_p13, release1_p10, release1_p11, release1_p16, release1_p17, release1_p14, release1_p15, release1_p8, release1_p9, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, checkpoint45b_post, release1_p7, checkpoint44g_post, release1_p13_pre, chkpt44a_pre, chkpt44a_post, release1-branch-end, chkpt44c_post, release1_final_v1, release1_chkpt44d_post, checkpoint44f_pre, checkpoint44b_pre, chkpt44c_pre, checkpoint45c_post, release1_p12_pre, checkpoint44e_pre, checkpoint43, release1_p6, checkpoint44, checkpoint45, checkpoint44f_post
Branch point for: release1, release1_coupled, release1-branch, release1_50yr, release1_final
update of aim-LatLon experiment following the recent changes in pkg/aim

1 C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_gt.F,v 1.36 2001/09/26 18:09:14 cnh Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: CALC_GT
8 C !INTERFACE:
9 SUBROUTINE CALC_GT(
10 I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
11 I xA,yA,uTrans,vTrans,rTrans,maskUp,
12 I KappaRT,
13 U fVerT,
14 I myTime,myIter,myThid )
15 C !DESCRIPTION: \bv
16 C *==========================================================*
17 C | SUBROUTINE CALC_GT
18 C | o Calculate the temperature tendency terms.
19 C *==========================================================*
20 C | A procedure called EXTERNAL_FORCING_T is called from
21 C | here. These procedures can be used to add per problem
22 C | heat 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 fVerT :: Flux of temperature (T) 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 fVerT (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 KappaRT(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 CEOP
79
80 #ifdef ALLOW_AUTODIFF_TAMC
81 C-- only the kUp part of fverT is set in this subroutine
82 C-- the kDown is still required
83 fVerT(1,1,kDown) = fVerT(1,1,kDown)
84 #endif
85
86 CALL GAD_CALC_RHS(
87 I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
88 I xA,yA,uTrans,vTrans,rTrans,maskUp,
89 I diffKhT, diffK4T, KappaRT, theta,
90 I GAD_TEMPERATURE, tempAdvScheme,
91 U fVerT, gT,
92 I myThid )
93
94 #ifdef INCLUDE_T_FORCING_CODE
95 C-- External thermal forcing term(s)
96 c CALL EXTERNAL_FORCING_T(
97 c I iMin,iMax,jMin,jMax,bi,bj,k,
98 c I myTime,myThid)
99 #endif /* INCLUDE_T_FORCING_CODE */
100
101 IF ( tempAdvScheme.EQ.ENUM_CENTERED_2ND
102 & .OR.tempAdvScheme.EQ.ENUM_UPWIND_3RD
103 & .OR.tempAdvScheme.EQ.ENUM_CENTERED_4TH ) THEN
104 CALL ADAMS_BASHFORTH2(
105 I bi, bj, K,
106 U gT, gTnm1,
107 I myIter, myThid )
108 ENDIF
109
110 #ifdef INCLUDE_T_FORCING_CODE
111 C-- External thermal forcing term(s)
112 CALL EXTERNAL_FORCING_T(
113 I iMin,iMax,jMin,jMax,bi,bj,k,
114 I myTime,myThid)
115 #endif /* INCLUDE_T_FORCING_CODE */
116
117 #ifdef NONLIN_FRSURF
118 IF (nonlinFreeSurf.GT.0) THEN
119 CALL FREESURF_RESCALE_G(
120 I bi, bj, K,
121 U gT,
122 I myThid )
123 ENDIF
124 #endif /* NONLIN_FRSURF */
125
126 RETURN
127 END

  ViewVC Help
Powered by ViewVC 1.1.22