[Met_help] PBL retrieval
Joe Eastman
Joe.Eastman at windlogics.com
Fri Mar 5 11:00:31 MST 2010
Hi John,
One more nuance I thought of. Does the data have to be in temporal sequential order? There are days when multiple flights are occurring at the same time. I guess I could write something to sort and combine these but it would be very easy if that was not necessary,
Take care,
Joe
Dr. Joseph L. Eastman
Senior Atmospheric Scientist
WindLogics Inc.
201 4th St NW
Grand Rapids, MN
55744
c - 410.279.9702
w - 651.556.4297
jeastman at windlogics.com
-----Original Message-----
From: John Halley Gotway [mailto:johnhg at ucar.edu]
Sent: Friday, March 05, 2010 10:48 AM
To: Joe Eastman
Cc: met_help at ucar.edu
Subject: Re: [Met_help] PBL retrieval
Joe,
OK, so you have observations of the "BOUNDARY LAYER DEPTH RETRIEVAL". What sort of comparison are you trying to do? Are you trying to compare it to GRIB code 221 (Planetary boundary layer height,
HPBL)? The units of HPBL are meters, and I'm guessing they correspond to above ground level (AGL), but I'm really not sure. You should figure out if you want to be comparing to AGL or MSL. I'm
really not sure.
So assuming you want to compare your forecast of HPBL to your lidar observations AGL, here's what I'd suggest doing...
- Reformat these ASCII files into an ASCII format for use by the ASCII2NC tool. That format consists of 10 columns:
(1) Message_Type - Set to "ADPSFC" to pretend that these are observations taken at the surface.
(2) Station_ID - Just put a dummy string in here. Doesn't really matter.
(3) Valid_Time - Put in the observation time in YYYYMMDD_HHMMSS format.
(4) Lat(Deg North) - Use the value in the LATITUDE column.
(5) Lon(Deg East) - Use the value in the LONGITUDE column.
(6) Elevation(msl) - Just put "-9999" since this isn't actually used in Point-Stat.
(7) Grib_Code - Set to 221 to mean that this observation should compared to Planetary boundary layer height forecasts.
(8) Level - Just put "-9999" since we don't have a pressure value for this observation. We're just pretending they were taken at the surface.
(9) Height(msl) - Just put "-9999" since this isn't actually used in Point-Stat.
(10) Observation_Value - Use the value in the ZI_AGL column. Be sure to use "-9999" for missing data.
I'd try reformatting the ASCII as described above, run it though ASCII2NC, and then use Point-Stat to compare your forecast of HPBL to these observations. So in the Point-Stat config file, set:
fcst_field[] = [ "HPBL/Z0" ];
message_type[] = [ "ADPSFC" ];
Give that a shot and let us know if that works for you.
Thanks,
John
Joe Eastman wrote:
> Hello,
>
> What is the best way to create an obs input file for this profiler PBL retrieval data? It is not a fixed station so I am not sure how to attack this.
>
> 198
> BOUNDARY LAYER DEPTH RETRIEVAL
> IHOP 2002 EXPERIMENT
> File: DLR0519leg01azi.txt
> Created: Tue Apr 12 11:54:16 EDT 2005
> Platform: DLR Lidar
> Date: 05/19/2002
> Time: 16:40:54 - 17:06:30
> Number of Data Points: 198
> Average Aircraft Speed (m/s): 142.279
> Starting Aircraft Altitude (m MSL): 4566.41
> Ending Aircraft Altitude (m MSL): 4851.70
> PI: Dr. Ken Davis (davis at met.psu.edu). Contact: Ken Craig (kcraig at met.psu.edu).
> Department of Meteorology, The Pennsylvania State University
> 503 Walker Building, University Park, PA 16802-5013
>
> TIME LATITUDE LONGITUDE ZI_MSL ZI_AGL
> 16:40:54 35.49359 -100.59423 -999.00 -999.00
> 16:41:00 35.49520 -100.60099 -999.00 -999.00
> 16:41:06 35.49749 -100.60725 -999.00 -999.00
> 16:41:12 35.50061 -100.61314 -999.00 -999.00
> 16:41:18 35.50465 -100.61868 -999.00 -999.00
> 16:41:24 35.50946 -100.62359 -999.00 -999.00
> 16:41:30 35.51509 -100.62767 -999.00 -999.00
> 16:41:36 35.52119 -100.63056 -999.00 -999.00
> 16:41:42 35.52773 -100.63193 -999.00 -999.00
> 16:41:48 35.53448 -100.63148 -999.00 -999.00
> 16:41:54 35.54134 -100.62993 -999.00 -999.00
> 16:42:00 35.54831 -100.62796 1755.00 840.00
> 16:42:06 35.55552 -100.62576 1740.00 825.00
> 16:42:18 35.57038 -100.62107 1755.00 840.00
> 16:42:24 35.57814 -100.61858 1755.00 825.00
> 16:42:30 35.58604 -100.61600 1740.00 810.00
> 16:42:36 35.59404 -100.61341 1755.00 810.00
> 16:42:42 35.60226 -100.61072 1740.00 795.00
> 16:42:48 35.61045 -100.60802 1755.00 810.00
> 16:42:54 35.61875 -100.60521 1770.00 825.00
> 16:43:00 35.62681 -100.60246 1740.00 795.00
> 16:43:06 35.63487 -100.59979 1695.00 765.00
> 16:43:12 35.64290 -100.59781 -999.00 -999.00
> 16:43:18 35.65094 -100.59665 -999.00 -999.00
> 16:43:24 35.65900 -100.59617 1695.00 780.00
> 16:43:30 35.66690 -100.59621 1680.00 765.00
> 16:43:36 35.67475 -100.59667 1650.00 735.00
> 16:43:42 35.68251 -100.59745 1695.00 795.00
> 16:43:54 35.69807 -100.59966 1725.00 825.00
> 16:44:00 35.70585 -100.60093 1680.00 825.00
> 16:44:06 35.71367 -100.60224 1725.00 870.00
> 16:44:12 35.72134 -100.60350 1740.00 900.00
> 16:44:18 35.72899 -100.60470 1740.00 915.00
> 16:44:24 35.73666 -100.60580 1740.00 885.00
> 16:44:30 35.74437 -100.60681 1740.00 900.00
> 16:44:36 35.75196 -100.60769 1665.00 810.00
> 16:44:42 35.75962 -100.60845 1710.00 855.00
> 16:44:48 35.76728 -100.60909 1785.00 915.00
> 16:44:54 35.77497 -100.60958 1815.00 945.00
> 16:45:00 35.78257 -100.60995 1755.00 870.00
> 16:45:06 35.79018 -100.61021 1680.00 780.00
> 16:45:12 35.79769 -100.61038 1725.00 810.00
> 16:45:18 35.80531 -100.61048 1680.00 780.00
> 16:45:30 35.82038 -100.61048 1635.00 750.00
> 16:45:36 35.82796 -100.61040 1650.00 795.00
> 16:45:42 35.83546 -100.61028 1620.00 750.00
> 16:45:48 35.84317 -100.61012 1650.00 780.00
> 16:45:54 35.85064 -100.60995 1680.00 825.00
> 16:46:00 35.85815 -100.60976 1695.00 840.00
> 16:46:06 35.86560 -100.60958 1650.00 810.00
> 16:46:12 35.87310 -100.60940 1665.00 795.00
> 16:46:18 35.88054 -100.60923 1680.00 825.00
> 16:46:24 35.88811 -100.60909 1620.00 750.00
> 16:46:30 35.89556 -100.60895 1620.00 780.00
> 16:46:36 35.90316 -100.60885 1650.00 840.00
> Etc
>
>
> Take care,
> Joe
>
> Dr. Joseph L. Eastman
> Senior Atmospheric Scientist
> WindLogics Inc.
> 201 4th St NW
> Grand Rapids, MN
> 55744
>
> c - 410.279.9702
> w - 651.556.4297
> jeastman at windlogics.com<mailto:jeastman at windlogics.com>
>
>
>
> ________________________________
> CONFIDENTIALITY NOTICE: This e-mail message is for the sole use of the intended recipient(s) and may contain confidential and privileged information. Any unauthorized review, use, disclosure or distribution of any kind is strictly prohibited. If you are not the intended recipient, please contact the sender via reply e-mail and destroy all copies of the original message. Thank you.
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> Met_help mailing list
> Met_help at mailman.ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/met_help
More information about the Met_help
mailing list