[Thread Prev][Thread Next][Index]

[las_users] Re: [ferret_users] bug with Ferret v6 ?



Hi Jaison and Jerome,
That's indeed what's happening. Interestingly, this does not always occur with variables whose missing_value flag is not in an input file; in this case a variable is are generally assigned the default missing-value flag of -1.e34. This is was something we test in our benchmark scripts before a Ferret release, but something's wrong which causes this one to get the wrong flag. I'll track this down and we'll get the fix into our upcoming bug-fix release of Ferret.
Ansley

Jaison Kurian wrote:
Hi Ansley & Jerome,
It seems like Ferret_V6 assigns "0" as the
missing value if the "variable in NetCDF file" does not have an
explicit "missing_value" attribute. (It was "-1.E+34" with earlier
versions of Ferret). Here the first value of trange is "0" and hence it will be set to "missing value" in Ferret_V6.

[jaison@mercury jerome]$ f
NOAA/PMEL TMAP
FERRET v6
Linux(g77) 2.4.21-32 - 08/23/06
30-Sep-06 01:32

yes? use 259da1dac3f7750fe1fc83930e5a24b6.nc
yes? say `trange,r=bad`
!-> MESSAGE/CONTINUE 0
0
yes? say `r_temp,r=bad`
!-> MESSAGE/CONTINUE -9999
-9999
yes? say `numprofs,r=bad`
!-> MESSAGE/CONTINUE 0
0
yes? say `numobs,r=bad`
!-> MESSAGE/CONTINUE 0
0
yes?

Variables "trange", "numprofs" and "numobs" does not have any explicit missing_value attribute, hence Ferret_V6 assigns "0"
as the missing value. For those variables (like r_temp) which have a valid missing_value attribute and value, Ferret_V6 behaves
properly.

This problem with missing_value appers to be there only for the
variables readed in from a NetCDF file. If we define a new variable
in Ferret_V6, the missing value is "-1.E+34" as with earlier versions.

Hope this helps,

Jaison

On Fri, 29 Sep 2006, Ansley Manke wrote:


Hi Jerome,
Thanks for sending the file. Indeed, running Ferret from the command line, with Ferret v5.81 the variable TRANGE is read correctly, but with Ferret v6.0

NOAA/PMEL TMAP
FERRET v6 Linux(g77) 2.4.21-32 - 08/23/06
29-Sep-06 11:34
yes? use 259da1dac3f7750fe1fc83930e5a24b6.nc
yes? list trange
VARIABLE : (hours)
FILENAME : 259da1dac3f7750fe1fc83930e5a24b6.nc
SUBSET : 2 points (TIME)
JAN-2001 / 1: ....
DEC-2002 / 2: 17496.


we'll track down why this is.

Ansley

Jerome King wrote:

Hi,

I am having an issue with using LAS with Ferret v6. I am not having this problem using LAS with Ferret 5.81. Both LAS are exactly the same (6.5.1), same xml files, same query.
So I am thinking it's a bug with Ferret v6 but I am not sure.

The issue happens after LAS creates a NetCDF file with the insitu package and asks Ferret to make a plot.
There seems to be an issue with Ferret reading "trange".
I have attached the NetCDF file created by LAS. and I have copied below the debug ferret part for both Ferret version.
I hope this helps find out what the issue is.
Hope this helps,
Jerome.

With Ferret 5.81, I get a nice plot:
-------------------------------------------------------------------------------------------------------------------------
yes? cancel mode verify
yes? cancel mode interp
yes? set mem/size=80
Cached data cleared from memory
yes? define symbol memsize = 80
yes? canc data/all
yes? go std_initialize "output/259da1dac3f7750fe1fc83930e5a24b6.nc" "1" "1" "R_TEMP"
yes? set region/x="-126.8":"-117.1"/y="34.1":"38.0"/z="0":"0"/t="01-Jan-2001 00:00:00":"31-Dec-2002 00:00:00"
yes? GO "jnls/insitu/insitu_setup.jnl" "output/259da1dac3f7750fe1fc83930e5a24b6.nc" "TAX" "CalCOFI Report Data" "Temperature"
yes? go insitu_meta_xy "output/259da1dac3f7750fe1fc83930e5a24b6.nc" TAX -126.8:-117.1 34.1:38.0 0:0
yes? FRAME/FORMAT=gif/FILE="output/3ccada4923f5f69e8e049265df5437b1.gif"
yes? quit
Waiting for program to exit...
program exited...
LAS::Server::Handler::run(): 0.003233 secs to perform all 'postExecute' activities (none typically)
-----------------------------------------------------------------------------------------------------------------------------------------

With Ferret v6, I get an error:
-------------------------------------------------------------------------------------------------------------------------
yes? cancel mode verify
yes? cancel mode interp
yes? set mem/size=80
Cached data cleared from memory
yes? define symbol memsize = 80
yes? canc data/all
yes? go std_initialize "output/259da1dac3f7750fe1fc83930e5a24b6.nc" "1" "1" "R_TEMP"
yes? set region/x="-126.8":"-117.1"/y="34.1":"38.0"/z="0":"0"/t="01-Jan-2001 00:00:00":"31-Dec-2002 00:00:00"
yes? GO "jnls/insitu/insitu_setup.jnl" "output/259da1dac3f7750fe1fc83930e5a24b6.nc" "TAX" "CalCOFI Report Data" "Temperature"
**ERROR: command syntax: bad
define axis/t=bad:17496/npoints=100/unit=hours/t0="01-JAN-2001 00:00:00" tfull
Command file, command group, or REPEAT execution aborted

killing...Waiting for program to exit...
program exited...
killing...Waiting for program to exit...
program exited...
------------------------------------------------------------------------------------------------------------------------------------------

--

~> ~> ~> ~> ~> ~> ~> ~> ~> ~> ~>
Ansley Manke, NOAA/PMEL ansley.b.manke@noaa.gov
7600 Sand Point Way NE, Seattle WA
Phone 206-526-6246, FAX 206-526-6744



[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement