/[MITgcm]/MITgcm/model/src/ini_theta.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_theta.F

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

revision 1.17 by dimitri, Thu Oct 23 07:14:49 2003 UTC revision 1.18 by edhill, Fri Oct 24 05:29:35 2003 UTC
# Line 51  C     I,J,K Line 51  C     I,J,K
51        INTEGER iG, jG        INTEGER iG, jG
52        INTEGER bi, bj        INTEGER bi, bj
53        INTEGER I, J, K, localWarnings        INTEGER I, J, K, localWarnings
54        _RL     rad, rD, Tfreezing        _RL     rad, rD
55        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
56  CEOP  CEOP
57    
# Line 75  C--   Initialise temperature field to th Line 75  C--   Initialise temperature field to th
75         _EXCH_XYZ_R8(theta,myThid)         _EXCH_XYZ_R8(theta,myThid)
76        ENDIF        ENDIF
77    
78  C--   Apply mask and test consistancy  C     Apply mask and test consistancy
79        localWarnings=0        localWarnings=0
80        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
81         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 100  C--   Apply mask and test consistancy Line 100  C--   Apply mask and test consistancy
100         STOP 'ABNORMAL END: S/R INI_THETA'         STOP 'ABNORMAL END: S/R INI_THETA'
101        ENDIF        ENDIF
102    
 C--   Check that there are no values of temperature below freezing point.  
       Tfreezing=-1.9  
       IF ( allowFreezing ) THEN  
        DO bj = myByLo(myThid), myByHi(myThid)  
         DO bi = myBxLo(myThid), myBxHi(myThid)  
          DO K=1,Nr  
           DO J=1-Oly,sNy+Oly  
            DO I=1-Olx,sNx+Olx  
             IF (theta(I,J,k,bi,bj) .LT. Tfreezing) THEN  
                theta(I,J,K,bi,bj) = Tfreezing  
             ENDIF  
            ENDDO  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDIF  
   
103        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
104       &                       Nr, 1, myThid )       &                       Nr, 1, myThid )
105    
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22