Re: [las_users] Lat/Lon problems

[Laurent Demagistri <laurent.demagistri@xxxxxxxxxx>]

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


--
_________________________________________________
Laurent DEMAGISTRI
Maison de la Télédétection - Pôle IRD
500 rue J.F. Breton - 34093 Montpellier Cedex 05 - France
Tél : 33 (0)4 67 54 87 01
Fax : 33 (0)4 67 54 87 00
°º¤ø,¸¸,ø¤º°`°º¤ø,¸¸,ø¤º°`°º¤ø,¸¸,ø¤º°`°º¤ø,¸¸,ø¤º°

begin:vcard
fn:Laurent DEMAGISTRI
n:DEMAGISTRI;Laurent
org;quoted-printable:Institut de Recherche pour le D=C3=A9veloppement
adr;quoted-printable;quoted-printable;dom:500 rue Jean-Fran=C3=A7ois Breton;;Maison de la T=C3=A9l=C3=A9d=C3=A9tection - P=C3=B4le IRD;Montpellier;;34093
email;internet:laurent.demagistri@xxxxxxxxxx
tel;work:0467548701
tel;fax:0467548700
x-mozilla-html:TRUE
version:2.1
end:vcard