Full cube for observations on 20080616, 20080617, 20080619 & 20081007. makecube in="a2008*_flaged.sdf" autogrid=true crota=325 pixsize=7.2761 sparse='nearest' specbounds=\"270.0 1040.0\" out=ngc3198_2008n #Projection parameters used: # CRPIX1 = 21.9102040827274 # CRPIX2 = 33.067346945405 # CRVAL1 = 154.979166666667 ( RA = 10:19:55.0 ) # CRVAL2 = 45.5497222222222 ( Dec = 45:32:59 ) # CDELT1 = -0.00202113893296984 ( -7.3 arcsec ) # CDELT2 = 0.00202113893296984 ( 7.3 arcsec ) # CROTA2 = 325 # Output cube pixel bounds: ( -21:22, -31:34, -907:905 ) # Output cube WCS bounds: # Right ascension: 10:19:29.8 -> 10:20:21.0 # Declination: 45:28:11 -> 45:37:48 # Optical velocity (USB): 269.6215 -> 1040.409 km/s mfittrend ngc3198_2008n.sdf ranges=\"270 480 840 1040\" order=3 axis=3 out=ngc3198_2008n_baseline # Fitting NDF Axis 3, using ranges in current co-ordinate Frame (pixels): # 270.0456 : 479.6758 km/s (-907 : -413) # 839.785 : 1039.982 km/s (434 : 904) sub ngc3198_2008n.sdf ngc3198_2008n_baseline.sdf ngc3198_2008n_sub compave ngc3198_2008n_sub.sdf ngc3198_2008n_sub_20kms compress="[1,1,47]" stats 'ngc3198_2008n_sub_20kms.sdf(~15,~15,:)' | grep Standard # Standard deviation : 0.0218001 stats 'ngc3198_2008n_sub_20kms.sdf(~15,~15,270.0:480.0)' | grep Standard # Standard deviation : 0.0174634 stats 'ngc3198_2008n_sub_20kms.sdf(~15,~15,840.0:1040.0)' | grep Standard # Standard deviation : 0.013771 ---------------------------------------------------------------------------------------- SINCSINC makecube in="a*" autogrid=true crota=325 pixsize=7.2761 sparse=false spread='sincsinc' params='[1,2]' specbounds=\"270.0 1040.0\" out=ngc3198_2008n_sincsinc # Projection parameters used: # CRPIX1 = 21.9102040827274 # CRPIX2 = 33.067346945405 # CRVAL1 = 154.979166666667 ( RA = 10:19:55.0 ) # CRVAL2 = 45.5497222222222 ( Dec = 45:32:59 ) # CDELT1 = -0.00202113893296984 ( -7.3 arcsec ) # CDELT2 = 0.00202113893296984 ( 7.3 arcsec ) # CROTA2 = 325 # # Output cube pixel bounds: ( -21:22, -31:34, -907:905 ) # Output cube WCS bounds: # Right ascension: 10:19:29.8 -> 10:20:21.0 # Declination: 45:28:11 -> 45:37:48 # Optical velocity (USB): 269.6215 -> 1040.409 km/s # # #WARNING: Cannot create output Tsys values. ndfcopy in=ngc3198_2008n_sincsinc'(-6:12,-19:25,)' out=ngc3198_2008n_sincsinc_trim mfittrend ngc3198_2008n_sincsinc_trim ranges=\"270 480 840 1040\" order=3 axis=3 out=ngc3198_2008n_sincsinc_trim_base3 # Fitting NDF Axis 3, using ranges in current co-ordinate Frame (pixels): # 270.0456 : 479.6758 km/s (-907 : -413) # 839.785 : 1039.982 km/s (434 : 904) sub ngc3198_2008n_sincsinc_trim ngc3198_2008n_sincsinc_trim_base3.sdf ngc3198_2008n_sincsinc_trim_sub3 compave ngc3198_2008n_sincsinc_trim_sub3.sdf ngc3198_2008n_sincsinc_trim_sub3_20kms compress="[1,1,47]" stats 'ngc3198_2008n_sincsinc_trim_sub3_20kms(~15,~15,:)' | grep Standard # Standard deviation : 0.0193093 stats 'ngc3198_2008n_sincsinc_trim_sub3_20kms(~15,~15,270.0:480.0)' | grep Standard # Standard deviation : 0.014679 stats 'ngc3198_2008n_sincsinc_trim_sub3_20kms(~15,~15,840.0:1040.0)' | grep Standard # Standard deviation : 0.0114971 --------------------------------------------------------------------------------------- BASELINE MASK PRODUCTION block in=ngc3198_2008n_sincsinc_trim.sdf out=ngc3198_2008n_sincsinc_trim_vsmooth 'axes=[1,3] box=[1,25]' block in=ngc3198_2008n_sincsinc_trim_vsmooth out=ngc3198_2008n_sincsinc_trim_smooth 'axes=[1,2] box=[3,3]' mfittrend in=ngc3198_2008n_sincsinc_trim_smooth axis=3 order=3 out=ngc3198_2008n_sincsinc_trim_smooth_base3 auto mask=ngc3198_2008n_sincsinc_trim_smooth_basemsk3 method=single variance=false subtract=false 'clip=[2,3]' mfittrend in=ngc3198_2008n_sincsinc_trim_smooth axis=3 order=3 out=ngc3198_2008n_sincsinc_trim_smooth_base3g auto mask=ngc3198_2008n_sincsinc_trim_smooth_basemsk3g method=global variance=false subtract=false 'clip=[2,2.5,3]' add ngc3198_2008n_sincsinc_trim_smooth_basemsk3 ngc3198_2008n_sincsinc_trim_smooth_basemsk3g ngc3198_2008n_trim_basemask0 chpix in=ngc3198_2008n_trim_basemask0 out=ngc3198_2008n_trim_basemask0_end1 section=\":,:,-907:-413\" newval=0 chpix in=ngc3198_2008n_trim_basemask0_end1 out=ngc3198_2008n_trim_basemask0_end2 section=\":,:,434:904\" newval=0 cadd ngc3198_2008n_trim_basemask0_end2 1 ngc3198_2008n_trim_basemask1 mult ngc3198_2008n_sincsinc_trim ngc3198_2008n_trim_basemask1 ngc3198_2008n_sincsinc_trim_masked mfittrend ngc3198_2008n_sincsinc_trim_masked ranges=\"270 1040\" order=3 axis=3 out=ngc3198_2008n_sincsinc_trim_masked_base3 #Fitting NDF Axis 3, using ranges in current co-ordinate Frame (pixels): # 270.0456 : 1039.982 km/s (-907 : 904) sub ngc3198_2008n_sincsinc_trim ngc3198_2008n_sincsinc_trim_masked_base3 ngc3198_2008n_sincsinc_trim_submask3 compave ngc3198_2008n_sincsinc_trim_submask3 ngc3198_2008n_sincsinc_trim_submask3_20kms compress="[1,1,47]" stats 'ngc3198_2008n_sincsinc_trim_submask3_20kms(~15,~15,:)' | grep Standard # Standard deviation : 0.0150699 stats 'ngc3198_2008n_sincsinc_trim_submask3_20kms(~15,~15,240.0:480.0)' | grep Standard # Standard deviation : 0.0155192 stats 'ngc3198_2008n_sincsinc_trim_submask3_20kms(~15,~15,840.0:1040.0)' | grep Standard # Standard deviation : 0.0135034 --------------------------------------------------------------------------------------- MAKE MOMENT MAPS mkdir Moment_sincsinc cp ngc3198_2008n_sincsinc_trim_submask3.sdf Moment_sincsinc/ngc3198_2008n_sincsinc_trim_submask3_rembsl.sdf cp ngc3198_2008n_trim_basemask1.sdf Moment_sincsinc/ngc3198_2008n_sincsinc_trim_submask3_msk.sdf cp ~/DATA/REDUCE_SCIENCE_GRADIENT Moment_sincsinc/ cd Moment_sincsinc *Edit RSG Script: emacs REDUCE_SCIENCE_GRADIENT & set object = "ngc3198_2008n_sincsinc_trim_submask3" set vrange = "low=480 high=840" # Quick jumps to repeat parts of the reduction #goto STEP1a #goto STEP1b #goto STEP2 goto STEP3 foreach rmslev (2.0 2.5 3.0 3.5 4.0 4.5 5.0) **Remove trim function: # Multiply baselined cube with clump mask ${kap}/mult ${sdf} ${sdf}_findclump${rmslab} ${sdf}_clmsk${rmslab} # Trim the noisy edges # ${kap}/ndfcopy temp'(-26:26,-36:40,)' ${sdf}_clmsk${rmslab} \rm -f ${sdf}_mask_temp1.sdf ${sdf}_mask_temp2.sdf >& /dev/null ./REDUCE_SCIENCE_GRADIENT #RMS MASKED SMOOTHED CUBE(-2:8,-11:17,-852:850): 0.00822879500607012 Rotate intvel maps: rotate ngc3198_2008n_sincsinc_trim_submask3_intvel20 ngc3198_2008n_sincsinc_trim_submask3_intvel20_rot angle=! rotate ngc3198_2008n_sincsinc_trim_submask3_intvel25 ngc3198_2008n_sincsinc_trim_submask3_intvel25_rot angle=! rotate ngc3198_2008n_sincsinc_trim_submask3_intvel30 ngc3198_2008n_sincsinc_trim_submask3_intvel30_rot angle=! rotate ngc3198_2008n_sincsinc_trim_submask3_intvel35 ngc3198_2008n_sincsinc_trim_submask3_intvel35_rot angle=! rotate ngc3198_2008n_sincsinc_trim_submask3_intvel40 ngc3198_2008n_sincsinc_trim_submask3_intvel40_rot angle=! rotate ngc3198_2008n_sincsinc_trim_submask3_intvel45 ngc3198_2008n_sincsinc_trim_submask3_intvel45_rot angle=! rotate ngc3198_2008n_sincsinc_trim_submask3_intvel50 ngc3198_2008n_sincsinc_trim_submask3_intvel50_rot angle=! (contoured totint maps and velocity maps then made using gaia) --------------------------------------------------------------------------------------------- MAKE CUBES DISPLAYING S/N RATIO: collapse in=ngc3198_2008n_sincsinc_trim_submask3.sdf low=270 high=480 axis=vopt out=sigma_low estimator=sigma reset collapse in=ngc3198_2008n_sincsinc_trim_submask3.sdf low=840 high=1040 axis=vopt out=sigma_high estimator=sigma reset block ngc3198_2008n_sincsinc_trim_submask3_20kms.sdf ngc3198_2008n_sincsinc_trim_submask3_20kms_bin5.sdf axes="[2,3]" box="[1,5]" estimator=mean block ngc3198_2008n_sincsinc_trim_submask3_20kms.sdf ngc3198_2008n_sincsinc_trim_submask3_20kms_bin9.sdf axes="[2,3]" box="[1,9]" estimator=mean manic in=sigma_low out=sigma_low_3d axes="[1,2,0]" lbound=-19 ubound=40 reset wcscopy ndf=sigma_low_3d like=ngc3198_2008n_sincsinc_trim_submask3_20kms_bin5.sdf confirm=false reset cdiv sigma_low_3d 6.85565 sigma_low_20kms div ngc3198_2008n_sincsinc_trim_submask3_20kms_bin5.sdf sigma_low_20kms.sdf ngc3198_2008n_sincsinc_trim_submask3_20kms_bin5_SN.sdf div ngc3198_2008n_sincsinc_trim_submask3_20kms_bin9.sdf sigma_low_20kms.sdf ngc3198_2008n_sincsinc_trim_submask3_20kms_bin9_SN.sdf ---------------------------------------------------------------------------------------------- MAKE TOTINT20 NOISE and SN MAPS: mkdir noise_maps cd noise_maps cp ../sigma_high.sdf . cp ../sigma_low.sdf . cp ../Moment_sincsinc/ngc3198_2008n_sincsinc_trim_submask3_rembsl_findclump20.sdf . cp ../Moment_sincsinc/ngc3198_2008n_sincsinc_trim_submask3_totint20.sdf . add sigma_high sigma_low out=sigma_sum cdiv sigma_sum.sdf 2 sigma_avg collapse in= ngc3198_2008n_sincsinc_trim_submask3_rembsl_findclump20.sdf out=test0 axis=vopt estimator=sum wlim=0.0 (must make mask to multiply by collapsed image so values are real, not integer, otherwise sqrt function below will round to nearest integer) cmult sigma_avg 0 test1 cadd test1 1.0 test2 mult test2 test0 test3 (now have same collapsed image with real, not integer values) maths exp="sqrt(IA)" out=test4 IA=test3 cmult test4 0.42 test5 mult sigma_avg test5 ngc3198_2008n_totint20_noise div ngc3198_2008n_sincsinc_trim_submask3_totint20.sdf ngc3198_2008n_totint20_noise ngc3198_2008n_totint20_SN (contoured SN maps then made using gaia) ----------------------------------------------------------------------------------------------- TRIM OUT NOISE FROM TOTINT20 MAP: *Try trimming regions with SN < 2.0: thresh ngc3198_2008n_totint20_SN.sdf thrlo=2 newlo=bad thrhi=500 newhi=bad out=ngc3198_2008n_totint20_SN_thresh2 #There were 431 elements changed in the DATA array below the threshold. #There were 0 elements changed in the DATA array above the threshold. cmult ngc3198_2008n_totint20_SN_thresh2.sdf 0 ngc3198_2008n_zeroall_2 cadd ngc3198_2008n_zeroall_2 1 ngc3198_2008n_msk_2 mult ngc3198_2008n_msk_2.sdf ngc3198_2008n_sincsinc_trim_submask3_totint20.sdf ngc3198_2008n_totint20_noisetrim2 *Try trimming regions with SN < 1.5: thresh ngc3198_2008n_totint20_SN.sdf thrlo=1.5 newlo=bad thrhi=500 newhi=bad out=ngc3198_2008n_totint20_SN_thresh15 #There were 370 elements changed in the DATA array below the threshold. #There were 0 elements changed in the DATA array above the threshold. cmult ngc3198_2008n_totint20_SN_thresh15.sdf 0 ngc3198_2008n_zeroall_15 cadd ngc3198_2008n_zeroall_15.sdf 1 ngc3198_2008n_msk_15 mult ngc3198_2008n_msk_15.sdf ngc3198_2008n_sincsinc_trim_submask3_totint20.sdf ngc3198_2008n_totint20_noisetrim_15 ---------------------------------------------------------------------------------------------- MAKE TOTINT30 NOISE and SN MAPS: cp ../Moment_sincsinc/ngc3198_2008n_sincsinc_trim_submask3_rembsl_findclump30.sdf . cp ../Moment_sincsinc/ngc3198_2008n_sincsinc_trim_submask3_totint30.sdf . collapse in= ngc3198_2008n_sincsinc_trim_submask3_rembsl_findclump30.sdf out=test6 axis=vopt estimator=sum wlim=0.0 mult test2 test6 test7 maths exp="sqrt(IA)" out=test8 IA=test7 cmult test8 0.42 test9 mult sigma_avg test9 ngc3198_2008n_totint30_noise div ngc3198_2008n_sincsinc_trim_submask3_totint30.sdf ngc3198_2008n_totint30_noise ngc3198_2008n_totint30_SN (contoured SN maps then made using gaia)