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