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 |
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 |
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 |