pkg/ebm/ebm_ini_vars.F

C $Header: /u/gcmpack/MITgcm/pkg/ebm/ebm_ini_vars.F,v 1.2 2004/05/19 20:38:53 heimbach Exp $
C $Name:  $

#include "EBM_OPTIONS.h"

      SUBROUTINE EBM_INI_VARS( myThid )
C     |==========================================================|
C     | S/R INITIALISE VARIABLES FOR EBM                         |
C     |==========================================================|

      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "GRID.h"
#include "FFIELDS.h"
#ifdef ALLOW_EBM
# include "EBM.h"
#endif

C     === Routine arguments ===
C     myThid - Instance number for this innvocation 
      INTEGER myThid

CEndOfInterface
C     == Local variables ==

#ifdef ALLOW_EBM

C     I, J, K - Loop counters
      INTEGER bi,bj,i,j,k
      INTEGER iMin, iMax, jMin, jMax

C--   initially all forcing fields are zero
      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        DO j=1-OLy,sNy+OLy
         DO i=1-OLx,sNx+OLx
          Run   (i,j,bi,bj) = 0.
          winPert(i,j,bi,bj) = 0.
         ENDDO
        ENDDO
       ENDDO
      ENDDO

      _EXCH_XY_RS(Run,   myThid )
      _EXCH_XY_RS(winPert,  myThid )

      CALL EBM_LOAD_CLIMATOLOGY( myThid )

      CALL EBM_DRIVER( startTime, niter0, myThid )

      print *, TmlS - 273.15, TtS - 273.15
      print *, TmlN - 273.15, TtN - 273.15

      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        iMin =   1-olx
        iMax = sNx+olx
        jMin =   1-oly
        jMax = sNy+oly

       CALL EBM_FORCING_SURF(
     I       bi, bj, iMin, iMax, jMin, jMax,
     I       startTime, niter0, myThid )
        ENDDO
       ENDDO

#endif /* ALLOW_EBM */

      RETURN
      END


1	jmc	1.3	C $Header: /u/gcmpack/MITgcm/pkg/ebm/ebm_ini_vars.F,v 1.2 2004/05/19 20:38:53 heimbach Exp $
2	heimbach	1.1	C $Name: $
3
4			#include "EBM_OPTIONS.h"
5
6			SUBROUTINE EBM_INI_VARS( myThid )
7			C \|==========================================================\|
8			C \| S/R INITIALISE VARIABLES FOR EBM \|
9			C \|==========================================================\|
10
11			IMPLICIT NONE
12
13			C === Global variables ===
14			#include "SIZE.h"
15			#include "EEPARAMS.h"
16			#include "PARAMS.h"
17			#include "DYNVARS.h"
18			#include "GRID.h"
19			#include "FFIELDS.h"
20			#ifdef ALLOW_EBM
21			# include "EBM.h"
22			#endif
23
24			C === Routine arguments ===
25			C myThid - Instance number for this innvocation
26			INTEGER myThid
27
28			CEndOfInterface
29			C == Local variables ==
30
31			#ifdef ALLOW_EBM
32
33			C I, J, K - Loop counters
34			INTEGER bi,bj,i,j,k
35			INTEGER iMin, iMax, jMin, jMax
36
37			C-- initially all forcing fields are zero
38			DO bj = myByLo(myThid), myByHi(myThid)
39			DO bi = myBxLo(myThid), myBxHi(myThid)
40			DO j=1-OLy,sNy+OLy
41			DO i=1-OLx,sNx+OLx
42			Run (i,j,bi,bj) = 0.
43			winPert(i,j,bi,bj) = 0.
44			ENDDO
45			ENDDO
46			ENDDO
47			ENDDO
48
49	jmc	1.3	_EXCH_XY_RS(Run, myThid )
50			_EXCH_XY_RS(winPert, myThid )
51	heimbach	1.1
52			CALL EBM_LOAD_CLIMATOLOGY( myThid )
53
54			CALL EBM_DRIVER( startTime, niter0, myThid )
55
56			print *, TmlS - 273.15, TtS - 273.15
57			print *, TmlN - 273.15, TtN - 273.15
58
59			DO bj = myByLo(myThid), myByHi(myThid)
60			DO bi = myBxLo(myThid), myBxHi(myThid)
61	heimbach	1.2	iMin = 1-olx
62			iMax = sNx+olx
63			jMin = 1-oly
64			jMax = sNy+oly
65	heimbach	1.1
66			CALL EBM_FORCING_SURF(
67			I bi, bj, iMin, iMax, jMin, jMax,
68			I startTime, niter0, myThid )
69			ENDDO
70			ENDDO
71
72			#endif /* ALLOW_EBM */
73
74			RETURN
75			END
76
77
78