25 |
#include "PARAMS.h" |
#include "PARAMS.h" |
26 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
27 |
|
|
28 |
|
#ifdef EXACT_CONSERV |
29 |
|
COMMON /TEMP_SURF_POS/ etaStar |
30 |
|
_RL etaStar(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy) |
31 |
|
#endif |
32 |
|
|
33 |
C == Routine arguments == |
C == Routine arguments == |
34 |
C myTime - Current time in simulation |
C myTime - Current time in simulation |
35 |
C myIter - Current iteration number in simulation |
C myIter - Current iteration number in simulation |
119 |
C-- End DO K=1,Nr |
C-- End DO K=1,Nr |
120 |
ENDDO |
ENDDO |
121 |
|
|
122 |
|
#ifdef EXACT_CONSERV |
123 |
|
IF (exactConserv) THEN |
124 |
|
C-- Compute again "eta" to satisfy exactly the total Volume Conservation : |
125 |
|
CALL CALC_EXACT_ETA( bi,bj, etaStar, phiSurfX,phiSurfY, |
126 |
|
I myTime, myIter, myThid ) |
127 |
|
ENDIF |
128 |
|
#endif /* EXACT_CONSERV */ |
129 |
|
|
130 |
C-- Convectively adjust new fields to be statically stable |
C-- Convectively adjust new fields to be statically stable |
131 |
CALL CONVECTIVE_ADJUSTMENT( |
CALL CONVECTIVE_ADJUSTMENT( |
132 |
I bi, bj, iMin, iMax, jMin, jMax, |
I bi, bj, iMin, iMax, jMin, jMax, |
141 |
ENDDO |
ENDDO |
142 |
ENDDO |
ENDDO |
143 |
|
|
144 |
|
#ifdef EXACT_CONSERV |
145 |
|
IF (exactConserv) _EXCH_XY_R8(etaN, myThid ) |
146 |
|
#endif /* EXACT_CONSERV */ |
147 |
|
|
148 |
RETURN |
RETURN |
149 |
END |
END |