NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
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.
thredds archives: