! Started logfile: /d2/scr/pima/rdv96_uvs/J1939-1525_X_uvs.log on Fri Dec 20 23:46:41 2013 obs /d2/scr/pima/rdv96_uvs/J1939-1525_X_uva.fits ! Reading UV FITS file: /d2/scr/pima/rdv96_uvs/J1939-1525_X_uva.fits ! AN table 1: 11 integrations on 55 of 55 possible baselines. ! AN table 2: 2 integrations on 28 of 28 possible baselines. ! AN table 3: 48 integrations on 3 of 3 possible baselines. ! AN table 4: 3 integrations on 3 of 3 possible baselines. ! AN table 5: 19 integrations on 6 of 6 possible baselines. ! AN table 6: 4 integrations on 15 of 15 possible baselines. ! AN table 7: 19 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.525773 visibilities/baseline/integration-bin. ! Found source: J1939-1525 ! ! 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 3876 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/J1939-1525_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:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HH self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 1:WF self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 3:HH self-cal status=correctable weight=1 ! 3:MA self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WF self-cal status=correctable weight=1 ! 5:HH self-cal status=correctable weight=1 ! 5:MA self-cal status=correctable weight=1 ! 5:WF self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:NL self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 6:WF self-cal status=correctable weight=1 ! 7:FD self-cal status=correctable weight=1 ! 7:HH self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:KP self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:MA self-cal status=correctable weight=1 ! 7:NL self-cal status=correctable weight=1 ! 7:SC self-cal status=correctable weight=1 ! 7:WF self-cal status=correctable weight=1 ! 7: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 67 telescope corrections were flagged in sub-array 1. ! A total of 42 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 5. ! A total of 78 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 67 telescope corrections were flagged in sub-array 1. ! A total of 42 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 5. ! A total of 78 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 67 telescope corrections were flagged in sub-array 1. ! A total of 42 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 5. ! A total of 78 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 67 telescope corrections were flagged in sub-array 1. ! A total of 42 telescope corrections were flagged in sub-array 3. ! A total of 13 telescope corrections were flagged in sub-array 5. ! A total of 78 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=1.235532Jy sigma=8.894492 ! Fit after self-cal, rms=0.414388Jy sigma=1.945517 ! 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 flag ! No telescope name matches 'flag'. ! Error occured in command: flag ![Exited script file: /tmp/ff.ff] ![Exited script file: /home/lpetrov/bin/rdv_muppet] ! wmap "/d2/scr/pima/rdv96_uvs/J1939-1525_X\_map.fits" ! Inverting map and beam ! Your choice of large map pixels excluded 8.19% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.5627 mas, bmaj=0.8952 mas, bpa=31.51 degrees ! Estimated noise=3.67692 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5627 x 0.8952 at 31.51 degrees (North through East) ! No model to restore with. ! Error occured in command: wmap ! float peak_flux peak_flux = peak(flux,max) invert ! uvinvert(): Neither beam nor a map requested. ! Error occured in command: invert ! 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.6694 rms= 0.03074 Jy quit ! Quitting program ! Log file /d2/scr/pima/rdv96_uvs/J1939-1525_X_uvs.log closed on Fri Dec 20 23:46:41 2013