1 |
|
C $Header$ |
2 |
|
C $Name$ |
3 |
|
|
4 |
#include "CPP_OPTIONS.h" |
#include "AUTODIFF_OPTIONS.h" |
5 |
|
|
6 |
CBOP |
CBOP |
7 |
C !ROUTINE: INI_AUTODIFF |
C !ROUTINE: INI_AUTODIFF |
21 |
#include "SIZE.h" |
#include "SIZE.h" |
22 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
23 |
#include "PARAMS.h" |
#include "PARAMS.h" |
24 |
|
#include "FFIELDS.h" |
25 |
#ifdef EXACT_CONSERV |
#ifdef EXACT_CONSERV |
26 |
# include "SURFACE.h" |
# include "SURFACE.h" |
27 |
#endif |
#endif |
|
#ifdef INCLUDE_EXTERNAL_FORCING_PACKAGE |
|
|
# include "FFIELDS.h" |
|
|
#endif |
|
28 |
|
|
29 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
30 |
C == Routine arguments == |
C == Routine arguments == |
38 |
C bi,bj - Loop counters |
C bi,bj - Loop counters |
39 |
C I,J,K |
C I,J,K |
40 |
INTEGER bi, bj |
INTEGER bi, bj |
41 |
INTEGER I, J, K |
INTEGER I, J, K |
42 |
CEOP |
CEOP |
43 |
|
|
44 |
C-- Over all tiles |
C-- Over all tiles |
46 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
47 |
|
|
48 |
C- 3D arrays |
C- 3D arrays |
49 |
|
DO K=1,Nr |
50 |
|
DO J=1-Oly,sNy+Oly |
51 |
|
DO I=1-Olx,sNx+Olx |
52 |
|
#ifdef ALLOW_EP_FLUX |
53 |
|
EfluxY(i,j,k,bi,bj)=0. |
54 |
|
EfluxP(i,j,k,bi,bj)=0. |
55 |
|
#endif |
56 |
|
ENDDO |
57 |
|
ENDDO |
58 |
|
ENDDO |
59 |
|
|
60 |
C- 2D arrays |
C- 2D arrays |
61 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
62 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
63 |
#ifdef EXACT_CONSERV |
#ifdef EXACT_CONSERV |
64 |
hDivFlow(I,J,bi,bj)=0. _d 0 |
dEtaHdt(i,j,bi,bj) = 0. _d 0 |
65 |
|
PmEpR(i,j,bi,bj) = 0. _d 0 |
66 |
|
cph hDivFlow(I,J,bi,bj)=0. _d 0 |
67 |
#endif |
#endif |
68 |
#ifdef INCLUDE_EXTERNAL_FORCING_PACKAGE |
#ifdef ALLOW_EXF |
69 |
EmPmR(I,J,bi,bj) = 0. _d 0 |
EmPmR(I,J,bi,bj) = 0. _d 0 |
70 |
#endif |
#endif |
71 |
ENDDO |
ENDDO |