[Met_help] [rt.rap.ucar.edu #67410] History for Wind direction statistics by station
John Halley Gotway via RT
met_help at ucar.edu
Tue Jun 3 11:20:05 MDT 2014
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear MET help team
hope you are doing fine at your end
I am Victor Almanza and am using MET to evaluate WRF model runs.
I would like to ask for your help with the following questions:
1) is it possible to obtain aggregated statistics of surface wind direction (10 m) by station with "stat_analysis".
I am interested in the model performance for wind direction at several monitoring stations. Specifically, to obtain aggregated "AGGR_WDIR" lines at each station over a defined time period (day,week, month) similar to when running for MPR lines (Temp for instance) with the "-by" flag.
I don't have it clear if this is possible, or if am overlooking something obvious. If so, could you provide please any hints/advice on that ?
2) I am interested in obtaining the BIAS of wind direction.
Is it possible to do it with MET ?.. if not, could you provide directions of which subroutines to look at in the code ?
Thank you in advance for your time
Best
Victor
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #67410] Wind direction statistics by station
From: John Halley Gotway
Time: Mon Jun 02 10:23:40 2014
Victor,
Unfortunately, Stat-Analysis does not currently have the capability of
computing wind direction errors from the MPR line type.
It can only compute wind direction errors from the VL1L2 line type.
There is a way to do the analysis you're after, but it's a bit clumsy.
When you run Point-Stat, you'd have to define each of the
stations in which you're interested as a separate masking region.
Then you'd end up with 1 VL1L2 line type for each station. It'd
basically just have 1 point in it. The name of the masking region
would just be the name of the station. Then you could run the
following type of Stat-Analysis job:
-job aggregate_stat -line_type VL1L2 -out_line_type WDIR -by
VX_MASK
That'd compute wind direction error for each unique entry in the
VX_MASK column - which would be the stations for which you've computed
VL1L2 output lines.
That would do it, but I realize that this is clunky to have to define
separate masking regions for each station of interest.
It occurs to me that we probably could create a new Stat-Analysis job
type that would convert MPR lines to wind direction errors. It would
search for pairs of UGRD and VGRD matched pairs and use them
to compute wind direction errors. I'll document this idea as an area
for future development, but for now you're stuck with the VL1L2 to
WDIR conversion.
Hope that helps clarify.
As for your question about bias, here's the info I received from our
resident statistician...
The aggregated wdir mean error is the bias. It will be most
meaningful when conditional biases do not cancel each other out, so
small areas (or single stations) with similar prevailing wind
directions (e.g. don't mix seasons, day and night, etc. if they are
very different).
Thanks,
John
On 05/30/2014 08:03 PM, Victor Almanza via RT wrote:
>
> Fri May 30 20:03:10 2014: Request 67410 was acted upon.
> Transaction: Ticket created by halmanset at yahoo.com
> Queue: met_help
> Subject: Wind direction statistics by station
> Owner: Nobody
> Requestors: halmanset at yahoo.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=67410 >
>
>
> Dear MET help team
>
> hope you are doing fine at your end
>
> I am Victor Almanza and am using MET to evaluate WRF model runs.
>
> I would like to ask for your help with the following questions:
>
> 1) is it possible to obtain aggregated statistics of surface wind
direction (10 m) by station with "stat_analysis".
>
> I am interested in the model performance for wind direction at
several monitoring stations. Specifically, to obtain aggregated
"AGGR_WDIR" lines at each station over a defined time period
(day,week, month) similar to when running for MPR lines (Temp for
instance) with the "-by" flag.
>
> I don't have it clear if this is possible, or if am overlooking
something obvious. If so, could you provide please any hints/advice
on that ?
>
>
> 2) I am interested in obtaining the BIAS of wind direction.
>
> Is it possible to do it with MET ?.. if not, could you provide
directions of which subroutines to look at in the code ?
>
>
> Thank you in advance for your time
> Best
>
> Victor
>
------------------------------------------------
Subject: Wind direction statistics by station
From: Victor Almanza
Time: Mon Jun 02 11:13:18 2014
Dear John,
It did clarify my questions.
That hint is helpful indeed.
I have a small region in complex terrain with rather few stations. I
think this helps in evaluating WRF results by station (T, Wind Speed,
RH, Wind Direction).
Thank you very much for your time
Best !
Victor
________________________________
De: John Halley Gotway via RT <met_help at ucar.edu>
Para: halmanset at yahoo.com
Enviado: Lunes, 2 de junio, 2014 9:23 A.M.
Asunto: Re: [rt.rap.ucar.edu #67410] Wind direction statistics by
station
Victor,
Unfortunately, Stat-Analysis does not currently have the capability of
computing wind direction errors from the MPR line type.
It can only compute wind direction errors from the VL1L2 line type.
There is a way to do the analysis you're after, but it's a bit
clumsy. When you run Point-Stat, you'd have to define each of the
stations in which you're interested as a separate masking region.
Then you'd end up with 1 VL1L2 line type for each station. It'd
basically just have 1 point in it. The name of the masking region
would just be the name of the station. Then you could run the
following type of Stat-Analysis job:
-job aggregate_stat -line_type VL1L2 -out_line_type WDIR -by
VX_MASK
That'd compute wind direction error for each unique entry in the
VX_MASK column - which would be the stations for which you've computed
VL1L2 output lines.
That would do it, but I realize that this is clunky to have to define
separate masking regions for each station of interest.
It occurs to me that we probably could create a new Stat-Analysis job
type that would convert MPR lines to wind direction errors. It would
search for pairs of UGRD and VGRD matched pairs and use them
to compute wind direction errors. I'll document this idea as an area
for future development, but for now you're stuck with the VL1L2 to
WDIR conversion.
Hope that helps clarify.
As for your question about bias, here's the info I received from our
resident statistician...
The aggregated wdir mean error is the bias. It will be most
meaningful when conditional biases do not cancel each other out, so
small areas (or single stations) with similar prevailing wind
directions (e.g. don't mix seasons, day and night, etc. if they are
very different).
Thanks,
John
On 05/30/2014 08:03 PM, Victor Almanza via RT wrote:
>
> Fri May 30 20:03:10 2014: Request 67410 was acted upon.
> Transaction: Ticket created by halmanset at yahoo.com
> Queue: met_help
> Subject: Wind direction statistics by station
> Owner: Nobody
> Requestors: halmanset at yahoo.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=67410 >
>
>
> Dear MET help team
>
> hope you are doing fine at your end
>
> I am Victor Almanza and am using MET to evaluate WRF model runs.
>
> I would like to ask for your help with the following questions:
>
> 1) is it possible to obtain aggregated statistics of surface wind
direction (10 m) by station with "stat_analysis".
>
> I am interested in the model performance for wind direction at
several monitoring stations. Specifically, to obtain aggregated
"AGGR_WDIR" lines at each station over a defined time period
(day,week, month) similar to when running for MPR lines (Temp for
instance) with the "-by" flag.
>
> I don't have it clear if this is possible, or if am overlooking
something obvious. If so, could you provide please any hints/advice
on that ?
>
>
> 2) I am interested in obtaining the BIAS of wind direction.
>
> Is it possible to do it with MET ?.. if not, could you provide
directions of which subroutines to look at in the code ?
>
>
> Thank you in advance for your time
> Best
>
> Victor
>
------------------------------------------------
Subject: Wind direction statistics by station
From: Victor Almanza
Time: Mon Jun 02 21:16:12 2014
Dear John
I followed your advice and think am having the WDIR by station.
I just want to make sure it is being doing right:
- I defined in the pointStat configuration file the stations ID (for
now it is just a test, so I didn't include all of them) am interested
in, within the "sid" field: a comma-separated array of stations ID
strings
- then ran pointStat for a small period...
- then ran statAnalysis
I obtained the following "out" file.
Thank you in advance for your time and help
Best
Victor
JOB_LIST: -job aggregate_stat -fcst_lev Z10 -interp_mthd UW_MEAN
-line_type VL1L2 -by VX_MASK -dump_row out/statAnalysis/aggr-stat_job-
mil-wdir-by-stn.stat -out_line_type WDIR
COL_NAME: VX_MASK TOTAL FBAR OBAR ME MAE
ROW_MEAN_WDIR: 76003 24 161.75494 175.26312 -10.22146 21.24104
AGGR_WDIR: 76003 24 175.05710 184.59730 -9.54019 NA
ROW_MEAN_WDIR: 76004 30 163.43686 211.67779 -50.11654 50.11654
AGGR_WDIR: 76004 30 176.53086 219.06307 -42.53221 NA
ROW_MEAN_WDIR: 76006 31 155.83850 173.87529 -17.68601 18.40499
AGGR_WDIR: 76006 31 165.33436 177.76244 -12.42808 NA
ROW_MEAN_WDIR: 76007 33 177.74966 186.75497 -10.33529 16.57426
AGGR_WDIR: 76007 33 179.14361 183.98678 -4.84317 NA
ROW_MEAN_WDIR: 76008 27 153.67546 154.30870 8.96001 19.77858
AGGR_WDIR: 76008 27 163.05044 158.08629 4.96415 NA
ROW_MEAN_WDIR: 76011 19 172.47562 195.93538 14.50380 42.45971
AGGR_WDIR: 76011 19 173.52899 161.25917 12.26982 NA
________________________________
De: John Halley Gotway via RT <met_help at ucar.edu>
Para: halmanset at yahoo.com
Enviado: Lunes, 2 de junio, 2014 10:13 A.M.
Asunto: [rt.rap.ucar.edu #67410] Resolved: Wind direction statistics
by station
According to our records, your request has been resolved. If you have
any
further questions or concerns, please respond to this message.
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #67410] Resolved: Wind direction statistics by station
From: John Halley Gotway
Time: Tue Jun 03 09:43:05 2014
Victor,
Yes, that all look correct to me. I'll go ahead and resolve this
ticket. If more questions arise, feel free to write again.
Thanks,
John
On 06/02/2014 09:16 PM, Victor Almanza via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=67410 >
>
> Dear John
>
> I followed your advice and think am having the WDIR by station.
>
> I just want to make sure it is being doing right:
>
> - I defined in the pointStat configuration file the stations ID (for
now it is just a test, so I didn't include all of them) am interested
in, within the "sid" field: a comma-separated array of stations ID
strings
> - then ran pointStat for a small period...
> - then ran statAnalysis
>
>
> I obtained the following "out" file.
>
>
> Thank you in advance for your time and help
>
> Best
>
> Victor
>
>
> JOB_LIST: -job aggregate_stat -fcst_lev Z10 -interp_mthd
UW_MEAN -line_type VL1L2 -by VX_MASK -dump_row out/statAnalysis/aggr-
stat_job-mil-wdir-by-stn.stat -out_line_type WDIR
> COL_NAME: VX_MASK TOTAL FBAR OBAR ME MAE
> ROW_MEAN_WDIR: 76003 24 161.75494 175.26312 -10.22146 21.24104
> AGGR_WDIR: 76003 24 175.05710 184.59730 -9.54019 NA
> ROW_MEAN_WDIR: 76004 30 163.43686 211.67779 -50.11654 50.11654
> AGGR_WDIR: 76004 30 176.53086 219.06307 -42.53221 NA
> ROW_MEAN_WDIR: 76006 31 155.83850 173.87529 -17.68601 18.40499
> AGGR_WDIR: 76006 31 165.33436 177.76244 -12.42808 NA
> ROW_MEAN_WDIR: 76007 33 177.74966 186.75497 -10.33529 16.57426
> AGGR_WDIR: 76007 33 179.14361 183.98678 -4.84317 NA
> ROW_MEAN_WDIR: 76008 27 153.67546 154.30870 8.96001 19.77858
> AGGR_WDIR: 76008 27 163.05044 158.08629 4.96415 NA
> ROW_MEAN_WDIR: 76011 19 172.47562 195.93538 14.50380 42.45971
> AGGR_WDIR: 76011 19 173.52899 161.25917 12.26982 NA
>
>
>
> ________________________________
> De: John Halley Gotway via RT <met_help at ucar.edu>
> Para: halmanset at yahoo.com
> Enviado: Lunes, 2 de junio, 2014 10:13 A.M.
> Asunto: [rt.rap.ucar.edu #67410] Resolved: Wind direction statistics
by station
>
>
> According to our records, your request has been resolved. If you
have any
> further questions or concerns, please respond to this message.
>
------------------------------------------------
Subject: Wind direction statistics by station
From: Victor Almanza
Time: Tue Jun 03 10:58:43 2014
Dear John
Once again, thank you very much for your hints/guidance on this
Best
Victor
________________________________
De: John Halley Gotway via RT <met_help at ucar.edu>
Para: halmanset at yahoo.com
Enviado: Martes, 3 de junio, 2014 8:43 A.M.
Asunto: [rt.rap.ucar.edu #67410] Resolved: Wind direction statistics
by station
According to our records, your request has been resolved. If you have
any
further questions or concerns, please respond to this message.
------------------------------------------------
More information about the Met_help
mailing list