/[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.4 by jmc, Fri Jul 6 22:13:37 2001 UTC revision 1.5 by molod, Mon Jul 19 15:13:07 2004 UTC
# Line 35  CEndOfInterface Line 35  CEndOfInterface
35    
36  C     == Functions ==  C     == Functions ==
37        Real*8  PORT_RAND        Real*8  PORT_RAND
38          real *8 seed
39    
40  C     == Local variables ==  C     == Local variables ==
41  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
# Line 47  C     I,J,K Line 48  C     I,J,K
48    
49        J = 99+myBxLo(myThid)+nPx*myByLo(myThid)        J = 99+myBxLo(myThid)+nPx*myByLo(myThid)
50  c     CALL SRAND( J )  c     CALL SRAND( J )
51    c     seed = j
52    
53        IF ( hydrogThetaFile .EQ. ' ' ) THEN        IF ( hydrogThetaFile .EQ. ' ' ) THEN
54  C--    Initialise temperature field to Held & Suarez equilibrium theta  C--    Initialise temperature field to Held & Suarez equilibrium theta
# Line 62  C--    Initialise temperature field to H Line 64  C--    Initialise temperature field to H
64             thetaEq=315. _d 0-term1-term2             thetaEq=315. _d 0-term1-term2
65              theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq )              theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq )
66  c    &                          + 0.01*(RAND()-0.5)  c    &                          + 0.01*(RAND()-0.5)
67  c    &                          + 0.01*(PORT_RAND()-0.5)  c    &                          + 0.01*(PORT_RAND(seed)-0.5)
68  c           theta(I,J,K,bi,bj) = tRef(K)  c           theta(I,J,K,bi,bj) = tRef(K)
69             ENDDO             ENDDO
70            ENDDO            ENDDO
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22