runvelfitss07 is a NEMO frontend to velfitss07, and fits a non-parametric non-axisymmetric flow to a planar disk galaxy velocity field. Radial (m=0) flow can be fitted by setting flags(4) , whereas bisymmetric (m=2) or lobsided (m=1) flow with flags(5) (see flags= below). The code is usually not well converging when fitting both. For more details and a full explanation of the parameters, see Spekkens and Sellwood (2007ApJ...664..204S) and the Users Guide on the VELFIT website.
runvelfitss07 examp.vels examp1.m examp1.p vsys=500 pa=80 inc=60 radii=2.5:50:2.5of which you can find the results in Figure 7. and the next 3 examples from the Users Guide are run as follows:
runvelfitss07 examp.vels examp2.m examp2.p vsys=500 pa=80 inc=60 radii=2.5:50:2.5 flags=1,1,1,1,0,0 rmsism=2.0 runvelfitss07 examp.vels examp3.m examp3.p vsys=500 pa=80 inc=60 radii=2.5:50:2.5 flags=1,1,1,1,0,1 rmsism=2.0 runvelfitss07 examp.vels examp4.m examp4.p vsys=500 pa=80 inc=60 radii=2.5:50:2.5 flags=1,1,1,0,1,1 rmsism=2.0which you can find the output .p parameter files in Figure 9, Figure 11 and Figure 3 resp.
1 1.29 2 4.86 3 - 4 -
2007ApJ...664..204S (Spekkens & Sellwood, 2007) WEBSITE: http://www.physics.rutgers.edu/~spekkens/velfit/ WEBSITE: http://www.physics.rutgers.edu/~spekkens/diskfit/
NEMO/usr/spekkens/velfit - original SS07 code/examples (see their WEBSITE) NEMOBIN/velfitss07 - SS07 executable (see also exe=) NEMO/src/image/rotcur - NEMO frontend
#! /bin/csh -f # # create an example input file for velfitss07 from a MIRIAD dataset # 1) select the inner 40x40 arcsec # 2) select points where flux mom0 map > 1.0 set region="arcsec,box(-20,-20,20,20)" set map=mom0a set vel=mom1a set sig=mom2a imtab in=$map region="$region" > tmp.map1.tab imtab in=$vel region="$region" > tmp.vel1.tab tail -n +3 tmp.map1.tab | grep -v Found > tmp.map2.tab tail -n +3 tmp.vel1.tab | grep -v Found > tmp.vel2.tab # here’s a simple example of constant errors and only use mom0 and mom1 echo "# input file for velfitss07: $0" echo "# X Y VEL EVEL" echo "# asec asec km/s km/s" echo "#----------------------------------------" paste tmp.map2.tab tmp.vel2.tab | awk ’{ if ($3 > 1.0) printf("%10.2f%10.2f%10.2f%10.2f0,$1,$2,$6,5.0)}’
27-oct-2008 V1.0 : frontend for velfitss07 V1.0 (oct 2008) PJT ??-???-???? V1.1 : should be compatible with their V1.1 (aug 2009) PJT ??-???-???? V2.0 : still not implemented PJT 14-sep-2012 noting this is now deprecated in favor of diskfit PJT 7-aug-2014 V1_0, V1_1, V2_0 now in a local CVS with NUMREC PJT