Opened 3 years ago

Closed 3 years ago

#459 closed defect (fixed)

interpweight crash when reading PFT map

Reported by: mmcgrath Owned by:
Priority: minor Milestone:
Component: Physical processes Version:
Keywords: PFT map, species, zoomed Cc:


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.

Attachments (1)

edit_netcdf.f90 (1.6 KB) - added by mmcgrath 3 years ago.
Code to modify a PFTmap by adding an attribute and changing NaN values

Download all attachments as: .zip

Change History (3)

Changed 3 years ago by mmcgrath

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

comment:1 Changed 3 years ago by 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.

comment:2 Changed 3 years ago by mmcgrath

  • Resolution set to fixed
  • Status changed from new to closed
Note: See TracTickets for help on using tickets.