/[MITgcm]/MITgcm/verification/advect_xz/code/ini_vel.F
ViewVC logotype

Diff of /MITgcm/verification/advect_xz/code/ini_vel.F

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

revision 1.1 by adcroft, Fri Sep 28 02:28:10 2001 UTC revision 1.2 by adcroft, Thu May 15 15:20:26 2003 UTC
# Line 25  C     bi,bj,i,j,k - Loop counters Line 25  C     bi,bj,i,j,k - Loop counters
25    
26        _RL psi,Lfac,Rpi        _RL psi,Lfac,Rpi
27        psi(i,j,k,bi,bj)=Lfac*sin(XG(i,j,bi,bj)*Rpi/Lfac)        psi(i,j,k,bi,bj)=Lfac*sin(XG(i,j,bi,bj)*Rpi/Lfac)
28       &    *sin( max(RF(k),R_low(i,j,bi,bj))*Rpi*recip_Rcol(i,j,bi,bj))       &    *sin( max(RF(k),
29         &       max(R_low(i,j,bi,bj),R_low(i-1,j,bi,bj))
30         &                   )*Rpi*recip_Rcol(i,j,bi,bj))
31    
32        Rpi=3.14159265358979        Rpi=3.14159265358979 _d 0
33        Lfac=(DELX(1)*float(Nx))        Lfac=(DELX(1)*float(Nx))
34    
35    
# Line 42  C Line 44  C
44         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
45          DO k=1,Nr          DO k=1,Nr
46           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
47            DO i=1-Olx,sNx+Olx            DO i=1-Olx+1,sNx+Olx
48             uVel (I,J,K,bi,bj) = 0. _d 0             uVel (I,J,K,bi,bj) = 0. _d 0
49       &         + (psi(I,J,K,bi,bj)-psi(I,J,K+1,bi,bj))*recip_dRF(k)       &         + (psi(I,J,K,bi,bj)-psi(I,J,K+1,bi,bj))*recip_dRF(k)
50       &         *recip_hFacW(i,j,k,bi,bj)       &         *recip_hFacW(i,j,k,bi,bj)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22