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.
It should be noted that velfit has been superseded by diskfit.
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