[Thread Prev][Thread Next][Index]
Re: ferret -gui gives errors with some longitudes (while non-guiferretdoes'nt)
Hi Emanuele,
I agree that the Ferret GUI should work reliably on files that the Ferret command line
works on. Unfortunately, we lost financial support for the Ferret GUI about five years
ago and no further work has been done on it since then. (The miracle is that it still
functions mostly after many Ferret upgrades!)
In older versions of Ferret the use of negative longitude values was a problem. There
is a high probability that this is the source of the problem in the GUI that you
reported. If your data file is not too large you can use command line Ferret to make a
new file that has positive longitude encodings. Something close to these commands:
yes? use topot.nc
yes? define symbol axname `topot,return=xaxis`
yes? define axis/x/units=degrees ($axname) = 360 + X[g=topot]
yes? save/file=new_name.nc topot
Hope this helps - steve
emanuele.lombardi@casaccia.enea.it wrote:
> Dear Ferreters,
>
> I'm having troubles using Ferret -gui with some netcdf files.
> I have same netcdf file which can be happily plotted using the non-GUI ferret,
> while they result in error while plotted or accessed by the ferret -gui
>
> Ferret version is 5.22 and it runs on COMPAQ tru64 v5.1A O.S.
>
> I have a topography netcdf file which longitudes go from -13 to 36 degrees_east.
> Ferret -gui returns the following error whenever I ask for
> "View-->Variable Values" or I try to plot the data,
>
> **ERROR: illegal limits: "TOPOT[D=TOPOT.NC]" \
> is not in the range X=13W:36E(396) \
> Axis extremes are X=13.1W(-13.1):36.1E
>
> The strange is that non-gui Ferret properly works with this file.
>
> If I produce another netcdf file with the same data but the longitudes, (adding 360
> to them so that they range from 347 to 396 degrees_east) then both
> Ferret -gui and non-gui Ferret work properly.
>
> So I'm here asking for your help to understand where I do make mistakes.
> After my signature I write some details to better show the problem.
>
> Thank you very much from Italy,
> Emanuele Lombardi
>
> $$$ mail: AMB-GEM-CLIM ENEA Casaccia
> $$$ I-00060 S.M. di Galeria (RM) ITALY
> $$$ mailto:emanuele.lombardi@casaccia.enea.it
> $$$ tel +39 06 30483366 fax +39 06 30484264
> $$$
> $$$* Contrary to popular belief, UNIX is user friendly.
> $$$ It's just very particular about who it makes friends with.
> $$$* Computers are not intelligent, but they think they are.
> $$$* True programmers never die, they just branch to an odd address
> $$$* THIS TRANSMISSION WAS MADE POSSIBLE BY 100% RECYCLED ELECTRONS
>
> --------- DETAILS -----------------------
>
> The "guilty" netcdf file is named topot.nc and it has the following ncdump
>
> netcdf topot {
> dimensions:
> nxt = 197 ;
> nyt = 57 ;
> variables:
> float nxt(nxt) ;
> nxt:var_code = 5001 ;
> nxt:long_name = "longitude" ;
> nxt:short_name = "lon" ;
> nxt:units = "degrees_east" ;
> nxt:coordinates = "nxt" ;
> nxt:axis = "X" ;
> float nyt(nyt) ;
> nyt:var_code = 5000 ;
> nyt:long_name = "latitude" ;
> nyt:short_name = "lat" ;
> nyt:units = "degrees_north" ;
> nyt:coordinates = "nyt" ;
> nyt:axis = "Y" ;
> float TOPOT(nyt, nxt) ;
> TOPOT:var_code = 6000 ;
> TOPOT:long_name = "topography" ;
> TOPOT:short_name = "topo" ;
> TOPOT:units = "m" ;
> TOPOT:coordinates = "nxt nyt" ;
> TOPOT:_FillValue = 9.96921e+36f ;
>
> // global attributes:
> :institution = "ENEA-CLIMA" ;
> :filetype = "type7" ;
> :model_code = 2 ;
> :run_date = "2001-12-06 17:26:06" ;
> :user_name = "lele" ;
> data:
>
> nxt = -13, -12.75, -12.5, -12.25, -12, -11.75, -11.5, -11.25, -11, -10.75,
> -10.5, -10.25, -10, -9.75, -9.5, -9.25, -9, -8.75, -8.5, -8.25, -8,
> -7.75, -7.5, -7.25, -7, -6.75, -6.5, -6.25, -6, -5.75, -5.5, -5.25, -5,
> -4.75, -4.5, -4.25, -4, -3.75, -3.5, -3.25, -3, -2.75, -2.5, -2.25, -2,
> -1.75, -1.5, -1.25, -1, -0.75, -0.5, -0.25, 0, 0.25, 0.5, 0.75, 1, 1.25,
> 1.5, 1.75, 2, 2.25, 2.5, 2.75, 3, 3.25, 3.5, 3.75, 4, 4.25, 4.5, 4.75, 5,
> 5.25, 5.5, 5.75, 6, 6.25, 6.5, 6.75, 7, 7.25, 7.5, 7.75, 8, 8.25, 8.5,
> 8.75, 9, 9.25, 9.5, 9.75, 10, 10.25, 10.5, 10.75, 11, 11.25, 11.5, 11.75,
> 12, 12.25, 12.5, 12.75, 13, 13.25, 13.5, 13.75, 14, 14.25, 14.5, 14.75,
> 15, 15.25, 15.5, 15.75, 16, 16.25, 16.5, 16.75, 17, 17.25, 17.5, 17.75,
> 18, 18.25, 18.5, 18.75, 19, 19.25, 19.5, 19.75, 20, 20.25, 20.5, 20.75,
> 21, 21.25, 21.5, 21.75, 22, 22.25, 22.5, 22.75, 23, 23.25, 23.5, 23.75,
> 24, 24.25, 24.5, 24.75, 25, 25.25, 25.5, 25.75, 26, 26.25, 26.5, 26.75,
> 27, 27.25, 27.5, 27.75, 28, 28.25, 28.5, 28.75, 29, 29.25, 29.5, 29.75,
> 30, 30.25, 30.5, 30.75, 31, 31.25, 31.5, 31.75, 32, 32.25, 32.5, 32.75,
> 33, 33.25, 33.5, 33.75, 34, 34.25, 34.5, 34.75, 35, 35.25, 35.5, 35.75, 36 ;
> }
>
> The "OK" file is named topot_ORIG.nc and it differs from the "guilty" one
> only for the longitudes values.
>
> netcdf topot_ORIG {
> dimensions:
> nxt = 197 ;
> nyt = 57 ;
> variables:
> float nxt(nxt) ;
> nxt:var_code = 5001 ;
> nxt:long_name = "longitude" ;
> nxt:short_name = "lon" ;
> nxt:units = "degrees_east" ;
> nxt:coordinates = "nxt" ;
> nxt:axis = "X" ;
> float nyt(nyt) ;
> nyt:var_code = 5000 ;
> nyt:long_name = "latitude" ;
> nyt:short_name = "lat" ;
> nyt:units = "degrees_north" ;
> nyt:coordinates = "nyt" ;
> nyt:axis = "Y" ;
> float TOPOT(nyt, nxt) ;
> TOPOT:var_code = 6000 ;
> TOPOT:long_name = "topography" ;
> TOPOT:short_name = "topo" ;
> TOPOT:units = "m" ;
> TOPOT:coordinates = "nxt nyt" ;
> TOPOT:_FillValue = 9.96921e+36f ;
>
> // global attributes:
> :institution = "ENEA-CLIMA" ;
> :filetype = "type7" ;
> :run_date = "2001-11-29 11:07:06" ;
> :user_name = "Rupolo" ;
> :model_code = 1 ;
> data:
>
> nxt = 347, 347.25, 347.5, 347.75, 348, 348.25, 348.5, 348.75, 349, 349.25,
> 349.5, 349.75, 350, 350.25, 350.5, 350.75, 351, 351.25, 351.5, 351.75,
> 352, 352.25, 352.5, 352.75, 353, 353.25, 353.5, 353.75, 354, 354.25,
> 354.5, 354.75, 355, 355.25, 355.5, 355.75, 356, 356.25, 356.5, 356.75,
> 357, 357.25, 357.5, 357.75, 358, 358.25, 358.5, 358.75, 359, 359.25,
> 359.5, 359.75, 360, 360.25, 360.5, 360.75, 361, 361.25, 361.5, 361.75,
> 362, 362.25, 362.5, 362.75, 363, 363.25, 363.5, 363.75, 364, 364.25,
> 364.5, 364.75, 365, 365.25, 365.5, 365.75, 366, 366.25, 366.5, 366.75,
> 367, 367.25, 367.5, 367.75, 368, 368.25, 368.5, 368.75, 369, 369.25,
> 369.5, 369.75, 370, 370.25, 370.5, 370.75, 371, 371.25, 371.5, 371.75,
> 372, 372.25, 372.5, 372.75, 373, 373.25, 373.5, 373.75, 374, 374.25,
> 374.5, 374.75, 375, 375.25, 375.5, 375.75, 376, 376.25, 376.5, 376.75,
> 377, 377.25, 377.5, 377.75, 378, 378.25, 378.5, 378.75, 379, 379.25,
> 379.5, 379.75, 380, 380.25, 380.5, 380.75, 381, 381.25, 381.5, 381.75,
> 382, 382.25, 382.5, 382.75, 383, 383.25, 383.5, 383.75, 384, 384.25,
> 384.5, 384.75, 385, 385.25, 385.5, 385.75, 386, 386.25, 386.5, 386.75,
> 387, 387.25, 387.5, 387.75, 388, 388.25, 388.5, 388.75, 389, 389.25,
> 389.5, 389.75, 390, 390.25, 390.5, 390.75, 391, 391.25, 391.5, 391.75,
> 392, 392.25, 392.5, 392.75, 393, 393.25, 393.5, 393.75, 394, 394.25,
> 394.5, 394.75, 395, 395.25, 395.5, 395.75, 396 ;
> }
>
> Nothing bad happens and ferret -gui succesfully plots it.
> The same it happens to non-gui Ferret wich properly reads and plot the file
> (as supposed, of course!)
>
> I don't attach the netcdf files to this mail, but if you want to have a
> look, please get them from anonymous ftp:/canova.casaccia.enea.it
>
> the "guilty" file is named topot.nc (46864 bytes)
> the "OK" file is named topot_ORIG.nc (46960 bytes)
--
Steve Hankin
NOAA/PMEL, 7600 Sand Point Way NE, Seattle, WA 98115-0070
ph. (206) 526-6080 -- FAX (206) 526-6744
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement