/[MITgcm]/MITgcm/verification/hs94.cs-32x32x5/code/ini_theta.F
ViewVC logotype

Diff of /MITgcm/verification/hs94.cs-32x32x5/code/ini_theta.F

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

revision 1.2 by adcroft, Tue May 29 14:01:58 2001 UTC revision 1.3 by adcroft, Wed Jun 6 19:46:43 2001 UTC
# Line 54  C--    Initialise temperature field to H Line 54  C--    Initialise temperature field to H
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)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22