1 |
C $Header$ |
C $Header$ |
2 |
C $Name$ |
C $Name$ |
3 |
|
#include "CPP_OPTIONS.h" |
4 |
#include "GMREDI_OPTIONS.h" |
#include "GMREDI_OPTIONS.h" |
5 |
|
|
6 |
C !ROUTINE: GMREDI_K3D |
C !ROUTINE: GMREDI_K3D |
267 |
C Find the zonal velocity at the cell centre |
C Find the zonal velocity at the cell centre |
268 |
C The logicals here are, in order: 1/ go from grid to north/east directions |
C The logicals here are, in order: 1/ go from grid to north/east directions |
269 |
C 2/ go from C to A grid and 3/ apply the mask |
C 2/ go from C to A grid and 3/ apply the mask |
270 |
CALL rotate_uv2en_rl(uVel, vVel, ubar, vbar, .TRUE., .TRUE., |
#ifdef ALLOW_EDDYPSI |
271 |
& .TRUE.,Nr,mythid) |
IF (GM_InMomAsStress) THEN |
272 |
|
CALL rotate_uv2en_rl(uMean, vMean, ubar, vbar, .TRUE., .TRUE., |
273 |
|
& .TRUE.,Nr,mythid) |
274 |
|
ELSE |
275 |
|
#endif |
276 |
|
CALL rotate_uv2en_rl(uVel, vVel, ubar, vbar, .TRUE., .TRUE., |
277 |
|
& .TRUE.,Nr,mythid) |
278 |
|
#ifdef ALLOW_EDDYPSI |
279 |
|
ENDIF |
280 |
|
#endif |
281 |
|
|
282 |
C Square of the buoyancy frequency at the top of a grid cell |
C Square of the buoyancy frequency at the top of a grid cell |
283 |
DO k=2,Nr |
DO k=2,Nr |