/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_dst3_adv_x.F
ViewVC logotype

Annotation of /MITgcm/pkg/generic_advdiff/gad_dst3_adv_x.F

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

Revision 1.15 - (hide annotations) (download)
Wed Jul 4 20:22:26 2012 UTC (11 years, 11 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint64, checkpoint65, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint65o, checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, HEAD
Changes since 1.14: +11 -11 lines 
comment out test on inAdMode (has no effect + OLD_DST3_FORMULATION no longer used)
1 jmc 1.15 C $Header: /u/gcmpack/MITgcm/pkg/generic_advdiff/gad_dst3_adv_x.F,v 1.14 2011/03/29 15:47:19 jmc Exp $
2 adcroft 1.2 C $Name: $
3 adcroft 1.1
4     #include "GAD_OPTIONS.h"
5    
6 jmc 1.10 CBOP
7     C !ROUTINE: GAD_DST3_ADV_X
8    
9     C !INTERFACE: ==========================================================
10 jmc 1.8 SUBROUTINE GAD_DST3_ADV_X(
11 jmc 1.12 I bi,bj,k, calcCFL, deltaTloc,
12 jmc 1.8 I uTrans, uFld,
13 jmc 1.4 I maskLocW, tracer,
14 adcroft 1.1 O uT,
15     I myThid )
16 jmc 1.10
17     C !DESCRIPTION:
18     C Calculates the area integrated zonal flux due to advection of a
19     C tracer using 3rd-order Direct Space and Time (DST-3) Advection Scheme
20    
21     C !USES: ===============================================================
22 adcroft 1.1 IMPLICIT NONE
23    
24     C == GLobal variables ==
25     #include "SIZE.h"
26 jmc 1.14 #ifdef OLD_DST3_FORMULATION
27 heimbach 1.5 #include "EEPARAMS.h"
28     #include "PARAMS.h"
29 jmc 1.14 #endif
30     #include "GRID.h"
31 adcroft 1.1 #include "GAD.h"
32    
33     C == Routine arguments ==
34 jmc 1.10 C !INPUT PARAMETERS: ===================================================
35     C bi,bj :: tile indices
36     C k :: vertical level
37 jmc 1.12 C calcCFL :: =T: calculate CFL number ; =F: take uFld as CFL.
38 jmc 1.10 C deltaTloc :: local time-step (s)
39     C uTrans :: zonal volume transport
40 jmc 1.12 C uFld :: zonal flow / CFL number
41 jmc 1.10 C tracer :: tracer field
42     C myThid :: thread number
43 adcroft 1.1 INTEGER bi,bj,k
44 jmc 1.12 LOGICAL calcCFL
45 heimbach 1.5 _RL deltaTloc
46 adcroft 1.1 _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
47 jmc 1.8 _RL uFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
48 jmc 1.4 _RS maskLocW(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
49 adcroft 1.1 _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
50 jmc 1.10 INTEGER myThid
51    
52     C !OUTPUT PARAMETERS: ==================================================
53     C uT :: zonal advective flux
54 adcroft 1.1 _RL uT (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
55    
56     C == Local variables ==
57 jmc 1.10 C !LOCAL VARIABLES: ====================================================
58     C i,j :: loop indices
59 jmc 1.11 C uCFL :: Courant-Friedrich-Levy number
60 adcroft 1.1 INTEGER i,j
61 jmc 1.11 _RL Rjm,Rj,Rjp,uCFL,d0,d1
62 jmc 1.10 #ifdef OLD_DST3_FORMULATION
63 adcroft 1.2 _RL psiP,psiM,thetaP,thetaM
64 heimbach 1.5 _RL smallNo
65 jmc 1.6 c _RL Rjjm,Rjjp
66 heimbach 1.5
67 jmc 1.15 c IF (inAdMode) THEN
68     c smallNo = 1.0D-20
69     c ELSE
70 heimbach 1.5 smallNo = 1.0D-20
71 jmc 1.15 c ENDIF
72 jmc 1.10 #endif
73 adcroft 1.1
74 jmc 1.15 DO j=1-OLy,sNy+OLy
75     uT(1-OLx,j)=0.
76     uT(2-OLx,j)=0.
77     uT(sNx+OLx,j)=0.
78 mlosch 1.13 ENDDO
79 jmc 1.15 DO j=1-OLy,sNy+OLy
80     DO i=1-OLx+2,sNx+OLx-1
81 jmc 1.4 Rjp=(tracer(i+1,j)-tracer( i ,j))*maskLocW(i+1,j)
82     Rj =(tracer( i ,j)-tracer(i-1,j))*maskLocW( i ,j)
83     Rjm=(tracer(i-1,j)-tracer(i-2,j))*maskLocW(i-1,j)
84 adcroft 1.1
85 jmc 1.12 uCFL = uFld(i,j)
86     IF ( calcCFL ) uCFL = ABS( uFld(i,j)*deltaTloc
87 jmc 1.11 & *recip_dxC(i,j,bi,bj)*recip_deepFacC(k) )
88     d0=(2.-uCFL)*(1.-uCFL)*oneSixth
89     d1=(1.-uCFL*uCFL)*oneSixth
90 jmc 1.10 #ifdef OLD_DST3_FORMULATION
91 heimbach 1.5 IF ( ABS(Rj).LT.smallNo .OR.
92     & ABS(Rjm).LT.smallNo ) THEN
93     thetaP=0.
94     psiP=0.
95     ELSE
96     thetaP=(Rjm+smallNo)/(smallNo+Rj)
97     psiP=d0+d1*thetaP
98     ENDIF
99     IF ( ABS(Rj).LT.smallNo .OR.
100     & ABS(Rjp).LT.smallNo ) THEN
101     thetaM=0.
102     psiM=0.
103     ELSE
104     thetaM=(Rjp+smallNo)/(smallNo+Rj)
105     psiM=d0+d1*thetaM
106     ENDIF
107 adcroft 1.1 uT(i,j)=
108 jmc 1.10 & 0.5*(uTrans(i,j)+ABS(uTrans(i,j)))
109 adcroft 1.2 & *( Tracer(i-1,j) + psiP*Rj )
110 jmc 1.10 & +0.5*(uTrans(i,j)-ABS(uTrans(i,j)))
111 adcroft 1.2 & *( Tracer( i ,j) - psiM*Rj )
112 jmc 1.10 #else /* OLD_DST3_FORMULATION */
113     uT(i,j)=
114     & 0.5*(uTrans(i,j)+ABS(uTrans(i,j)))
115     & *( Tracer(i-1,j) + (d0*Rj+d1*Rjm) )
116     & +0.5*(uTrans(i,j)-ABS(uTrans(i,j)))
117     & *( Tracer( i ,j) - (d0*Rj+d1*Rjp) )
118     #endif /* OLD_DST3_FORMULATION */
119 adcroft 1.1
120     ENDDO
121     ENDDO
122    
123     RETURN
124     END
  ViewVC Help
Powered by ViewVC 1.1.22