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

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

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

Revision 1.2 - (hide annotations) (download)
Tue May 29 14:01:39 2001 UTC (23 years ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint40pre3, checkpoint40pre1, checkpoint40pre7, checkpoint40pre6, checkpoint40pre9, checkpoint40pre8, release1_b1, checkpoint43, checkpoint40pre2, checkpoint40pre4, checkpoint40pre5, ecco-branch-mod1, release1_beta1, checkpoint42, checkpoint40, checkpoint41
Branch point for: release1, ecco-branch, release1_coupled
Changes since 1.1: +78 -0 lines 
Merge from branch pre38:
 o essential mods for cubed sphere
 o debugged atmosphere, dynamcis + physics (aim)
 o new packages (mom_vecinv, mom_fluxform, ...)
1 adcroft 1.2 C $Header: /u/gcmpack/models/MITgcmUV/pkg/mom_vecinv/Attic/mom_vi_del2uv.F,v 1.1.2.2 2001/05/07 16:19:17 adcroft Exp $
2     C $Name: pre38-close $
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
  ViewVC Help
Powered by ViewVC 1.1.22