[gradsusr] R: I: Why cannot contour grid for this file?

Wesley Ebisuzaki Wesley.Ebisuzaki at noaa.gov
Thu May 13 09:18:51 EDT 2010


Diego,

When I have missing vertical levels, I use the grib2ctl/g2ctl option 
-no_profile. Then
each field will have it's own name such as HGT100mb. Another option 
besides editing
the files.

Wesley

Diego Morelli wrote:
> Hello Wesley, how are you?.
>
> Yes, yes, I understood. You and Jennifer have been very clear.
>
> I just tried to edit the descriptor file ctl as you suggest, and it works
> fine for level 100.
>
> Ok, it means that, if necessary, I will split my file so that I will ensure
> that there are only, in the produced files, parameters on the same list of
> levels.
>
> Thank you vey much to you all.
>
> Diego
>
>
>
> -----Messaggio originale-----
> Da: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] Per
> conto di Wesley Ebisuzaki
> Inviato: mercoledì 12 maggio 2010 20.54
> A: GrADS Users Forum
> Oggetto: Re: [gradsusr] I: Why cannot contour grid for this file?
>
> Diego,
>
> grib2ctl.pl doesn't handle missing levels very well. grib2ctl found 3 
> pressure levels as expected.
>
>   
>> * z has 3 levels, for prs
>> zdef 3 levels
>> 850 150 100
>>     
> grib2ctl found UGRD on two levels and wrote the line
>   
>> UGRDprs 2 33,100,0 ** (profile) u wind [m/s]
>>     
> With NCEP files, the missing levels occur at the smallest pressures, so 
> "2" is ok. With
> your file, the missing level was at 850, so you have to edit the ctl 
> file and change the "2" into a "3".
>   
>> UGRDprs 3 33,100,0 ** (profile) u wind [m/s]
>>     
> For what it is worth, a grib2 version of your data wouldn't have had a 
> problem.
>
> Wesley Ebisuzaki
>
>
>
>
>
> Diego Morelli wrote:
>   
>> Hello.
>>
>> I don’t know if my question is clear, so let me repeat.
>>
>> Using in GrADS the attached file, the ugrd can’t be displaied for the 
>> level 100.
>>
>> Look at this pic in order to have clear the situation I have in GrADS
>>
>> What is wrong in the attached file that prevent GrADS to display the 
>> data?
>>
>> In the warning message I don’t see any abnormal situation: In my 
>> opinion, the “grid limits of request” appears correct in front of the 
>> “grid limits of file”, then what this message means?
>>
>> Regards,
>>
>> Diego
>>
>> * Da: * Diego Morelli [mailto:diego.morelli at trs.it]
>> *Inviato:* mercoledì 12 maggio 2010 16.34
>> *A:* ' GrADS Users Forum '
>> *Oggetto:* Why cannot contour grid for this file?
>>
>> Hello.
>>
>> I have this file in the attachment.
>>
>> It get 3 contains this parameters: ugrdprs, for levels 150/100, 
>> vgrdprs for levels 850/150/100.
>>
>> If I use this file in grads, it can display all the 3 levels for vgrdprs.
>>
>> For ugrdprs the display command produces a correct map for level 150 
>> and an empty map for level 100 (Entire grid is undefined) and this 
>> message:
>>
>> Why this behaviour? Can you say me what doesn’t like to GrADS for 
>> ugrdprs on level 100?
>>
>> The products I’m using are:
>>
>> grib2ctl v0.9.12.5p45
>>
>> wgrib v1.8.0.12
>>
>> GrADS Ver. 2.0.a7.1
>>
>> The command I’m using before grads are:
>>
>> perl grib2ctl.pl 
>> LFPW00190281.20090718151333_P4400PT7758_BPMG12.101562_1.b > 
>> LFPW00190281.20090718151333_P4400PT7758_BPMG12.101562_1.ctl
>>
>> ./gribmap – i 
>> LFPW00190281.20090718151333_P4400PT7758_BPMG12.101562_1.ctl – 0.
>>
>> Thank you in advance.
>>
>> Diego
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>   
>>     
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>   




More information about the gradsusr mailing list