Re: GFS precipitaion data

Hi Dave,
That is a great suggestion; I wasn't aware of that option. Offline Jean
determined that while reading in the different files in a loop he was
exactly specifying the variable A_PCP_3_SFC_acc_3h, which works for one
file but doesn't work for the next file, which has the variable
A_PCP_3_SFC_acc_6h in it.
Adam

On 07/27/2010 03:35 PM, David Brown wrote:
> Hi Jean,
>
> Although it would still be good for you to send a representative
> sample of your data files, you might be able to fix this problem
> using the file option "TimePeriodSuffix". Try
> putting the following line prior to your addfiles call:
>
> setfileoption("grb", "TimePeriodSuffix", False)
>
> This will remove the time period suffixes and make all instances of
> the variable have the spelling, A_PCP_3_SFC_acc.
> -dave
>
>
> On Jul 27, 2010, at 1:02 PM, Adam Phillips wrote:
>
>> Hi Jean,
>> Can you ftp the data file over so we can take a look?
>>
>> ftp ftp.cgd.ucar.edu
>> anonymous
>> Use your email address for the password
>> cd incoming
>> put ...
>> quit
>>
>> Then email me back the name of the file, as we cannot see what is
>> placed
>> in the incoming directory.. Thanks!
>> Adam
>>
>> On 07/27/2010 12:01 PM, Jean Peres wrote:
>>> Hello,
>>>
>>> I'm trying to open GFS data to plot the precipitation variable but
>>> when
>>> I put this command:
>>>
>>> precip(:,:,i) = f[i]->A_PCP_3_SFC_acc3h
>>>
>>> I get the following error:
>>>
>>> fatal:No valid instance of variable A_PCP_3_SFC_acc3h found in file
>>> list
>>>
>>> So I ran the ncl_filedump and the result was that:
>>>
>>> float A_PCP_3_SFC_acc3h ( lat_3, lon_3 )
>>> center : US National Weather Service - NCEP (WMC)
>>> long_name : Total precipitation
>>> units : kg/m^2
>>> _FillValue : 1e+20
>>> level_indicator : 1
>>> grid_number : 3
>>> parameter_table_version : 2
>>> parameter_number : 61
>>> model : Global Forecast System Model (formerly known as
>>> the
>>> Aviation), T382 - Forecast hours 00-180, T190 - Forecast hours
>>> 192-384
>>> forecast_time : 9
>>> forecast_time_units : hours
>>> initial_time : 07/27/2010 (00:00)
>>>
>>> however when I use the getfilevarnames command:
>>>
>>> vNames = getfilevarnames (f[0])
>>> print (vNames)
>>>
>>> it show me the following name to preciptation variable:
>>>
>>> A_PCP_3_SFC_acc6h'
>>>
>>> And both variable give me the same error.
>>>
>>> What do I do in this case?
>>>
>>> Thanks
>>>
>>> --
>>> Jean Peres
>>> Meteorologista
>>> Laboratório MASTER
>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>> --
>> __________________________________________________
>> Adam Phillips
>> asphilli@ucar.edu
>> National Center for Atmospheric Research tel: (303) 497-1726
>> Climate and Global Dynamics Division fax: (303) 497-1333
>> P.O. Box 3000
>> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
__________________________________________________
Adam Phillips 
asphilli@ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
Climate and Global Dynamics Division         fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk