pkg/seaice/seaice_init_fixed.F

C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_init_fixed.F,v 1.1 2006/10/22 01:11:02 heimbach 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 k, bi, bj

#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	jmc	1.2	C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_init_fixed.F,v 1.1 2006/10/22 01:11:02 heimbach Exp $
2	heimbach	1.1	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	jmc	1.2	INTEGER k, bi, bj
35	heimbach	1.1
36			#ifdef ALLOW_DIAGNOSTICS
37			IF ( useDiagnostics ) THEN
38			CALL SEAICE_DIAGNOSTICS_INIT( myThid )
39			ENDIF
40			#endif
41
42			#ifdef ALLOW_TIMEAVE
43			C Initialize averages to zero
44			DO bj = myByLo(myThid), myByHi(myThid)
45			DO bi = myBxLo(myThid), myBxHi(myThid)
46			CALL TIMEAVE_RESET(FUtave ,1,bi,bj,myThid)
47			CALL TIMEAVE_RESET(FVtave ,1,bi,bj,myThid)
48			CALL TIMEAVE_RESET(EmPmRtave,1,bi,bj,myThid)
49			CALL TIMEAVE_RESET(QNETtave ,1,bi,bj,myThid)
50			CALL TIMEAVE_RESET(QSWtave ,1,bi,bj,myThid)
51			CALL TIMEAVE_RESET(UICEtave ,1,bi,bj,myThid)
52			CALL TIMEAVE_RESET(VICEtave ,1,bi,bj,myThid)
53			CALL TIMEAVE_RESET(HEFFtave ,1,bi,bj,myThid)
54			CALL TIMEAVE_RESET(AREAtave ,1,bi,bj,myThid)
55			DO k=1,Nr
56			SEAICE_TimeAve(k,bi,bj)=ZERO
57			ENDDO
58			ENDDO
59			ENDDO
60			#endif /* ALLOW_TIMEAVE */
61
62			RETURN
63			END