model/src/freeze_surface.F

C $Header: /u/gcmpack/MITgcm/model/src/freeze_surface.F,v 1.2 2003/11/17 23:43:06 heimbach Exp $
C $Name:  $

#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: FREEZE_SURFACE
C     !INTERFACE:
      SUBROUTINE FREEZE_SURFACE( myTime, myIter, myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | S/R FREEZE_SURFACE                                                
C     | o Check water temperature and limit range of temperature  
C     | appropriately.                                            
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     == Global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
c #include "GRID.h"
#include "FFIELDS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine Arguments ==
C     myTime - Current time in simulation
C     myIter - Current iteration number in simulation
C     myThid :: Thread no. that called this routine.
      _RL myTime
      INTEGER myIter
      INTEGER myThid

C     !LOCAL VARIABLES:
C     == Local variables ==
C     Tfreezing :: Freezing threshold temperature.
      INTEGER bi,bj,i,j,k
      _RL Tfreezing
CEOP

      IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN
        k = Nr
      ELSE
        k= 1
      ENDIF

      Tfreezing = -1.9 _d 0

C     Check for water that should have frozen
      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        DO j=1-OLy,sNy+OLy
         DO i=1-OLx,sNx+OLx
          IF (theta(I,J,k,bi,bj) .LT. Tfreezing) THEN
             surfaceTendencyTice(I,J,bi,bj) =
     &            ( Tfreezing - theta(I,J,k,bi,bj) ) / deltaTtracer
             theta(I,J,k,bi,bj) = Tfreezing
          ELSE
             surfacetendencyTice(I,J,bi,bj) = 0. _d 0
          ENDIF
         ENDDO
        ENDDO
       ENDDO
      ENDDO

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/model/src/freeze_surface.F,v 1.2 2003/11/17 23:43:06 heimbach Exp $
2	C $Name: $
3
4	#include "PACKAGES_CONFIG.h"
5	#include "CPP_OPTIONS.h"
6
7	CBOP
8	C !ROUTINE: FREEZE_SURFACE
9	C !INTERFACE:
10	SUBROUTINE FREEZE_SURFACE( myTime, myIter, myThid )
11	C !DESCRIPTION: \bv
12	C ==========================================================
13	C \| S/R FREEZE_SURFACE
14	C \| o Check water temperature and limit range of temperature
15	C \| appropriately.
16	C ==========================================================
17	C \ev
18
19	C !USES:
20	IMPLICIT NONE
21	C == Global variables ==
22	#include "SIZE.h"
23	#include "EEPARAMS.h"
24	#include "PARAMS.h"
25	#include "DYNVARS.h"
26	c #include "GRID.h"
27	#include "FFIELDS.h"
28
29	C !INPUT/OUTPUT PARAMETERS:
30	C == Routine Arguments ==
31	C myTime - Current time in simulation
32	C myIter - Current iteration number in simulation
33	C myThid :: Thread no. that called this routine.
34	_RL myTime
35	INTEGER myIter
36	INTEGER myThid
37
38	C !LOCAL VARIABLES:
39	C == Local variables ==
40	C Tfreezing :: Freezing threshold temperature.
41	INTEGER bi,bj,i,j,k
42	_RL Tfreezing
43	CEOP
44
45	IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN
46	k = Nr
47	ELSE
48	k= 1
49	ENDIF
50
51	Tfreezing = -1.9 _d 0
52
53	C Check for water that should have frozen
54	DO bj=myByLo(myThid),myByHi(myThid)
55	DO bi=myBxLo(myThid),myBxHi(myThid)
56	DO j=1-OLy,sNy+OLy
57	DO i=1-OLx,sNx+OLx
58	IF (theta(I,J,k,bi,bj) .LT. Tfreezing) THEN
59	surfaceTendencyTice(I,J,bi,bj) =
60	& ( Tfreezing - theta(I,J,k,bi,bj) ) / deltaTtracer
61	theta(I,J,k,bi,bj) = Tfreezing
62	ELSE
63	surfacetendencyTice(I,J,bi,bj) = 0. _d 0
64	ENDIF
65	ENDDO
66	ENDDO
67	ENDDO
68	ENDDO
69
70	RETURN
71	END