model/src/freeze.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/freeze.F,v 1.16 1998/11/06 22:44:49 cnh Exp $

#include "CPP_OPTIONS.h"

C     /==========================================================\
C     | S/R FREEZE                                               |
C     | o Check water temperature and limit range of temperature |
C     | appropriately.                                           |
C     \==========================================================/
      SUBROUTINE FREEZE( bi, bj, iMin, iMax, jMin, jMax,
     I                     K,
     I                     myThid )
      implicit none
C Common
#include "SIZE.h"
#include "DYNVARS.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
C     == Routine Arguments ==
      INTEGER bi,bj,iMin,iMax,jMin,jMax
      INTEGER K
      INTEGER myThid
C     == Local variables ==
      INTEGER i,j
      _RL Tfreezing

      Tfreezing=-1.9

C     Check for water that should have frozen
       DO j=jMin,jMax
        DO i=iMin,iMax
         IF (gTNm1(i,j,k,bi,bj) .LT. Tfreezing) THEN
          gTNm1(i,j,k,bi,bj)=Tfreezing
         ENDIF
        ENDDO
       ENDDO

      RETURN
      END
1	adcroft	1.1	C $Header: /u/gcmpack/models/MITgcmUV/model/src/freeze.F,v 1.16 1998/11/06 22:44:49 cnh Exp $
2
3			#include "CPP_OPTIONS.h"
4
5			C /==========================================================\
6			C \| S/R FREEZE \|
7			C \| o Check water temperature and limit range of temperature \|
8			C \| appropriately. \|
9			C \==========================================================/
10			SUBROUTINE FREEZE( bi, bj, iMin, iMax, jMin, jMax,
11			I K,
12			I myThid )
13			implicit none
14			C Common
15			#include "SIZE.h"
16			#include "DYNVARS.h"
17			#include "EEPARAMS.h"
18			#include "PARAMS.h"
19			#include "GRID.h"
20			C == Routine Arguments ==
21			INTEGER bi,bj,iMin,iMax,jMin,jMax
22			INTEGER K
23			INTEGER myThid
24			C == Local variables ==
25			INTEGER i,j
26			_RL Tfreezing
27
28			Tfreezing=-1.9
29
30			C Check for water that should have frozen
31			DO j=jMin,jMax
32			DO i=iMin,iMax
33			IF (gTNm1(i,j,k,bi,bj) .LT. Tfreezing) THEN
34			gTNm1(i,j,k,bi,bj)=Tfreezing
35			ENDIF
36			ENDDO
37			ENDDO
38
39			RETURN
40			END