[Met_help] [rt.rap.ucar.edu #67713] History for interpolation using constant height fields
John Halley Gotway via RT
met_help at ucar.edu
Fri Jan 22 13:18:25 MST 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hi John,
When I use temperature (on constant pressure fields made by UPP) to
verify the vertical profile, the foretasted field gets interpolated to
the observation height. Everything looks very smooth (see
constant_pressure_field.png).
Now when I use temperature (on constant height fields made by UPP), my
foretasted field does not get interpolated (see
constant_height_field.png).
I am assuming this has something to do with the fact that the fields
with constant height are in the same grib file that has fields of
constant pressure (and the native coordinate in the UPP grib file is
pressure, not sure how to make it height). my wgrib is below. Any idea
on how to interpolate the forecast fields?
Thanks,
Travis
(note: I am using VTMP)
30:569006:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=50:TR=0:P1=0:P2=0:TimeU=1:50 mb:anl:NAve=0
46:853000:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=100:TR=0:P1=0:P2=0:TimeU=1:100 mb:anl:NAve=0
61:1130366:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=150:TR=0:P1=0:P2=0:TimeU=1:150 mb:anl:NAve=0
76:1427466:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=200:TR=0:P1=0:P2=0:TimeU=1:200 mb:anl:NAve=0
91:1737568:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=250:TR=0:P1=0:P2=0:TimeU=1:250 mb:anl:NAve=0
106:2041070:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=300:TR=0:P1=0:P2=0:TimeU=1:300 mb:anl:NAve=0
121:2357640:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=350:TR=0:P1=0:P2=0:TimeU=1:350 mb:anl:NAve=0
136:2667676:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=400:TR=0:P1=0:P2=0:TimeU=1:400 mb:anl:NAve=0
151:2975534:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=450:TR=0:P1=0:P2=0:TimeU=1:450 mb:anl:NAve=0
166:3276924:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=500:TR=0:P1=0:P2=0:TimeU=1:500 mb:anl:NAve=0
181:3573958:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=550:TR=0:P1=0:P2=0:TimeU=1:550 mb:anl:NAve=0
196:3868814:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=600:TR=0:P1=0:P2=0:TimeU=1:600 mb:anl:NAve=0
211:4165848:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=650:TR=0:P1=0:P2=0:TimeU=1:650 mb:anl:NAve=0
226:4465060:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=700:TR=0:P1=0:P2=0:TimeU=1:700 mb:anl:NAve=0
241:4903598:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=750:TR=0:P1=0:P2=0:TimeU=1:750 mb:anl:NAve=0
256:5344314:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=775:TR=0:P1=0:P2=0:TimeU=1:775 mb:anl:NAve=0
271:5778496:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=800:TR=0:P1=0:P2=0:TimeU=1:800 mb:anl:NAve=0
286:6210500:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=825:TR=0:P1=0:P2=0:TimeU=1:825 mb:anl:NAve=0
301:6642504:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=850:TR=0:P1=0:P2=0:TimeU=1:850 mb:anl:NAve=0
316:7074508:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=875:TR=0:P1=0:P2=0:TimeU=1:875 mb:anl:NAve=0
331:7506512:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=900:TR=0:P1=0:P2=0:TimeU=1:900 mb:anl:NAve=0
346:7945050:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=910:TR=0:P1=0:P2=0:TimeU=1:910 mb:anl:NAve=0
361:8379232:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=920:TR=0:P1=0:P2=0:TimeU=1:920 mb:anl:NAve=0
376:8813414:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=930:TR=0:P1=0:P2=0:TimeU=1:930 mb:anl:NAve=0
391:9247596:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=940:TR=0:P1=0:P2=0:TimeU=1:940 mb:anl:NAve=0
406:9686134:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=950:TR=0:P1=0:P2=0:TimeU=1:950 mb:anl:NAve=0
421:10126850:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=960:TR=0:P1=0:P2=0:TimeU=1:960 mb:anl:NAve=0
436:10565388:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=970:TR=0:P1=0:P2=0:TimeU=1:970 mb:anl:NAve=0
451:11003926:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=980:TR=0:P1=0:P2=0:TimeU=1:980 mb:anl:NAve=0
466:11444642:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=990:TR=0:P1=0:P2=0:TimeU=1:990 mb:anl:NAve=0
480:11852594:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=1000:TR=0:P1=0:P2=0:TimeU=1:1000 mb:anl:NAve=0
578:13222406:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=30:TR=0:P1=0:P2=0:TimeU=1:30 m above MSL:anl:NAve=0
579:13240112:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=80:TR=0:P1=0:P2=0:TimeU=1:80 m above MSL:anl:NAve=0
580:13258406:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=140:TR=0:P1=0:P2=0:TimeU=1:140 m above MSL:anl:NAve=0
581:13277194:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=200:TR=0:P1=0:P2=0:TimeU=1:200 m above MSL:anl:NAve=0
582:13296424:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=300:TR=0:P1=0:P2=0:TimeU=1:300 m above MSL:anl:NAve=0
583:13316380:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=400:TR=0:P1=0:P2=0:TimeU=1:400 m above MSL:anl:NAve=0
584:13336962:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=500:TR=0:P1=0:P2=0:TimeU=1:500 m above MSL:anl:NAve=0
585:13358144:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=600:TR=0:P1=0:P2=0:TimeU=1:600 m above MSL:anl:NAve=0
586:13380028:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=700:TR=0:P1=0:P2=0:TimeU=1:700 m above MSL:anl:NAve=0
587:13402522:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=800:TR=0:P1=0:P2=0:TimeU=1:800 m above MSL:anl:NAve=0
588:13425618:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=900:TR=0:P1=0:P2=0:TimeU=1:900 m above MSL:anl:NAve=0
589:13449322:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1000:TR=0:P1=0:P2=0:TimeU=1:1000 m above MSL:anl:NAve=0
590:13473584:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1200:TR=0:P1=0:P2=0:TimeU=1:1200 m above MSL:anl:NAve=0
591:13498886:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1400:TR=0:P1=0:P2=0:TimeU=1:1400 m above MSL:anl:NAve=0
592:13526228:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1600:TR=0:P1=0:P2=0:TimeU=1:1600 m above MSL:anl:NAve=0
593:13556146:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1800:TR=0:P1=0:P2=0:TimeU=1:1800 m above MSL:anl:NAve=0
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: interpolation using constant height fields
From: Randy Bullock
Time: Tue Jun 17 11:02:07 2014
Travis -
I had a look at the code and also a few old met_help questions in
order to
get
a handle on this. As far as I can tell, interpolation should be
happening. The only
times when vertical interpolation doesn't happen is when you're
dealing
with surface
fields, but it looks from your config file that that's not the case
here.
I'll talk this over with John when gets back from vacation next week,
and
see if he
has any insights on this.
Randy
On Mon, Jun 16, 2014 at 8:07 PM, Travis Wilson via RT
<met_help at ucar.edu>
wrote:
>
> Mon Jun 16 20:07:20 2014: Request 67713 was acted upon.
> Transaction: Ticket created by Wilson0028 at ucla.edu
> Queue: met_help
> Subject: interpolation using constant height fields
> Owner: Nobody
> Requestors: Wilson0028 at ucla.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=67713 >
>
>
> Hi John,
> When I use temperature (on constant pressure fields made by
UPP) to
> verify the vertical profile, the foretasted field gets interpolated
to
> the observation height. Everything looks very smooth (see
> constant_pressure_field.png).
>
> Now when I use temperature (on constant height fields made
by
> UPP), my
> foretasted field does not get interpolated (see
> constant_height_field.png).
>
> I am assuming this has something to do with the fact that the fields
> with constant height are in the same grib file that has fields of
> constant pressure (and the native coordinate in the UPP grib file is
> pressure, not sure how to make it height). my wgrib is below. Any
idea
> on how to interpolate the forecast fields?
>
> Thanks,
> Travis
>
> (note: I am using VTMP)
>
30:569006:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=50:TR=0:P1=0:P2=0:TimeU=1:50
> mb:anl:NAve=0
>
46:853000:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=100:TR=0:P1=0:P2=0:TimeU=1:100
> mb:anl:NAve=0
>
61:1130366:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=150:TR=0:P1=0:P2=0:TimeU=1:150
> mb:anl:NAve=0
>
76:1427466:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=200:TR=0:P1=0:P2=0:TimeU=1:200
> mb:anl:NAve=0
>
91:1737568:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=250:TR=0:P1=0:P2=0:TimeU=1:250
> mb:anl:NAve=0
>
106:2041070:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=300:TR=0:P1=0:P2=0:TimeU=1:300
> mb:anl:NAve=0
>
121:2357640:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=350:TR=0:P1=0:P2=0:TimeU=1:350
> mb:anl:NAve=0
>
136:2667676:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=400:TR=0:P1=0:P2=0:TimeU=1:400
> mb:anl:NAve=0
>
151:2975534:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=450:TR=0:P1=0:P2=0:TimeU=1:450
> mb:anl:NAve=0
>
166:3276924:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=500:TR=0:P1=0:P2=0:TimeU=1:500
> mb:anl:NAve=0
>
181:3573958:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=550:TR=0:P1=0:P2=0:TimeU=1:550
> mb:anl:NAve=0
>
196:3868814:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=600:TR=0:P1=0:P2=0:TimeU=1:600
> mb:anl:NAve=0
>
211:4165848:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=650:TR=0:P1=0:P2=0:TimeU=1:650
> mb:anl:NAve=0
>
226:4465060:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=700:TR=0:P1=0:P2=0:TimeU=1:700
> mb:anl:NAve=0
>
241:4903598:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=750:TR=0:P1=0:P2=0:TimeU=1:750
> mb:anl:NAve=0
>
256:5344314:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=775:TR=0:P1=0:P2=0:TimeU=1:775
> mb:anl:NAve=0
>
271:5778496:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=800:TR=0:P1=0:P2=0:TimeU=1:800
> mb:anl:NAve=0
>
286:6210500:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=825:TR=0:P1=0:P2=0:TimeU=1:825
> mb:anl:NAve=0
>
301:6642504:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=850:TR=0:P1=0:P2=0:TimeU=1:850
> mb:anl:NAve=0
>
316:7074508:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=875:TR=0:P1=0:P2=0:TimeU=1:875
> mb:anl:NAve=0
>
331:7506512:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=900:TR=0:P1=0:P2=0:TimeU=1:900
> mb:anl:NAve=0
>
346:7945050:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=910:TR=0:P1=0:P2=0:TimeU=1:910
> mb:anl:NAve=0
>
361:8379232:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=920:TR=0:P1=0:P2=0:TimeU=1:920
> mb:anl:NAve=0
>
376:8813414:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=930:TR=0:P1=0:P2=0:TimeU=1:930
> mb:anl:NAve=0
>
391:9247596:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=940:TR=0:P1=0:P2=0:TimeU=1:940
> mb:anl:NAve=0
>
406:9686134:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=950:TR=0:P1=0:P2=0:TimeU=1:950
> mb:anl:NAve=0
>
421:10126850:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=960:TR=0:P1=0:P2=0:TimeU=1:960
> mb:anl:NAve=0
>
436:10565388:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=970:TR=0:P1=0:P2=0:TimeU=1:970
> mb:anl:NAve=0
>
451:11003926:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=980:TR=0:P1=0:P2=0:TimeU=1:980
> mb:anl:NAve=0
>
466:11444642:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=990:TR=0:P1=0:P2=0:TimeU=1:990
> mb:anl:NAve=0
>
480:11852594:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=1000:TR=0:P1=0:P2=0:TimeU=1:1000
> mb:anl:NAve=0
>
578:13222406:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=30:TR=0:P1=0:P2=0:TimeU=1:30
> m above MSL:anl:NAve=0
>
579:13240112:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=80:TR=0:P1=0:P2=0:TimeU=1:80
> m above MSL:anl:NAve=0
>
580:13258406:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=140:TR=0:P1=0:P2=0:TimeU=1:140
> m above MSL:anl:NAve=0
>
581:13277194:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=200:TR=0:P1=0:P2=0:TimeU=1:200
> m above MSL:anl:NAve=0
>
582:13296424:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=300:TR=0:P1=0:P2=0:TimeU=1:300
> m above MSL:anl:NAve=0
>
583:13316380:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=400:TR=0:P1=0:P2=0:TimeU=1:400
> m above MSL:anl:NAve=0
>
584:13336962:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=500:TR=0:P1=0:P2=0:TimeU=1:500
> m above MSL:anl:NAve=0
>
585:13358144:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=600:TR=0:P1=0:P2=0:TimeU=1:600
> m above MSL:anl:NAve=0
>
586:13380028:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=700:TR=0:P1=0:P2=0:TimeU=1:700
> m above MSL:anl:NAve=0
>
587:13402522:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=800:TR=0:P1=0:P2=0:TimeU=1:800
> m above MSL:anl:NAve=0
>
588:13425618:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=900:TR=0:P1=0:P2=0:TimeU=1:900
> m above MSL:anl:NAve=0
>
589:13449322:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1000:TR=0:P1=0:P2=0:TimeU=1:1000
> m above MSL:anl:NAve=0
>
590:13473584:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1200:TR=0:P1=0:P2=0:TimeU=1:1200
> m above MSL:anl:NAve=0
>
591:13498886:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1400:TR=0:P1=0:P2=0:TimeU=1:1400
> m above MSL:anl:NAve=0
>
592:13526228:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1600:TR=0:P1=0:P2=0:TimeU=1:1600
> m above MSL:anl:NAve=0
>
593:13556146:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1800:TR=0:P1=0:P2=0:TimeU=1:1800
> m above MSL:anl:NAve=0
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #67713] interpolation using constant height fields
From: Travis Wilson
Time: Tue Jun 17 11:35:35 2014
Hi Randy,
I appreciate the help. I uploaded everything to
/incoming/irap/met_help/wilson_data so you can reproduce the problem
if need be. Thanks!
Travis Wilson
UCLA Atmospheric and Oceanic Sciences
model viewer: http://web.atmos.ucla.edu/~wilson28/pages/viewer.html
-----Original Message-----
From: Randy Bullock via RT [mailto:met_help at ucar.edu]
Sent: Tuesday, June 17, 2014 10:02 AM
To: Wilson0028 at ucla.edu
Subject: Re: [rt.rap.ucar.edu #67713] interpolation using constant
height fields
Travis -
I had a look at the code and also a few old met_help questions in
order to get a handle on this. As far as I can tell, interpolation
should be happening. The only times when vertical interpolation
doesn't happen is when you're dealing with surface fields, but it
looks from your config file that that's not the case here.
I'll talk this over with John when gets back from vacation next week,
and see if he has any insights on this.
Randy
On Mon, Jun 16, 2014 at 8:07 PM, Travis Wilson via RT
<met_help at ucar.edu>
wrote:
>
> Mon Jun 16 20:07:20 2014: Request 67713 was acted upon.
> Transaction: Ticket created by Wilson0028 at ucla.edu
> Queue: met_help
> Subject: interpolation using constant height fields
> Owner: Nobody
> Requestors: Wilson0028 at ucla.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=67713
> >
>
>
> Hi John,
> When I use temperature (on constant pressure fields made by
> UPP) to verify the vertical profile, the foretasted field gets
interpolated to
> the observation height. Everything looks very smooth (see
> constant_pressure_field.png).
>
> Now when I use temperature (on constant height fields made
by
> UPP), my foretasted field does not get interpolated (see
> constant_height_field.png).
>
> I am assuming this has something to do with the fact that the fields
> with constant height are in the same grib file that has fields of
> constant pressure (and the native coordinate in the UPP grib file is
> pressure, not sure how to make it height). my wgrib is below. Any
> idea on how to interpolate the forecast fields?
>
> Thanks,
> Travis
>
> (note: I am using VTMP)
>
30:569006:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=50:TR=0:P1=0:P2=0:T
> imeU=1:50
> mb:anl:NAve=0
>
46:853000:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=100:TR=0:P1=0:P2=0:
> TimeU=1:100
> mb:anl:NAve=0
>
61:1130366:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=150:TR=0:P1=0:P2=0
> :TimeU=1:150
> mb:anl:NAve=0
>
76:1427466:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=200:TR=0:P1=0:P2=0
> :TimeU=1:200
> mb:anl:NAve=0
>
91:1737568:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=250:TR=0:P1=0:P2=0
> :TimeU=1:250
> mb:anl:NAve=0
>
106:2041070:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=300:TR=0:P1=0:P2=
> 0:TimeU=1:300
> mb:anl:NAve=0
>
121:2357640:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=350:TR=0:P1=0:P2=
> 0:TimeU=1:350
> mb:anl:NAve=0
>
136:2667676:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=400:TR=0:P1=0:P2=
> 0:TimeU=1:400
> mb:anl:NAve=0
>
151:2975534:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=450:TR=0:P1=0:P2=
> 0:TimeU=1:450
> mb:anl:NAve=0
>
166:3276924:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=500:TR=0:P1=0:P2=
> 0:TimeU=1:500
> mb:anl:NAve=0
>
181:3573958:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=550:TR=0:P1=0:P2=
> 0:TimeU=1:550
> mb:anl:NAve=0
>
196:3868814:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=600:TR=0:P1=0:P2=
> 0:TimeU=1:600
> mb:anl:NAve=0
>
211:4165848:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=650:TR=0:P1=0:P2=
> 0:TimeU=1:650
> mb:anl:NAve=0
>
226:4465060:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=700:TR=0:P1=0:P2=
> 0:TimeU=1:700
> mb:anl:NAve=0
>
241:4903598:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=750:TR=0:P1=0:P2=
> 0:TimeU=1:750
> mb:anl:NAve=0
>
256:5344314:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=775:TR=0:P1=0:P2=
> 0:TimeU=1:775
> mb:anl:NAve=0
>
271:5778496:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=800:TR=0:P1=0:P2=
> 0:TimeU=1:800
> mb:anl:NAve=0
>
286:6210500:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=825:TR=0:P1=0:P2=
> 0:TimeU=1:825
> mb:anl:NAve=0
>
301:6642504:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=850:TR=0:P1=0:P2=
> 0:TimeU=1:850
> mb:anl:NAve=0
>
316:7074508:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=875:TR=0:P1=0:P2=
> 0:TimeU=1:875
> mb:anl:NAve=0
>
331:7506512:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=900:TR=0:P1=0:P2=
> 0:TimeU=1:900
> mb:anl:NAve=0
>
346:7945050:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=910:TR=0:P1=0:P2=
> 0:TimeU=1:910
> mb:anl:NAve=0
>
361:8379232:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=920:TR=0:P1=0:P2=
> 0:TimeU=1:920
> mb:anl:NAve=0
>
376:8813414:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=930:TR=0:P1=0:P2=
> 0:TimeU=1:930
> mb:anl:NAve=0
>
391:9247596:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=940:TR=0:P1=0:P2=
> 0:TimeU=1:940
> mb:anl:NAve=0
>
406:9686134:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=950:TR=0:P1=0:P2=
> 0:TimeU=1:950
> mb:anl:NAve=0
>
421:10126850:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=960:TR=0:P1=0:P2
> =0:TimeU=1:960
> mb:anl:NAve=0
>
436:10565388:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=970:TR=0:P1=0:P2
> =0:TimeU=1:970
> mb:anl:NAve=0
>
451:11003926:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=980:TR=0:P1=0:P2
> =0:TimeU=1:980
> mb:anl:NAve=0
>
466:11444642:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=990:TR=0:P1=0:P2
> =0:TimeU=1:990
> mb:anl:NAve=0
>
480:11852594:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=1000:TR=0:P1=0:P
> 2=0:TimeU=1:1000
> mb:anl:NAve=0
>
578:13222406:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=30:TR=0:P1=0:P2=
> 0:TimeU=1:30
> m above MSL:anl:NAve=0
>
579:13240112:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=80:TR=0:P1=0:P2=
> 0:TimeU=1:80
> m above MSL:anl:NAve=0
>
580:13258406:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=140:TR=0:P1=0:P2
> =0:TimeU=1:140
> m above MSL:anl:NAve=0
>
581:13277194:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=200:TR=0:P1=0:P2
> =0:TimeU=1:200
> m above MSL:anl:NAve=0
>
582:13296424:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=300:TR=0:P1=0:P2
> =0:TimeU=1:300
> m above MSL:anl:NAve=0
>
583:13316380:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=400:TR=0:P1=0:P2
> =0:TimeU=1:400
> m above MSL:anl:NAve=0
>
584:13336962:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=500:TR=0:P1=0:P2
> =0:TimeU=1:500
> m above MSL:anl:NAve=0
>
585:13358144:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=600:TR=0:P1=0:P2
> =0:TimeU=1:600
> m above MSL:anl:NAve=0
>
586:13380028:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=700:TR=0:P1=0:P2
> =0:TimeU=1:700
> m above MSL:anl:NAve=0
>
587:13402522:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=800:TR=0:P1=0:P2
> =0:TimeU=1:800
> m above MSL:anl:NAve=0
>
588:13425618:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=900:TR=0:P1=0:P2
> =0:TimeU=1:900
> m above MSL:anl:NAve=0
>
589:13449322:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1000:TR=0:P1=0:P
> 2=0:TimeU=1:1000
> m above MSL:anl:NAve=0
>
590:13473584:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1200:TR=0:P1=0:P
> 2=0:TimeU=1:1200
> m above MSL:anl:NAve=0
>
591:13498886:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1400:TR=0:P1=0:P
> 2=0:TimeU=1:1400
> m above MSL:anl:NAve=0
>
592:13526228:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1600:TR=0:P1=0:P
> 2=0:TimeU=1:1600
> m above MSL:anl:NAve=0
>
593:13556146:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1800:TR=0:P1=0:P
> 2=0:TimeU=1:1800
> m above MSL:anl:NAve=0
>
>
------------------------------------------------
Subject: interpolation using constant height fields
From: John Halley Gotway
Time: Mon Jun 23 14:58:08 2014
Travis,
Let me start by telling you how Point-Stat does vertical
interpolation.
For a range of pressure levels (specified with a letter P, like P250-
350),
the vertical interpolation is done linear in the log of pressure. For
a
range of any other level type (like Z250-350 for height), it's done
linearly by height.
I took a look at the data you sent and ran the same case for which you
sent
me output. I'm reproducing the numbers you're getting identically.
To
narrow things down, I reran for just one range:
level = [ "12/P250-350", "12/Z250-350" ]
Here I'm saying to verify virtual temperature between 250 and 350
millibars
(P250-350) and verify it between 250 and 350 meters above ground
(Z250-350). Running point_stat at verbosity level 2, I see the
following
output:
DEBUG 2: For VTMP/P350-250 found 3 forecast levels and 0
climatology
levels.
...
DEBUG 2: For VTMP/Z350-250 found 1 forecast levels and 0
climatology
levels.
For the range P250-350, it found data for pressure levels 250, 300,
and
350. But for the range Z250-350, it only found data for a single
level 300
meters above ground.
When doing vertical interpolation on pressure levels, since we found 3
forecast levels, there's always a pressure level above/below each
observation. So the vertical interpolation can be done smoothly.
But when doing vertical interpolation on height levels, there's only a
single level in that range. So all the observations are matched to
that
single level.
I think the problem here is a poor choice of levels for verification.
Based on your settings, it looks like you want to verify using
observations
in bands of 100 meters. I'd suggest setting up the level for the
forecast
and observation fields a little differently. You could try this:
fcst = {
wind_thresh = [ NA ];
field = [
{
name = "VTMP";
level = [ "Z200-400" ];
cat_thresh = [];
}
];
};
obs = {
field = [
{
name = "VTMP";
level = [ "Z250-350" ];
cat_thresh = [];
}
];
};
This tells point_stat to extract levels 200, 300, and 400 m from the
forecast file and use it for this verification task. But only use
observations falling between 250 and 350 m. That will give point_stat
forecast level data to use for the vertical interpolation above/below
each
observation value. But it'll only use the obs falling in your
requested
range. And that should enable the vertical interpolation to be done
more
smoothly.
I do realize this is confusing. Hopefully that helps clarify.
Thanks,
John
On Tue, Jun 17, 2014 at 11:35 AM, Travis Wilson via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=67713 >
>
> Hi Randy,
> I appreciate the help. I uploaded everything to
> /incoming/irap/met_help/wilson_data so you can reproduce the problem
if
> need be. Thanks!
>
> Travis Wilson
> UCLA Atmospheric and Oceanic Sciences
> model viewer: http://web.atmos.ucla.edu/~wilson28/pages/viewer.html
>
>
> -----Original Message-----
> From: Randy Bullock via RT [mailto:met_help at ucar.edu]
> Sent: Tuesday, June 17, 2014 10:02 AM
> To: Wilson0028 at ucla.edu
> Subject: Re: [rt.rap.ucar.edu #67713] interpolation using constant
height
> fields
>
> Travis -
>
> I had a look at the code and also a few old met_help questions in
order to
> get a handle on this. As far as I can tell, interpolation should be
> happening. The only times when vertical interpolation doesn't
happen is
> when you're dealing with surface fields, but it looks from your
config file
> that that's not the case here.
>
> I'll talk this over with John when gets back from vacation next
week, and
> see if he has any insights on this.
>
> Randy
>
>
> On Mon, Jun 16, 2014 at 8:07 PM, Travis Wilson via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Mon Jun 16 20:07:20 2014: Request 67713 was acted upon.
> > Transaction: Ticket created by Wilson0028 at ucla.edu
> > Queue: met_help
> > Subject: interpolation using constant height fields
> > Owner: Nobody
> > Requestors: Wilson0028 at ucla.edu
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=67713
> > >
> >
> >
> > Hi John,
> > When I use temperature (on constant pressure fields made
by
> > UPP) to verify the vertical profile, the foretasted field gets
> interpolated to
> > the observation height. Everything looks very smooth (see
> > constant_pressure_field.png).
> >
> > Now when I use temperature (on constant height fields made
by
> > UPP), my foretasted field does not get interpolated (see
> > constant_height_field.png).
> >
> > I am assuming this has something to do with the fact that the
fields
> > with constant height are in the same grib file that has fields of
> > constant pressure (and the native coordinate in the UPP grib file
is
> > pressure, not sure how to make it height). my wgrib is below.
Any
> > idea on how to interpolate the forecast fields?
> >
> > Thanks,
> > Travis
> >
> > (note: I am using VTMP)
> >
30:569006:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=50:TR=0:P1=0:P2=0:T
> > imeU=1:50
> > mb:anl:NAve=0
> >
46:853000:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=100:TR=0:P1=0:P2=0:
> > TimeU=1:100
> > mb:anl:NAve=0
> >
61:1130366:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=150:TR=0:P1=0:P2=0
> > :TimeU=1:150
> > mb:anl:NAve=0
> >
76:1427466:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=200:TR=0:P1=0:P2=0
> > :TimeU=1:200
> > mb:anl:NAve=0
> >
91:1737568:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=250:TR=0:P1=0:P2=0
> > :TimeU=1:250
> > mb:anl:NAve=0
> >
106:2041070:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=300:TR=0:P1=0:P2=
> > 0:TimeU=1:300
> > mb:anl:NAve=0
> >
121:2357640:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=350:TR=0:P1=0:P2=
> > 0:TimeU=1:350
> > mb:anl:NAve=0
> >
136:2667676:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=400:TR=0:P1=0:P2=
> > 0:TimeU=1:400
> > mb:anl:NAve=0
> >
151:2975534:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=450:TR=0:P1=0:P2=
> > 0:TimeU=1:450
> > mb:anl:NAve=0
> >
166:3276924:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=500:TR=0:P1=0:P2=
> > 0:TimeU=1:500
> > mb:anl:NAve=0
> >
181:3573958:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=550:TR=0:P1=0:P2=
> > 0:TimeU=1:550
> > mb:anl:NAve=0
> >
196:3868814:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=600:TR=0:P1=0:P2=
> > 0:TimeU=1:600
> > mb:anl:NAve=0
> >
211:4165848:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=650:TR=0:P1=0:P2=
> > 0:TimeU=1:650
> > mb:anl:NAve=0
> >
226:4465060:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=700:TR=0:P1=0:P2=
> > 0:TimeU=1:700
> > mb:anl:NAve=0
> >
241:4903598:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=750:TR=0:P1=0:P2=
> > 0:TimeU=1:750
> > mb:anl:NAve=0
> >
256:5344314:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=775:TR=0:P1=0:P2=
> > 0:TimeU=1:775
> > mb:anl:NAve=0
> >
271:5778496:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=800:TR=0:P1=0:P2=
> > 0:TimeU=1:800
> > mb:anl:NAve=0
> >
286:6210500:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=825:TR=0:P1=0:P2=
> > 0:TimeU=1:825
> > mb:anl:NAve=0
> >
301:6642504:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=850:TR=0:P1=0:P2=
> > 0:TimeU=1:850
> > mb:anl:NAve=0
> >
316:7074508:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=875:TR=0:P1=0:P2=
> > 0:TimeU=1:875
> > mb:anl:NAve=0
> >
331:7506512:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=900:TR=0:P1=0:P2=
> > 0:TimeU=1:900
> > mb:anl:NAve=0
> >
346:7945050:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=910:TR=0:P1=0:P2=
> > 0:TimeU=1:910
> > mb:anl:NAve=0
> >
361:8379232:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=920:TR=0:P1=0:P2=
> > 0:TimeU=1:920
> > mb:anl:NAve=0
> >
376:8813414:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=930:TR=0:P1=0:P2=
> > 0:TimeU=1:930
> > mb:anl:NAve=0
> >
391:9247596:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=940:TR=0:P1=0:P2=
> > 0:TimeU=1:940
> > mb:anl:NAve=0
> >
406:9686134:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=950:TR=0:P1=0:P2=
> > 0:TimeU=1:950
> > mb:anl:NAve=0
> >
421:10126850:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=960:TR=0:P1=0:P2
> > =0:TimeU=1:960
> > mb:anl:NAve=0
> >
436:10565388:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=970:TR=0:P1=0:P2
> > =0:TimeU=1:970
> > mb:anl:NAve=0
> >
451:11003926:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=980:TR=0:P1=0:P2
> > =0:TimeU=1:980
> > mb:anl:NAve=0
> >
466:11444642:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=990:TR=0:P1=0:P2
> > =0:TimeU=1:990
> > mb:anl:NAve=0
> >
480:11852594:d=11010400:VTMP:kpds5=12:kpds6=100:kpds7=1000:TR=0:P1=0:P
> > 2=0:TimeU=1:1000
> > mb:anl:NAve=0
> >
578:13222406:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=30:TR=0:P1=0:P2=
> > 0:TimeU=1:30
> > m above MSL:anl:NAve=0
> >
579:13240112:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=80:TR=0:P1=0:P2=
> > 0:TimeU=1:80
> > m above MSL:anl:NAve=0
> >
580:13258406:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=140:TR=0:P1=0:P2
> > =0:TimeU=1:140
> > m above MSL:anl:NAve=0
> >
581:13277194:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=200:TR=0:P1=0:P2
> > =0:TimeU=1:200
> > m above MSL:anl:NAve=0
> >
582:13296424:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=300:TR=0:P1=0:P2
> > =0:TimeU=1:300
> > m above MSL:anl:NAve=0
> >
583:13316380:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=400:TR=0:P1=0:P2
> > =0:TimeU=1:400
> > m above MSL:anl:NAve=0
> >
584:13336962:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=500:TR=0:P1=0:P2
> > =0:TimeU=1:500
> > m above MSL:anl:NAve=0
> >
585:13358144:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=600:TR=0:P1=0:P2
> > =0:TimeU=1:600
> > m above MSL:anl:NAve=0
> >
586:13380028:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=700:TR=0:P1=0:P2
> > =0:TimeU=1:700
> > m above MSL:anl:NAve=0
> >
587:13402522:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=800:TR=0:P1=0:P2
> > =0:TimeU=1:800
> > m above MSL:anl:NAve=0
> >
588:13425618:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=900:TR=0:P1=0:P2
> > =0:TimeU=1:900
> > m above MSL:anl:NAve=0
> >
589:13449322:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1000:TR=0:P1=0:P
> > 2=0:TimeU=1:1000
> > m above MSL:anl:NAve=0
> >
590:13473584:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1200:TR=0:P1=0:P
> > 2=0:TimeU=1:1200
> > m above MSL:anl:NAve=0
> >
591:13498886:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1400:TR=0:P1=0:P
> > 2=0:TimeU=1:1400
> > m above MSL:anl:NAve=0
> >
592:13526228:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1600:TR=0:P1=0:P
> > 2=0:TimeU=1:1600
> > m above MSL:anl:NAve=0
> >
593:13556146:d=11010400:VTMP:kpds5=12:kpds6=103:kpds7=1800:TR=0:P1=0:P
> > 2=0:TimeU=1:1800
> > m above MSL:anl:NAve=0
> >
> >
>
>
>
>
>
------------------------------------------------
More information about the Met_help
mailing list