floating invalid for 'lon' when mapping

[Xinjie Huang]

Hi all,

When I am using epa_anthro_emis to process the dataset, I have the following floating invalid error: 

srun: error: della-h12n16: task 0: Aborted (core dumped)
srun: Terminating StepId=63585872.0
forrtl: error (65): floating invalid
Image              PC                Routine            Line        Source            
anthro_emis        000000000056BEDB  Unknown               Unknown  Unknown
libpthread-2.28.s  00001523A84CD990  Unknown               Unknown  Unknown
anthro_emis        0000000000453B97  area_mapper_mp_ll         992  area_mapper.f90
anthro_emis        0000000000452BAD  area_mapper_mp_ll         760  area_mapper.f90
anthro_emis        0000000000452CD3  area_mapper_mp_ll         779  area_mapper.f90
anthro_emis        000000000048DD1E  stack_mp_stkfilei         165  stack.f90
anthro_emis        00000000004CF2DA  data_file_utils_m         620  data_file_utils.f90
anthro_emis        0000000000501990  MAIN__                    487  anthro_emis.f90
anthro_emis        000000000040C9A2  Unknown               Unknown  Unknown
libc-2.28.so       00001523A811F7E5  __libc_start_main     Unknown  Unknown
anthro_emis        000000000040C8AE  Unknown               Unknown  Unknown

The error pointed to the Line 992 in area_mapper_mp_ll:

 IF (deltalon .GT. +180.) deltalon = deltalon - 360.

So I printed the value of 'lon' out and found it became unstable with extremely large numbers (both positive and negative), and showed NAN in the end which finally caused this issue. I am attaching the .out file, .inp file, and my namelists. 

Any help would be so much appreciated. 

Best, 

Xinjie

[Gabriele Pfister]

Hi Xinje

this is very strange and I have never encountered this problem before. Have you tried running only d01 and does that work? 

Gabi

--

================================
Gabriele Pfister

Atmospheric Chemistry Observations & Modeling (ACOM) Lab

NSF National Center for Atmospheric Research (NCAR)

Email: pfi...@ucar.edu

Work Phone: +1 303 497 2915

Web: https://staff.ucar.edu/users/pfister

[Xinjie Huang]

Thanks Gabi for your reply.

This issue always occurs with stack emission files, i.e.,  running only d01 with stack emissions also failed. 

I tried to exclude stack emissions and it works for d01, but not for d02 and d03 for other reasons. 

I don't know if it is related to the compilation - I modified the line 68-71 of make_anthro as follows due to a library issue:

set ar_libs = "-lnetcdf -lnetcdff"

if( -e ${NETCDF_DIR}/lib/libnetcdff.a ) then

  set ar_libs = "${ar_libs} -lnetcdff"

endif

I am also attaching my Makefile and make_anthro for your review. Let me know if you have any thoughts! Thanks!

Xinjie