pkg/seaice/seaice_init_fixed.F

C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_init.F,v 1.37 2006/10/05 18:41:32 jmc Exp $
C $Name:  $

#include "SEAICE_OPTIONS.h"

CStartOfInterface
      SUBROUTINE SEAICE_INIT_FIXED( myThid )
C     /==========================================================\
C     | SUBROUTINE SEAICE_INIT_FIXED                             |
C     | o Initialization of sea ice model.                       |
C     |==========================================================|
C     \==========================================================/
      IMPLICIT NONE
 
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "SEAICE.h"
CML#include "SEAICE_GRID.h"
#include "SEAICE_DIAGS.h"
#include "SEAICE_PARAMS.h"
#include "FFIELDS.h"

C     === Routine arguments ===
C     myThid - Thread no. that called this routine.
      INTEGER myThid
CEndOfInterface
 
C     === Local variables ===
C     i,j,k,bi,bj - Loop counters

      INTEGER i, j, k, bi, bj
      INTEGER myIter, myTile

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

#ifdef ALLOW_TIMEAVE
C     Initialize averages to zero
      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        CALL TIMEAVE_RESET(FUtave   ,1,bi,bj,myThid)
        CALL TIMEAVE_RESET(FVtave   ,1,bi,bj,myThid)
        CALL TIMEAVE_RESET(EmPmRtave,1,bi,bj,myThid)
        CALL TIMEAVE_RESET(QNETtave ,1,bi,bj,myThid)
        CALL TIMEAVE_RESET(QSWtave  ,1,bi,bj,myThid)
        CALL TIMEAVE_RESET(UICEtave ,1,bi,bj,myThid)
        CALL TIMEAVE_RESET(VICEtave ,1,bi,bj,myThid)
        CALL TIMEAVE_RESET(HEFFtave ,1,bi,bj,myThid)
        CALL TIMEAVE_RESET(AREAtave ,1,bi,bj,myThid)
        DO k=1,Nr
         SEAICE_TimeAve(k,bi,bj)=ZERO
        ENDDO
       ENDDO
      ENDDO
#endif /* ALLOW_TIMEAVE */

      RETURN
      END
1	heimbach	1.1	C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_init.F,v 1.37 2006/10/05 18:41:32 jmc Exp $
2			C $Name: $
3
4			#include "SEAICE_OPTIONS.h"
5
6			CStartOfInterface
7			SUBROUTINE SEAICE_INIT_FIXED( myThid )
8			C /==========================================================\
9			C \| SUBROUTINE SEAICE_INIT_FIXED \|
10			C \| o Initialization of sea ice model. \|
11			C \|==========================================================\|
12			C \==========================================================/
13			IMPLICIT NONE
14
15			C === Global variables ===
16			#include "SIZE.h"
17			#include "EEPARAMS.h"
18			#include "PARAMS.h"
19			#include "GRID.h"
20			#include "SEAICE.h"
21			CML#include "SEAICE_GRID.h"
22			#include "SEAICE_DIAGS.h"
23			#include "SEAICE_PARAMS.h"
24			#include "FFIELDS.h"
25
26			C === Routine arguments ===
27			C myThid - Thread no. that called this routine.
28			INTEGER myThid
29			CEndOfInterface
30
31			C === Local variables ===
32			C i,j,k,bi,bj - Loop counters
33
34			INTEGER i, j, k, bi, bj
35			INTEGER myIter, myTile
36
37			#ifdef ALLOW_DIAGNOSTICS
38			IF ( useDiagnostics ) THEN
39			CALL SEAICE_DIAGNOSTICS_INIT( myThid )
40			ENDIF
41			#endif
42
43			#ifdef ALLOW_TIMEAVE
44			C Initialize averages to zero
45			DO bj = myByLo(myThid), myByHi(myThid)
46			DO bi = myBxLo(myThid), myBxHi(myThid)
47			CALL TIMEAVE_RESET(FUtave ,1,bi,bj,myThid)
48			CALL TIMEAVE_RESET(FVtave ,1,bi,bj,myThid)
49			CALL TIMEAVE_RESET(EmPmRtave,1,bi,bj,myThid)
50			CALL TIMEAVE_RESET(QNETtave ,1,bi,bj,myThid)
51			CALL TIMEAVE_RESET(QSWtave ,1,bi,bj,myThid)
52			CALL TIMEAVE_RESET(UICEtave ,1,bi,bj,myThid)
53			CALL TIMEAVE_RESET(VICEtave ,1,bi,bj,myThid)
54			CALL TIMEAVE_RESET(HEFFtave ,1,bi,bj,myThid)
55			CALL TIMEAVE_RESET(AREAtave ,1,bi,bj,myThid)
56			DO k=1,Nr
57			SEAICE_TimeAve(k,bi,bj)=ZERO
58			ENDDO
59			ENDDO
60			ENDDO
61			#endif /* ALLOW_TIMEAVE */
62
63			RETURN
64			END