pkg/ebm/ebm_ini_vars.F

C $Header: /u/gcmpack/MITgcm/pkg/ebm/ebm_ini_vars.F,v 1.1 2004/05/14 21:10:34 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_R4(Run,   myThid )
      _EXCH_XY_R4(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	C $Header: /u/gcmpack/MITgcm/pkg/ebm/ebm_ini_vars.F,v 1.1 2004/05/14 21:10:34 heimbach Exp $
2	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	_EXCH_XY_R4(Run, myThid )
50	_EXCH_XY_R4(winPert, myThid )
51
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	iMin = 1-olx
62	iMax = sNx+olx
63	jMin = 1-oly
64	jMax = sNy+oly
65
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