Tidal current analysis in ROMS

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
stvyng
Posts: 10
Joined: Tue Jan 31, 2012 10:35 pm
Location: NOAA

Tidal current analysis in ROMS

#1 Unread post by stvyng »

I got the following run error message using ROMS rev 748:

GET_NGFLD - tidal period
(Min = 4.47141644E+04 Max = 4.47141644E+04)
GET_2DFLD - tidal elevation amplitude
(Min = 0.00000000E+00 Max = 4.78616861E+00)
GET_2DFLD - tidal elevation phase angle
(Min = 1.64720212E+00 Max = 4.21640664E+00)
GET_2DFLD - tidal current inclination angle
(Min = 0.00000000E+00 Max = 6.28317406E+00)
GET_2DFLD - tidal current phase angle
(Min = 0.00000000E+00 Max = 6.28318268E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 0.00000000E+00 Max = 2.38517276E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -7.83491968E-01 Max = 2.64810100E-01)

DEF_VAR - unable to define variable: ubar_tide
in NetCDF file: ./TIDE_FORCING.nc

Here are what I used:
(1) main input file:
Aout(idUbar) == T !F ! ubar 2D U-velocity
Aout(idVbar) == T !F ! vbar 2D V-velocity
Aout(idu2dE) == T !F ! ubar_eastward 2D U-eastward at RHO-points
Aout(idv2dN) == T !F ! vbar_northward 2D V-northward at RHO-points
Aout(idFsur) == T !F ! zeta free-surface

Aout(idu3dD) == T !F ! u_detided detided 3D U-velocity
Aout(idv3dD) == T !F ! v_detided detided 3D V-velocity
Aout(idu2dD) == T !F ! ubar_detided detided 2D U-velocity
Aout(idv2dD) == T !F ! vbar_detided detided 2D V-velocity
Aout(idFsuD) == T !F ! zeta_detided detided free-surface


(2) cpp.h options:
# define RAMP_TIDES
# define SSH_TIDES
# define UV_TIDES
# define AVERAGES_DETIDE
# define WET_DRY
# define AVERAGES

However, when I specified "Aout(idu2dD) == F and Aout(idv2dD) == F" in the input file, the model ran well and created zeta_detided in average file and appended variables for zeta's harmonic analysis in the tidal forcing file.

Any hints about the issue? Thanks in advance for your help!

Post Reply