[Met_help] Troubleshooting test cases with MET 2.0

[John Halley Gotway]

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
>>>
>>>     
>>
>>   
>