[Met_help] A problem

John Halley Gotway johnhg at rap.ucar.edu
Tue Jun 30 07:16:28 MDT 2009


Hello,

You are correct.  The 2 most important parameters in the MODE configuration file are the convolution radius and the convolution threshold.  They determine how the objects are defined from the raw field.

And actually, there are 4 entries in the config file for these - 2 for the forecast field and 2 for the observation field:
fcst_conv_radius
obs_conv_radius
fcst_conv_thresh
obs_conv_thresh

But if you're comparing the same field between the fcst and obs, just set the obs settings to the same values as the fcst.

To address your question... unfortunately, there is no one answer based on the model resolution.  Instead, it's up to the user to decide what scale of objects to analyze.

But here are some guidelines:
- Setting the convolution radius lower will make the objects less smooth and more detailed.
- Setting the convolution radius higher will make the objects smoother.
- Setting the threshold lower will make the objects bigger.
- Setting the threshold higher will make the objects smaller.

It's up to you to decide how you'd like to define the objects based on what type of features you'd like to extract from the field.  For example, when running MODE on precip, you may define objects
that are pretty large and smooth to represent large systems or MCS's.  Or you may define them to be pretty small and detailed to capture convection.  You need to ask yourself the question, what
features am I trying to analyze in my output?

You'll just need to play around with running MODE to get a sense of how the settings work.  Choose a setting for the convolution radius (5, 10, 15, 20 grid squares?).  Just pick one and set it.  Then
set the convolution threshold to some value - for precip, just try "gt0.0".  Then run MODE.  And bring up the output PostScript file in a window.  Next, adjust the convolution threshold, rerun MODE,
and observe how the objects change.  Once you get a sense for how the convolution threshold works, try fixing that one and playing with the convolution radius, moving it up and down, and observing the
changes in the PostScript output.  At first, I'd suggest keeping radius or threshold fixed, while you adjust the other one.

Sorry I can't provide more help.  You'll need to play around with it to get a sense of what it's doing.  If you try to run MODE to define a lot of objects (more than about 15 in each field), it'll
take a lot longer to run since there's more calculations to perform for each object.

Once you do find some settings you like, feel free to send some sample data to us for advice on interpreting the output.  If you do send some output, be sure to send the MODE configuration file you
used and the PostScript plots.

Thanks,
John Halley Gotway
johnhg at ucar.edu

zhxubinchaoshan wrote:
> Hello,I am using the MODE tool in MET for verification. 
> And I have some problems about the parameters in it. In the WrfModeConfig, there are two important parameters.They are obs_conv_radius and obs_conv_thresh.I have read some referances, and I found that how to specify these two parameters is dependent on my model resolution. If my model resolution is 12km, how can I specify these two parameters? 
> Thank you very much!
> 
> 
> 
> ÔÚ2009-05-13£¬"John Halley Gotway" <johnhg at rap.ucar.edu> дµÀ£º
>> Great.  Glad that did the trick for you.
>>
>> John
>>
>> zhxubinchaoshan wrote:
>>> Thanks for your reply!It can work now! 
>>> The problem is the "fcst_var". By your suggestion, I change "APCP_6" to "APCP_06", and then it can work successfully! 
>>> Thank you! 
>>>
>>>
>>>
>>>
>>> ÔÚ2009-05-11£¬"John Halley Gotway" <johnhg at rap.ucar.edu> дµÀ£º
>>>> Hello,
>>>>
>>>> What this warning message is telling you is that the Stat-Analysis tool is
>>>> looking in the directory "out/rain" for *.stat files, but it is not
>>>> finding any lines in those files that EXACTLY match the job you've
>>>> specified.
>>>>
>>>> When I experience a problem like this, my usual approach is to take a step
>>>> back and under-specify the job just to make sure I'm matching some lines. 
>>>> Then I take a look at the STAT lines that went into the job and adjust
>>> >from there.  For example, you may try the following:
>>>> (1) Run this Stat-Analysis job on the command line:
>>>> ../bin/stat_analysis -lookin ./out/rain \
>>>> -job aggregate_stat -line_type FHO -out_line_type CTS \
>>>> -dump_row tmp.stat
>>>>
>>>> (2) Did this job run successfully? Did it find any input FHO lines to use?
>>>> If not, then maybe there aren't any FHO lines in your STAT files.  If it
>>>> did find FHO lines, and the job ran fine, proceed to the next step.
>>>>
>>>> (3) Take a look in the file "tmp.stat".  How would you like to filter this
>>>> data down more?  Do you see "APCP_6" in the column for "fcst_var"?  Or
>>>> maybe it's actually "APCP_06"?  Try adding that to the job:
>>>>
>>>> ../bin/stat_analysis -lookin ./out/rain \
>>>> -job aggregate_stat -line_type FHO -fcst_var APCP_6 -out_line_type CTS \
>>>> -dump_row tmp.stat
>>>>
>>>> (4) Did this job run successfully, or was there a problem with how you
>>>> specified the "fcst_var".  Now check the "tmp.stat" file again, and adjust
>>>> your job as necessary.  One hint though, when running Stat-Analysis jobs
>>>> like this on the command line, you'll usually need to put a backslash
>>>> before special characters like this: -fcst_thresh "\>0.000"
>>>>
>>>> So I'd suggest approaching the problem like that.  If you can't figure it
>>>> out, feel free to send along all of the STAT files in your "out/rain"
>>>> directory, the Stat-Analysis config file you're using, and the command
>>>> line you're using... and I'll try to figure out what's going on.
>>>>
>>>> Good luck.
>>>>
>>>> John Halley Gotway
>>>> johnhg at ucar.edu
>>>>
>>>>> When I run the Stat_analysis Tool, I get the error below:
>>>>> WARNING: do_job_aggr_stat() -> no matching STAT lines found for job: -job
>>>>> aggregate_stat -fcst_var APCP_6 -fcst_thresh >0.000 -line_type FHO
>>>>> -dump_row ./out/stat_analysis/job_aggregate_stat_FHO_CTS.stat
>>>>> -out_line_type CTS -out_alpha 0.050000
>>>>>
>>>>> And my command is :
>>>>>  ../bin/stat_analysis -config ./STATAnalysisConfig_07 -lookin ./out/rain
>>>>> -out ./out/stat_analysis/stat_analysis.out
>>>>>
>>>>> Can you tell me how to figure out the problem above?
>>>>> Thank you very much!
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> 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