Parent Directory
|
Revision Log
|
Revision Graph
Merge from branch pre38: o essential mods for cubed sphere o debugged atmosphere, dynamcis + physics (aim) o new packages (mom_vecinv, mom_fluxform, ...)
1 | C $Header: /u/gcmpack/models/MITgcmUV/pkg/mom_fluxform/Attic/mom_v_adv_uv.F,v 1.1.2.1 2001/03/28 19:51:14 adcroft Exp $ |
2 | C $Name: pre38-close $ |
3 | |
4 | #include "CPP_OPTIONS.h" |
5 | |
6 | SUBROUTINE MOM_V_ADV_UV( |
7 | I bi,bj,k, |
8 | I uTrans, vFld, |
9 | O AdvectFluxUV, |
10 | I myThid) |
11 | IMPLICIT NONE |
12 | C |
13 | C Calculate zonal advective flux of V using centered second order diff. |
14 | C - is not vector invariant ... |
15 | C |
16 | |
17 | C == Global variables == |
18 | #include "SIZE.h" |
19 | #include "EEPARAMS.h" |
20 | #include "PARAMS.h" |
21 | #include "GRID.h" |
22 | |
23 | C == Routine arguments == |
24 | INTEGER bi,bj,k |
25 | _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
26 | _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
27 | _RL AdvectFluxUV(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
28 | INTEGER myThid |
29 | |
30 | C == Local variables == |
31 | INTEGER I,J |
32 | |
33 | DO j=1-Oly+1,sNy+Oly |
34 | DO i=1-Olx+1,sNx+Olx |
35 | AdvectFluxUV(i,j) = |
36 | & 0.25*( uTrans(i,j) + uTrans(i,j-1) ) |
37 | & *( vFld(i,j) + vFld(i-1,j) ) |
38 | #ifdef OLD_ADV_BCS |
39 | & *_maskS(i,j,k,bi,bj) |
40 | & *_maskS(i-1,j,k,bi,bj) |
41 | #endif /* OLD_ADV_BCS */ |
42 | ENDDO |
43 | ENDDO |
44 | |
45 | RETURN |
46 | END |
ViewVC Help | |
Powered by ViewVC 1.1.22 |