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

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

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

Revision 1.14 - (hide annotations) (download)
Wed Jul 4 20:22:26 2012 UTC (11 years, 10 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.13: +11 -11 lines 
comment out test on inAdMode (has no effect + OLD_DST3_FORMULATION no longer used)
1 jmc 1.14 C $Header: /u/gcmpack/MITgcm/pkg/generic_advdiff/gad_dst3_adv_y.F,v 1.13 2011/03/29 15:47:19 jmc Exp $
2 jmc 1.3 C $Name: $
3 adcroft 1.1
4     #include "GAD_OPTIONS.h"
5    
6 jmc 1.10 CBOP
7     C !ROUTINE: GAD_DST3_ADV_Y
8    
9     C !INTERFACE: ==========================================================
10 jmc 1.8 SUBROUTINE GAD_DST3_ADV_Y(
11 jmc 1.12 I bi,bj,k, calcCFL, deltaTloc,
12 jmc 1.8 I vTrans, vFld,
13 jmc 1.4 I maskLocS, tracer,
14 adcroft 1.1 O vT,
15     I myThid )
16 jmc 1.10 C !DESCRIPTION:
17     C Calculates the area integrated Meridional flux due to advection of a
18     C tracer using 3rd-order Direct Space and Time (DST-3) Advection Scheme
19    
20     C !USES: ===============================================================
21 adcroft 1.1 IMPLICIT NONE
22    
23     C == GLobal variables ==
24     #include "SIZE.h"
25 jmc 1.13 #ifdef OLD_DST3_FORMULATION
26 heimbach 1.5 #include "EEPARAMS.h"
27     #include "PARAMS.h"
28 jmc 1.13 #endif
29     #include "GRID.h"
30 adcroft 1.1 #include "GAD.h"
31    
32     C == Routine arguments ==
33 jmc 1.10 C !INPUT PARAMETERS: ===================================================
34     C bi,bj :: tile indices
35     C k :: vertical level
36 jmc 1.12 C calcCFL :: =T: calculate CFL number ; =F: take vFld as CFL
37 jmc 1.10 C deltaTloc :: local time-step (s)
38     C vTrans :: meridional volume transport
39 jmc 1.12 C vFld :: meridional flow / CFL number
40 jmc 1.10 C tracer :: tracer field
41     C myThid :: thread number
42 adcroft 1.1 INTEGER bi,bj,k
43 jmc 1.12 LOGICAL calcCFL
44 heimbach 1.5 _RL deltaTloc
45 adcroft 1.1 _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
46 jmc 1.8 _RL vFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
47 jmc 1.4 _RS maskLocS(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
48 adcroft 1.1 _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
49 jmc 1.10 INTEGER myThid
50    
51     C !OUTPUT PARAMETERS: ==================================================
52     C vT :: meridional advective flux
53 adcroft 1.1 _RL vT (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
54    
55     C == Local variables ==
56 jmc 1.10 C !LOCAL VARIABLES: ====================================================
57     C i,j :: loop indices
58 jmc 1.11 C vCFL :: Courant-Friedrich-Levy number
59 adcroft 1.1 INTEGER i,j
60 jmc 1.11 _RL Rjm,Rj,Rjp,vCFL,d0,d1
61 jmc 1.10 #ifdef OLD_DST3_FORMULATION
62 adcroft 1.2 _RL psiP,psiM,thetaP,thetaM
63 heimbach 1.5 _RL smallNo
64 jmc 1.6 c _RL Rjjm,Rjjp
65 heimbach 1.5
66 jmc 1.14 c IF (inAdMode) THEN
67     c smallNo = 1.0D-20
68     c ELSE
69 heimbach 1.5 smallNo = 1.0D-20
70 jmc 1.14 c ENDIF
71 jmc 1.10 #endif
72 adcroft 1.1
73 jmc 1.14 DO i=1-OLx,sNx+OLx
74     vT(i,1-OLy)=0.
75     vT(i,2-OLy)=0.
76     vT(i,sNy+OLy)=0.
77 adcroft 1.1 ENDDO
78 jmc 1.14 DO j=1-OLy+2,sNy+OLy-1
79     DO i=1-OLx,sNx+OLx
80 jmc 1.4 Rjp=(tracer(i,j+1)-tracer(i, j ))*maskLocS(i,j+1)
81     Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j )
82     Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1)
83 adcroft 1.1
84 jmc 1.12 vCFL = vFld(i,j)
85     IF ( calcCFL ) vCFL = ABS( vFld(i,j)*deltaTloc
86 jmc 1.11 & *recip_dyC(i,j,bi,bj)*recip_deepFacC(k) )
87     d0=(2.-vCFL)*(1.-vCFL)*oneSixth
88     d1=(1.-vCFL*vCFL)*oneSixth
89 jmc 1.10 #ifdef OLD_DST3_FORMULATION
90 heimbach 1.5 IF ( ABS(Rj).LT.smallNo .OR.
91     & ABS(Rjm).LT.smallNo ) THEN
92     thetaP=0.
93     psiP=0.
94 jmc 1.8 ELSE
95 heimbach 1.5 thetaP=(Rjm+smallNo)/(smallNo+Rj)
96     psiP=d0+d1*thetaP
97     ENDIF
98     IF ( ABS(Rj).LT.smallNo .OR.
99     & ABS(Rjp).LT.smallNo ) THEN
100     thetaM=0.
101     psiM=0.
102     ELSE
103     thetaM=(Rjp+smallNo)/(smallNo+Rj)
104     psiM=d0+d1*thetaM
105     ENDIF
106 adcroft 1.1 vT(i,j)=
107 jmc 1.10 & 0.5*(vTrans(i,j)+ABS(vTrans(i,j)))
108 adcroft 1.2 & *( Tracer(i,j-1) + psiP*Rj )
109 jmc 1.10 & +0.5*(vTrans(i,j)-ABS(vTrans(i,j)))
110 adcroft 1.2 & *( Tracer(i, j ) - psiM*Rj )
111 jmc 1.10 #else /* OLD_DST3_FORMULATION */
112     vT(i,j)=
113     & 0.5*(vTrans(i,j)+ABS(vTrans(i,j)))
114     & *( Tracer(i,j-1) + (d0*Rj+d1*Rjm) )
115     & +0.5*(vTrans(i,j)-ABS(vTrans(i,j)))
116     & *( Tracer(i, j ) - (d0*Rj+d1*Rjp) )
117     #endif /* OLD_DST3_FORMULATION */
118 adcroft 1.1
119     ENDDO
120     ENDDO
121    
122     RETURN
123     END
  ViewVC Help
Powered by ViewVC 1.1.22