pkg/generic_advdiff/gad_init_varia.F

C $Header: /u/gcmpack/MITgcm/pkg/generic_advdiff/gad_init_varia.F,v 1.2 2008/05/09 21:43:16 jmc Exp $
C $Name:  $

#include "GAD_OPTIONS.h"

CBOP
C     !ROUTINE: GAD_INIT_VARIA
C     !INTERFACE:
      SUBROUTINE GAD_INIT_VARIA( myThid )

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE GAD_INIT_VARIA
C     | o Routine to initialize GAD variables
C     *==========================================================*
C     \ev
C     !USES:
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
c #include "GRID.h"
#include "GAD.h"
#include "GAD_SOM_VARS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     myThid ::  my Thread Id number
      INTEGER myThid
CEOP

#ifdef GAD_ALLOW_TS_SOM_ADV
C     !LOCAL VARIABLES:
C     === Local variables ===
      INTEGER i,j,k,bi,bj
      INTEGER n

      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)

C--   Initialize 1rst Order (3:x,y,z) and 2nd Order (6:xx,yy,zz,xy,xz,yz) moments
        DO n=1,nSOM
         DO k=1,Nr
          DO j=1-Oly,sNy+OLy
           DO i=1-Olx,sNx+Olx
C       pot.temp field:
            som_T(i,j,k,bi,bj,n) = 0. _d 0
C       salinity field:
            som_S(i,j,k,bi,bj,n) = 0. _d 0
           ENDDO
          ENDDO
         ENDDO
        ENDDO

C- end bi,bj loops
       ENDDO
      ENDDO

      IF ( tempSOM_Advection .OR. saltSOM_Advection ) THEN
       IF ( startTime.EQ.baseTime .AND. nIter0.EQ.0
     &                      .AND. pickupSuff.EQ.' ' ) THEN
C-     Need to synchronize here before doing master-thread IO
c        _BARRIER
C-     Read-in initial conditions:
C        --- not yet coded ---
c        CALL GAD_SOM_EXCHANGES( myThid )
       ELSE
C-     Read 2nd-O moments state from pickup files
         CALL GAD_READ_PICKUP( nIter0, myThid )
         CALL GAD_SOM_EXCHANGES( myThid )
       ENDIF
      ENDIF

#endif /* GAD_ALLOW_TS_SOM_ADV */

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/generic_advdiff/gad_init_varia.F,v 1.2 2008/05/09 21:43:16 jmc Exp $
2	C $Name: $
3
4	#include "GAD_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: GAD_INIT_VARIA
8	C !INTERFACE:
9	SUBROUTINE GAD_INIT_VARIA( myThid )
10
11	C !DESCRIPTION: \bv
12	C ==========================================================
13	C \| SUBROUTINE GAD_INIT_VARIA
14	C \| o Routine to initialize GAD variables
15	C ==========================================================
16	C \ev
17	C !USES:
18	IMPLICIT NONE
19
20	C === Global variables ===
21	#include "SIZE.h"
22	#include "EEPARAMS.h"
23	#include "PARAMS.h"
24	c #include "GRID.h"
25	#include "GAD.h"
26	#include "GAD_SOM_VARS.h"
27
28	C !INPUT/OUTPUT PARAMETERS:
29	C === Routine arguments ===
30	C myThid :: my Thread Id number
31	INTEGER myThid
32	CEOP
33
34	#ifdef GAD_ALLOW_TS_SOM_ADV
35	C !LOCAL VARIABLES:
36	C === Local variables ===
37	INTEGER i,j,k,bi,bj
38	INTEGER n
39
40	DO bj = myByLo(myThid), myByHi(myThid)
41	DO bi = myBxLo(myThid), myBxHi(myThid)
42
43	C-- Initialize 1rst Order (3:x,y,z) and 2nd Order (6:xx,yy,zz,xy,xz,yz) moments
44	DO n=1,nSOM
45	DO k=1,Nr
46	DO j=1-Oly,sNy+OLy
47	DO i=1-Olx,sNx+Olx
48	C pot.temp field:
49	som_T(i,j,k,bi,bj,n) = 0. _d 0
50	C salinity field:
51	som_S(i,j,k,bi,bj,n) = 0. _d 0
52	ENDDO
53	ENDDO
54	ENDDO
55	ENDDO
56
57	C- end bi,bj loops
58	ENDDO
59	ENDDO
60
61	IF ( tempSOM_Advection .OR. saltSOM_Advection ) THEN
62	IF ( startTime.EQ.baseTime .AND. nIter0.EQ.0
63	& .AND. pickupSuff.EQ.' ' ) THEN
64	C- Need to synchronize here before doing master-thread IO
65	c _BARRIER
66	C- Read-in initial conditions:
67	C --- not yet coded ---
68	c CALL GAD_SOM_EXCHANGES( myThid )
69	ELSE
70	C- Read 2nd-O moments state from pickup files
71	CALL GAD_READ_PICKUP( nIter0, myThid )
72	CALL GAD_SOM_EXCHANGES( myThid )
73	ENDIF
74	ENDIF
75
76	#endif /* GAD_ALLOW_TS_SOM_ADV */
77
78	RETURN
79	END