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

Contents of /MITgcm/pkg/generic_advdiff/gad_c2_adv_y.F

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

Revision 1.1 - (show annotations) (download)
Wed May 30 19:34:48 2001 UTC (22 years, 11 months ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint40pre3, checkpoint40pre1, checkpoint40pre7, checkpoint40pre6, checkpoint40pre9, checkpoint40pre8, checkpoint40pre2, checkpoint40pre4, checkpoint40pre5, checkpoint40 
Added "gad" package. Needs generalizing to allow selection
of advection schemes at run-time and different schemes for
each tracer.
1 C $Header: $
2 C $Name: $
3
4 #include "GAD_OPTIONS.h"
5
6 SUBROUTINE GAD_C2_ADV_Y(
7 I bi,bj,k,
8 I vTrans,
9 I tracer,
10 O vT,
11 I myThid )
12 IMPLICIT NONE
13
14 C == GLobal variables ==
15 #include "SIZE.h"
16 #include "GRID.h"
17
18 C == Routine arguments ==
19 INTEGER bi,bj,k
20 _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
21 _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
22 _RL vT (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
23 INTEGER myThid
24
25 C == Local variables ==
26 INTEGER i,j
27
28 DO i=1-Olx,sNx+Olx
29 vT(i,1-Oly)=0.
30 ENDDO
31 DO j=1-Oly+1,sNy+Oly
32 DO i=1-Olx,sNx+Olx
33 vT(i,j) =
34 & vTrans(i,j)*(Tracer(i,j)+Tracer(i,j-1))*0.5 _d 0
35 ENDDO
36 ENDDO
37
38 RETURN
39 END
  ViewVC Help
Powered by ViewVC 1.1.22