/[MITgcm]/MITgcm/pkg/mom_vecinv/mom_vi_del2uv.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_vecinv/mom_vi_del2uv.F

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

revision 1.1 by adcroft, Fri Mar 30 21:13:33 2001 UTC revision 1.2 by adcroft, Tue May 29 14:01:39 2001 UTC
# Line 0  Line 1 
1    C $Header$
2    C $Name$
3    
4    #include "CPP_OPTIONS.h"
5    
6          SUBROUTINE MOM_VI_DEL2UV(
7         I        bi,bj,k,
8         I        hDiv,vort3,hFacZ,
9         O        del2u,del2v,
10         I        myThid)
11          IMPLICIT NONE
12    C
13    C     Calculate del^2 of (u,v) in terms of hDiv and vort3
14    C
15    
16    C     == Global variables ==
17    #include "SIZE.h"
18    #include "GRID.h"
19    
20    C     == Routine arguments ==
21          INTEGER bi,bj,k
22          _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
23          _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
24          _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
25          _RL del2u(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
26          _RL del2v(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
27          INTEGER myThid
28    
29    C     == Local variables ==
30          INTEGER I,J
31          _RL Zip,Zij,Zpj,Dim,Dij,Dmj
32    
33    C     - Laplacian  and bi-harmonic terms
34          DO j=2-Oly,sNy+Oly-1
35           DO i=2-Olx,sNx+Olx-1
36    
37    c       Dim=dyF( i ,j-1,bi,bj)*hFacC( i ,j-1,k,bi,bj)*hDiv( i ,j-1)
38    c       Dij=dyF( i , j ,bi,bj)*hFacC( i , j ,k,bi,bj)*hDiv( i , j )
39    c       Dmj=dyF(i-1, j ,bi,bj)*hFacC(i-1, j ,k,bi,bj)*hDiv(i-1, j )
40    c       Dim=dyF( i ,j-1,bi,bj)*                       hDiv( i ,j-1)
41    c       Dij=dyF( i , j ,bi,bj)*                       hDiv( i , j )
42    c       Dmj=dyF(i-1, j ,bi,bj)*                       hDiv(i-1, j )
43            Dim=                                          hDiv( i ,j-1)
44            Dij=                                          hDiv( i , j )
45            Dmj=                                          hDiv(i-1, j )
46    
47    c       Zip=dxV( i ,j+1,bi,bj)*hFacZ( i ,j+1)*vort3( i ,j+1)
48    c       Zij=dxV( i , j ,bi,bj)*hFacZ( i , j )*vort3( i , j )
49    c       Zpj=dxV(i+1, j ,bi,bj)*hFacZ(i+1, j )*vort3(i+1, j )
50            Zip=                   hFacZ( i ,j+1)*vort3( i ,j+1)
51            Zij=                   hFacZ( i , j )*vort3( i , j )
52            Zpj=                   hFacZ(i+1, j )*vort3(i+1, j )
53    
54    c       del2u(i,j) = recip_rAw(i,j,bi,bj)*(
55    c    &   +recip_hFacW(i,j,k,bi,bj)*( Dij-Dmj )
56    c    &   -recip_hFacW(i,j,k,bi,bj)*( Zip-Zij ) )
57    c       del2u(i,j) = recip_rAw(i,j,bi,bj)*(
58    c    &   +                         ( Dij-Dmj )
59    c    &   -recip_hFacW(i,j,k,bi,bj)*( Zip-Zij ) )
60            del2u(i,j) =
61         &   +                         ( Dij-Dmj )*recip_DXC(i,j,bi,bj)
62         &   -recip_hFacW(i,j,k,bi,bj)*( Zip-Zij )*recip_DYU(i,j,bi,bj)
63    
64    c       del2v(i,j) = recip_rAs(i,j,bi,bj)*(
65    c    &    recip_hFacS(i,j,k,bi,bj)*( Zpj-Zij )
66    c    &   +recip_hFacS(i,j,k,bi,bj)*( Dij-Dim ) )
67    c       del2v(i,j) = recip_rAs(i,j,bi,bj)*(
68    c    &    recip_hFacS(i,j,k,bi,bj)*( Zpj-Zij )
69    c    &   +                         ( Dij-Dim ) )
70            del2v(i,j) =
71         &    recip_hFacS(i,j,k,bi,bj)*( Zpj-Zij )*recip_DXV(i,j,bi,bj)
72         &   +                         ( Dij-Dim )*recip_DYC(i,j,bi,bj)
73    
74           ENDDO
75          ENDDO
76    
77          RETURN
78          END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22