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 |
|
|
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) |