Grid is not recognized by ROMS

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
aniskiewicz.paulina

Grid is not recognized by ROMS

#1 Unread post by aniskiewicz.paulina »

Hi :)

I have a problem with my grid file. In mu output file I see that model doesn't recognized grid file. My input grid looks like this:
Image (see below if you don't see the images here)

Unfortunately in my output file I see something like this:
Image

Here are the information about grid file:

Code: Select all

netcdf norfjords_160m_grid {
dimensions:
	xi_rho = 901 ;
	eta_rho = 451 ;
	xi_u = 900 ;
	eta_u = 451 ;
	xi_v = 901 ;
	eta_v = 450 ;
	bath = UNLIMITED ; // (2 currently)
	eta_psi = 450 ;
	xi_psi = 900 ;
variables:
	double xi_rho(xi_rho) ;
		xi_rho:long_name = "X coordinate of RHO-points" ;
		xi_rho:standard_name = "projection_x_coordinate" ;
		xi_rho:units = "meter" ;
	double eta_rho(eta_rho) ;
		eta_rho:long_name = "Y coordinate of RHO-points" ;
		eta_rho:standard_name = "projection_y_coordinate" ;
		eta_rho:units = "meter" ;
	double xi_u(xi_u) ;
		xi_u:long_name = "X coordinate of U-points" ;
		xi_u:standard_name = "projection_x_coordinate" ;
		xi_u:units = "meter" ;
	double eta_u(eta_u) ;
		eta_u:long_name = "Y coordinate of U-points" ;
		eta_u:standard_name = "projection_y_coordinate" ;
		eta_u:units = "meter" ;
	double xi_v(xi_v) ;
		xi_v:long_name = "X coordinate of V-points" ;
		xi_v:standard_name = "projection_x_coordinate" ;
		xi_v:units = "meter" ;
	double eta_v(eta_v) ;
		eta_v:long_name = "Y coordinate of V-points" ;
		eta_v:standard_name = "projection_y_coordinate" ;
		eta_v:units = "meter" ;
	double lon_rho(eta_rho, xi_rho) ;
		lon_rho:long_name = "longitude of RHO-points" ;
		lon_rho:standard_name = "longitude" ;
		lon_rho:units = "degrees_east" ;
	double lat_rho(eta_rho, xi_rho) ;
		lat_rho:long_name = "latitude of RHO-points" ;
		lat_rho:standard_name = "latitude" ;
		lat_rho:units = "degrees_north" ;
	double lon_u(eta_u, xi_u) ;
		lon_u:long_name = "longitude of U-points" ;
		lon_u:standard_name = "longitude" ;
		lon_u:units = "degrees_east" ;
	double lat_u(eta_u, xi_u) ;
		lat_u:long_name = "latitude of U-points" ;
		lat_u:standard_name = "latitude" ;
		lat_u:units = "degrees_north" ;
	double lon_v(eta_v, xi_v) ;
		lon_v:long_name = "longitude of V-points" ;
		lon_v:standard_name = "longitude" ;
		lon_v:units = "degrees_east" ;
	double lat_v(eta_v, xi_v) ;
		lat_v:long_name = "latitude of V-points" ;
		lat_v:standard_name = "latitude" ;
		lat_v:units = "degrees_north" ;
	double angle(eta_rho, xi_rho) ;
		angle:long_name = "angle between xi axis and east" ;
		angle:standard_name = "angle_of_rotation_from_east_to_x" ;
		angle:units = "radian" ;
		angle:coordinates = "lon_rho lat_rho" ;
		angle:grid_mapping = "grid_mapping" ;
	char spherical ;
		spherical:long_name = "grid type logical switch" ;
		spherical:flag_values = "T, F" ;
		spherical:flag_meanings = "spherical Cartesian" ;
	double pm(eta_rho, xi_rho) ;
		pm:long_name = "curvilinear coordinate metric in XI" ;
		pm:units = "meter-1" ;
		pm:field = "pm, scalar" ;
		pm:coordinates = "lon_rho lat_rho" ;
		pm:grid_mapping = "grid_mapping" ;
	double pn(eta_rho, xi_rho) ;
		pn:long_name = "curvilinear coordinate metric in ETA" ;
		pn:units = "meter-1" ;
		pn:field = "pn, scalar" ;
		pn:coordinates = "lon_rho lat_rho" ;
		pn:grid_mapping = "grid_mapping" ;
	double dndx(eta_rho, xi_rho) ;
		dndx:long_name = "xi derivative of inverse metric factor pn" ;
		dndx:units = "meter" ;
		dndx:field = "dndx, scalar" ;
		dndx:coordinates = "lon_rho lat_rho" ;
		dndx:grid_mapping = "grid_mapping" ;
	double dmde(eta_rho, xi_rho) ;
		dmde:long_name = "eta derivative of inverse metric factor pm" ;
		dmde:units = "meter" ;
		dmde:field = "pn, scalar" ;
		dmde:coordinates = "lon_rho lat_rho" ;
		dmde:grid_mapping = "grid_mapping" ;
	int grid_mapping ;
		grid_mapping:long_name = "grid mapping" ;
		grid_mapping:straight_vertical_longitude_from_pole = 70 ;
		grid_mapping:false_easting = 1480640. ;
		grid_mapping:standard_parallel = 60. ;
		grid_mapping:semi_major_axis = 6378137. ;
		grid_mapping:false_northing = 1496320. ;
		grid_mapping:grid_mapping_name = "polar_stereographic" ;
		grid_mapping:inverse_flattening = 298.257223563 ;
		grid_mapping:latitude_of_projection_origin = 90. ;
		grid_mapping:ellipsoid = "WGS84" ;
		grid_mapping:dx = 160 ;
		grid_mapping:proj4string = "+proj=stere +ellps=WGS84 +lat_0=90.0 +lat_ts=60.0 +x_0=1480640 +y_0=1496320 +lon_0=70" ;
	double hraw(bath, eta_rho, xi_rho) ;
		hraw:long_name = "Working bathymetry at RHO-points" ;
		hraw:standard_name = "sea_floor_depth" ;
		hraw:units = "meter" ;
		hraw:field = "bath, scalar" ;
		hraw:coordinates = "lon_rho lat_rho" ;
		hraw:grid_mapping = "grid_mapping" ;
	double h(eta_rho, xi_rho) ;
		h:long_name = "Final bathymetry at RHO-points" ;
		h:standard_name = "sea_floor_depth" ;
		h:units = "meter" ;
		h:field = "bath, scalar" ;
		h:coordinates = "lon_rho lat_rho" ;
		h:grid_mapping = "grid_mapping" ;
	double f(eta_rho, xi_rho) ;
		f:long_name = "Coriolis parameter at RHO-points" ;
		f:standard_name = "coriolis_parameter" ;
		f:units = "second-1" ;
		f:field = "Coriolis, scalar" ;
		f:coordinates = "lon_rho lat_rho" ;
		f:grid_mapping = "grid_mapping" ;
	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" ;
		mask_rho:coordinates = "lon_rho lat_rho" ;
		mask_rho:grid_mapping = "grid_mapping" ;
	double mask_u(eta_u, xi_u) ;
		mask_u:long_name = "mask on U-points" ;
		mask_u:option_0 = "land" ;
		mask_u:option_1 = "water" ;
		mask_u:coordinates = "lon_u lat_u" ;
		mask_u:grid_mapping = "grid_mapping" ;
	double mask_v(eta_v, xi_v) ;
		mask_v:long_name = "mask on V-points" ;
		mask_v:option_0 = "land" ;
		mask_v:option_1 = "water" ;
		mask_v:coordinates = "lon_v lat_v" ;
		mask_v:grid_mapping = "grid_mapping" ;
	double mask_psi(eta_psi, xi_psi) ;
		mask_psi:long_name = "mask on PSI-points" ;
		mask_psi:option_0 = "land" ;
		mask_psi:option_1 = "water" ;
	double xl ;
		xl:long_name = "domain length in the XI-direction" ;
		xl:units = "meter" ;
	double el ;
		el:long_name = "domain length in the ETA-direction" ;
		el:units = "meter" ;

I've tried to use this assumption:
https://code.mpimet.mpg.de/boards/1/topics/913

I used ncks function:
ncks -v h,lon_rho,lat_rho norfjords_160m_grid.nc norfjords_160m_grid6.nc

And my result was smaller grid file:

Code: Select all

 
dimensions:
	eta_rho = 451 ;
	xi_rho = 901 ;
variables:
	double eta_rho(eta_rho) ;
		eta_rho:long_name = "Y coordinate of RHO-points" ;
		eta_rho:standard_name = "projection_y_coordinate" ;
		eta_rho:units = "meter" ;
	double h(eta_rho, xi_rho) ;
		h:long_name = "Final bathymetry at RHO-points" ;
		h:standard_name = "sea_floor_depth" ;
		h:units = "meter" ;
		h:field = "bath, scalar" ;
		h:coordinates = "lon_rho lat_rho" ;
		h:grid_mapping = "grid_mapping" ;
	double lat_rho(eta_rho, xi_rho) ;
		lat_rho:long_name = "latitude of RHO-points" ;
		lat_rho:standard_name = "latitude" ;
		lat_rho:units = "degrees_north" ;
	double lon_rho(eta_rho, xi_rho) ;
		lon_rho:long_name = "longitude of RHO-points" ;
		lon_rho:standard_name = "longitude" ;
		lon_rho:units = "degrees_east" ;
	double xi_rho(xi_rho) ;
		xi_rho:long_name = "X coordinate of RHO-points" ;
		xi_rho:standard_name = "projection_x_coordinate" ;
		xi_rho:units = "meter" ;
Nevertheless, model still can not see the grid.
Did anybody have this problem?

Thank you for your response :)
Attachments
ROMS input grid
ROMS input grid
ROMS output grid
ROMS output grid
Top
User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Grid is not recognized by ROMS

#2 Unread post by kate »

I can't see what the error is from ROMS. What does it say?
Top
User avatar
wilkin
Posts: 879
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:
Contact wilkin

Re: Grid is not recognized by ROMS

#3 Unread post by wilkin »

When you say ROMS does not recognize the grid, do you mean it did not read it at all or that it gave an error. If it did not even read perhaps you have #define ANA_GRID by mistake.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu
Top
Post Reply

Return to “ROMS Problems”