shelfice_remeshing/code/shelfice_remeshing.F

C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/remeshing/code/shelfice_remeshing.F,v 1.2 2017/03/06 10:29:19 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_REMESHING
C     | o If SHELFICERemeshFrequencey is a multiple of mytime,
C     |   checks if etaN under shelfice triggers a remesh
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
#ifdef ALLOW_SHELFICE_REMESHING
          IF ((myTime/ShelficeRemeshFrequency) .EQ.
     &     INT (myTime/ShelficeRemeshFrequency) .AND. myTime
     &      .NE. 0.0        ) THEN
         CALL INI_MASKS_REMESH (myThid)
         CALL SHELFICE_INIT_FIXED_REMESH( myThid )
#ifdef ALLOW_OBCS
         if (useobcs) then
         CALL OBCS_INIT_FIXED( myThid )
         endif
#endif
         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, myTime, niter0, myThid)
         ENDIF 
#endif /* ALLOW_SHELFICE_REMESHING */
#endif /* ALLOW_SHELFICE */

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/remeshing/code/shelfice_remeshing.F,v 1.2 2017/03/06 10:29:19 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_REMESHING
16	C \| o If SHELFICERemeshFrequencey is a multiple of mytime,
17	C \| checks if etaN under shelfice triggers a remesh
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
40
41
42
43
44
45	#ifdef ALLOW_SHELFICE
46	#ifdef ALLOW_SHELFICE_REMESHING
47	IF ((myTime/ShelficeRemeshFrequency) .EQ.
48	& INT (myTime/ShelficeRemeshFrequency) .AND. myTime
49	& .NE. 0.0 ) THEN
50	CALL INI_MASKS_REMESH (myThid)
51	CALL SHELFICE_INIT_FIXED_REMESH( myThid )
52	#ifdef ALLOW_OBCS
53	if (useobcs) then
54	CALL OBCS_INIT_FIXED( myThid )
55	endif
56	#endif
57	CALL INI_LINEAR_PHISURF(myThid)
58	CALL INI_CG2D (myThid)
59	CALL CALC_SURF_DR( etaH, myTime, myIter, myThid )
60	CALL CALC_SURF_DR( etaH, myTime, myIter, myThid )
61	CALL UPDATE_SURF_DR(.FALSE., myTime, myIter, myThid)
62	CALL UPDATE_CG2D(myTime,myIter, myThid)
63	CALL INTEGR_CONTINUITY (uVel ,vVel, myTime, niter0, myThid)
64	ENDIF
65	#endif /* ALLOW_SHELFICE_REMESHING */
66	#endif /* ALLOW_SHELFICE */
67
68	RETURN
69	END