0925 Arrive. Going to work on debugging that glish script first... 1040 After reading up on coordsys/image tools/modules/etc, there's a more obvious and simpler way to do what I was trying to do. Instead of tweaking with coordinate systems, it seems easiest to simply re-assign (by .setcoordsys()) the old coord sys onto the sub-cube on the velocity cut. Duh. 1050 subv.g --> subv() seems to work. documented in the code. also making subv2() (in subv.g), to do velocity cuts across the whole cube, rather than just one chunk at a time (to save time loading images) 1200 found a bug in subv2.g (and subv.g, too)... i only calculated change in velocity, and always scaled from 0 plane (-250 km/s). this has been fixed in both now. 1220 latest bug: when doing a range from 50-100, for example, the velocity it thinks you start from is still -250... gotta fix the coordinate reassigment. 1255 SEEMS TO WORK! Hooray. there is chance for slight errors, in that I'm resetting the reference value/pixel on the velocities, without actually reading the velocities off the cube, but instead by assigning whole numbers (like 250) to them. eh? also updated subv.g with the fix 1305 finally i can get the 4 sub-cubes I started trying for yesterday: g_v_-250_50.fits g_v_51_124.fits g_v_125_200.fits g_v_200_250.fits AND gsmall_v_-250_50.fits gsmall_v_51_124.fits gsmall_v_125_200.fits gsmall_v_200_250.fits those names are pretty explanatory 1325 reformatted some logs from earlier in the week. going to find a way to do the moment maps... perhaps by adding to subv2.g, or maybe making a subv3.g ? first, going to take care of a bill 1335 back - paid for deposit ($135), may rent ($104.40), & june rent($270), $509.40, check 0098 from First Citizens acct. now going to use the GUI to make moment maps... 1350 Yes, so this works, but making it scripted/automated would be vastly superior! time for subvmom.g! 1440 testing subvmom.g... errors of about 0.001 between doing it by hand (errors from velocity assignments mentioned above) 1505 generating moment maps for selected range of g and gsmall... [integrated values (0)] format is: subcube name(spatially)_ v (for velocities)_ lower cutoff velocity (km/s)_ upper cutoff velocity (km/s)_ moment map computation_ number of moment(0=integration,-1=mean) .fits 1550 g_v_120_150_mom_0.fits g_v_125_200_mom_0.fits g_v_200_250_mom_0.fits g_v_-250_-50_mom_0.fits g_v_-250_50_mom_0.fits g_v_-25_25_mom_0.fits g_v_25_50_mom_0.fits g_v_50_125_mom_0.fits g_v_-50_-25_mom_0.fits g_v_75_125_mom_0.fits ^^these were just picked out to highlight interesting clouds/features going to work on some other sub-cubes(spatially) now... 1620 ran it through for sub-cube h, results are in /h/ folder. interesting results for certain velocity ranges... some striping/gridding patterns... in -50,-25 km/s range, so I'm going to go through that part again, by hand, to check it. --nothing obvious when scanning through those ranges in the big cube --manually, going to pull velocity planes 243-274 using these commands: - im := image('h.img') - im2 := im.getchunk([1,1,243],[x,x,274]) - im3 := imagefromarray(pixels=im2) - im3mom := im3.moments(moments=0) - ok := im3mom.tofits('h_v_-50_-25_mom_0_doublecheck.fits') can look at this image ^^^ to see the grid pattern, still. **interesting side note:** in that doublecheck image, the overall max and min are a bit different than in the original (script-made) one max min full filename dblchk: 1.18184 0.0733172 h_v_-50_-25_mom_0_doublecheck.fits script: 4.13246 0.216722 h_v_-050_-025_mom_0.fits 1700 TO DO: ask about this next week. either the grid pattern is left over from part of the reduction, or I've bungled something. TO DO: also ask about the difference in numbers, and if that's ok 1715 Toney Minter came down and cleared up some of our questions about Hannah house stuff: -leave trash out on M/Th mornings, by driveway -house doesn't have an address, have mail sent to NRAO: PO Box 2, Green Bank, WV 4944944949494 -he doesn't know who the other 3 people are in the house, but thinks Danielle is coming Sunday, and Katie on Monday... -mouse in kitchen - use trap, or buy them and save receipts for reimbursement -fire outside fine whenever, no burning house down -well is on electric. no power, no water 1527 leaving w/ Toney and his wife for Pizza and Dairy King... wahoo! 1945 back from pizza/ice cream, hooray! They treated us out to a fun meal/dessert and conversation. She's sort of a surrogate mother with lots of "inside info" and stories to tell of previous summers/etc. she also mentioned Sunday 5pm, 2 houses down from theirs in the rabbit patch, "Larry" is having a party since he got engaged. going to head back now, taking book with me for the weekend to read.