Loading sources/data_types.f90 +1 −1 Original line number Diff line number Diff line Loading @@ -1025,7 +1025,7 @@ CONTAINS ! real(kind=DBL) :: xtmp ! tmp variable real(kind=DBL) :: tau, qabs real(kind=DBL) :: ud_err2 !!real(kind=DBL) :: ud_err2 !! (paz), unused variable !!!>>> preliminary hint >>> !!! if(present(ps_background)) then Loading sources/drspine.f90 +1 −1 Original line number Diff line number Diff line Loading @@ -3577,7 +3577,7 @@ ntl: do i=0, nt if(verbose) write(*,*)"TEST: ntl ", j, i, ii, tau(ii), fqt(ii) write(*,*)"TEST NTL: ",j,i,ii, tau(ii), fqt(ii), fqt_alt(ii) , cdX%tq_bin(i,j)%sqt, cdX%tq_bin(i,j)%updown !!TEST!! defintion/role of updown?? !!write(*,*)"TEST NTL: ",j,i,ii, tau(ii), fqt(ii), fqt_alt(ii) , cdX%tq_bin(i,j)%sqt, cdX%tq_bin(i,j)%updown !!TEST!! defintion/role of updown?? end do ntl ntx = ii Loading sources/matching.f90 +2 −2 Original line number Diff line number Diff line Loading @@ -181,7 +181,7 @@ contains integer :: i, j character(len=MAX_NAME_LENGTH) :: cmsg integer :: l !>>test !! integer :: l !>>test !! (paz) unused variables source_loop: do i=1, source_data%number_of_points ! >>new>> mm Loading Loading @@ -251,7 +251,7 @@ contains integer, intent(in) :: match_role ! integer :: i, j integer :: l,lp ! TEST !!integer :: l,lp ! TEST !! (paz), unused variables ! source_loop: do i=1, scan_size if ( .not. is_valid_scan(scan_data(i)) ) cycle source_loop Loading sources/reduce_data.f90 +13 −13 Original line number Diff line number Diff line Loading @@ -109,15 +109,15 @@ dix: do ix=1, sig_point%no_xpix diy: do iy=1, sig_point%no_ypix if ( sig_point%pixelbin(it, ix, iy)%status /= PIXEL_OK ) cycle if(ix==4 .and. iy==4) write(*,*) & !!TEST!! "T combine ref",res_point%id,sig_point%matching_ref_arr(1)%ptr%id,& !!TEST!! associated(sig_point%matching_ref, sig_point%matching_ref_arr(1)%ptr) !!if(ix==4 .and. iy==4) write(*,*) & !!TEST!! !! "T combine ref",res_point%id,sig_point%matching_ref_arr(1)%ptr%id,& !!TEST!! !! associated(sig_point%matching_ref, sig_point%matching_ref_arr(1)%ptr) call combine_runs(sig_point%matching_ref_arr, it, ix, iy, resol ) ! later use resol%fqt if ( res_point%pixelbin(it, ix, iy)%status /= PIXEL_OK ) cycle if ( bgr_flag > 0 .and. has_bgr & .and. .not. sig_point%pixelbin(it, ix, iy)%background_subtracted ) then if(ix==4 .and. iy==4) write(*,*)"T combine bgr" !!TEST!! !!if(ix==4 .and. iy==4) write(*,*)"T combine bgr" !!TEST!! call combine_runs(sig_point%matching_bgr_arr, it, ix, iy, bgr ) ! for later use instead of bgr_point%pixelbin if (bgr_point%pixelbin(it, ix, iy)%status /= PIXEL_OK ) then sig_point%pixelbin(it, ix, iy)%status = & Loading @@ -132,11 +132,11 @@ if(ix==4 .and. iy==4) write(*,*)"T combine bgr" !!TEST!! !!TEST!! fqt = get_fqt_bgr(sig_point%pixelbin(it, ix, iy), bgr_point%pixelbin(it, ix, iy), tfac, volume_fraction) !!TEST!! if(ix==4 .and. iy==4) write(*,*)"T fqt1_sngl(old)",fqt !!TEST!! !!if(ix==4 .and. iy==4) write(*,*)"T fqt1_sngl(old)",fqt !!TEST!! fqt = get_fqt_bgr(sig_point%pixelbin(it, ix, iy), bgr, tfac, volume_fraction) if(ix==4 .and. iy==4) write(*,*)"T fqt2_mult(new)",fqt !!TEST!! !!if(ix==4 .and. iy==4) write(*,*)"T fqt2_mult(new)",fqt !!TEST!! else fqt = sig_point%pixelbin(it, ix, iy)%fqt endif Loading Loading @@ -202,7 +202,7 @@ if(ix==4 .and. iy==4) write(*,*)"T fqt2_mult(new)",fqt !!TEST!! endif stat = PIXEL_OK if(ix==4 .and. iy==4) write(*,'(a,i4,i9,2(i3,2f18.9))') "Tcombine_runs i amp fqt: ",i,run_pointer_array(i)%ptr%id, amp,fqt !!TEST!! !if(ix==4 .and. iy==4) write(*,'(a,i4,i9,2(i3,2f18.9))') "Tcombine_runs i amp fqt: ",i,run_pointer_array(i)%ptr%id, amp,fqt !!TEST!! !!TEST!! if(ix==4 .and. iy==4) write(*,'(a,i4,i9,2(i3,2f18.9))') "Tcombine_runs o amp fqt: ",i,run_pointer_array(i)%ptr%id,&!!TEST!! !!TEST!! run_pointer_array(i)%ptr%pixelbin(it, ix, iy)%amplitude , run_pointer_array(i)%ptr%pixelbin(it, ix, iy)%fqt !!TEST!! !!TEST!! if(ix==4 .and. iy==4) write(*,*)"FQT ", 2*amp%value/(up%value-down%value), run_pointer_array(i)%ptr%pixelbin(it, ix, iy)%fqt !!TEST!! Loading Loading @@ -239,7 +239,7 @@ if(ix==4 .and. iy==4) write(*,'(a,i4,i9,2(i3,2f18.9))') "Tcombine_runs i amp fqt type(value_struct), dimension(:), allocatable :: bgr_pixel_bin character(len=MAX_NAME_LENGTH) :: cmsg !character(len=MAX_NAME_LENGTH) :: cmsg, (paz), unused variable sample_loop: do i=1, sample_data%number_of_points Loading Loading @@ -453,7 +453,7 @@ bgr_pixel_bin(l) integer , intent(in) :: scan_size real(kind=DBL) , intent(in) :: transmission ! integer :: i, j integer :: i !!, j !! (paz) unused variable ! source_loop: do i=1, scan_size if ( .not. is_valid_scan(scan_data(i)) ) cycle source_loop Loading Loading
sources/data_types.f90 +1 −1 Original line number Diff line number Diff line Loading @@ -1025,7 +1025,7 @@ CONTAINS ! real(kind=DBL) :: xtmp ! tmp variable real(kind=DBL) :: tau, qabs real(kind=DBL) :: ud_err2 !!real(kind=DBL) :: ud_err2 !! (paz), unused variable !!!>>> preliminary hint >>> !!! if(present(ps_background)) then Loading
sources/drspine.f90 +1 −1 Original line number Diff line number Diff line Loading @@ -3577,7 +3577,7 @@ ntl: do i=0, nt if(verbose) write(*,*)"TEST: ntl ", j, i, ii, tau(ii), fqt(ii) write(*,*)"TEST NTL: ",j,i,ii, tau(ii), fqt(ii), fqt_alt(ii) , cdX%tq_bin(i,j)%sqt, cdX%tq_bin(i,j)%updown !!TEST!! defintion/role of updown?? !!write(*,*)"TEST NTL: ",j,i,ii, tau(ii), fqt(ii), fqt_alt(ii) , cdX%tq_bin(i,j)%sqt, cdX%tq_bin(i,j)%updown !!TEST!! defintion/role of updown?? end do ntl ntx = ii Loading
sources/matching.f90 +2 −2 Original line number Diff line number Diff line Loading @@ -181,7 +181,7 @@ contains integer :: i, j character(len=MAX_NAME_LENGTH) :: cmsg integer :: l !>>test !! integer :: l !>>test !! (paz) unused variables source_loop: do i=1, source_data%number_of_points ! >>new>> mm Loading Loading @@ -251,7 +251,7 @@ contains integer, intent(in) :: match_role ! integer :: i, j integer :: l,lp ! TEST !!integer :: l,lp ! TEST !! (paz), unused variables ! source_loop: do i=1, scan_size if ( .not. is_valid_scan(scan_data(i)) ) cycle source_loop Loading
sources/reduce_data.f90 +13 −13 Original line number Diff line number Diff line Loading @@ -109,15 +109,15 @@ dix: do ix=1, sig_point%no_xpix diy: do iy=1, sig_point%no_ypix if ( sig_point%pixelbin(it, ix, iy)%status /= PIXEL_OK ) cycle if(ix==4 .and. iy==4) write(*,*) & !!TEST!! "T combine ref",res_point%id,sig_point%matching_ref_arr(1)%ptr%id,& !!TEST!! associated(sig_point%matching_ref, sig_point%matching_ref_arr(1)%ptr) !!if(ix==4 .and. iy==4) write(*,*) & !!TEST!! !! "T combine ref",res_point%id,sig_point%matching_ref_arr(1)%ptr%id,& !!TEST!! !! associated(sig_point%matching_ref, sig_point%matching_ref_arr(1)%ptr) call combine_runs(sig_point%matching_ref_arr, it, ix, iy, resol ) ! later use resol%fqt if ( res_point%pixelbin(it, ix, iy)%status /= PIXEL_OK ) cycle if ( bgr_flag > 0 .and. has_bgr & .and. .not. sig_point%pixelbin(it, ix, iy)%background_subtracted ) then if(ix==4 .and. iy==4) write(*,*)"T combine bgr" !!TEST!! !!if(ix==4 .and. iy==4) write(*,*)"T combine bgr" !!TEST!! call combine_runs(sig_point%matching_bgr_arr, it, ix, iy, bgr ) ! for later use instead of bgr_point%pixelbin if (bgr_point%pixelbin(it, ix, iy)%status /= PIXEL_OK ) then sig_point%pixelbin(it, ix, iy)%status = & Loading @@ -132,11 +132,11 @@ if(ix==4 .and. iy==4) write(*,*)"T combine bgr" !!TEST!! !!TEST!! fqt = get_fqt_bgr(sig_point%pixelbin(it, ix, iy), bgr_point%pixelbin(it, ix, iy), tfac, volume_fraction) !!TEST!! if(ix==4 .and. iy==4) write(*,*)"T fqt1_sngl(old)",fqt !!TEST!! !!if(ix==4 .and. iy==4) write(*,*)"T fqt1_sngl(old)",fqt !!TEST!! fqt = get_fqt_bgr(sig_point%pixelbin(it, ix, iy), bgr, tfac, volume_fraction) if(ix==4 .and. iy==4) write(*,*)"T fqt2_mult(new)",fqt !!TEST!! !!if(ix==4 .and. iy==4) write(*,*)"T fqt2_mult(new)",fqt !!TEST!! else fqt = sig_point%pixelbin(it, ix, iy)%fqt endif Loading Loading @@ -202,7 +202,7 @@ if(ix==4 .and. iy==4) write(*,*)"T fqt2_mult(new)",fqt !!TEST!! endif stat = PIXEL_OK if(ix==4 .and. iy==4) write(*,'(a,i4,i9,2(i3,2f18.9))') "Tcombine_runs i amp fqt: ",i,run_pointer_array(i)%ptr%id, amp,fqt !!TEST!! !if(ix==4 .and. iy==4) write(*,'(a,i4,i9,2(i3,2f18.9))') "Tcombine_runs i amp fqt: ",i,run_pointer_array(i)%ptr%id, amp,fqt !!TEST!! !!TEST!! if(ix==4 .and. iy==4) write(*,'(a,i4,i9,2(i3,2f18.9))') "Tcombine_runs o amp fqt: ",i,run_pointer_array(i)%ptr%id,&!!TEST!! !!TEST!! run_pointer_array(i)%ptr%pixelbin(it, ix, iy)%amplitude , run_pointer_array(i)%ptr%pixelbin(it, ix, iy)%fqt !!TEST!! !!TEST!! if(ix==4 .and. iy==4) write(*,*)"FQT ", 2*amp%value/(up%value-down%value), run_pointer_array(i)%ptr%pixelbin(it, ix, iy)%fqt !!TEST!! Loading Loading @@ -239,7 +239,7 @@ if(ix==4 .and. iy==4) write(*,'(a,i4,i9,2(i3,2f18.9))') "Tcombine_runs i amp fqt type(value_struct), dimension(:), allocatable :: bgr_pixel_bin character(len=MAX_NAME_LENGTH) :: cmsg !character(len=MAX_NAME_LENGTH) :: cmsg, (paz), unused variable sample_loop: do i=1, sample_data%number_of_points Loading Loading @@ -453,7 +453,7 @@ bgr_pixel_bin(l) integer , intent(in) :: scan_size real(kind=DBL) , intent(in) :: transmission ! integer :: i, j integer :: i !!, j !! (paz) unused variable ! source_loop: do i=1, scan_size if ( .not. is_valid_scan(scan_data(i)) ) cycle source_loop Loading