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

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

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

revision 1.3 by adcroft, Tue Jul 20 17:46:38 2004 UTC revision 1.4 by heimbach, Fri Dec 10 20:15:10 2004 UTC
# Line 65  c    &               +vFld(i-1,j)+vFld(i Line 65  c    &               +vFld(i-1,j)+vFld(i
65            vort3u=0.5*(omega3(i,j)*r_hFacZ(i,j)            vort3u=0.5*(omega3(i,j)*r_hFacZ(i,j)
66       &               +omega3(i,j+1)*r_hFacZ(i,j+1))       &               +omega3(i,j+1)*r_hFacZ(i,j+1))
67           ENDIF           ENDIF
68    
69    cph(
70    cph The following block with 'interlaced' comments
71    cph is bad for TAMC and will be replaced
72    cph)
73    #undef THIS_IS_BAD_FOR_TAMC
74    #ifdef THIS_IS_BAD_FOR_TAMC
75           uCoriolisTerm(i,j)=           uCoriolisTerm(i,j)=
76  C high order vorticity advection term  C high order vorticity advection term
77       &   +vort3u*vBarXY*recip_dxc(i,j,bi,bj)       &   +vort3u*vBarXY*recip_dxc(i,j,bi,bj)
# Line 78  c    &          fCori(I-1,J,bi,bj)*vBarY Line 85  c    &          fCori(I-1,J,bi,bj)*vBarY
85  C original form of Coriolis term (copied from calc_mom_rhs)  C original form of Coriolis term (copied from calc_mom_rhs)
86  c    &   +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY  c    &   +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY
87       &   *_maskW(I,J,K,bi,bj)       &   *_maskW(I,J,K,bi,bj)
88    #else
89             uCoriolisTerm(i,j)=
90         &   +vort3u*vBarXY*recip_dxc(i,j,bi,bj)
91         &   *_maskW(I,J,K,bi,bj)
92    #endif
93    
94          ENDDO          ENDDO
95         ENDDO         ENDDO
96    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22