36 |
C o Constant F case |
C o Constant F case |
37 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
38 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
39 |
DO K=1,Nz |
DO K=1,Nr |
40 |
DO J=1,sNy |
DO J=1,sNy |
41 |
DO I=1,sNx |
DO I=1,sNx |
42 |
fCori(i,j,bi,bj)=f0 |
fCori(i,j,bi,bj)=f0 |
51 |
IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere |
IF ( usingSphericalPolarGrid ) facGrid = deg2rad*rSphere |
52 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
53 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
54 |
DO K=1,Nz |
DO K=1,Nr |
55 |
DO J=1,sNy |
DO J=1,sNy |
56 |
DO I=1,sNx |
DO I=1,sNx |
57 |
fCori(i,j,bi,bj)=f0+beta*_yC(i,j,bi,bj)*facGrid |
fCori(i,j,bi,bj)=f0+beta*_yC(i,j,bi,bj)*facGrid |
65 |
C Note in this case we assume yC is in degrees. |
C Note in this case we assume yC is in degrees. |
66 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
67 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
68 |
DO K=1,Nz |
DO K=1,Nr |
69 |
DO J=1,sNy |
DO J=1,sNy |
70 |
DO I=1,sNx |
DO I=1,sNx |
71 |
fCori(i,j,bi,bj)=2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad) |
fCori(i,j,bi,bj)=2. _d 0*omega*sin(_yC(i,j,bi,bj)*deg2rad) |
78 |
C o Special custom form |
C o Special custom form |
79 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
80 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
81 |
DO K=1,Nz |
DO K=1,Nr |
82 |
DO J=1,sNy |
DO J=1,sNy |
83 |
DO I=1,sNx |
DO I=1,sNx |
84 |
fCori(i,j,bi,bj)=0. |
fCori(i,j,bi,bj)=0. |