Opened 4 years ago
Closed 4 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: |
Description
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)
Change History (3)
Changed 4 years ago by mmcgrath
comment:1 Changed 4 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 4 years ago by mmcgrath
- Resolution set to fixed
- Status changed from new to closed
Code to modify a PFTmap by adding an attribute and changing NaN values