/[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.16 by cnh, Wed Sep 26 18:09:15 2001 UTC revision 1.17 by dimitri, Thu Oct 23 07:14:49 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        _RL     rad, rD, Tfreezing
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    
103    C--   Check that there are no values of temperature below freezing point.
104          Tfreezing=-1.9
105          IF ( allowFreezing ) THEN
106           DO bj = myByLo(myThid), myByHi(myThid)
107            DO bi = myBxLo(myThid), myBxHi(myThid)
108             DO K=1,Nr
109              DO J=1-Oly,sNy+Oly
110               DO I=1-Olx,sNx+Olx
111                IF (theta(I,J,k,bi,bj) .LT. Tfreezing) THEN
112                   theta(I,J,K,bi,bj) = Tfreezing
113                ENDIF
114               ENDDO
115              ENDDO
116             ENDDO
117            ENDDO
118           ENDDO
119          ENDIF
120    
121        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
122       &                       Nr, 1, myThid )       &                       Nr, 1, myThid )
123    
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22