dgoldberg/streamice/streamice_init_fixed.F

C $Header: /u/gcmpack/MITgcm/pkg/streamice/streamice_init_fixed.F,v 1.2 2006/08/07 21:13:52 dng Exp $
C $Name:  $

#include "STREAMICE_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP 0
C !ROUTINE: STREAMICE_INIT_FIXED

C !INTERFACE:
      SUBROUTINE STREAMICE_INIT_FIXED( myThid )

C     !DESCRIPTION:
C     Initialize STREAMICE variables that are kept fixed during the run.

C     !USES:
      IMPLICIT NONE
#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
#include "STREAMICE.h"
#include "STREAMICE_CG.h"
#include "STREAMICE_BDRY.h"
#include "GRID.h"

C     myThid ::  my Thread Id number
      INTEGER myThid
CEOP

C     !LOCAL VARIABLES:
C     === Local variables ===
      INTEGER bi, bj, i, j, Gi, Gj
      _RL x, y
      CHARACTER*(MAX_LEN_MBUF) msgBuf

#ifdef ALLOW_STREAMICE

#ifdef ALLOW_MNC

      IF ( useMNC .AND. (streamice_tave_mnc.OR.streamice_dump_mnc)
     &    ) THEN
C        CALL STREAMICE_MNC_INIT( myThid )
      ENDIF
#endif /* ALLOW_MNC */

#ifdef ALLOW_DIAGNOSTICS
      IF ( useDiagnostics ) THEN
        CALL STREAMICE_INIT_DIAGNOSTICS( myThid )
      ENDIF
#endif

      IF ( STREAMICE_calve_to_mask ) THEN
       IF ( STREAMICEcalveMaskFile .NE. ' ') THEN
        _BARRIER
C       The 0 is the "iteration" argument. The ' ' is an empty suffix
        CALL READ_FLD_XY_RL( STREAMICEcalveMaskFile, ' ', 
     &   STREAMICE_calve_mask, 0, myThid )
       ELSE
        WRITE(msgBuf,'(A)') 'INIT CALVE MASK - NOT IMPLENTED'
        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)
       ENDIF
      ENDIF

C- fill in the overlap (+ BARRIER):
      _EXCH_XY_RL(STREAMICE_calve_mask, myThid )

      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        DO j=1,sNy
         DO i=1,sNx

C         INIT VALUES FOR METRIC TERMS

          k1AtC_str(i,j,bi,bj) = recip_rA(i,j,bi,bj) * 
     &     (dxG(i+1,j,bi,bj)-dxG(i,j,bi,bj))
          k2AtC_str(i,j,bi,bj) = recip_rA(i,j,bi,bj) * 
     &     (dyG(i,j+1,bi,bj)-dyG(i,j,bi,bj))

C         INIT BDRY CONDITIONS

          Gi = (myXGlobalLo-1)+(bi-1)*sNx+i
          Gj = (myYGlobalLo-1)+(bj-1)*sNy+j
          x = xC(i,j,bi,bj)
          y = yC(i,j,bi,bj)

          STREAMICE_ufacemask_bdry (i,j,bi,bj) = -1.0
          STREAMICE_vfacemask_bdry (i,j,bi,bj) = -1.0

          IF (Gi .EQ. 1) THEN
C         INITIALIZE BOUNDARY CONDS AT WEST
           IF ((y .gt. min_y_noflow_WEST) .AND.
     &         (y .le. max_y_noflow_WEST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 0.0
           ENDIF
           IF ((y .gt. min_y_noStress_WEST) .AND.
     &         (y .le. max_y_noStress_WEST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 1.0
           ENDIF
           IF ((y .gt. min_y_FluxBdry_WEST) .AND.
     &         (y .le. max_y_FluxBdry_WEST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 4.0
            u_flux_bdry_SI (i,j,bi,bj) = flux_bdry_val_WEST
           ENDIF
           IF ((y .gt. min_y_Dirich_WEST) .AND.
     &         (y .le. max_y_Dirich_WEST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 3.0
            u_bdry_values_SI (i,j,bi,bj) = 0.0
           ENDIF
           IF ((y .gt. min_y_CFBC_WEST) .AND.
     &         (y .le. max_y_CFBC_WEST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 2.0
           ENDIF
          ENDIF

          IF (Gi .EQ. Nx) THEN
C         INITIALIZE BOUNDARY CONDS AT EAST
           IF ((y .gt. min_y_noflow_EAST) .AND.
     &         (y .le. max_y_noflow_EAST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 0.0
           ENDIF
           IF ((y .gt. min_y_noStress_EAST) .AND.
     &         (y .le. max_y_noStress_EAST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 1.0
           ENDIF
           IF ((y .gt. min_y_FluxBdry_EAST) .AND.
     &         (y .le. max_y_FluxBdry_EAST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 4.0 
            u_flux_bdry_SI (i,j,bi,bj) = flux_bdry_val_EAST
           ENDIF
           IF ((y .gt. min_y_Dirich_EAST) .AND.
     &         (y .le. max_y_Dirich_EAST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 3.0
            u_bdry_values_SI (i,j,bi,bj) = 0. _d 0
           ENDIF
           IF ((y .gt. min_y_CFBC_EAST) .AND.
     &         (y .le. max_y_CFBC_EAST)) THEN
            STREAMICE_ufacemask_bdry (i,j,bi,bj) = 2.0
           ENDIF
          ENDIF

          IF (Gj .EQ. 1) THEN
C         INITIALIZE BOUNDARY CONDS AT SOUTH
           IF ((x .gt. min_x_noflow_SOUTH) .AND.
     &         (x .le. max_x_noflow_SOUTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 0.0
           ENDIF
           IF ((x .gt. min_x_noStress_SOUTH) .AND.
     &         (x .le. max_x_noStress_SOUTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 1.0
           ENDIF
           IF ((x .gt. min_x_FluxBdry_SOUTH) .AND.
     &         (x .le. max_x_FluxBdry_SOUTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 4.0
            u_flux_bdry_SI (i,j,bi,bj) = flux_bdry_val_SOUTH
           ENDIF
           IF ((x .gt. min_x_Dirich_SOUTH) .AND.
     &         (x .le. max_x_Dirich_SOUTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 3.0
            u_bdry_values_SI (i,j,bi,bj) = 0. _d 0
           ENDIF
           IF ((x .gt. min_x_CFBC_SOUTH) .AND.
     &         (x .le. max_x_CFBC_SOUTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 2.0
           ENDIF
          ENDIF

          IF (Gj .EQ. Ny) THEN
C         INITIALIZE BOUNDARY CONDS AT NORTH
           IF ((x .gt. min_x_noflow_NORTH) .AND.
     &         (x .le. max_x_noflow_NORTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 0. _d 0
           ENDIF
           IF ((x .gt. min_x_noStress_NORTH) .AND.
     &         (x .le. max_x_noStress_NORTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 1.0
           ENDIF
           IF ((x .gt. min_x_FluxBdry_NORTH) .AND.
     &         (x .le. max_x_FluxBdry_NORTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 4.0
            v_flux_bdry_SI (i,j,bi,bj) = flux_bdry_val_NORTH
           ENDIF
           IF ((x .gt. min_x_Dirich_NORTH) .AND.
     &         (x .le. max_x_Dirich_NORTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 3.0
            v_bdry_values_SI (i,j,bi,bj) = 0. _d 0
           ENDIF
           IF ((x .gt. min_x_CFBC_NORTH) .AND.
     &         (x .le. max_x_CFBC_NORTH)) THEN
            STREAMICE_vfacemask_bdry (i,j,bi,bj) = 2.0
           ENDIF
          ENDIF
                   
         ENDDO
        ENDDO
       ENDDO
      ENDDO

      _EXCH_XY_RL(k1AtC_str, myThid )
      _EXCH_XY_RL(k2AtC_str, myThid )
      _EXCH_XY_RL(STREAMICE_ufacemask_bdry, myThid )
      _EXCH_XY_RL(STREAMICE_vfacemask_bdry, myThid )
      _EXCH_XY_RL(u_bdry_values_SI, myThid )
      _EXCH_XY_RL(v_bdry_values_SI, myThid )
      _EXCH_XY_RL(u_flux_bdry_SI, myThid )
      _EXCH_XY_RL(v_flux_bdry_SI, myThid )

      Xquad (1) = .5 * (1.-1./sqrt(3.))
      Xquad (2) = .5 * (1.+1./sqrt(3.))

      CALL STREAMICE_INIT_PHI( myThid )

      

#endif

      RETURN
      END

1	heimbach	1.1	C $Header: /u/gcmpack/MITgcm/pkg/streamice/streamice_init_fixed.F,v 1.2 2006/08/07 21:13:52 dng Exp $
2			C $Name: $
3
4			#include "STREAMICE_OPTIONS.h"
5
6			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7			CBOP 0
8			C !ROUTINE: STREAMICE_INIT_FIXED
9
10			C !INTERFACE:
11			SUBROUTINE STREAMICE_INIT_FIXED( myThid )
12
13			C !DESCRIPTION:
14			C Initialize STREAMICE variables that are kept fixed during the run.
15
16			C !USES:
17			IMPLICIT NONE
18			#include "EEPARAMS.h"
19			#include "SIZE.h"
20			#include "PARAMS.h"
21			#include "STREAMICE.h"
22			#include "STREAMICE_CG.h"
23			#include "STREAMICE_BDRY.h"
24			#include "GRID.h"
25
26			C myThid :: my Thread Id number
27			INTEGER myThid
28			CEOP
29
30			C !LOCAL VARIABLES:
31			C === Local variables ===
32			INTEGER bi, bj, i, j, Gi, Gj
33			_RL x, y
34			CHARACTER*(MAX_LEN_MBUF) msgBuf
35
36			#ifdef ALLOW_STREAMICE
37
38			#ifdef ALLOW_MNC
39
40			IF ( useMNC .AND. (streamice_tave_mnc.OR.streamice_dump_mnc)
41			& ) THEN
42			C CALL STREAMICE_MNC_INIT( myThid )
43			ENDIF
44			#endif /* ALLOW_MNC */
45
46			#ifdef ALLOW_DIAGNOSTICS
47			IF ( useDiagnostics ) THEN
48			CALL STREAMICE_INIT_DIAGNOSTICS( myThid )
49			ENDIF
50			#endif
51
52			IF ( STREAMICE_calve_to_mask ) THEN
53			IF ( STREAMICEcalveMaskFile .NE. ' ') THEN
54			_BARRIER
55			C The 0 is the "iteration" argument. The ' ' is an empty suffix
56			CALL READ_FLD_XY_RL( STREAMICEcalveMaskFile, ' ',
57			& STREAMICE_calve_mask, 0, myThid )
58			ELSE
59			WRITE(msgBuf,'(A)') 'INIT CALVE MASK - NOT IMPLENTED'
60			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
61			& SQUEEZE_RIGHT , 1)
62			ENDIF
63			ENDIF
64
65			C- fill in the overlap (+ BARRIER):
66			_EXCH_XY_RL(STREAMICE_calve_mask, myThid )
67
68			DO bj = myByLo(myThid), myByHi(myThid)
69			DO bi = myBxLo(myThid), myBxHi(myThid)
70			DO j=1,sNy
71			DO i=1,sNx
72
73			C INIT VALUES FOR METRIC TERMS
74
75			k1AtC_str(i,j,bi,bj) = recip_rA(i,j,bi,bj) *
76			& (dxG(i+1,j,bi,bj)-dxG(i,j,bi,bj))
77			k2AtC_str(i,j,bi,bj) = recip_rA(i,j,bi,bj) *
78			& (dyG(i,j+1,bi,bj)-dyG(i,j,bi,bj))
79
80			C INIT BDRY CONDITIONS
81
82			Gi = (myXGlobalLo-1)+(bi-1)*sNx+i
83			Gj = (myYGlobalLo-1)+(bj-1)*sNy+j
84			x = xC(i,j,bi,bj)
85			y = yC(i,j,bi,bj)
86
87			STREAMICE_ufacemask_bdry (i,j,bi,bj) = -1.0
88			STREAMICE_vfacemask_bdry (i,j,bi,bj) = -1.0
89
90			IF (Gi .EQ. 1) THEN
91			C INITIALIZE BOUNDARY CONDS AT WEST
92			IF ((y .gt. min_y_noflow_WEST) .AND.
93			& (y .le. max_y_noflow_WEST)) THEN
94			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 0.0
95			ENDIF
96			IF ((y .gt. min_y_noStress_WEST) .AND.
97			& (y .le. max_y_noStress_WEST)) THEN
98			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 1.0
99			ENDIF
100			IF ((y .gt. min_y_FluxBdry_WEST) .AND.
101			& (y .le. max_y_FluxBdry_WEST)) THEN
102			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 4.0
103			u_flux_bdry_SI (i,j,bi,bj) = flux_bdry_val_WEST
104			ENDIF
105			IF ((y .gt. min_y_Dirich_WEST) .AND.
106			& (y .le. max_y_Dirich_WEST)) THEN
107			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 3.0
108			u_bdry_values_SI (i,j,bi,bj) = 0.0
109			ENDIF
110			IF ((y .gt. min_y_CFBC_WEST) .AND.
111			& (y .le. max_y_CFBC_WEST)) THEN
112			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 2.0
113			ENDIF
114			ENDIF
115
116			IF (Gi .EQ. Nx) THEN
117			C INITIALIZE BOUNDARY CONDS AT EAST
118			IF ((y .gt. min_y_noflow_EAST) .AND.
119			& (y .le. max_y_noflow_EAST)) THEN
120			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 0.0
121			ENDIF
122			IF ((y .gt. min_y_noStress_EAST) .AND.
123			& (y .le. max_y_noStress_EAST)) THEN
124			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 1.0
125			ENDIF
126			IF ((y .gt. min_y_FluxBdry_EAST) .AND.
127			& (y .le. max_y_FluxBdry_EAST)) THEN
128			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 4.0
129			u_flux_bdry_SI (i,j,bi,bj) = flux_bdry_val_EAST
130			ENDIF
131			IF ((y .gt. min_y_Dirich_EAST) .AND.
132			& (y .le. max_y_Dirich_EAST)) THEN
133			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 3.0
134			u_bdry_values_SI (i,j,bi,bj) = 0. _d 0
135			ENDIF
136			IF ((y .gt. min_y_CFBC_EAST) .AND.
137			& (y .le. max_y_CFBC_EAST)) THEN
138			STREAMICE_ufacemask_bdry (i,j,bi,bj) = 2.0
139			ENDIF
140			ENDIF
141
142			IF (Gj .EQ. 1) THEN
143			C INITIALIZE BOUNDARY CONDS AT SOUTH
144			IF ((x .gt. min_x_noflow_SOUTH) .AND.
145			& (x .le. max_x_noflow_SOUTH)) THEN
146			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 0.0
147			ENDIF
148			IF ((x .gt. min_x_noStress_SOUTH) .AND.
149			& (x .le. max_x_noStress_SOUTH)) THEN
150			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 1.0
151			ENDIF
152			IF ((x .gt. min_x_FluxBdry_SOUTH) .AND.
153			& (x .le. max_x_FluxBdry_SOUTH)) THEN
154			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 4.0
155			u_flux_bdry_SI (i,j,bi,bj) = flux_bdry_val_SOUTH
156			ENDIF
157			IF ((x .gt. min_x_Dirich_SOUTH) .AND.
158			& (x .le. max_x_Dirich_SOUTH)) THEN
159			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 3.0
160			u_bdry_values_SI (i,j,bi,bj) = 0. _d 0
161			ENDIF
162			IF ((x .gt. min_x_CFBC_SOUTH) .AND.
163			& (x .le. max_x_CFBC_SOUTH)) THEN
164			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 2.0
165			ENDIF
166			ENDIF
167
168			IF (Gj .EQ. Ny) THEN
169			C INITIALIZE BOUNDARY CONDS AT NORTH
170			IF ((x .gt. min_x_noflow_NORTH) .AND.
171			& (x .le. max_x_noflow_NORTH)) THEN
172			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 0. _d 0
173			ENDIF
174			IF ((x .gt. min_x_noStress_NORTH) .AND.
175			& (x .le. max_x_noStress_NORTH)) THEN
176			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 1.0
177			ENDIF
178			IF ((x .gt. min_x_FluxBdry_NORTH) .AND.
179			& (x .le. max_x_FluxBdry_NORTH)) THEN
180			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 4.0
181			v_flux_bdry_SI (i,j,bi,bj) = flux_bdry_val_NORTH
182			ENDIF
183			IF ((x .gt. min_x_Dirich_NORTH) .AND.
184			& (x .le. max_x_Dirich_NORTH)) THEN
185			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 3.0
186			v_bdry_values_SI (i,j,bi,bj) = 0. _d 0
187			ENDIF
188			IF ((x .gt. min_x_CFBC_NORTH) .AND.
189			& (x .le. max_x_CFBC_NORTH)) THEN
190			STREAMICE_vfacemask_bdry (i,j,bi,bj) = 2.0
191			ENDIF
192			ENDIF
193
194			ENDDO
195			ENDDO
196			ENDDO
197			ENDDO
198
199			_EXCH_XY_RL(k1AtC_str, myThid )
200			_EXCH_XY_RL(k2AtC_str, myThid )
201			_EXCH_XY_RL(STREAMICE_ufacemask_bdry, myThid )
202			_EXCH_XY_RL(STREAMICE_vfacemask_bdry, myThid )
203			_EXCH_XY_RL(u_bdry_values_SI, myThid )
204			_EXCH_XY_RL(v_bdry_values_SI, myThid )
205			_EXCH_XY_RL(u_flux_bdry_SI, myThid )
206			_EXCH_XY_RL(v_flux_bdry_SI, myThid )
207
208			Xquad (1) = .5 * (1.-1./sqrt(3.))
209			Xquad (2) = .5 * (1.+1./sqrt(3.))
210
211			CALL STREAMICE_INIT_PHI( myThid )
212
213
214
215			#endif
216
217			RETURN
218			END
219