Well I scratched my head about that code and then realized it is the code in
TDS 5.0. I tried the TDS 5 beta on my local machine and it did work on our
netCDF. It was pretty tolerant on a lot of things. Looks like we will work the
new version.
From: Sean Arms <sarms@xxxxxxxx>
Sent: Monday, 10 August 2020 23:36
To: Phil Scadden <P.Scadden@xxxxxxxxxx>
Subject: Re: [thredds] Limit on diminsion when using Godiva?
Oops! Never mind, I see you replied to John (it was just deeper down in my
inbox). Ultimately, I don't think the WMS server will work with non-spatial
dimensions like realization. That would be a quite useful feature, and the
underlying netCDF-Java library does support the ensemble dimension, so it might
be possible. I believe the code that would need to be modified lives here:
https://github.com/Reading-eScience-Centre/edal-java
Cheers,
Sean
On Mon, Aug 10, 2020 at 5:31 AM Sean Arms
<sarms@xxxxxxxx<mailto:sarms@xxxxxxxx>> wrote:
Greetings Phil,
Not sure if the message from Jon Blower came through or not, but just in case:
https://www.unidata.ucar.edu/mailing_lists/archives/thredds/2020/msg00088.html
He is the lead developer of ncWMS, which the TDS uses for its WMS service. In
his reply he mentions that the WMS server will only work with spatial
dimensions, so the ensemble dimension in your datasets isn't recognized.
Cheers,
Sean
On Sun, Aug 9, 2020 at 11:21 PM Phil Scadden
<P.Scadden@xxxxxxxxxx<mailto:P.Scadden@xxxxxxxxxx>> wrote:
We are making slow progress here. 4 dimensions are fine. Eg
netcdf hpm_outputs {
dimensions:
time = 76 ;
lay = 8 ;
x = 501 ;
y = 302 ;
etc
However, this is output from a stochastic model and we would really like the
data from all realizations. The moment we put in realization, eg:
netcdf hpm_outputs {
dimensions:
time = 76 ;
lay = 8 ;
x = 501 ;
y = 302 ;
realn = 10 ;
then, no go. Godiva only shows us latitude, longitude variables. Full header
below. Has anyone succeeded with 5D data?
netcdf hpm_outputs {
dimensions:
time = 76 ;
lay = 8 ;
x = 501 ;
y = 302 ;
realn = 10 ;
variables:
int time(time) ;
time:standard_name = "time" ;
time:axis = "T" ;
time:long_name = "time" ;
time:units = "days since 1940-07-01 00:00:00" ;
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]]" ;
double lon(y, x) ;
lon:units = "degrees_east" ;
lon:long_name = "longitude coordinate" ;
lon:standard_name = "longitude" ;
double lat(y, x) ;
lat:units = "degrees_north" ;
lat:long_name = "latitude coordinate" ;
lat:standard_name = "latitude" ;
int lay(lay) ;
lay:standard_name = "model_level_number" ;
lay:long_name = "Model layer" ;
lay:positive = "down" ;
lay:axis = "Z" ;
lay:units = "1" ;
double x(x) ;
x:standard_name = "projection_x_coordinate" ;
x:long_name = "x coordinate of projection" ;
x:axis = "X" ;
x:units = "m" ;
double y(y) ;
y:standard_name = "projection_y_coordinate" ;
y:long_name = "y coordinate of projection" ;
y:axis = "Y" ;
y:units = "m" ;
int realn(realn) ;
realn:long_name = "Realisation" ;
realn:standard_name = "realization" ;
realn:units = "1" ;
float head(realn, time, 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" ;
head:coordinates = "lat lon" ;
float s_flow(realn, time, 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" ;
s_flow:coordinates = "lat lon" ;
float s_flux(realn, time, 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" ;
s_flux:coordinates = "lat lon" ;
// global attributes:
:title = "hpm model outputs" ;
:description = "ensemble simulated outputs" ;
:file_creation_time = "2020-08-10
17:03:38.054010" ;
:Conventions = "CF-1.6" ;
:institution = "GNS Science" ;
:source = "PEST++/pyEMU" ;
}
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.
_______________________________________________
NOTE: All exchanges posted to Unidata maintained email lists are
recorded in the Unidata inquiry tracking system and made publicly
available through the web. Users who post to any of the lists we
maintain are reminded to remove any personal information that they
do not want to be made public.
thredds mailing list
thredds@xxxxxxxxxxxxxxxx<mailto:thredds@xxxxxxxxxxxxxxxx>
For list information or to unsubscribe, visit:
https://www.unidata.ucar.edu/mailing_lists/
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.
