[MITgcm-support] Run-time errors on Archer

[Martin Losch]

Dan,

maybe you can find out what’s happening in your line 3331 in print.f (that should be within PRINT_LIST_I).
There are a couple of internal writes of the type
 WRITE(fmt1,’(A,I1,A)’) ’some string’,someInteger,’some string’

For debugging I would identify the line, and print to screen whatever is written to the character variable. Maybe there is a type mismatch that hasn’t been noticed before, because previous compilers have more mercy than yours? This is unlikely, but you never know. If this is not the problem, this debugging procedure may give you a hint of what may have gone wrong (maybe some type mismatch in the name list that wasn’t caught, etc.).

Martin


On Mar 7, 2014, at 2:16 PM, Dan Jones <dcjones.work at gmail.com> wrote:

> Correction:
> 
> The snippet of source code listed above actually comes from yet *another* error produced if the obcs package is included in packages.conf:
> 
> forrtl: error (63): output conversion error, unit -5, file Internal Formatted Write 
> Image              PC                Routine            Line        Source             
> mitgcmuv           000000000302A3EA  Unknown               Unknown  Unknown 
> 
> mitgcmuv           00000000023B5F85  print_list_i_            3331  print.f 
> 
> mitgcmuv           0000000001A302E7  obcs_readparms_          2425  obcs_readparms.f 
> mitgcmuv           000000000297EBF0  packages_readparm        1907  packages_readparms.f 
> 
> mitgcmuv           0000000002932168  initialise_fixed_        1874  initialise_fixed.f 
> 
> mitgcmuv           0000000002B03175  the_model_main_          3052  the_model_main.f 
> mitgcmuv           00000000023A38F1  MAIN__                   4407  main.f 
> 
> mitgcmuv           0000000000400F06  Unknown               Unknown  Unknown 
> 
> mitgcmuv           00000000030A8124  Unknown               Unknown  Unknown 
> mitgcmuv           0000000000400DD1  Unknown               Unknown  Unknown 
> 
> 
>       IF ( debugLevel.GE.debLevA ) THEN 
>         CALL PRINT_LIST_I( OB_Jnorth, 1, OBNS_Nx, INDEX_I, 
>     &                    .FALSE., .TRUE., standardMessageUnit ) 
> It does *not* come from ini_depths, but from obcs_readparms.  Sorry about that.
> 
> Dan
> 
> 
> On Fri, Mar 7, 2014 at 1:12 PM, Dan Jones <dcjones.work at gmail.com> wrote:
> Greetings:
> 
> I am having trouble getting MITgcm to run on Archer.  I am using the Intel compiler (14.0.1.106) with the following defines/flags in the build options file:
> 
> DEFINES='-DALLOW_USE_MPI -DALWAYS_USE_MPI -D_BYTESWAPIO -DWORDLENGTH=4 -DHAVE_FLUSH'
> LIBS='-L${CRAY_MPICH2_DIR}/lib -L${HDF5_DIR}/lib -L$NETCDF_DIR/lib -lnetcdf -lnetcdff  -lhdf5 -lhdf5_hl'
> INCLUDES='-I${CRAY_MPICH2_DIR}/include -I${HDF5_DIR}/include -I${NETCDF_DIR}/include -I${HDF5_INCLUDE_OPTS}'
> FFLAGS='-h byteswapio -assume byterecl -convert big_endian -heap-arrays -O2 -g -traceback'
> 
> The code compiles with no errors, but it does not run.  The code crashes with the error:
> 
> ABNORMAL END: S/R INI_THETA
> 
> with no other information.  The initial theta file is fine and has been used successfully in other MITgcm model setups.  When I turn on the debugger (i.e. set debugMode=.TRUE. in input/eedata and set the debugLevel=4 in input/data), I get a *different* error that appears to occur in an *earlier* part of the code.  The code crashes as it tries to read in the bathymetry file:
> 
> forrtl: error (63): output conversion error, unit -5, file Internal Formatted Write 
> 
> Image              PC                Routine            Line        Source             
> mitgcmuv    00000000023CCCBC  print_maprs_         4981  print.f
> mitgcmuv    000000000297FD81   plot_field_xyrs_       1841  plot_field.f
> mitgcmuv    000000000276CDC2  ini_depths_           3271  ini_depths.f
> mitgcmuv    0000000002932628   initialise_fixed_      1908  initialise_fixed.f
> mitgcmuv    0000000002B03175   the_model_main       3052  the_model_main.f
> mitgcmuv    00000000023A38F1   MAIN__           4407  main.f
> 
> Again, the bathymetry file is fine and has been used successfully before.  The problem indicated above in ini_depths.f happens in this function:
>         CALL PRINT_LIST_I( OB_Jnorth, 1, OBNS_Nx, INDEX_I, 
>     &                    .FALSE., .TRUE., standardMessageUnit ) 
> I can suppress the output conversion error by re-compiling with a -check nooutput_conversion flag, but the code quickly produces a segmentation fault at about the same place (ini_depths.f and the functions that call it):
> 
> forrtl: severe (194): SIGSEGV, segmentation fault occurred
> mitgcmuv           00000000006FFAE7  print_maprs_             4982  print.f 
> mitgcmuv           00000000007950D3  plot_field_xyrs_         1841  plot_field.f 
> 
> 
> mitgcmuv           0000000000760811  ini_depths_              3271  ini_depths.f 
> 
> mitgcmuv           000000000078699A  initialise_fixed_        1908  initialise_fixed.f 
> mitgcmuv           00000000007AFA93  the_model_main_          3052  the_model_main.f 
> 
> 
> mitgcmuv           00000000006F5C41  MAIN__                   4407  main.f 
> 
> The fact that turning on the debugger produces an error *earlier* in the code is the most interesting/distressing bit here.  Is this an I/O issue?  Has anyone else run into something like this?  I have contacted the Archer support team, but I thought it would be worth asking around here as well.
> 
> Many thanks,
> Dan
> 
> -- 
> *************************************************
> 
> Dr Dan Jones
> Open Oceans Group
> British Antarctic Survey
> Cambridge, UK
> 
> Phone: +44 (0)1223 221505
> Fax: +44 (0)1223 362616
> 
> *************************************************
> 
> 
> 
> -- 
> *************************************************
> 
> Dr Dan Jones
> Open Oceans Group
> British Antarctic Survey
> Cambridge, UK
> 
> Phone: +44 (0)1223 221505
> Fax: +44 (0)1223 362616
> 
> *************************************************
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support