[MITgcm-support] Initialization question

[William Aiken]

I have been investigating why I get

>STOP: NORMAL END
>>  Note: IEEE floating-point exception flags raised:
>>     Inexact;  Invalid Operation;

messages when I build and run the exp2 model and verification case
on my Solaris x64 system.

I followed the logic through:

main
 the_main_model
  initialize_varia
   ini_fields
    ini_theta
     read_fld_xyz_rl
      mdsio_readfield
       mdsio_seg4torl

I'm opening an unnamed global file in mdsio_readfield, reading each
record into r4seg and passing this array to MDS_SEG4toRL, where the data
is copied to array arr. The values in r4seg seem to be random, differing
by 30 orders of magnitude or more, but at least one of the values is
-NaN. The above IEEE warnings are triggered when this value is copied to
arr.

I'd like to sanitize this part of MITgcm to detect NaN's when reading
the input, but am an MITgcm newbie and unsure where the best place in
the code is to do these checks.

Is mdsio_readfield the best routine to change in this case, or is there
another routine which writes the random data to a file which is then
read by mdsio_readfield?

Any suggestions?

Thanks,

Bill

-- 
William Aiken          Sun Microsystems, Inc.
Phone: 781.442.3312    Market Development Engineering
FAX: 781.442.1678      One Network Drive UBUR01-204
                       Burlington, Massachusetts 01803