11 |
SUBROUTINE MOM_INIT_FIXED( myThid ) |
SUBROUTINE MOM_INIT_FIXED( myThid ) |
12 |
|
|
13 |
C !DESCRIPTION: |
C !DESCRIPTION: |
14 |
C Initialize fixed quantities |
C Initialize fixed quantities |
15 |
C for momentum (common to fluxform & vecinv) packages |
C for momentum (common to fluxform & vecinv) packages |
16 |
|
|
17 |
C !USES: |
C !USES: |
19 |
#include "SIZE.h" |
#include "SIZE.h" |
20 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
21 |
#include "PARAMS.h" |
#include "PARAMS.h" |
22 |
|
#include "GRID.h" |
23 |
|
#include "SURFACE.h" |
24 |
|
|
25 |
C !INPUT PARAMETERS: |
C !INPUT PARAMETERS: |
26 |
C myThid :: thread number |
C myThid :: thread number |
27 |
INTEGER myThid |
INTEGER myThid |
28 |
CEOP |
CEOP |
29 |
|
|
30 |
|
C !LOCAL VARIABLES: |
31 |
|
#ifdef NONLIN_FRSURF |
32 |
|
C i,j,k,bi,bj :: loop counter |
33 |
|
INTEGER i,j,k,bi,bj |
34 |
|
#endif |
35 |
|
|
36 |
|
#ifdef NONLIN_FRSURF |
37 |
|
C- used for side-drag: |
38 |
|
DO bj=myByLo(myThid), myByHi(myThid) |
39 |
|
DO bi=myBxLo(myThid), myBxHi(myThid) |
40 |
|
DO k=1,Nr |
41 |
|
DO j=1-Oly,sNy+Oly |
42 |
|
DO i=1-Olx,sNx+Olx |
43 |
|
h0FacC(i,j,k,bi,bj) = hFacC(i,j,k,bi,bj) |
44 |
|
h0FacW(i,j,k,bi,bj) = hFacW(i,j,k,bi,bj) |
45 |
|
h0FacS(i,j,k,bi,bj) = hFacS(i,j,k,bi,bj) |
46 |
|
ENDDO |
47 |
|
ENDDO |
48 |
|
ENDDO |
49 |
|
ENDDO |
50 |
|
ENDDO |
51 |
|
#endif /* NONLIN_FRSURF */ |
52 |
|
|
53 |
|
|
54 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
55 |
IF ( useDiagnostics ) THEN |
IF ( useDiagnostics ) THEN |
56 |
CALL MOM_DIAGNOSTICS_INIT( myThid ) |
CALL MOM_DIAGNOSTICS_INIT( myThid ) |