Can't add the tidal forcing!

[backkom]

Hi,everyone!
i want to run a test 2D case only forced by tidal harmonics! the grid,bathy and tidal harmonics are given by external input files while others are described by analitical.but the model results are all zero which indicates the tidal forcing is not imposed,but why???
the header file is given as follows:
/* Basic physics options */
#define UV_ADV
#define UV_COR
#define UV_QDRAG
#define TS_FIXED

/* Basic numerics options */
#define UV_SADVECTION
#define CURVGRID
#define MASKING

/* Outputs */
#define AVERAGES
#define AVERAGES_QUADRATIC
#define AVERAGES_FLUXES
#define DIAGNOSTICS_UV
#define DIAGNOSTICS_TS
#define STATIONS

/* All forcing except tidal will be difined by analitical*/
#define ANA_SMFLUX /* All set to zero */
#define ANA_INITIAL /* All set to zero */
#define ANA_BSFLUX /* All set to zero */
#define ANA_BTFLUX /* All set to zero */
#define ANA_FSOBC /* All set to zero */
#define ANA_M2OBC /* All set to zero */



/*define boundary conditon */

# define SSH_TIDES
# define UV_TIDES
# define RAMP_TIDES
# define ADD_FSOBC
# define ADD_M2OBC

# define EAST_FSGRADIENT
# define EAST_M2GRADIENT

# define SOUTH_FSGRADIENT
# define SOUTH_M2GRADIENT

# define WEST_FSGRADIENT
# define WEST_M2GRADIENT

# define NORTH_FSGRADIENT
# define NORTH_M2GRADIENT


Thanks in advance

Yao

[kate]

Do you have a tides file? Is the model reading it during initialization? I get lines like:

Code: Select all

    GET_NGFLD   - tidal period
                   (Min =  4.30820446E+04 Max =  9.67260841E+04)
    GET_2DFLD   - tidal elevation amplitude
                   (Min = -9.99900000E+03 Max =  1.52504425E+01)
    GET_2DFLD   - tidal elevation phase angle
                   (Min = -1.74515472E+02 Max =  6.28274642E+00)
    GET_2DFLD   - tidal current inclination angle
                   (Min = -1.74515472E+02 Max =  3.14159265E+00)

and so on. So you can see the min and max are not zero, though that -9999 minimum amplitude is worrisome. If yours looks like that too, then I suggest a debugger.

[backkom]

thanks kate for your kind reply!
when runs,the output is as followd:

INITIAL: Configuring and initializing forward nonlinear model ...


Minimum X-grid spacing, DXmin = 6.93178755E+00 km
Maximum X-grid spacing, DXmax = 8.64595367E+00 km
Minimum Y-grid spacing, DYmin = 6.94181911E+00 km
Maximum Y-grid spacing, DYmax = 8.64369654E+00 km

Minimum barotropic Courant Number = 3.74115727E-02
Maximum barotropic Courant Number = 1.31520820E+00
Maximum Coriolis Courant Number = 2.89353685E-03

GET_NGFLD - tidal period
(Min = 4.47141598E+04 Max = 4.47141598E+04)
GET_2DFLD - tidal elevation amplitude
(Min = 3.01320199E-04 Max = 1.76557302E+00)
GET_2DFLD - tidal elevation phase angle
(Min = 1.45915342E-04 Max = 6.28262751E+00)
GET_2DFLD - tidal current inclination angle
(Min = 6.05481596E-05 Max = 6.28317587E+00)
GET_2DFLD - tidal current phase angle
(Min = 1.75422949E-05 Max = 6.28316549E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 3.65175881E-04 Max = 1.44343474E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -3.40958586E-01 Max = 3.59562539E-01)

Initial basin volumes: TotVolume = 6.04932098207068E+15 m3
MinVolume = 2.74901520081810E+08 m3
MaxVolume = 4.90707078019366E+11 m3
Max/Min = 1.78502860905728E+03

which should indicates that reading and initialization are ok,right?
but the output is still all zeros.
as you say, i should try debugger? but how?
thanks in advance !

[kate]

Do you have gdb? Last I checked, ddd didn't handle Fortran. Here, we have a totalview license so that is what I use. The more you learn about your available debuggers this time, the easier they will be to use next time. Then again, you can always try print statements. Search on TIDES to see how they ought to be coming in and examine various values there.