/Virgo/FeldPNe/PNe8alpha/subs.sm
OR
/Virgo/FeldPNe/vat/PNe8alpha/subs.sm

subs
    #subsampling PNe catalog
    data 8.PNe.3mm.fits.alpha2523.RADEC.dat
    read { RA3 1 DEC3 2 }
   
    data 8.PNe.4mm.fits.alpha2523.RADEC.dat
    read { RA4 1 DEC4 2 }
   
    data 8.PNe.7mm.fits.alpha2523.RADEC.dat
    read { RA7 1 DEC7 2 }
   
    data 8.PNe.Submm.fits.alpha2523.RADEC.dat
    read { RASub 1 DECSub 2 }
   
    #set limits of subsample
    define num3 75
    define num4 11
    define num7 141
    define numSub 34
   
    define w 0
    set dimen(PNeRA) = $($num3+$num4+$num7+$numSub)
    set dimen(PNeDEC) = dimen(PNeRA)
   
    #do the subsampling
    #for 3
    set select3 = random($(dimen(RA3)))
    sort {select3 RA3 DEC3}
    do i = 0, $($num3-1) {
        set PNeRA[$w] = $(RA3[$i])
        set PNeDEC[$w] = $(DEC3[$i])
        define w $($w+1)
    }
    #for 4
    set select4 = random($(dimen(RA4)))
    sort {select4 RA4 DEC4}
    do i = 0, $($num4-1) {
        set PNeRA[$w] = $(RA4[$i])
        set PNeDEC[$w] = $(DEC4[$i])
        define w $($w+1)
    }
    #for 7
    set select7 = random($(dimen(RA7)))
    sort {select7 RA7 DEC7}
    do i = 0, $($num7-1) {
        set PNeRA[$w] = $(RA7[$i])
        set PNeDEC[$w] = $(DEC7[$i])
        define w $($w+1)
    }
    #for Sub
    set selectSub = random($(dimen(RASub)))
    sort {selectSub RASub DECSub}
    do i = 0, $($numSub-1) {
        set PNeRA[$w] = $(RASub[$i])
        set PNeDEC[$w] = $(DECSub[$i])
        define w $($w+1)
    }
    do i = 0, $(dimen(PNeRA)-1) {
        write PNeFAKE.dat $(PNeRA[$i]) $(PNeDEC[$i])
    }