25 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
26 |
#include "PARAMS.h" |
#include "PARAMS.h" |
27 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
28 |
|
#ifdef ALLOW_SOLVE4_PS_AND_DRAG |
29 |
|
#include "GRID.h" |
30 |
|
#endif /* ALLOW_SOLVE4_PS_AND_DRAG */ |
31 |
|
|
32 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
33 |
C == Routine arguments == |
C == Routine arguments == |
81 |
totPhiHyd (i,j,k,bi,bj) = 0. _d 0 |
totPhiHyd (i,j,k,bi,bj) = 0. _d 0 |
82 |
rhoInSitu (i,j,k,bi,bj) = 0. _d 0 |
rhoInSitu (i,j,k,bi,bj) = 0. _d 0 |
83 |
IVDConvCount(i,j,k,bi,bj) = 0. _d 0 |
IVDConvCount(i,j,k,bi,bj) = 0. _d 0 |
84 |
|
#ifdef ALLOW_SOLVE4_PS_AND_DRAG |
85 |
|
dU_psFacX (i,j,k,bi,bj) = 0. _d 0 |
86 |
|
dV_psFacY (i,j,k,bi,bj) = 0. _d 0 |
87 |
|
#endif /* ALLOW_SOLVE4_PS_AND_DRAG */ |
88 |
ENDDO |
ENDDO |
89 |
ENDDO |
ENDDO |
90 |
ENDDO |
ENDDO |
99 |
ENDDO |
ENDDO |
100 |
ENDDO |
ENDDO |
101 |
|
|
102 |
|
#ifdef ALLOW_SOLVE4_PS_AND_DRAG |
103 |
|
C-- Put some constant values (if not part of S/R INI_FIELDS): |
104 |
|
DO k=1,Nr |
105 |
|
DO j=1-OLy,sNy+OLy |
106 |
|
DO i=1-OLx,sNx+OLx |
107 |
|
dU_psFacX(i,j,k,bi,bj) = maskW(i,j,k,bi,bj) |
108 |
|
& *recip_deepFacC(k)*recip_rhoFacC(k) |
109 |
|
dV_psFacY(i,j,k,bi,bj) = maskS(i,j,k,bi,bj) |
110 |
|
& *recip_deepFacC(k)*recip_rhoFacC(k) |
111 |
|
ENDDO |
112 |
|
ENDDO |
113 |
|
ENDDO |
114 |
|
#endif /* ALLOW_SOLVE4_PS_AND_DRAG */ |
115 |
|
|
116 |
ENDDO |
ENDDO |
117 |
ENDDO |
ENDDO |
118 |
|
|