mksurfdata_map resulting segmentation fault occurred

[evasinha]

I am trying to generate a surface data set using the `mksurfdata_map` tool but keep getting this error even when I have tried using multiple nodes for running the tool:

(gridmap_map_read) reading mapping matrix data...
(gridmap_map_read) * file name : /global/cfs/cdirs/m2702/liliyao/inputdata/cesm_inputdata/lnd/clm2/mappingdata/maps/0.125x0.125/map_0.25x0.25_MODIS_to_0.125x0.125_nomask_aave_da_c170321.nc
* matrix dimensions rows x cols : 1036800 x 4147200
* number of non-zero elements: 2184582
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
libpthread-2.31.s 00007FE5AC811910 Unknown Unknown Unknown
mksurfdata_map 00000000004571B7 mkpftmod_mp_mkpft 563 mkpftMod.F90
mksurfdata_map 000000000047768C MAIN__ 568 mksurfdat.F90

How can I resolve this error?

 

[slevis]

As you suggested, the PE layout that you requested may explain the failure. You may have to continue experimenting with different PE layouts, until you find one that works. If someone in your group has run this tool successfully before on the same machine (liliyao?) you may get good guidance from them.

 

[evasinha]

@slevis Thank you for response and suggestion. I tried few more PE layouts and even tried to generate a coarse resolution file (0.9x1.25) but was running in to the same error. I touched base with Lili (@liliyao ) but she is using `mksurfdata_esmf` as opposed to `mksurfdata_map` tool. In the end I tried using CESM mksurfdata_map on a different machine (PNNL's compy machine as opposed to NERSC perlmutter) and was able to generate the surface dataset and land use timeseries at 0.25x0.25 resolution on compy.