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

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

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

revision 1.2 by adcroft, Tue May 29 14:01:39 2001 UTC revision 1.3 by heimbach, Wed Sep 5 17:46:03 2001 UTC
# Line 67  C     Special stuff for Cubed Sphere Line 67  C     Special stuff for Cubed Sphere
67           vort3(I,J)=           vort3(I,J)=
68       &     +recip_rAz(I,J,bi,bj)*(       &     +recip_rAz(I,J,bi,bj)*(
69       &      vFld(I,J)*dyc(I,J,bi,bj)       &      vFld(I,J)*dyc(I,J,bi,bj)
 c    &     -vFld(I-1,J)*dyc(I-1,J,bi,bj)  
70       &     -uFld(I,J)*dxc(I,J,bi,bj)       &     -uFld(I,J)*dxc(I,J,bi,bj)
71       &     +uFld(I,J-1)*dxc(I,J-1,bi,bj)       &     +uFld(I,J-1)*dxc(I,J-1,bi,bj)
72       &     )       &     )
73    cph    &     -vFld(I-1,J)*dyc(I-1,J,bi,bj)
74           I=sNx+1           I=sNx+1
75           J=1           J=1
76           vort3(I,J)=           vort3(I,J)=
77       &     +recip_rAz(I,J,bi,bj)*(       &     +recip_rAz(I,J,bi,bj)*(
 c    &      vFld(I,J)*dyc(I,J,bi,bj)  
78       &     -vFld(I-1,J)*dyc(I-1,J,bi,bj)       &     -vFld(I-1,J)*dyc(I-1,J,bi,bj)
79       &     -uFld(I,J)*dxc(I,J,bi,bj)       &     -uFld(I,J)*dxc(I,J,bi,bj)
80       &     +uFld(I,J-1)*dxc(I,J-1,bi,bj)       &     +uFld(I,J-1)*dxc(I,J-1,bi,bj)
81       &     )       &     )
82    cph    &      vFld(I,J)*dyc(I,J,bi,bj)
83           I=1           I=1
84           J=sNy+1           J=sNy+1
85           vort3(I,J)=           vort3(I,J)=
86       &     +recip_rAz(I,J,bi,bj)*(       &     +recip_rAz(I,J,bi,bj)*(
87       &      vFld(I,J)*dyc(I,J,bi,bj)       &      vFld(I,J)*dyc(I,J,bi,bj)
 c    &     -vFld(I-1,J)*dyc(I-1,J,bi,bj)  
88       &     -uFld(I,J)*dxc(I,J,bi,bj)       &     -uFld(I,J)*dxc(I,J,bi,bj)
89       &     +uFld(I,J-1)*dxc(I,J-1,bi,bj)       &     +uFld(I,J-1)*dxc(I,J-1,bi,bj)
90       &     )       &     )
91    cph    &     -vFld(I-1,J)*dyc(I-1,J,bi,bj)
92           I=sNx+1           I=sNx+1
93           J=sNy+1           J=sNy+1
94           vort3(I,J)=           vort3(I,J)=
95       &     +recip_rAz(I,J,bi,bj)*(       &     +recip_rAz(I,J,bi,bj)*(
 c    &      vFld(I,J)*dyc(I,J,bi,bj)  
96       &     -vFld(I-1,J)*dyc(I-1,J,bi,bj)       &     -vFld(I-1,J)*dyc(I-1,J,bi,bj)
97       &     -uFld(I,J)*dxc(I,J,bi,bj)       &     -uFld(I,J)*dxc(I,J,bi,bj)
98       &     +uFld(I,J-1)*dxc(I,J-1,bi,bj)       &     +uFld(I,J-1)*dxc(I,J-1,bi,bj)
99       &     )       &     )
100    cph    &      vFld(I,J)*dyc(I,J,bi,bj)
101         ENDIF         ENDIF
102    
103        RETURN        RETURN
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22