54 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
55 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
56 |
DO K=1,Nr |
DO K=1,Nr |
57 |
Ro_SeaLevel=1.E5 |
Ro_SeaLevel=1. _d 5 |
58 |
thKappa = 2./7. |
thKappa = 2. _d 0/7. _d 0 |
59 |
thetaLim = 200. / ((rC(K)/Ro_SeaLevel)**thKappa) |
thetaLim = 200. _d 0/((rC(K)/Ro_SeaLevel)**thKappa) |
60 |
DO J=1,sNy |
DO J=1,sNy |
61 |
DO I=1,sNx |
DO I=1,sNx |
62 |
term1=60.*(sin(yC(I,J,bi,bj)*deg2rad)**2) |
term1=60. _d 0*(sin(yC(I,J,bi,bj)*deg2rad)**2) |
63 |
term2=10.*log((rC(K)/Ro_SeaLevel)) |
term2=10. _d 0*log((rC(K)/Ro_SeaLevel)) |
64 |
& *(cos(yC(I,J,bi,bj)*deg2rad)**2) |
& *(cos(yC(I,J,bi,bj)*deg2rad)**2) |
65 |
thetaEq=315.-term1-term2 |
thetaEq=315. _d 0-term1-term2 |
66 |
theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq ) |
theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq ) |
67 |
c & + 0.01*(RAND()-0.5) |
c & + 0.01*(RAND()-0.5) |
68 |
c & + 0.01*(PORT_RAND()-0.5) |
c & + 0.01*(PORT_RAND()-0.5) |