[Met_help] Troubleshooting test cases with MET 2.0

Remi Montroty remi.montroty at mfi.fr
Mon Apr 20 08:13:01 MDT 2009


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

-- 

MFI <http://www.mfi.fr> 	Rémi MONTROTY
Research Engineer -
Numerical Modeling and Data Assimilation
9, impasse Michel Labrousse -
31100 Toulouse - FRANCE
Phone: +33 (0)5.61.43.29.75 - Fax: +33 (0)5.61.43.29.41
Web: http://www.mfi.fr




More information about the Met_help mailing list