AVERAGES_DETIDES result problem

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
jkoceans11
Posts: 4
Joined: Thu Jan 06, 2011 2:52 pm

AVERAGES_DETIDES result problem

#1 Unread post by jkoceans11 »

Hi every one

I performed detide test using the ROMS/TOMS 3.7 version.
In my cpp option, AVERAGES_DETIDE, AVERAGES, UV_TIDES, SSH_TIDES was applied.
Compiling and model run very well. To confirm detide result, average file output is hourly mean.
But result is that zeta_detided is too high from 1 to 55 hour. Initial max zeta_detided height appear higher sea level than 1000000 m. But since 60 hour, zeta_detided is stable and appear removed tide signal.

My problem is that why too high zeta_detided such as 1000000 m appear and how solved this phenomenon (or problem ?). Below figure is my result.

thanks.
Attachments
detide_test.png

prakash22
Posts: 24
Joined: Thu Jan 08, 2015 11:16 pm
Location: iit delhi

Re: AVERAGES_DETIDES result problem

#2 Unread post by prakash22 »

Dear All,

I am getting the following run error message while simulating ROMS model

GET_NGFLD - tidal period
(Min = 2.19720215E+04 Max = 2.38071489E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 0.00000000E+00 Max = 1.64313322E+00)
GET_2DFLD - tidal elevation phase angle
(Min = -2.86857319E+00 Max = 1.06235184E+01)
GET_2DFLD - tidal current inclination angle
(Min = 0.00000000E+00 Max = 6.28318242E+00)
GET_2DFLD - tidal current phase angle
(Min = 0.00000000E+00 Max = 6.28318375E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 0.00000000E+00 Max = 7.56411946E-01)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -2.30937586E-01 Max = 2.50106657E-01)

DEF_VAR - unable to define variable: zeta_tide
in NetCDF file: ./Projects/RAVI_EXP_8/tide_ravi_exp_8.nc

Here are what I used:

Aout(idUbar) == T ! ubar 2D U-velocity
Aout(idVbar) == T ! vbar 2D V-velocity
Aout(idu2dE) == T ! ubar_eastward 2D U-eastward at RHO-points
Aout(idv2dN) == T ! vbar_northward 2D V-northward at RHO-points
Aout(idFsur) == T ! zeta free-surface

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


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

But once I removed # define AVERAGES_DETIDE model is running well
Could you help regarding this issue? Thanks in advance for your help!

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: AVERAGES_DETIDES result problem

#3 Unread post by wilkin »

Do you have write permission on NetCDF file: ./Projects/RAVI_EXP_8/tide_ravi_exp_8.nc ?

AVERAGES_DETIDE is going to modify that file. ROMS reports it read successfully, so it must be in the correct path. But opening it to write there is a problem.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

prakash22
Posts: 24
Joined: Thu Jan 08, 2015 11:16 pm
Location: iit delhi

Re: AVERAGES_DETIDES result problem

#4 Unread post by prakash22 »

Thank you so much for your kind reply.

Could you tell me how to make write permission in existing netcdf of tide forcing file.

and how to check it whether it has writing permission or not.

prakash22
Posts: 24
Joined: Thu Jan 08, 2015 11:16 pm
Location: iit delhi

Re: AVERAGES_DETIDES result problem

#5 Unread post by prakash22 »

Now, I am getting following error

GET_NGFLD - tidal period
(Min = 2.19720215E+04 Max = 2.38071489E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 0.00000000E+00 Max = 1.64313322E+00)
GET_2DFLD - tidal elevation phase angle
(Min = -2.86857319E+00 Max = 1.06235184E+01)
GET_2DFLD - tidal current inclination angle
(Min = 0.00000000E+00 Max = 6.28318242E+00)
GET_2DFLD - tidal current phase angle
(Min = 0.00000000E+00 Max = 6.28318375E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 0.00000000E+00 Max = 7.56411946E-01)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -2.30937586E-01 Max = 2.50106657E-01)

DEF_VAR - unable to define variable: Hcount
in NetCDF file: ./Projects/TACLOBAN/tide_frc.nc

Please suggest me regarding this issue

kleinr8913
Posts: 1
Joined: Mon Nov 28, 2016 5:11 pm
Location: National Marine Environmental Monitoring Center

Re: AVERAGES_DETIDES result problem

#6 Unread post by kleinr8913 »

prakash22 wrote:Now, I am getting following error

GET_NGFLD - tidal period
(Min = 2.19720215E+04 Max = 2.38071489E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 0.00000000E+00 Max = 1.64313322E+00)
GET_2DFLD - tidal elevation phase angle
(Min = -2.86857319E+00 Max = 1.06235184E+01)
GET_2DFLD - tidal current inclination angle
(Min = 0.00000000E+00 Max = 6.28318242E+00)
GET_2DFLD - tidal current phase angle
(Min = 0.00000000E+00 Max = 6.28318375E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 0.00000000E+00 Max = 7.56411946E-01)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -2.30937586E-01 Max = 2.50106657E-01)

DEF_VAR - unable to define variable: Hcount
in NetCDF file: ./Projects/TACLOBAN/tide_frc.nc

Please suggest me regarding this issue
Try cover the current tide forcing file with the original one which you have used when you tried to run your case the first time.

prakash22
Posts: 24
Joined: Thu Jan 08, 2015 11:16 pm
Location: iit delhi

Re: AVERAGES_DETIDES result problem

#7 Unread post by prakash22 »

hi jkoceans11,


I am also facing similar problem like you (Zeta_detide value is too large).

Did you resolve or understand the problem? If yes than Could you tell me how to resolve this issue.


With regards
-Ravi

User avatar
susonic
Posts: 167
Joined: Tue Aug 21, 2007 5:44 pm
Location: UST21 / Korea
Contact:

Re: AVERAGES_DETIDES result problem

#8 Unread post by susonic »

Regarding this issue,
please give definition to your tide file with all available point when you build a tide file(say u,v,rho,psi dimension).

-JH
Joonho Lee

fabeobd
Posts: 13
Joined: Thu Jun 09, 2011 2:10 pm
Location: CCRC - UNSW
Contact:

Re: AVERAGES_DETIDES result problem

#9 Unread post by fabeobd »

Hi,

I'm trying to output the detide variables on my runs but I can't get it to work.

I have a tidal file:

Code: Select all

   
   TIDENAME ==/g/data3/hh5/tmp/access-om/fbd581/ROMS/waom10_frc/waom10extend_tds_nodal_2007.nc
that contains the following:

Code: Select all

dimensions:
	xi_rho = 630 ;
	eta_rho = 560 ;
	tide_period = UNLIMITED ; // (13 currently)
	string = 4 ;
variables:
	double lat_rho(eta_rho, xi_rho) ;
		lat_rho:long_name = "latitude of RHO-points" ;
		lat_rho:units = "degree_north" ;
		lat_rho:field = "lat_rho, scalar" ;
	double lon_rho(eta_rho, xi_rho) ;
		lon_rho:long_name = "longitude of RHO-points" ;
		lon_rho:units = "degree_east" ;
		lon_rho:field = "lon_rho, scalar" ;
	double mask_rho(eta_rho, xi_rho) ;
		mask_rho:long_name = "mask on RHO-points" ;
		mask_rho:option_0 = "land" ;
		mask_rho:option_1 = "water" ;
	double tide_period(tide_period) ;
		tide_period:long_name = "tide angular period" ;
		tide_period:units = "hours" ;
		tide_period:field = "tide_period, scalar " ;
	double tide_Eamp(tide_period, eta_rho, xi_rho) ;
		tide_Eamp:long_name = "tidal elevation amplitude" ;
		tide_Eamp:units = "meter" ;
		tide_Eamp:field = "tide_Eamp, scalar" ;
	double tide_Ephase(tide_period, eta_rho, xi_rho) ;
		tide_Ephase:long_name = "tidal elevation phase angle" ;
		tide_Ephase:units = "degrees, time of maximum elevation with respect to chosen time origin" ;
		tide_Ephase:field = "tide_Ephase, scalar " ;
	double tide_Cphase(tide_period, eta_rho, xi_rho) ;
		tide_Cphase:long_name = "tidal current phase angle" ;
		tide_Cphase:units = "degrees, time of maximum velocity with respect chosen time origin" ;
		tide_Cphase:field = "tide_Cphase, scalar " ;
	double tide_Cangle(tide_period, eta_rho, xi_rho) ;
		tide_Cangle:long_name = "tidal current inclination angle" ;
		tide_Cangle:units = "degrees between semi-major axis and East" ;
		tide_Cangle:field = "tide_Cangle, scalar " ;
	double tide_Cmin(tide_period, eta_rho, xi_rho) ;
		tide_Cmin:long_name = "minimum tidal current, ellipse semi-minor axis" ;
		tide_Cmin:units = "meter second-1" ;
		tide_Cmin:field = "tide_Cmin, scalar " ;
	double tide_Cmax(tide_period, eta_rho, xi_rho) ;
		tide_Cmax:long_name = "maximum tidal current, ellipse semi-major axis" ;
		tide_Cmax:units = "meter second-1" ;
		tide_Cmax:field = "tide_Cmax, scalar " ;
	char tidal_constituents(string, tide_period) ;
I see that the AVERAGES_DETIDE wants an harmonic file (default HARNAME == ocean_har.nc), but it crashes either away if

(1) if a just give a dummy name (should I expect the model to create this file?)

Code: Select all

 DEF_TIDES_NF90   - inquiring harmonics file,       Grid 01: /scratch/gi0/fbd581/waom10extend_testcode/output_21yr/ocean_har.nc
 Found Error: 02   Line: 8961     Source: ROMS/Modules/mod_netcdf.F, netcdf_open

 NETCDF_OPEN - unable to open existing NetCDF file:
               /scratch/gi0/fbd581/waom10extend_testcode/output_21yr/ocean_har.nc
               call from:  ROMS/Utility/def_tides.F, def_tides_nf90
               No such file or directory
 Found Error: 03   Line: 525      Source: ROMS/Utility/def_tides.F, def_tides_nf90

 DEF_TIDES_NF90 - unable to open harmonics NetCDF file: /scratch/gi0/fbd581/waom10extend_testcode/output_21yr/ocean_har.nc
or (2) I give my tidal file for the HARNAME)

Code: Select all

 
 DEF_TIDES_NF90   - inquiring harmonics file,       Grid 01: /g/data3/hh5/tmp/access-om/fbd581/ROMS/waom10_frc/waom10extend_tds_nodal_2007.nc

 DEF_TIDES_NF90 - unable to find variable: ocean_time   in detide harmonics NetCDF file: /g/data3/hh5/tmp/access-om/fbd581/ROMS/waom10_frc/waom10extend_tds_nodal_2007.nc
There are any instructions on what this harmonic file should look like somewhere?

Thanks in advance!
Fabio Boeira Dias, PhD
Postdoctoral researcher
Climate Change Research Centre (CCRC)
University of New South Wales

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: AVERAGES_DETIDES result problem

#10 Unread post by wilkin »

ROMS should create the file defined by HARNAME. When we first introduced this option, ROMS would add to the forcing harmonics file but that was not an option once nesting was introduced. See ticket https://www.myroms.org/projects/src/ticket/736

To get to that line in the def_tides.F code where you have the error, def_tides must have been called from get_idata.F with LdefTIDE = false. The only place I can see that flag being set to FALSE is in read_phypar.F

Code: Select all

# if defined AVERAGES  && defined AVERAGES_DETIDE && \
    (defined SSH_TIDES || defined UV_TIDES)
!
!  If restart, turn off definition of detide harmonics NetCDF file
!  since we only need to update time-averaged accumulated harmonics.
!
        IF (nrrec(ng).ne.0) THEN
          LdefTIDE(ng)=.FALSE.
        END IF
# endif

Are you trying to do a hot restart with the detide option newly introduced? That won't work. You need to treat this as starting from initial conditions with nrrec=0, then the template detide file should be created anew, and the detiding will start.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

fabeobd
Posts: 13
Joined: Thu Jun 09, 2011 2:10 pm
Location: CCRC - UNSW
Contact:

Re: AVERAGES_DETIDES result problem

#11 Unread post by fabeobd »

Hi John,

thanks a lot for your help. Indeed I was trying to do a restart from a previous run and turning on the detide option. Setting NRREC == 0 solved the problem! Thanks for sending the ticket too, I should have searched for things there as well.

Best,
Fabio
Fabio Boeira Dias, PhD
Postdoctoral researcher
Climate Change Research Centre (CCRC)
University of New South Wales

Post Reply