Project

General

Profile

Segmentation fault when using ECMWF analysis data

Added by Fanni Dora Kelemen about 6 years ago

Dear All,

I have problem using int2lm with ECMWF analysis data, which are converted into netcdf format.

The program stops at the vertical interpolation:

0: in vert_interpol: first pressure adaptation
0: forrtl: severe (174): SIGSEGV, segmentation fault occurred

I attached the namelist, and the error message.

Do you have any idea what could be wrong?

The header of one of the netcdf input files is :

netcdf cas2014060900 {
dimensions:
bnds = 2 ;
lon = 269 ;
lat = 268 ;
level = 137 ;
level1 = 138 ;
soil = 3 ;
soil1 = 4 ;
time = UNLIMITED ; // (1 currently)
pressure = 1 ;
variables:
float lon(lon) ;
lon:axis = “X” ;
lon:standard_name = “longitude” ;
lon:long_name = “longitude” ;
lon:units = “degrees_east” ;
float lat(lat) ;
lat:axis = “Y” ;
lat:standard_name = “latitude” ;
lat:long_name = “latitude” ;
lat:units = “degrees_north” ;
int level(level) ;
level:standard_name = “atmosphere_hybrid_sigma_pressure_coordinate” ;
level:long_name = “mid-layer hybrid coordinate (sigma-pressure)” ;
level:positive = “down” ;
level:formula_terms = “ap: akm b: bkm ps: PS” ;
float akm(level) ;
akm:long_name = “hybrid coordinate (sigma-pressure)” ;
akm:units = “Pa” ;
float bkm(level) ;
bkm:long_name = “hybrid coordinate (sigma-pressure)” ;
bkm:units = “1” ;
int level1(level1) ;
level1:standard_name = “atmosphere_hybrid_sigma_pressure_coordinate” ;
level1:long_name = “hybrid coordinate (sigma-pressure)” ;
level1:positive = “down” ;
level1:formula_terms = “ap: ak b: bk ps: PS” ;
float ak(level1) ;
ak:long_name = “hybrid coordinate (sigma-pressure)” ;
ak:units = “Pa” ;
float bk(level1) ;
bk:long_name = “hybrid coordinate (sigma-pressure)” ;
bk:units = “1” ;
float soil1(soil1) ;
soil1:standard_name = “depth” ;
soil1:long_name = “depth of soil layers” ;
soil1:units = “m” ;
soil1:positive = “down” ;
soil1:bounds = “soil1_bnds” ;
float soil1_bnds(soil1, bnds) ;
double time(time) ;
time:units = “seconds since 2014-06-09 00:00:00” ;
time:standard_name = “time” ;
time:long_name = “time” ;
time:calendar = “standard” ;
time:bounds = “time_bnds” ;
double time_bnds(time, bnds) ;
time_bnds:long_name = “time bounds” ;
float pressure(pressure) ;
float T(time, level, lat, lon) ;
float U(time, level, lat, lon) ;
float V(time, level, lat, lon) ;
float QV ;
float QC ;
float QI ;
float PS ;
float FIS ;
float FR_LAND(time, lat, lon) ;
float W_SNOW(time, lat, lon) ;
float T_SNOW(time, lat, lon) ;
float T_SKIN(time, lat, lon) ;
float W_I(time, lat, lon) ;
float FR_SEA_ICE(time, lat, lon) ;
float T_CL(time, lat, lon) ;
float W_SO_REL(time, soil1, lat, lon) ;
float FI ;

// global attributes:
:institution = “ECMWF” ;
:Conventions = “CF-1.4” ;
:conventionsURL = “http://www.cfconventions.org” ;
:source = “analysis/reanalysis” ;
:creation_date = “2016-03-22 17:08:46”

Thank you for your help!
Best regards,
Fanni