Dear Laurent,Thank you for your reply. I have just talked with my fellow researchers and have been discussing the issues regarding our NetCDF files. I have come to know some other details, hence my next set of question.
Is there a way for LAS to display non-geographic data? (i.e. data that is not based on lon,lat coordinates but is based on arbitrary x,y coordinates.)
We can plot the data using Ferret. (See attached plot.) But we're wondering how to get the same plot from LAS. The plot does show a map, but its North is not along the Y axis but about 45degrees off, as you can see in the plot.
Any pointers and advice would be welcome. Regards, Anton Yoel Buenavista Laurent Demagistri wrote:
Dear Anton, I am new to LAS too and had to face that problem too. Here is part of my data conversion script in IDL. In my memory, I had to add the "..'positive', 'down'.." attribute to the NetCDF file "Z"-dimension for addXML to be able to recognize it the right way. .... zid = NCDF_VARDEF(id, 'zlev', [zlev_id], /FLOAT) NCDF_ATTPUT, id, zid, 'units', 'meters' NCDF_ATTPUT, id, zid, 'long_name', 'Depth' NCDF_ATTPUT, id, zid, *'positive', 'down'* .... Hope this helps. Regards. Laurent. Anton Yoel Buenavista a écrit , Le 18.02.2010 09:26:Hello everyone, I have just recently installed LAS v7.0.3. Everything works perfectly and the sample COADS NetCDF files provided displays perfectly as well. The NetCDF files provided by my colleagues seem to follow a convention of their own, yet the files display perfectly in Ferret. Thus, when I tried to use the addXML.sh tool, the following error was returned: "File parsed. No Lat/Lon grids found." Can someone guide me as to how I can go about this? (i.e. modify parts of the NetCDF file, use another scripting tool, etc.) Below is the dump of one of the files: ==START== netcdf 2D_001 { dimensions: two = 2 ; three = 3 ; time = UNLIMITED ; // (84 currently) xaxis = 1320 ; yaxis = 1482 ; zaxis = 1 ; zip2 = 1149078 ; zip3 = 1099600 ; variables: double time(time) ; time:long_name = "coordinate variables of time" ; time:units = "second since 2008-01-22" ; time:figure = "na" ; float h(yaxis, xaxis) ; h:long_name = "distance from still water level to sea bed" ; h:units = "meter" ; h:figure = "scalar_h" ; h:_FillValue = 9.96921e+36f ; float lon(yaxis, xaxis) ; lon:long_name = "longitude(east)" ; lon:units = "degree_east" ; lon:figure = "na" ; lon:_FillValue = 9.96921e+36f ; float lat(yaxis, xaxis) ; lat:long_name = "latitude(north)" ; lat:units = "degree_north" ; lat:figure = "na" ; lat:_FillValue = 9.96921e+36f ; float mask(yaxis, xaxis) ; mask:long_name = "computational mask. 1 is a computational cell, 0 is a land cell." ; mask:units = "none" ; mask:figure = "mask" ; mask:_FillValue = 9.96921e+36f ; float e(time, yaxis, xaxis) ; e:long_name = "water surface elevation from still water level" ; e:units = "meter" ; e:figure = "scalar" ; e:_FillValue = 9.96921e+36f ; float va(time, yaxis, xaxis) ; va:long_name = "depth averaged current speed along y-axis" ; va:units = "meter/second" ; va:figure = "vector_y" ; va:_FillValue = 9.96921e+36f ; float ua(time, yaxis, xaxis) ; ua:long_name = "depth averaged water current speed along x-axis" ; ua:units = "meter/second" ; ua:figure = "vector_x" ; ua:_FillValue = 9.96921e+36f ; // global attributes: :title = "standard output file number 1" ; :institution = "university" ; :source = "hydro program" ; :comment = "" ; :format = "version 1.1" ; :reference_point = "none" ; :previous_file = "na" ; :next_file = "2D_002.nc" ; :created = "2010-02-17 12:35:04" ; :first_record = "2008-01-23 00:00:00" ; :last_record = "2008-01-23 20:45:00" ; } ==END== -- Regards, Anton Yoel Buenavista
Attachment:
TMH.gif
Description: GIF image