/[MITgcm]/MITgcm/pkg/ebm/ebm_ini_vars.F
ViewVC logotype

Contents of /MITgcm/pkg/ebm/ebm_ini_vars.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.3 - (show annotations) (download)
Tue Apr 28 18:11:51 2009 UTC (15 years, 1 month ago) by jmc
Branch: MAIN
CVS Tags: checkpoint63a, checkpoint62, checkpoint63, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint61n, checkpoint61o, checkpoint61m, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.2: +3 -3 lines
change macros (EXCH & GLOBAL_SUM/MAX) sufix _R4/_R8 to _RS/_RL
 when applied to _RS/_RL variable

1 C $Header: /u/gcmpack/MITgcm/pkg/ebm/ebm_ini_vars.F,v 1.2 2004/05/19 20:38:53 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_RS(Run, myThid )
50 _EXCH_XY_RS(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

  ViewVC Help
Powered by ViewVC 1.1.22