/[MITgcm]/MITgcm/verification/hs94.1x64x5/code/ini_theta.F
ViewVC logotype

Diff of /MITgcm/verification/hs94.1x64x5/code/ini_theta.F

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

revision 1.3 by cnh, Sun Feb 4 14:38:53 2001 UTC revision 1.4 by adcroft, Mon Jun 4 20:20:18 2001 UTC
# Line 49  C--    Initialise temperature field to H Line 49  C--    Initialise temperature field to H
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)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22