Commit 0562bab4 authored by Zolnierczuk, Piotr's avatar Zolnierczuk, Piotr
Browse files

reporting error 

do not save status variables between calls, e.g.
this implies save
 real(kind=DBL) :: dphase_deg_max = 0
and this does not
 real(kind=DBL) :: dphase_deg_max
 dphase_deg_max = 0
parent 8c2adadb

Makefile.version

+1 −1
Original line number Diff line number Diff line
@@ -3,7 +3,7 @@ export PROJECT=drspine 


export VERSION_MAJOR=1

export VERSION_MINOR=3

export VERSION_RELEASE=1

export VERSION_RELEASE=2



export PROJLIB=lib$(PROJECT).a

export PROJARCH=$(PROJECT)-$(VERSION_MAJOR).$(VERSION_MINOR)

sources/drspine.f90

+12 −6
Original line number Diff line number Diff line
@@ -2435,7 +2435,7 @@ CONTAINS 


    double precision    :: scpar(5), last_scpar(5) = 0d0

    real(kind=DBL)      :: dj_average, dj_variance, dj_value, dphase_deg

    real(kind=DBL)      :: dphase_deg_max = 0

    real(kind=DBL)      :: dphase_deg_max ! = 0 this implies (save)

    logical             :: phase_warning



!! -- data for bvariation checking ---
@@ -2451,11 +2451,11 @@ CONTAINS 
!    real(kind=DBL), save           :: bfield_max_diff = 1d-2

!    real(kind=DBL), save           :: PHASE_OFFSET_WARN_LIMIT = 15.0d0



    real(kind=DBL)                 :: chi2_max        = 0

    real(kind=DBL)                 :: bdiff_max       = 0

    logical                        :: bvalue_is_strange = .FALSE.

    logical                        :: bvalue_differs_from_ref = .FALSE.

    integer                        :: plot_start_number = 0

    real(kind=DBL)                 :: chi2_max                != 0

    real(kind=DBL)                 :: bdiff_max               != 0

    logical                        :: bvalue_is_strange       != .FALSE.

    logical                        :: bvalue_differs_from_ref != .FALSE.

    integer                        :: plot_start_number       != 0

    integer                        :: ipl, ips



    integer                        :: MAX_PIC_FIG
@@ -2473,6 +2473,12 @@ CONTAINS 
    !

    character(len=MAX_PATH_LENGTH)        :: report_name ! to allow to override



    dphase_deg_max = 0

    chi2_max                = 0

    bdiff_max               = 0

    bvalue_is_strange       = .FALSE.

    bvalue_differs_from_ref = .FALSE.

    plot_start_number       = 0



    verbose  = get_loglevel()>LOG_INFO

sources/read_data.f90

+1 −1
Original line number Diff line number Diff line
@@ -240,7 +240,7 @@ contains 
             nphase = nphase - 1

             if ((number_of_phases-nphase)<=max_number_of_missing_phases) then ! FIXME: (paz) HACK - accept if at most two up points at the end of a scan

                call msg_warning('read_echo_data', &

                  trim(msg_fmt("('scan pointt:', i0)", nscan))//' is missing '//&

                  trim(msg_fmt("('scan point:', i0)", nscan))//' is missing '//&

                  trim(msg_fmt("(i0, ' phase point(s), accepting the data')", number_of_phases-nphase)))

                exit scan_point_do

             end if