/[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.7 by jmc, Tue Apr 28 18:06:14 2009 UTC revision 1.8 by jmc, Sun Sep 27 23:49:28 2009 UTC
# Line 10  C     !INTERFACE: Line 10  C     !INTERFACE:
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    
# Line 41  C     myThid -  Number of this instance Line 41  C     myThid -  Number of this instance
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 ==
# Line 54  C     I,J,K Line 54  C     I,J,K
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
# Line 83  c           theta(I,J,K,bi,bj) = tRef(K) Line 81  c           theta(I,J,K,bi,bj) = tRef(K)
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
# Line 127  C--   Apply mask and test consistency Line 124  C--   Apply mask and test consistency
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
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22