interpweight crash when reading PFT map

[mmcgrath]

I get the following error at the beginning of an ORCHIDEE run in full debug mode on Obelix.

orchidee_ol 0000000002A71A95 interpweight_mp_i 1146 interpweight.f90
orchidee_ol 0000000000D94556 slowproc_mp_slowp 3879 slowproc.f90
orchidee_ol 0000000000C302A1 slowproc_mp_slowp 1680 slowproc.f90
orchidee_ol 0000000000BDD324 slowproc_mp_slowp 213 slowproc.f90
orchidee_ol 00000000009BF6A0 sechiba_mp_sechib 613 sechiba.f90
orchidee_ol 00000000005E0E13 intersurf_mp_inte 388 intersurf.f90
orchidee_ol 00000000004DD198 MAIN 1259 dim2_driver.f90

it happens even for a single processor. I've traced it back to the reading in of the PFT map. In particular, it's the combined species/PFT map that is zoomed over Europe (a grid of 128x119x28), and it doesn't happen with a 13 PFT map at half degree (720x360x13).

After some more tracing, the map comes back all NaNs? from NF90_GET_VAR on line 1636 of flincom.f90, but it doesn't seem to throw an error then.

[mmcgrath]

Code to modify a PFTmap by adding an attribute and changing NaN values

[mmcgrath]

It appears to have been a problem with the map itself. I wrote a code to add an attribute and change all the NaN in the non-data points to the Fillvalue/missing_value (attached). Now the code runs past interpweight and crashes in sechiba. I will close this ticket and look into that to see if it is related. If so, I will reopen this ticket.

I compile the code with

mpif90 edit_netcdf.f90 -I/usr/local/install/netcdf-4.3.2p/include -L/usr/local/install/netcdf-4.3.2p/lib -lnetcdff -lnetcdf

on obelix.