[Met_help] Troubleshooting test cases with MET 2.0

John Halley Gotway johnhg at rap.ucar.edu
Wed Apr 22 12:39:12 MDT 2009


Remi,

I spoke with the other MET developer about this and spent some time looking into it.  Our initial guess is that it may be some floating point issue particular to your machine/OS/compiler.  However,
that'd be pretty difficult to debug remotely.  I'd suggest continuing to run MODE on your local WRF output, and if you encounter any more errors similar to these, please let us know.  You could send
us your input data and configuration files, and we could try to reproduce the error here.  We'd be happy to help.

I agree that this issue is likely somewhere in vx_wrfdata.cc around line 1045.  In particular, I'm wondering about lines 1038 to 1040.  Is the angle value being computed actually <= 360.0 or is
something odd going on?
   angle = deg_per_rad*atan2(beta, alpha);
   angle -= 360.0*floor(angle/360.0);

Please let us know if you encounter any more errors or have questions about running MET.

Thanks,
John Halley Gotway

John Halley Gotway wrote:
> Remi,
> 
> The developer who really needs to look at this is out sick today.  But I'll talk to him about it tomorrow.
> 
> I'm very glad to hear that you're able to run MODE on your local WRF data without any errors.  Ultimately, that's more important than the test scripts.  But that makes the errors you're seeing even
> more puzzling.  Why would you see this error message with the test data, but not with your local GRIB output?
> 
> Please let me know if you see this or any other errors when running MET on your local WRF output.
> 
> Thanks,
> John
> 
> Remi Montroty wrote:
>> John,
>>
>> Yes, we were looking at that portion of the code and it would seem that
>> in vx_wrfdata, the angle is never set properly. Thus j stays at -1 and
>> we get that error message.
>> (around line 1045 of vx_wrfdata.cc). Given how the test that is done, we
>> don't quite understand why it is not set to k (j=k in the loop)...
>>
>> Regarding the GNU compilers, we did not manage to compile partly because
>> netcdf had been compiled with intel compilers first and that caused
>> problems.
>> However, even after recompiling netcdf with GNU compilers, we still had
>> an error message and I must say we gladly gave up....  So in other
>> terms, we never got that message about convex hulls since we never
>> managed to compile! I could have another look at it if need be.
>>
>> Any debugging file welcome. FYI, on the suggestion of a colleague, we
>> tested the MODE binary on two of our GRIB1 files (India domain) and made
>> it work smoothly. It would seem the test_mode.sh is thus causing more
>> problems than our local WRF run output files.
>>
>> Thanks,
>>
>> Remi
>>
>>
>> John Halley Gotway wrote:
>>> Remi,
>>>
>>> Thanks for registering.
>>>
>>> The error message you're seeing is rather odd.  This is the first time
>>> I've seen this one.  I really don't think that it has anything to do
>>> with NetCDF because you're getting the same error message
>>> from the last of the 3 MODE test cases which uses only GRIB input files.
>>>
>>> Just to clarify, when you say that you tried the GNU compilers to no
>>> avail, do you mean that you got the exact same error message about the
>>> convex hulls?  Or were there other problems?
>>>
>>> I'll talk to the developer who wrote this section of code.  We may
>>> want to send you a debug version of a file or two, have you recompile
>>> MET, and rerun the MODE test case.  But I'll let you know.
>>>
>>> Thanks,
>>> John
>>>
>>> Remi Montroty wrote:
>>>  
>>>> Hi,
>>>>
>>>> I am trying to make sure MET v2.0 installed properly.
>>>>
>>>> I seem to have a problem with test_mode.sh  as follows:
>>>> [met at cipsproc3 scripts]$ ./test_mode.sh
>>>>
>>>> *** Running MODE on APCP using netCDF input for both forecast and
>>>> observation ***
>>>> Forecast File: ../out/pcp_combine/sample_fcst_12L_2005080712V_12A.nc
>>>> Observation File: ../out/pcp_combine/sample_obs_2005080712V_12A.nc
>>>> Match Config File: config/WrfModeConfig_APCP_12
>>>> Merge Config File: config/WrfModeConfig_APCP_12
>>>> Forecast Field: APCP_12 at SFC
>>>> Observation Field: APCP_12 at SFC
>>>> Identifying objects in the forecast and observation fields...
>>>> Computing contingency table statistics...
>>>> Identified: 7 forecast objects and 4 observation objects.
>>>> Performing merging (threshold merging) in the forecast field.
>>>> Performing merging (threshold merging) in the observation field.
>>>> Remaining: 7 forecast objects and 4 observation objects.
>>>> Performing matching between the forecast and observation fields.
>>>>
>>>>
>>>> ERROR: convex_hull(Polyline &) -> can't find next hull point
>>>>
>>>>
>>>> *** Running MODE on APCP using a GRIB forecast and netCDF observation
>>>> ***
>>>> Forecast File: ../data/sample_fcst/2005080700/wrfprs_ruc13_24.tm00_G212
>>>> Observation File: ../out/pcp_combine/sample_obs_2005080800V_24A.nc
>>>> Match Config File: config/WrfModeConfig_APCP_24
>>>> Merge Config File: config/WrfModeConfig_APCP_24
>>>> Searching GRIB forecast file for GRIB code 61 with accumulation/level
>>>> indicator of A24.
>>>> Forecast Field: APCP_24 at SFC
>>>> Observation Field: APCP_24 at SFC
>>>> Identifying objects in the forecast and observation fields...
>>>> Computing contingency table statistics...
>>>> Identified: 4 forecast objects and 6 observation objects.
>>>> Performing merging (no merging) in the forecast field.
>>>> Performing merging (no merging) in the observation field.
>>>> Remaining: 4 forecast objects and 6 observation objects.
>>>> Performing matching between the forecast and observation fields.
>>>>
>>>>
>>>> ERROR: convex_hull(Polyline &) -> can't find next hull point
>>>>
>>>>
>>>> *** Running MODE on RH at 500mb using GRIB input for both forecast
>>>> and observation ***
>>>> Forecast File: ../data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212
>>>> Observation File:
>>>> ../data/sample_fcst/2005080712/wrfprs_ruc13_00.tm00_G212
>>>> Match Config File: config/WrfModeConfig_RH
>>>> Merge Config File: config/WrfModeConfig_RH
>>>> Searching GRIB forecast file for GRIB code 52 with accumulation/level
>>>> indicator of P500.
>>>> Searching GRIB observation file for GRIB code 52 with
>>>> accumulation/level indicator of P500.
>>>> Forecast Field: RH at ISBL_500
>>>> Observation Field: RH at ISBL_500
>>>> Identifying objects in the forecast and observation fields...
>>>> Computing contingency table statistics...
>>>> Identified: 8 forecast objects and 12 observation objects.
>>>> Performing merging (no merging) in the forecast field.
>>>> Performing merging (no merging) in the observation field.
>>>> Remaining: 8 forecast objects and 12 observation objects.
>>>> Performing matching between the forecast and observation fields.
>>>>
>>>> ERROR: convex_hull(Polyline &) -> can't find next hull point
>>>>
>>>> The fact that the convex_hull function fails seems problematic. What
>>>> should I look into?
>>>>
>>>> Thanks!
>>>>
>>>> Remi
>>>>
>>>>     
>>>   
> _______________________________________________
> 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