10 |
SUBROUTINE INI_THETA( myThid ) |
SUBROUTINE INI_THETA( myThid ) |
11 |
C !DESCRIPTION: \bv |
C !DESCRIPTION: \bv |
12 |
C *==========================================================* |
C *==========================================================* |
13 |
C | SUBROUTINE INI_THETA |
C | SUBROUTINE INI_THETA |
14 |
C | o Set model initial temperature field. |
C | o Set model initial temperature field. |
15 |
C *==========================================================* |
C *==========================================================* |
16 |
C | There are several options for setting the initial |
C | There are several options for setting the initial |
17 |
C | temperature file |
C | temperature file |
18 |
C | 1. Inline code |
C | 1. Inline code |
19 |
C | 2. Vertical profile ( uniform T in X and Y ) |
C | 2. Vertical profile ( uniform T in X and Y ) |
20 |
C | 3. Three-dimensional data from a file. For example from |
C | 3. Three-dimensional data from a file. For example from |
21 |
C | Levitus or from a checkpoint file from a previous |
C | Levitus or from a checkpoint file from a previous |
22 |
C | integration. |
C | integration. |
23 |
C | In addition to setting the temperature field we also |
C | In addition to setting the temperature field we also |
24 |
C | set the initial temperature tendency term here. |
C | set the initial temperature tendency term here. |
25 |
C *==========================================================* |
C *==========================================================* |
26 |
C \ev |
C \ev |
27 |
|
|
41 |
INTEGER myThid |
INTEGER myThid |
42 |
|
|
43 |
C == Functions == |
C == Functions == |
44 |
real*8 PORT_RAND |
c real*8 PORT_RAND |
45 |
real*8 seed |
c real*8 seed |
46 |
|
|
47 |
C !LOCAL VARIABLES: |
C !LOCAL VARIABLES: |
48 |
C == Local variables == |
C == Local variables == |
54 |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
55 |
CEOP |
CEOP |
56 |
|
|
|
_BARRIER |
|
|
|
|
57 |
J = 99+myBxLo(myThid)+nPx*myByLo(myThid) |
J = 99+myBxLo(myThid)+nPx*myByLo(myThid) |
58 |
c CALL SRAND( J ) |
c CALL SRAND( J ) |
59 |
c seed = j |
c seed = j |
81 |
C-- Zonal FFT filter initial conditions |
C-- Zonal FFT filter initial conditions |
82 |
IF (useZONAL_FILT) THEN |
IF (useZONAL_FILT) THEN |
83 |
CALL ZONAL_FILTER( |
CALL ZONAL_FILTER( |
84 |
U theta, hFacC, |
U theta(1-OLx,1-OLy,1,bi,bj), |
85 |
I 1, sNy, 1, Nr, bi, bj, 1, myThid) |
I hFacC(1-OLx,1-OLy,1,bi,bj), |
86 |
|
I 1, sNy, Nr, bi, bj, 1, myThid ) |
87 |
ENDIF |
ENDIF |
88 |
#endif /* ALLOW_ZONAL_FILT */ |
#endif /* ALLOW_ZONAL_FILT */ |
89 |
ENDDO |
ENDDO |
90 |
ENDDO |
ENDDO |
91 |
ELSE |
ELSE |
|
_BEGIN_MASTER( myThid ) |
|
92 |
CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid ) |
CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid ) |
|
_END_MASTER(myThid) |
|
93 |
ENDIF |
ENDIF |
94 |
C-- Apply mask and test consistency |
C-- Apply mask and test consistency |
95 |
localWarnings=0 |
localWarnings=0 |
124 |
|
|
125 |
_EXCH_XYZ_RL(theta , myThid ) |
_EXCH_XYZ_RL(theta , myThid ) |
126 |
|
|
127 |
CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' , |
IF (debugMode) THEN |
128 |
& Nr, 1, myThid ) |
CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' , |
129 |
|
& Nr, 1, myThid ) |
130 |
|
ENDIF |
131 |
|
|
132 |
RETURN |
RETURN |
133 |
END |
END |