MANUAL/code/shelfice_remeshing.F

C $Header: /u/gcmpack/MITgcm_contrib/shelfice_remeshing/code/shelfice_remeshing.F,v 1.2 2015/08/07 10:35:32 dgoldberg Exp $
C $Name:  $
#include "SHELFICE_OPTIONS.h"
#ifdef ALLOW_CTRL
# include "CTRL_OPTIONS.h"
#endif

CBOP
C     !ROUTINE: SHELFICE_REMESHING
C     !INTERFACE:
      SUBROUTINE SHELFICE_REMESHING(
     I    myTime, myIter, myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE SHELFICE_UPDATE_MASKS
C     | o modify topography factor hFacC according to ice shelf
C     |   topography
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "SURFACE.h"
#include "GRID.h"
#include "CG2D.h" 
#ifdef ALLOW_SHELFICE
# include "SHELFICE.h"
#endif /* ALLOW_SHELFICE */

      INTEGER myThid
      INTEGER myIter
      _RL myTime

#ifdef ALLOW_SHELFICE
          IF ((myTime/ShelficeRemeshFrequency) .EQ.
     &     INT (myTime/ShelficeRemeshFrequency) .AND. myTime
     &      .NE. 0.0        ) THEN


       CALL INI_MASKS_ETC_JJ (myThid)

       CALL   OBCS_INIT_FIXED( myThid )

       CALL INI_LINEAR_PHISURF(myThid)

       CALL INI_CG2D (myThid) 

       CALL CALC_SURF_DR( etaH, myTime, myIter, myThid )
       CALL CALC_SURF_DR( etaH, myTime, myIter, myThid )

       CALL UPDATE_SURF_DR(.FALSE., myTime, myIter, myThid)
    
       CALL UPDATE_CG2D(myTime,myIter, myThid)

       CALL INTEGR_CONTINUITY (uVel ,vVel , starttime, myIter, myThid)
   

         ENDIF 
#endif /* ALLOW_SHELFICE */

      RETURN
      END
1	dgoldberg	1.1	C $Header: /u/gcmpack/MITgcm_contrib/shelfice_remeshing/code/shelfice_remeshing.F,v 1.2 2015/08/07 10:35:32 dgoldberg Exp $
2			C $Name: $
3			#include "SHELFICE_OPTIONS.h"
4			#ifdef ALLOW_CTRL
5			# include "CTRL_OPTIONS.h"
6			#endif
7
8			CBOP
9			C !ROUTINE: SHELFICE_REMESHING
10			C !INTERFACE:
11			SUBROUTINE SHELFICE_REMESHING(
12			I myTime, myIter, myThid )
13			C !DESCRIPTION: \bv
14			C ==========================================================
15			C \| SUBROUTINE SHELFICE_UPDATE_MASKS
16			C \| o modify topography factor hFacC according to ice shelf
17			C \| topography
18			C ==========================================================
19			C \ev
20
21			C !USES:
22			IMPLICIT NONE
23			C === Global variables ===
24			#include "SIZE.h"
25			#include "EEPARAMS.h"
26			#include "PARAMS.h"
27			#include "DYNVARS.h"
28			#include "SURFACE.h"
29			#include "GRID.h"
30			#include "CG2D.h"
31			#ifdef ALLOW_SHELFICE
32			# include "SHELFICE.h"
33			#endif /* ALLOW_SHELFICE */
34
35			INTEGER myThid
36			INTEGER myIter
37			_RL myTime
38
39			#ifdef ALLOW_SHELFICE
40			IF ((myTime/ShelficeRemeshFrequency) .EQ.
41			& INT (myTime/ShelficeRemeshFrequency) .AND. myTime
42			& .NE. 0.0 ) THEN
43
44
45			CALL INI_MASKS_ETC_JJ (myThid)
46
47			CALL OBCS_INIT_FIXED( myThid )
48
49			CALL INI_LINEAR_PHISURF(myThid)
50
51			CALL INI_CG2D (myThid)
52
53			CALL CALC_SURF_DR( etaH, myTime, myIter, myThid )
54			CALL CALC_SURF_DR( etaH, myTime, myIter, myThid )
55
56			CALL UPDATE_SURF_DR(.FALSE., myTime, myIter, myThid)
57
58			CALL UPDATE_CG2D(myTime,myIter, myThid)
59
60			CALL INTEGR_CONTINUITY (uVel ,vVel , starttime, myIter, myThid)
61
62
63
64
65			ENDIF
66			#endif /* ALLOW_SHELFICE */
67
68			RETURN
69			END