Segmentation fault when using ECMWF analysis data – in #10: INT2LM
in #10: INT2LM
Cookies disclaimer
Our site saves small pieces of text information (cookies) on your
device in order to verify your login. These cookies are essential
to provide access to resources on this website and it will not
work properly without.
Learn more
<p>
Dear All,
</p>
<p>
I have problem using int2lm with
<span class="caps">
ECMWF
</span>
analysis data, which are converted into netcdf format.
</p>
<p>
The program stops at the vertical interpolation:
</p>
<p>
0: in vert_interpol: first pressure adaptation
<br/>
0: forrtl: severe (174):
<span class="caps">
SIGSEGV
</span>
, segmentation fault occurred
</p>
<p>
I attached the namelist, and the error message.
</p>
<p>
Do you have any idea what could be wrong?
</p>
<p>
The header of one of the netcdf input files is :
</p>
<p>
netcdf cas2014060900 {
<br/>
dimensions:
bnds = 2 ;
lon = 269 ;
lat = 268 ;
level = 137 ;
level1 = 138 ;
soil = 3 ;
soil1 = 4 ;
time =
<span class="caps">
UNLIMITED
</span>
; // (1 currently)
pressure = 1 ;
<br/>
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(time, level, lat, lon) ;
float QC(time, level, lat, lon) ;
float QI(time, level, lat, lon) ;
float PS(time, lat, lon) ;
float
<acronym title="time, lat, lon">
<span class="caps">
FIS
</span>
</acronym>
;
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(time, lat, lon) ;
</p>
<p>
// global attributes:
:institution = “
<span class="caps">
ECMWF
</span>
” ;
:Conventions = “CF-1.4” ;
:conventionsURL = “http://www.cfconventions.org” ;
:source = “analysis/reanalysis” ;
:creation_date = “2016-03-22 17:08:46”
</p>
<p>
Thank you for your help!
<br/>
Best regards,
<br/>
Fanni
</p>
Segmentation fault when using ECMWF analysis data
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(time, level, lat, lon) ; float QC(time, level, lat, lon) ; float QI(time, level, lat, lon) ; float PS(time, lat, lon) ; 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(time, lat, lon) ;
// 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