Segmentation fault when using ECMWF analysis data – in #10: INT2LM

  @fannidorakelemen in #0c17036

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