! Started logfile: /d2/scr/pima/rdv96_uvs/J0038+4137_X_uvs.log on Fri Dec 20 23:43:16 2013 obs /d2/scr/pima/rdv96_uvs/J0038+4137_X_uva.fits ! Reading UV FITS file: /d2/scr/pima/rdv96_uvs/J0038+4137_X_uva.fits ! AN table 1: 10 integrations on 3 of 3 possible baselines. ! AN table 2: 29 integrations on 66 of 66 possible baselines. ! AN table 3: 7 integrations on 45 of 45 possible baselines. ! AN table 4: 2 integrations on 3 of 3 possible baselines. ! AN table 5: 3 integrations on 28 of 28 possible baselines. ! AN table 6: 14 integrations on 1 of 1 possible baselines. ! Apparent sampling: 0.570038 visibilities/baseline/integration-bin. ! Found source: J0038+4137 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40599e+09 8e+06 1 8e+06 ! 02 2 8.47599e+09 8e+06 1 8e+06 ! 03 3 8.79099e+09 8e+06 1 8e+06 ! 04 4 8.89599e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 5388 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/rdv_muppet /d2/scr/pima/rdv96_uvs/J0038+4137_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:HH self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WF self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:HH self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.200x0.200 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 153 telescope corrections were flagged in sub-array 2. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 153 telescope corrections were flagged in sub-array 2. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 193 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 193 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.051625Jy sigma=4.363776 ! Fit after self-cal, rms=0.718593Jy sigma=2.688495 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 1:HH ! Flagging baselines of antenna 1:HH in the currently selected channels. flag 2:HH ! Flagging baselines of antenna 2:HH in the currently selected channels. flag 2:NY ! Flagging baselines of antenna 2:NY in the currently selected channels. flag 2:WF ! Flagging baselines of antenna 2:WF in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 3:NY ! Flagging baselines of antenna 3:NY in the currently selected channels. flag 6:HH ! Flagging baselines of antenna 6:HH in the currently selected channels. flag 6:WZ ! Flagging baselines of antenna 6:WZ in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.122 mas, bmaj=1.532 mas, bpa=-49.07 degrees ! Estimated noise=4.2144 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.260331 Jy ! Component: 100 - total flux cleaned = 0.327255 Jy ! Total flux subtracted in 100 components = 0.327255 Jy ! Clean residual min=-0.009716 max=0.024644 Jy/beam ! Clean residual mean=0.000008 rms=0.002802 Jy/beam ! Combined flux in latest and established models = 0.327255 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 9 components and 0.327255 Jy ! ! Correcting IF 1. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 88 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 88 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 48 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 48 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.078388Jy sigma=0.412056 ! Fit after self-cal, rms=0.078256Jy sigma=0.413052 ! Inverting map and beam ! Estimated beam: bmin=1.124 mas, bmaj=1.527 mas, bpa=-48.56 degrees ! Estimated noise=4.23526 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.025742 Jy ! Component: 100 - total flux cleaned = 0.038255 Jy ! Total flux subtracted in 100 components = 0.038255 Jy ! Clean residual min=-0.007384 max=0.011936 Jy/beam ! Clean residual mean=0.000002 rms=0.001965 Jy/beam ! Combined flux in latest and established models = 0.36551 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 16 components and 0.36551 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.066802Jy sigma=0.338266 ! Fit after self-cal, rms=0.066866Jy sigma=0.338015 ! Inverting map ! Added new window around map position (3.8, -0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0116443 Jy ! Component: 100 - total flux cleaned = 0.0186456 Jy ! Total flux subtracted in 100 components = 0.0186456 Jy ! Clean residual min=-0.007451 max=0.007113 Jy/beam ! Clean residual mean=0.000002 rms=0.001744 Jy/beam ! Combined flux in latest and established models = 0.384156 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 24 components and 0.384156 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.063885Jy sigma=0.318537 ! Fit after self-cal, rms=0.063826Jy sigma=0.318073 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.397 mas, bmaj=1.873 mas, bpa=-57.23 degrees ! Estimated noise=2.65678 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0103826 Jy ! Component: 100 - total flux cleaned = 0.0166545 Jy ! Component: 150 - total flux cleaned = 0.0208315 Jy ! Component: 200 - total flux cleaned = 0.0236327 Jy ! Total flux subtracted in 200 components = 0.0236327 Jy ! Clean residual min=-0.007644 max=0.005834 Jy/beam ! Clean residual mean=0.000007 rms=0.001567 Jy/beam ! Combined flux in latest and established models = 0.407789 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 36 components and 0.407789 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.062127Jy sigma=0.306809 ! Fit after self-cal, rms=0.062042Jy sigma=0.306440 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* HN 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.96 FD 1.03 HH 1.00* HN 0.99 ! KP 1.02 LA 1.02 NL 0.96 NY 1.00* ! PT 0.99 SC 0.99 WF 1.00* WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00 FD 1.01 HN 1.02 KP 1.04 ! LA 1.00 MK 1.01 NL 1.00 NY 1.00* ! OV 0.96 PT 0.96 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.01 FD 1.05 HN 1.08 KP 1.01 ! LA 1.03 NL 0.97 OV 1.06 PT 1.03 ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* WZ 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* HN 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.96 FD 1.01 HH 1.00* HN 0.99 ! KP 0.96 LA 0.98 NL 0.97 NY 1.00* ! PT 0.99 SC 1.02 WF 1.00* WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00 FD 0.99 HN 1.01 KP 0.96 ! LA 0.99 MK 1.01 NL 0.99 NY 1.00* ! OV 0.97 PT 1.00 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.01 FD 1.00 HN 1.08 KP 0.98 ! LA 1.03 NL 0.97 OV 0.99 PT 1.04 ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* WZ 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* HN 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.00 HH 1.00* HN 1.02 ! KP 0.88 LA 1.05 NL 1.02 NY 1.00* ! PT 0.99 SC 0.95 WF 1.00* WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.01 FD 1.00 HN 1.07 KP 0.85 ! LA 1.05 MK 1.02 NL 1.01 NY 1.00* ! OV 1.02 PT 0.96 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.05 FD 1.02 HN 1.13 KP 0.85 ! LA 1.09 NL 1.04 OV 1.06 PT 1.05 ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* WZ 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* HN 0.00* SC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.06 HH 1.00* HN 1.02 ! KP 0.90 LA 1.03 NL 1.01 NY 1.00* ! PT 1.02 SC 1.00 WF 1.00* WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.05 FD 1.06 HN 1.05 KP 0.90 ! LA 1.01 MK 0.98 NL 1.04 NY 1.00* ! OV 1.02 PT 1.04 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.00 FD 1.12 HN 1.10 KP 0.93 ! LA 1.08 NL 1.01 OV 1.01 PT 1.10 ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* WZ 0.00* ! ! ! Fit before self-cal, rms=0.062042Jy sigma=0.306440 ! Fit after self-cal, rms=0.058031Jy sigma=0.280176 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.41 mas, bmaj=1.873 mas, bpa=-58.35 degrees ! Estimated noise=2.66508 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000496852 Jy ! Component: 100 - total flux cleaned = 0.000499266 Jy ! Component: 150 - total flux cleaned = 0.000498497 Jy ! Component: 200 - total flux cleaned = 0.000590663 Jy ! Total flux subtracted in 200 components = 0.000590663 Jy ! Clean residual min=-0.003448 max=0.004713 Jy/beam ! Clean residual mean=0.000008 rms=0.000961 Jy/beam ! Combined flux in latest and established models = 0.408379 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 51 components and 0.408379 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.057894Jy sigma=0.279457 ! Fit after self-cal, rms=0.057883Jy sigma=0.279331 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.057883Jy sigma=0.279331 ! Fit after self-cal, rms=0.057887Jy sigma=0.279301 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.411 mas, bmaj=1.874 mas, bpa=-58.34 degrees ! Estimated noise=2.66345 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000439487 Jy ! Component: 100 - total flux cleaned = -0.000295582 Jy ! Component: 150 - total flux cleaned = -0.000167483 Jy ! Component: 200 - total flux cleaned = -0.000126588 Jy ! Total flux subtracted in 200 components = -0.000126588 Jy ! Clean residual min=-0.003416 max=0.004680 Jy/beam ! Clean residual mean=0.000008 rms=0.000949 Jy/beam ! Combined flux in latest and established models = 0.408253 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 55 components and 0.408253 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.057795Jy sigma=0.278956 ! Fit after self-cal, rms=0.057760Jy sigma=0.278898 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-0.4, -36.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00420998 Jy ! Component: 100 - total flux cleaned = 0.00624059 Jy ! Component: 150 - total flux cleaned = 0.00688079 Jy ! Component: 200 - total flux cleaned = 0.00714482 Jy ! Total flux subtracted in 200 components = 0.00714482 Jy ! Clean residual min=-0.003289 max=0.003432 Jy/beam ! Clean residual mean=0.000003 rms=0.000818 Jy/beam ! Combined flux in latest and established models = 0.415397 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 72 components and 0.415397 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.057400Jy sigma=0.277200 ! Fit after self-cal, rms=0.057374Jy sigma=0.277119 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.057374Jy sigma=0.277119 ! Fit after self-cal, rms=0.057363Jy sigma=0.276993 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.412 mas, bmaj=1.876 mas, bpa=-58.33 degrees ! Estimated noise=2.66218 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000737034 Jy ! Component: 100 - total flux cleaned = 0.00108202 Jy ! Component: 150 - total flux cleaned = 0.0012137 Jy ! Component: 200 - total flux cleaned = 0.00133522 Jy ! Total flux subtracted in 200 components = 0.00133522 Jy ! Clean residual min=-0.002975 max=0.003500 Jy/beam ! Clean residual mean=0.000001 rms=0.000789 Jy/beam ! Combined flux in latest and established models = 0.416733 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 81 components and 0.416733 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.057269Jy sigma=0.276645 ! Fit after self-cal, rms=0.057251Jy sigma=0.276621 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-38.8, 46.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00352185 Jy ! Component: 100 - total flux cleaned = 0.00469685 Jy ! Component: 150 - total flux cleaned = 0.00548951 Jy ! Component: 200 - total flux cleaned = 0.00617613 Jy ! Total flux subtracted in 200 components = 0.00617613 Jy ! Clean residual min=-0.003103 max=0.003089 Jy/beam ! Clean residual mean=-0.000000 rms=0.000709 Jy/beam ! Combined flux in latest and established models = 0.422909 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 98 components and 0.422909 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.057001Jy sigma=0.275359 ! Fit after self-cal, rms=0.056989Jy sigma=0.275307 ! Inverting map ! Added new window around map position (34.4, 27.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00283713 Jy ! Component: 100 - total flux cleaned = 0.00407932 Jy ! Component: 150 - total flux cleaned = 0.00478143 Jy ! Component: 200 - total flux cleaned = 0.00500257 Jy ! Total flux subtracted in 200 components = 0.00500257 Jy ! Clean residual min=-0.002571 max=0.002729 Jy/beam ! Clean residual mean=0.000005 rms=0.000641 Jy/beam ! Combined flux in latest and established models = 0.427911 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 114 components and 0.427911 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056804Jy sigma=0.274503 ! Fit after self-cal, rms=0.056799Jy sigma=0.274469 ! Inverting map ! Added new window around map position (-3.8, 6.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00144829 Jy ! Component: 100 - total flux cleaned = 0.00120133 Jy ! Component: 150 - total flux cleaned = 0.0010374 Jy ! Component: 200 - total flux cleaned = 0.00118558 Jy ! Total flux subtracted in 200 components = 0.00118558 Jy ! Clean residual min=-0.002557 max=0.002524 Jy/beam ! Clean residual mean=0.000003 rms=0.000630 Jy/beam ! Combined flux in latest and established models = 0.429097 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 127 components and 0.429097 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056598Jy sigma=0.273878 ! Fit after self-cal, rms=0.056596Jy sigma=0.273866 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056596Jy sigma=0.273866 ! Fit after self-cal, rms=0.056549Jy sigma=0.273616 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.673 mas, bmaj=3.199 mas, bpa=-76.83 degrees ! Estimated noise=2.85883 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000786711 Jy ! Component: 100 - total flux cleaned = 0.00130638 Jy ! Component: 150 - total flux cleaned = 0.00168768 Jy ! Component: 200 - total flux cleaned = 0.00205619 Jy ! Total flux subtracted in 200 components = 0.00205619 Jy ! Clean residual min=-0.002143 max=0.002411 Jy/beam ! Clean residual mean=0.000004 rms=0.000646 Jy/beam ! Combined flux in latest and established models = 0.431153 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 153 components and 0.431153 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056526Jy sigma=0.273407 ! Fit after self-cal, rms=0.056520Jy sigma=0.273392 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 9.02498e-05 Jy ! Component: 100 - total flux cleaned = 0.00028803 Jy ! Component: 150 - total flux cleaned = 0.000415174 Jy ! Component: 200 - total flux cleaned = 0.000462424 Jy ! Total flux subtracted in 200 components = 0.000462424 Jy ! Clean residual min=-0.002074 max=0.002348 Jy/beam ! Clean residual mean=0.000005 rms=0.000640 Jy/beam ! Combined flux in latest and established models = 0.431616 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 162 components and 0.431615 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056509Jy sigma=0.273290 ! Fit after self-cal, rms=0.056504Jy sigma=0.273280 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] flag 1:HH ! Flagging baselines of antenna 1:HH in the currently selected channels. flag 2:HH ! Flagging baselines of antenna 2:HH in the currently selected channels. ! Applying 2560 buffered edits. flag 2:NY ! Flagging baselines of antenna 2:NY in the currently selected channels. flag 2:WF ! Flagging baselines of antenna 2:WF in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 3:NY ! Flagging baselines of antenna 3:NY in the currently selected channels. flag 6:HH ! Flagging baselines of antenna 6:HH in the currently selected channels. flag 6:WZ ! Flagging baselines of antenna 6:WZ in the currently selected channels. ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 241 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 241 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 241 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 241 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.056535Jy sigma=0.272245 ! Fit after self-cal, rms=0.056526Jy sigma=0.273253 selfant "",false ! 1:HH self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WF self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:HH self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056526Jy sigma=0.273253 ! Fit after self-cal, rms=0.056526Jy sigma=0.273252 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.149 mas, bmaj=1.504 mas, bpa=-48.45 degrees ! Estimated noise=4.30132 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (13.8, 34.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00220537 Jy ! Component: 100 - total flux cleaned = 0.00208443 Jy ! Component: 150 - total flux cleaned = 0.00228749 Jy ! Component: 200 - total flux cleaned = 0.00276162 Jy ! Total flux subtracted in 200 components = 0.00276162 Jy ! Clean residual min=-0.003408 max=0.003581 Jy/beam ! Clean residual mean=0.000000 rms=0.000814 Jy/beam ! Combined flux in latest and established models = 0.434377 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 186 components and 0.434377 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056375Jy sigma=0.272842 ! Fit after self-cal, rms=0.056378Jy sigma=0.272788 ! Inverting map ! Added new window around map position (-19.6, -12.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00202364 Jy ! Component: 100 - total flux cleaned = 0.00229461 Jy ! Component: 150 - total flux cleaned = 0.00246002 Jy ! Component: 200 - total flux cleaned = 0.0024605 Jy ! Total flux subtracted in 200 components = 0.0024605 Jy ! Clean residual min=-0.003055 max=0.003251 Jy/beam ! Clean residual mean=-0.000000 rms=0.000772 Jy/beam ! Combined flux in latest and established models = 0.436838 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 203 components and 0.436838 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056313Jy sigma=0.272702 ! Fit after self-cal, rms=0.056314Jy sigma=0.272667 ! Inverting map ! Added new window around map position (-27, 17.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0014357 Jy ! Component: 100 - total flux cleaned = 0.00192216 Jy ! Component: 150 - total flux cleaned = 0.00245915 Jy ! Component: 200 - total flux cleaned = 0.00260457 Jy ! Total flux subtracted in 200 components = 0.00260457 Jy ! Clean residual min=-0.002955 max=0.003155 Jy/beam ! Clean residual mean=0.000000 rms=0.000731 Jy/beam ! Combined flux in latest and established models = 0.439442 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 217 components and 0.439442 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056228Jy sigma=0.272488 ! Fit after self-cal, rms=0.056227Jy sigma=0.272462 ! Inverting map ! Added new window around map position (-47.4, -32.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0019675 Jy ! Component: 100 - total flux cleaned = 0.00181171 Jy ! Component: 150 - total flux cleaned = 0.00181342 Jy ! Component: 200 - total flux cleaned = 0.00188159 Jy ! Total flux subtracted in 200 components = 0.00188159 Jy ! Clean residual min=-0.002955 max=0.002891 Jy/beam ! Clean residual mean=-0.000000 rms=0.000703 Jy/beam ! Combined flux in latest and established models = 0.441324 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 241 components and 0.441324 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.056173Jy sigma=0.272288 ! Fit after self-cal, rms=0.056172Jy sigma=0.272267 ! Inverting map uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.414 mas, bmaj=1.877 mas, bpa=-58.04 degrees ! Estimated noise=2.6605 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00167093 Jy ! Component: 100 - total flux cleaned = -0.00241922 Jy ! Component: 150 - total flux cleaned = -0.0029837 Jy ! Component: 200 - total flux cleaned = -0.00341083 Jy ! Component: 250 - total flux cleaned = -0.00376573 Jy ! Component: 300 - total flux cleaned = -0.0040564 Jy ! Component: 350 - total flux cleaned = -0.00421517 Jy ! Component: 400 - total flux cleaned = -0.00436679 Jy ! Total flux subtracted in 400 components = -0.00436679 Jy ! Clean residual min=-0.002238 max=0.002171 Jy/beam ! Clean residual mean=0.000003 rms=0.000542 Jy/beam ! Combined flux in latest and established models = 0.436957 Jy ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 282 components and 0.436957 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.055992Jy sigma=0.271575 ! Fit after self-cal, rms=0.055986Jy sigma=0.271535 ! Inverting map dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000682742 Jy ! Component: 100 - total flux cleaned = -0.000720205 Jy ! Component: 150 - total flux cleaned = -0.000791608 Jy ! Component: 200 - total flux cleaned = -0.000825366 Jy ! Component: 250 - total flux cleaned = -0.00082531 Jy ! Component: 300 - total flux cleaned = -0.000857763 Jy ! Component: 350 - total flux cleaned = -0.000888329 Jy ! Component: 400 - total flux cleaned = -0.000858281 Jy ! Total flux subtracted in 400 components = -0.000858281 Jy ! Clean residual min=-0.002107 max=0.002044 Jy/beam ! Clean residual mean=0.000003 rms=0.000520 Jy/beam ! Combined flux in latest and established models = 0.436099 Jy ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 294 components and 0.436099 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.055896Jy sigma=0.271219 ! Fit after self-cal, rms=0.055895Jy sigma=0.271209 ! Inverting map print "=== Finished natural weighting clean for al sites ===" ! === Finished natural weighting clean for al sites === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 177 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 177 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 196 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 196 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.119378Jy sigma=0.265978 ! Fit after self-cal, rms=0.106928Jy sigma=0.246796 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! Estimated beam: bmin=1.008 mas, bmaj=1.271 mas, bpa=49.42 degrees ! Estimated noise=2.78942 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00187159 Jy ! Component: 100 - total flux cleaned = 0.00212118 Jy ! Component: 150 - total flux cleaned = 0.00198286 Jy ! Component: 200 - total flux cleaned = 0.00191411 Jy ! Component: 250 - total flux cleaned = 0.00172126 Jy ! Component: 300 - total flux cleaned = 0.00153871 Jy ! Component: 350 - total flux cleaned = 0.00136261 Jy ! Component: 400 - total flux cleaned = 0.00124988 Jy ! Total flux subtracted in 400 components = 0.00124988 Jy ! Clean residual min=-0.002726 max=0.002907 Jy/beam ! Clean residual mean=0.000071 rms=0.000582 Jy/beam ! Combined flux in latest and established models = 0.437349 Jy selfcal ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 339 components and 0.437349 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106575Jy sigma=0.246113 ! Fit after self-cal, rms=0.106013Jy sigma=0.245221 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (12.2, -9). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00280119 Jy ! Component: 100 - total flux cleaned = 0.00275651 Jy ! Component: 150 - total flux cleaned = 0.00245556 Jy ! Component: 200 - total flux cleaned = 0.00222694 Jy ! Component: 250 - total flux cleaned = 0.00200766 Jy ! Component: 300 - total flux cleaned = 0.00190216 Jy ! Component: 350 - total flux cleaned = 0.00190211 Jy ! Component: 400 - total flux cleaned = 0.00185299 Jy ! Total flux subtracted in 400 components = 0.00185299 Jy ! Clean residual min=-0.002535 max=0.002415 Jy/beam ! Clean residual mean=0.000064 rms=0.000533 Jy/beam ! Combined flux in latest and established models = 0.439202 Jy ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 379 components and 0.439202 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.105728Jy sigma=0.244652 ! Fit after self-cal, rms=0.105700Jy sigma=0.244613 ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.105700Jy sigma=0.244613 ! Fit after self-cal, rms=0.105365Jy sigma=0.244483 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! Estimated beam: bmin=1.007 mas, bmaj=1.27 mas, bpa=49.4 degrees ! Estimated noise=2.78916 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000907292 Jy ! Component: 100 - total flux cleaned = 0.000907584 Jy ! Component: 150 - total flux cleaned = 0.000855848 Jy ! Component: 200 - total flux cleaned = 0.000760459 Jy ! Component: 250 - total flux cleaned = 0.000807346 Jy ! Component: 300 - total flux cleaned = 0.000715716 Jy ! Component: 350 - total flux cleaned = 0.000805042 Jy ! Component: 400 - total flux cleaned = 0.000762344 Jy ! Total flux subtracted in 400 components = 0.000762344 Jy ! Clean residual min=-0.002434 max=0.002335 Jy/beam ! Clean residual mean=0.000051 rms=0.000508 Jy/beam ! Combined flux in latest and established models = 0.439964 Jy selfcal ! Performing phase self-cal ! Adding 77 model components to the UV plane model. ! The established model now contains 398 components and 0.439964 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.105175Jy sigma=0.244189 ! Fit after self-cal, rms=0.105151Jy sigma=0.244177 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 4.74999e-05 Jy ! Component: 100 - total flux cleaned = 4.73192e-05 Jy ! Component: 150 - total flux cleaned = 4.67313e-05 Jy ! Component: 200 - total flux cleaned = 4.7036e-05 Jy ! Component: 250 - total flux cleaned = 4.7206e-05 Jy ! Component: 300 - total flux cleaned = 7.17003e-06 Jy ! Component: 350 - total flux cleaned = 7.68824e-06 Jy ! Component: 400 - total flux cleaned = 7.81655e-06 Jy ! Total flux subtracted in 400 components = 7.81655e-06 Jy ! Clean residual min=-0.002407 max=0.002313 Jy/beam ! Clean residual mean=0.000050 rms=0.000499 Jy/beam ! Combined flux in latest and established models = 0.439972 Jy ! Performing phase self-cal ! Adding 92 model components to the UV plane model. ! The established model now contains 415 components and 0.439972 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.104999Jy sigma=0.243980 ! Fit after self-cal, rms=0.104976Jy sigma=0.243975 ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime3 ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.104976Jy sigma=0.243975 ! Fit after self-cal, rms=0.104999Jy sigma=0.240543 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! Estimated beam: bmin=2.561 mas, bmaj=3.132 mas, bpa=87.22 degrees ! Estimated noise=2.76045 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000405737 Jy ! Component: 100 - total flux cleaned = 0.000315158 Jy ! Component: 150 - total flux cleaned = 0.000205957 Jy ! Component: 200 - total flux cleaned = 5.46659e-05 Jy ! Component: 250 - total flux cleaned = -1.57015e-05 Jy ! Component: 300 - total flux cleaned = -8.34078e-05 Jy ! Component: 350 - total flux cleaned = -0.000189829 Jy ! Component: 400 - total flux cleaned = -0.000210384 Jy ! Total flux subtracted in 400 components = -0.000210384 Jy ! Clean residual min=-0.002032 max=0.002046 Jy/beam ! Clean residual mean=0.000077 rms=0.000516 Jy/beam ! Combined flux in latest and established models = 0.439762 Jy selfcal ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 446 components and 0.439762 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.104983Jy sigma=0.240447 ! Fit after self-cal, rms=0.104994Jy sigma=0.240435 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-51.2, -11). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00191947 Jy ! Component: 100 - total flux cleaned = 0.00266774 Jy ! Component: 150 - total flux cleaned = 0.00308559 Jy ! Component: 200 - total flux cleaned = 0.00336282 Jy ! Component: 250 - total flux cleaned = 0.00358965 Jy ! Component: 300 - total flux cleaned = 0.00382726 Jy ! Component: 350 - total flux cleaned = 0.00407904 Jy ! Component: 400 - total flux cleaned = 0.00436704 Jy ! Total flux subtracted in 400 components = 0.00436704 Jy ! Clean residual min=-0.001884 max=0.001600 Jy/beam ! Clean residual mean=0.000045 rms=0.000475 Jy/beam ! Combined flux in latest and established models = 0.444129 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 487 components and 0.444129 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.104931Jy sigma=0.240118 ! Fit after self-cal, rms=0.104930Jy sigma=0.240104 ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for all stations ===" ! === Finished amplitude self-cal for all stations === restore ! Inverting map and beam ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! Estimated beam: bmin=1.007 mas, bmaj=1.27 mas, bpa=49.38 degrees ! Estimated noise=2.79122 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.007 x 1.27 at 49.38 degrees (North through East) ! Clean map min=-0.0022535 max=0.31399 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if print "=== CLearing model and starting over ===" ! === CLearing model and starting over === clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! Estimated beam: bmin=0.7733 mas, bmaj=1.025 mas, bpa=46.52 degrees ! Estimated noise=4.94202 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.257361 Jy ! Component: 100 - total flux cleaned = 0.327919 Jy ! Component: 150 - total flux cleaned = 0.354549 Jy ! Component: 200 - total flux cleaned = 0.369172 Jy ! Component: 250 - total flux cleaned = 0.379214 Jy ! Component: 300 - total flux cleaned = 0.386876 Jy ! Component: 350 - total flux cleaned = 0.393255 Jy ! Component: 400 - total flux cleaned = 0.39821 Jy ! Total flux subtracted in 400 components = 0.39821 Jy ! Clean residual min=-0.004528 max=0.006590 Jy/beam ! Clean residual mean=0.000032 rms=0.001140 Jy/beam ! Combined flux in latest and established models = 0.39821 Jy uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Adding 52 model components to the UV plane model. ! The established model now contains 52 components and 0.39821 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! Estimated beam: bmin=1.007 mas, bmaj=1.27 mas, bpa=49.38 degrees ! Estimated noise=2.79122 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (1.8, -0.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00805688 Jy ! Component: 100 - total flux cleaned = 0.0131607 Jy ! Component: 150 - total flux cleaned = 0.0171828 Jy ! Component: 200 - total flux cleaned = 0.0206286 Jy ! Component: 250 - total flux cleaned = 0.0236689 Jy ! Component: 300 - total flux cleaned = 0.026425 Jy ! Component: 350 - total flux cleaned = 0.0286641 Jy ! Component: 400 - total flux cleaned = 0.0301948 Jy ! Total flux subtracted in 400 components = 0.0301948 Jy ! Clean residual min=-0.002949 max=0.002926 Jy/beam ! Clean residual mean=0.000050 rms=0.000640 Jy/beam ! Combined flux in latest and established models = 0.428405 Jy ! Adding 69 model components to the UV plane model. ! The established model now contains 112 components and 0.428405 Jy ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! Estimated beam: bmin=2.561 mas, bmaj=3.132 mas, bpa=87.22 degrees ! Estimated noise=2.76045 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00258346 Jy ! Component: 100 - total flux cleaned = 0.00446728 Jy ! Component: 150 - total flux cleaned = 0.00596667 Jy ! Component: 200 - total flux cleaned = 0.00719783 Jy ! Component: 250 - total flux cleaned = 0.00821852 Jy ! Component: 300 - total flux cleaned = 0.00917561 Jy ! Component: 350 - total flux cleaned = 0.0100245 Jy ! Component: 400 - total flux cleaned = 0.010781 Jy ! Total flux subtracted in 400 components = 0.010781 Jy ! Clean residual min=-0.002211 max=0.002566 Jy/beam ! Clean residual mean=0.000058 rms=0.000580 Jy/beam ! Combined flux in latest and established models = 0.439186 Jy flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (10.4, -4.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00253452 Jy ! Component: 100 - total flux cleaned = 0.00334592 Jy ! Component: 150 - total flux cleaned = 0.00411149 Jy ! Component: 200 - total flux cleaned = 0.0046248 Jy ! Component: 250 - total flux cleaned = 0.00505994 Jy ! Component: 300 - total flux cleaned = 0.00534895 Jy ! Component: 350 - total flux cleaned = 0.00566308 Jy ! Component: 400 - total flux cleaned = 0.00588987 Jy ! Total flux subtracted in 400 components = 0.00588987 Jy ! Clean residual min=-0.001930 max=0.001749 Jy/beam ! Clean residual mean=0.000019 rms=0.000459 Jy/beam ! Combined flux in latest and established models = 0.445076 Jy ! Adding 141 model components to the UV plane model. ! The established model now contains 248 components and 0.445076 Jy ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! Estimated beam: bmin=1.007 mas, bmaj=1.27 mas, bpa=49.38 degrees ! Estimated noise=2.79122 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000862039 Jy ! Component: 100 - total flux cleaned = 0.00160152 Jy ! Component: 150 - total flux cleaned = 0.00208963 Jy ! Component: 200 - total flux cleaned = 0.00242491 Jy ! Component: 250 - total flux cleaned = 0.00268541 Jy ! Component: 300 - total flux cleaned = 0.00287469 Jy ! Component: 350 - total flux cleaned = 0.00305928 Jy ! Component: 400 - total flux cleaned = 0.00317899 Jy ! Total flux subtracted in 400 components = 0.00317899 Jy ! Clean residual min=-0.002262 max=0.002205 Jy/beam ! Clean residual mean=0.000010 rms=0.000474 Jy/beam ! Combined flux in latest and established models = 0.448255 Jy selfcal ! Performing phase self-cal ! Adding 60 model components to the UV plane model. ! The established model now contains 294 components and 0.448255 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.105092Jy sigma=0.240634 ! Fit after self-cal, rms=0.104538Jy sigma=0.240319 wmodel /d2/scr/pima/rdv96_uvs/J0038+4137_X_map.mod ! Writing 294 model components to file: /d2/scr/pima/rdv96_uvs/J0038+4137_X_map.mod wobs /d2/scr/pima/rdv96_uvs/J0038+4137_X_uvs.fits ! Applying 2488 buffered edits. ! Writing UV FITS file: /d2/scr/pima/rdv96_uvs/J0038+4137_X_uvs.fits wwins /d2/scr/pima/rdv96_uvs/J0038+4137_X_map.win ! wwins: Wrote 14 windows to /d2/scr/pima/rdv96_uvs/J0038+4137_X_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000530278 Jy ! Component: 100 - total flux cleaned = -0.000744209 Jy ! Component: 150 - total flux cleaned = -0.00101689 Jy ! Component: 200 - total flux cleaned = -0.00124793 Jy ! Component: 250 - total flux cleaned = -0.00144155 Jy ! Component: 300 - total flux cleaned = -0.00159855 Jy ! Component: 350 - total flux cleaned = -0.00159901 Jy ! Component: 400 - total flux cleaned = -0.00162967 Jy ! Component: 450 - total flux cleaned = -0.00163014 Jy ! Component: 500 - total flux cleaned = -0.00160136 Jy ! Component: 550 - total flux cleaned = -0.00154389 Jy ! Component: 600 - total flux cleaned = -0.00137415 Jy ! Component: 650 - total flux cleaned = -0.00131822 Jy ! Component: 700 - total flux cleaned = -0.001235 Jy ! Component: 750 - total flux cleaned = -0.00109831 Jy ! Component: 800 - total flux cleaned = -0.00101737 Jy ! Component: 850 - total flux cleaned = -0.000856544 Jy ! Component: 900 - total flux cleaned = -0.000750142 Jy ! Component: 950 - total flux cleaned = -0.000592639 Jy ! Component: 1000 - total flux cleaned = -0.000566306 Jy ! Component: 1050 - total flux cleaned = -0.000411393 Jy ! Component: 1100 - total flux cleaned = -0.000257869 Jy ! Component: 1150 - total flux cleaned = -0.000232428 Jy ! Component: 1200 - total flux cleaned = -8.11476e-05 Jy ! Component: 1250 - total flux cleaned = 4.38211e-05 Jy ! Component: 1300 - total flux cleaned = 0.000143264 Jy ! Component: 1350 - total flux cleaned = 0.000192626 Jy ! Component: 1400 - total flux cleaned = 0.000217057 Jy ! Component: 1450 - total flux cleaned = 0.000265783 Jy ! Component: 1500 - total flux cleaned = 0.000386424 Jy ! Component: 1550 - total flux cleaned = 0.000458525 Jy ! Component: 1600 - total flux cleaned = 0.000529878 Jy ! Component: 1650 - total flux cleaned = 0.000529875 Jy ! Component: 1700 - total flux cleaned = 0.000600491 Jy ! Component: 1750 - total flux cleaned = 0.000647167 Jy ! Component: 1800 - total flux cleaned = 0.00071686 Jy ! Component: 1850 - total flux cleaned = 0.000763246 Jy ! Component: 1900 - total flux cleaned = 0.000809126 Jy ! Component: 1950 - total flux cleaned = 0.000809061 Jy ! Component: 2000 - total flux cleaned = 0.000831818 Jy ! Component: 2050 - total flux cleaned = 0.000990062 Jy ! Component: 2100 - total flux cleaned = 0.000967728 Jy ! Component: 2150 - total flux cleaned = 0.00110177 Jy ! Component: 2200 - total flux cleaned = 0.00105735 Jy ! Component: 2250 - total flux cleaned = 0.00119007 Jy ! Component: 2300 - total flux cleaned = 0.00114602 Jy ! Component: 2350 - total flux cleaned = 0.0012555 Jy ! Component: 2400 - total flux cleaned = 0.00127751 Jy ! Component: 2450 - total flux cleaned = 0.00149452 Jy ! Component: 2500 - total flux cleaned = 0.00147311 Jy ! Component: 2550 - total flux cleaned = 0.00160209 Jy ! Component: 2600 - total flux cleaned = 0.00173031 Jy ! Component: 2650 - total flux cleaned = 0.00170899 Jy ! Component: 2700 - total flux cleaned = 0.00173025 Jy ! Component: 2750 - total flux cleaned = 0.00179369 Jy ! Component: 2800 - total flux cleaned = 0.00187774 Jy ! Component: 2850 - total flux cleaned = 0.00187779 Jy ! Component: 2900 - total flux cleaned = 0.00198189 Jy ! Component: 2950 - total flux cleaned = 0.00202342 Jy ! Component: 3000 - total flux cleaned = 0.00200267 Jy ! Component: 3050 - total flux cleaned = 0.0020644 Jy ! Component: 3100 - total flux cleaned = 0.0021054 Jy ! Component: 3150 - total flux cleaned = 0.00212591 Jy ! Component: 3200 - total flux cleaned = 0.00214617 Jy ! Component: 3250 - total flux cleaned = 0.00210566 Jy ! Component: 3300 - total flux cleaned = 0.0021258 Jy ! Component: 3350 - total flux cleaned = 0.00224609 Jy ! Component: 3400 - total flux cleaned = 0.00218603 Jy ! Component: 3450 - total flux cleaned = 0.00216616 Jy ! Component: 3500 - total flux cleaned = 0.00218614 Jy ! Component: 3550 - total flux cleaned = 0.00222568 Jy ! Component: 3600 - total flux cleaned = 0.00222567 Jy ! Component: 3650 - total flux cleaned = 0.00224527 Jy ! Component: 3700 - total flux cleaned = 0.00224545 Jy ! Component: 3750 - total flux cleaned = 0.00230386 Jy ! Component: 3800 - total flux cleaned = 0.00236216 Jy ! Component: 3850 - total flux cleaned = 0.00232355 Jy ! Component: 3900 - total flux cleaned = 0.00232361 Jy ! Component: 3950 - total flux cleaned = 0.00238106 Jy ! Component: 4000 - total flux cleaned = 0.00238105 Jy ! Component: 4050 - total flux cleaned = 0.0024383 Jy ! Total flux subtracted in 4096 components = 0.00245735 Jy ! Clean residual min=-0.001121 max=0.001148 Jy/beam ! Clean residual mean=0.000011 rms=0.000303 Jy/beam ! Combined flux in latest and established models = 0.450712 Jy keep ! Adding 277 model components to the UV plane model. ! The established model now contains 569 components and 0.450712 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/d2/scr/pima/rdv96_uvs/J0038+4137_X\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.007 x 1.27 at 49.38 degrees (North through East) ! Clean map min=-0.0017118 max=0.31377 Jy/beam ! Writing clean map to FITS file: /d2/scr/pima/rdv96_uvs/J0038+4137_X_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map ! Your choice of large map pixels excluded 2.33% of the data. ! The x-axis pixel size should ideally be below 0.1932 milli-arcsec ! The y-axis pixel size should ideally be below 0.1774 milli-arcsec float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 0.3137 rms= 0.0003 Jy quit ! Quitting program ! Log file /d2/scr/pima/rdv96_uvs/J0038+4137_X_uvs.log closed on Fri Dec 20 23:43:26 2013