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:
Unfortunately in my output file I see something like this:
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" ;
Did anybody have this problem?
Thank you for your response
