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.
