free surface

[asujrpv]

To whom it may concern:

I wonder if the *h file that I have could be modified to give me free-surface. It is driven by barotropic tides in the yellow sea but even though I apply tides, the results from the fortran plot programs do not seem to provide a variation of vater surface elevation. Below is my cdw.h ... I would appreciate your help.

Rafael

#define SOLVE3D
#define SALINITY
#define UV_ADV
#define DJ_GRADPS
#define UV_COR
#define UV_QDRAG
#define NONLIN_EOS
#define SPLINES
#define GLS_MIXING /* use if Generic Length-Scale mixing */
#define N2S2_HORAVG
#define KANTHA_CLAYSON
#define MIX_S_UV /* use if mixing along constant S-surfaces */
#define MIX_S_TS /* use if mixing along constant S-surfaces */
#define RAMP_TIDES
/*
#define MASKING
*/


#define SSH_TIDES
#define UV_TIDES

#define NORTH_FSCHAPMAN
#define NORTH_M2FLATHER
#define SOUTH_FSCHAPMAN
#define SOUTH_M2FLATHER
#define EAST_FSCHAPMAN
#define EAST_M2FLATHER
#define WEST_FSCHAPMAN
#define WEST_M2FLATHER

#define ANA_INITIAL
#define ANA_FSOBC
#define ANA_M2OBC
#define ANA_SMFLUX
#define ANA_STFLUX
#define ANA_SSFLUX
#define ANA_SRFLUX
#define ANA_BSFLUX
#define ANA_BTFLUX


/* Lateral boundary conditions*/
/* #define EAST_FSGRADIENT /* use if free-surface radiation condition */
#define EAST_M3GRADIENT /* use if 3D momentum radiation condition */
#define EAST_KGRADIENT /* use if TKE fields radiation condition */
/* #define WEST_FSGRADIENT /* use if free-surface radiation condition */
#define WEST_M3GRADIENT /* use if 3D momentum radiation condition */
#define WEST_KGRADIENT /* use if TKE fields radiation condition */
/* #define NORTH_FSGRADIENT /* use if free-surface radiation condition */
#define NORTH_M3GRADIENT /* use if 3D momentum radiation condition */
#define NORTH_KGRADIENT /* use if TKE fields radiation condition */
/* #define SOUTH_FSGRADIENT /* use if free-surface radiation condition */
#define SOUTH_M3GRADIENT /* use if 3D momentum radiation condition */
#define SOUTH_KGRADIENT /* use if TKE fields radiation condition */
~

[kate]

I don't immediately see why you wouldn't get tides from this. Have you tried turning off the TIDES_RAMP? Where are you getting your tides from? Have you checked that file?

[asujrpv]

Kate,
I do provide tides and they seem to work fine, but when I plot (using the fortran plot package or any other) I do not see variation of my water surface elevation, I see the elevation set at -0.5 m, which is the non-zero initial water surface that I was told ROMS accepts (do not take zero).

Is this a problem of plotting?

thanks
Rafael

[kate]

I use ncview for a quick and dirty look at what's happening in the output.

ROMS can have a zero zeta as initial conditions. What can't be zero is zeta+depth (total water column depth).