Initialize the NaN_rk constant
Must be called before using NaN_rk. Automatically called by other functions.
subroutine init_nan() if (.not. nan_initialized) then NaN_rk = ieee_value(0.0_rk, ieee_quiet_nan) nan_initialized = .true. end if end subroutine init_nan