We have a 5D netCDF which works fine in Panoply, but when we try to use it in
thredds, we don’t get a list of variable. Digging deeper, I see that call to
/thredds/wms/gwAll/hpm_outputs.nc?item=menu&menu=&request=GetMetadata returns:
{"label":"GNS Thredds","children":[{"label":"hpm model
outputs","children":[]}]} I don’t see any errors reports in the
threddsServlet.log.
A very similar file in the same datascan but with one less dimension is working
ok.
The header for the failing file is:
netcdf hpm_outputs {
dimensions:
lay = 8 ;
x = 501 ;
y = 302 ;
realn = 91 ;
datetime = 76 ;
variables:
char transverse_mercator ;
transverse_mercator:grid_mapping_name =
"transverse_mercator" ;
transverse_mercator:longitude_of_central_meridian = 173. ;
transverse_mercator:false_easting = 1600000. ;
transverse_mercator:false_northing = 10000000. ;
transverse_mercator:latitude_of_projection_origin = 0. ;
transverse_mercator:scale_factor_at_central_meridian = 0.9996 ;
transverse_mercator:long_name = "CRS
definition" ;
transverse_mercator:longitude_of_prime_meridian
= 0. ;
transverse_mercator:semi_major_axis = 6378137. ;
transverse_mercator:semi_minor_axis =
6356752.31414036 ;
transverse_mercator:reference_ellipsoid_name =
"GRS 1980" ;
transverse_mercator:prime_meridian_name =
"Greenwich" ;
transverse_mercator:geographic_crs_name =
"NZGD2000" ;
transverse_mercator:horizontal_datum_name =
"New Zealand Geodetic Datum 2000" ;
transverse_mercator:projected_crs_name =
"NZGD2000 / New Zealand Transverse Mercator 2000" ;
transverse_mercator:inverse_flattening =
298.257222101 ;
transverse_mercator:spatial_ref =
"PROJCRS[\"NZGD2000 / New Zealand Transverse Mercator
2000\",BASEGEOGCRS[\"NZGD2000\",DATUM[\"New Zealand Geodetic Datum
2000\",ELLIPSOID[\"GRS
1980\",6378137,298.257222101,LENGTHUNIT[\"metre\",1]]],PRIMEM[\"Greenwich\",0,ANGLEUNIT[\"degree\",0.0174532925199433]],
ID[\"EPSG\",4167]],CONVERSION[\"New Zealand Transverse Mercator
2000\",METHOD[\"Transverse Mercator\", ID[\"EPSG\",9807]],PARAMETER[\"Latitude
of natural origin\",0,ANGLEUNIT[\"degree\",0.0174532925199433],
ID[\"EPSG\",8801]],PARAMETER[\"Longitude of natural
origin\",173,ANGLEUNIT[\"degree\",0.0174532925199433],
ID[\"EPSG\",8802]],PARAMETER[\"Scale factor at natural
origin\",0.9996,SCALEUNIT[\"unity\",1],ID[\"EPSG\",8805]],PARAMETER[\"False
easting\",1600000,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8806]],PARAMETER[\"False
northing\",10000000,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8807]]],
CS[Cartesian,2], AXIS[\"northing (N)\",north,ORDER[1],LENGTHUNIT[\"metre\",1]],
AXIS[\"easting
(E)\",east,ORDER[2],LENGTHUNIT[\"metre\",1]],USAGE[SCOPE[\"unknown\"],AREA[\"New
Zealand - onshore\"],BBOX[-47.33,166.37,-34.1,178.63]],ID[\"EPSG\",2193]]" ;
int lay(lay) ;
lay:standard_name = "(no standard name)" ;
lay:long_name = "Model layer" ;
double x(x) ;
x:standard_name = "projection_x_coordinate" ;
x:long_name = "x coordinate of projection" ;
x:units = "m" ;
double y(y) ;
y:standard_name = "projection_y_coordinate" ;
y:long_name = "y coordinate of projection" ;
y:units = "m" ;
string realn(realn) ;
realn:long_name = "Realisation" ;
realn:standard_name = "realization" ;
int64 datetime(datetime) ;
datetime:standard_name = "(no standard name)" ;
datetime:long_name = "Datetime" ;
datetime:units = "YYYYMMDD" ;
float head(realn, datetime, lay, y, x) ;
head:_FillValue = -1.e+30f ;
head:grid_mapping = "transverse_mercator" ;
head:long_name = "Simulated GW Head" ;
head:standard_name = "(no standard name)" ;
head:units = "m" ;
float s_flow(realn, datetime, y, x) ;
s_flow:_FillValue = -1.e+30f ;
s_flow:grid_mapping = "transverse_mercator" ;
s_flow:long_name = "Simulated Stream Flow" ;
s_flow:standard_name = "(no standard name)" ;
s_flow:units = "m3/d" ;
float s_flux(realn, datetime, y, x) ;
s_flux:_FillValue = -1.e+30f ;
s_flux:grid_mapping = "transverse_mercator" ;
s_flux:long_name = "Simulated Stream Flux to
GW" ;
s_flux:standard_name = "(no standard name)" ;
s_flux:units = "m3/d" ;
// global attributes:
:title = "hpm model outputs" ;
:description = "ensemble simulated outputs" ;
:file_creation_time = "2020-08-03
15:03:00.336462" ;
:Conventions = "CF-1.7" ;
:institution = "GNS Science" ;
:source = "PEST++/pyEMU" ;
}
And for the good file is:
netcdf hpm_inputs {
dimensions:
lay = 8 ;
x = 501 ;
y = 302 ;
realn = 100 ;
variables:
char transverse_mercator ;
transverse_mercator:grid_mapping_name =
"transverse_mercator" ;
transverse_mercator:longitude_of_central_meridian = 173. ;
transverse_mercator:false_easting = 1600000. ;
transverse_mercator:false_northing = 10000000. ;
transverse_mercator:latitude_of_projection_origin = 0. ;
transverse_mercator:scale_factor_at_central_meridian = 0.9996 ;
transverse_mercator:long_name = "CRS
definition" ;
transverse_mercator:longitude_of_prime_meridian
= 0. ;
transverse_mercator:semi_major_axis = 6378137. ;
transverse_mercator:semi_minor_axis =
6356752.31414036 ;
transverse_mercator:reference_ellipsoid_name =
"GRS 1980" ;
transverse_mercator:prime_meridian_name =
"Greenwich" ;
transverse_mercator:geographic_crs_name =
"NZGD2000" ;
transverse_mercator:horizontal_datum_name =
"New Zealand Geodetic Datum 2000" ;
transverse_mercator:projected_crs_name =
"NZGD2000 / New Zealand Transverse Mercator 2000" ;
transverse_mercator:inverse_flattening =
298.257222101 ;
transverse_mercator:spatial_ref =
"PROJCRS[\"NZGD2000 / New Zealand Transverse Mercator
2000\",BASEGEOGCRS[\"NZGD2000\",DATUM[\"New Zealand Geodetic Datum
2000\",ELLIPSOID[\"GRS
1980\",6378137,298.257222101,LENGTHUNIT[\"metre\",1]]],PRIMEM[\"Greenwich\",0,ANGLEUNIT[\"degree\",0.0174532925199433]],
ID[\"EPSG\",4167]],CONVERSION[\"New Zealand Transverse Mercator
2000\",METHOD[\"Transverse Mercator\", ID[\"EPSG\",9807]],PARAMETER[\"Latitude
of natural origin\",0,ANGLEUNIT[\"degree\",0.0174532925199433],
ID[\"EPSG\",8801]],PARAMETER[\"Longitude of natural
origin\",173,ANGLEUNIT[\"degree\",0.0174532925199433],
ID[\"EPSG\",8802]],PARAMETER[\"Scale factor at natural
origin\",0.9996,SCALEUNIT[\"unity\",1],ID[\"EPSG\",8805]],PARAMETER[\"False
easting\",1600000,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8806]],PARAMETER[\"False
northing\",10000000,LENGTHUNIT[\"metre\",1],ID[\"EPSG\",8807]]],
CS[Cartesian,2], AXIS[\"northing (N)\",north,ORDER[1],LENGTHUNIT[\"metre\",1]],
AXIS[\"easting
(E)\",east,ORDER[2],LENGTHUNIT[\"metre\",1]],USAGE[SCOPE[\"unknown\"],AREA[\"New
Zealand - onshore\"],BBOX[-47.33,166.37,-34.1,178.63]],ID[\"EPSG\",2193]]" ;
int lay(lay) ;
lay:standard_name = "(no standard name)" ;
lay:long_name = "Model layer" ;
double x(x) ;
x:standard_name = "projection_x_coordinate" ;
x:long_name = "x coordinate of projection" ;
x:units = "m" ;
double y(y) ;
y:standard_name = "projection_y_coordinate" ;
y:long_name = "y coordinate of projection" ;
y:units = "m" ;
string realn(realn) ;
realn:long_name = "Realisation" ;
realn:standard_name = "realization" ;
float K(realn, lay, y, x) ;
K:_FillValue = -1.e+30f ;
K:grid_mapping = "transverse_mercator" ;
K:long_name = "Hydraulic conductivity" ;
K:standard_name = "(no standard name)" ;
K:units = "m day-1" ;
// global attributes:
:title = "hpm model pars" ;
:description = "ensemble parameters for
MODFLOW" ;
:file_creation_time = "2020-08-03
13:14:55.609150" ;
:Conventions = "CF-1.7" ;
:institution = "GNS Science" ;
:source = "PEST++/pyEMU" ;
}
https://data-dev.gns.cri.nz/thredds/wms/gwAll/hpm_outputs1.nc?item=menu&menu=&request=GetMetadata
________________________________________________
Ngā mihi, Nā Phil Scadden
Te Aroturuki Matepā Aronuku me te Pūtaiao Raraunga
GNS Science Ltd 764 Cumberland St, Private Bag 1930,
Dunedin, New Zealand Ph +64 3 4799663, 027 3463185
“Whāia te iti kahurangi ki te tūohu koe me he maunga teitei”
Notice: This email and any attachments are confidential and may not be used,
published or redistributed without the prior written consent of the Institute
of Geological and Nuclear Sciences Limited (GNS Science). If received in error
please destroy and immediately notify GNS Science. Do not copy or disclose the
contents.
