ungrib.exe error with 1/12 degree SST from NCEP

[Egsar]

We are using the 0.5 degree resolution from NCEP ( NCEP Data Products SST ) as constant SST field for our WRF/real simulations.

To achieve this, we use ungrib.exe binary from WPS package with Vtable.SST also provided by WPS.

This works realy well, but now we plan to use the 1/12 degree resolution provided in the same place.

Unfortunately this time ungrib.exe fails with error:

forrtl: severe (174): SIGSEGV, segmentation fault occurred

Searching the internet I could not find a solution for the problem and I'm clueless here so any advice would be helpful.

 

[Ming Chen]

Have you compared 1/12 degree data with 0.5 degree SST data? I am not sure whether they have the same GRIB parameter, etc. If not, it is possible that ungrib may extract wrong information from the 1/12 degree data.

 

[Egsar]

Thank you for your comment. Yes we have checked both files (1/2 and 1/12 degree resolution)


* This is the output of "g2print.exe rtgssthr_grb_0.5.grib2":

ungrib - grib edition num 2
reading from grib file = rtgssthr_grb_0.5.grib2

NCEP SST Analysis
---------------------------------------------------------------------------------------
rec Prod Cat Param Lvl Lvl Lvl Prod Name Time Fcst
num Disc num code one two Templ hour
---------------------------------------------------------------------------------------
1 0 0 0 1 0 0 0 TMP 2024-04-22_00:00:00 00


Successful completion of g2print


* And this is the output from "g2print.exe rtgssthr_grb_0.083.grib2":

ungrib - grib edition num 2
reading from grib file = rtgssthr_grb_0.083.grib2

NCEP SST Analysis
---------------------------------------------------------------------------------------
rec Prod Cat Param Lvl Lvl Lvl Prod Name Time Fcst
num Disc num code one two Templ hour
---------------------------------------------------------------------------------------
1 0 0 0 1 0 0 0 TMP 2024-04-22_00:00:00 00


Successful completion of g2print


* On the other hand "wgrib2 -V rtgssthr_grb_0.5.grib2" reports as follows:

1:0:vt=2024042200:surface:anl:TMP Temperature [K]:
ndata=259200:undef=0:mean=285.125:min=271.35:max=307.94
grid_template=0:winds(N/S):
lat-lon grid720 x 360) units 1e-06 input WE:NS output WE:SN res 48
lat 89.750000 to -89.750000 by 0.500000
lon 0.250000 to 359.750000 by 0.500000 #points=259200

* Finally "wgrib2 -V rtgssthr_grb_0.083.grib2" reports:

1:0:vt=2024042200:surface:anl:TMP Temperature [K]:
ndata=9331200:undef=0:mean=285.125:min=271.35:max=309.47
grid_template=0:winds(N/S):
lat-lon grid4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48
lat 89.958000 to -89.958000 by 0.083000
lon 0.042000 to 359.958000 by 0.083000 #points=9331200


Honestly I cannot see any significant difference that justify why "ungrib.exe" works with 1/2 degree resolution but not with 1/12.

I welcome any comments
Thank you