[Met_help] grib-Stat output

[John Halley Gotway]

Walford,

Good news.  I figured out what's going on, and it's an easy fix.

The reason why Grid-Stat isn't writing all the data in the output ASCII files is that there's a typo in the GridStatConfig file you're using.

Right now your config file contains the following line:
  mask_grids[] = [ "FULL", "FULL" ];

So Grid-Stat is applying 2 identical masking regions, both named "FULL".  When it goes to create the output NetCDF matched pairs file, it names the NetCDF variables using the name of the masking
region that was applied.  It writes the first one fine, but when it goes to write the second one, NetCDF errors out because we're trying to define a new variable using a variable name that's already
in use - since the names for the masking regions are identical.

Actually, if you modify the "output_flag" in your config file to turn OFF the NetCDF output file, you'll see that Grid-Stat runs fine without any errors.

But really, you probably just want to set mask_grids as follows:
  mask_grids[] = [ "FULL" ];
I can't see a reason why you'd want to apply the same masking region twice.  So it was probably just a typo.

John

walford02 at hotmail.com wrote:
> I am running at ncar on the lightning system.
> 
> my forecast and observations files are output from pcp_combine -
> 0_to_18.nc (forcast) and obs_0_to_18.nc (obs)
> 
> --------------------------------------------------
> From: "John Halley Gotway" <johnhg at rap.ucar.edu>
> Sent: Tuesday, October 28, 2008 2:50 PM
> To: <walford02 at hotmail.com>
> Cc: <met_help at mailman.ucar.edu>
> Subject: Re: [Met_help] grib-Stat output
> 
>> Seagayle,
>>
>> I'm not sure what's going on but I'd be happy to take a look at it.
>>
>> Please send me the following:
>> - a forecast file
>> - a corresponding observation file
>> - the GridStatConfig file you're using
>> - any masking polyline regions you're applying
>>
>> And I'll try running it here to see if I can figure out what's going on.
>>
>> Thanks,
>> John
>>
>> walford02 at hotmail.com wrote:
>>> Hello ,
>>>
>>> I am running the grid-stat tool and obtaining output, however, many
>>> of my output files contain no information. For example I only get
>>> data from
>>> two files CTN and CNT. All the other output just has the header
>>> information. The output flag[] is set to 2 for all fields. Any ideas
>>> as to why I am not getting information in all the files?
>>>
>>> Segayle
>>>
>>>
>>> ------------------------------------------------------------------------
>>>
>>> _______________________________________________
>>> Met_help mailing list
>>> Met_help at mailman.ucar.edu
>>> http://mailman.ucar.edu/mailman/listinfo/met_help
>>