49 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
50 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
51 |
DO K=1,Nr |
DO K=1,Nr |
52 |
rSurf=1.E5 |
rSurf=1. _d 5 |
53 |
thKappa = 2./7. |
thKappa = 2. _d 0/7. _d 0 |
54 |
thetaLim = 200. / ((rC(K)/rSurf)**thKappa) |
thetaLim = 200. _d 0 / ((rC(K)/rSurf)**thKappa) |
55 |
C thetaLim = 170. / ((rC(K)/rSurf)**thKappa) |
C thetaLim = 170. _d 0 / ((rC(K)/rSurf)**thKappa) |
56 |
DO J=1,sNy |
DO J=1,sNy |
57 |
term1=60.*(sin(yC(1,J,bi,bj)*deg2rad)**2) |
term1=60. _d 0*(sin(yC(1,J,bi,bj)*deg2rad)**2) |
58 |
term2=10.*log((rC(K)/rSurf)) |
term2=10. _d 0*log((rC(K)/rSurf)) |
59 |
& *(cos(yC(1,J,bi,bj)*deg2rad)**2) |
& *(cos(yC(1,J,bi,bj)*deg2rad)**2) |
60 |
thetaEq=315.-term1-term2 |
thetaEq=315. _d 0-term1-term2 |
61 |
DO I=1,sNx |
DO I=1,sNx |
62 |
theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq ) |
theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq ) |
63 |
C theta(I,J,K,bi,bj) = tRef(K) |
C theta(I,J,K,bi,bj) = tRef(K) |