People often make the mistake of putting " in place of `` in LaTeX documents. To repair this, the only easy solution is something like:
s/\(\s\)"/\1``/g
in VIM or sed.
I've seen alternate solutions posted but they're all too complicated. One recommendation was to switch to XeTeX, which sounds like overkill, and others all stated that a complicated perl script is required. The latter is technically true, but why isn't such a thing readily available?
Any other ideas?
Saturday, December 27, 2008
Sunday, December 21, 2008
Thursday, December 11, 2008
skill set
is it a marketable skill that I can bring down almost any computer system I'm given access too entirely unintentionally?
Friday, December 5, 2008
post rate
my post rate is too low I need to up it so that my front page doesn't show the same posts any more
Wednesday, December 3, 2008
scary
I think this is my worst nightmare:
http://antwrp.gsfc.nasa.gov/apod/ap081203.html
LA. Not the smiley.
http://antwrp.gsfc.nasa.gov/apod/ap081203.html
LA. Not the smiley.
Monday, November 24, 2008
.inputrc
New unix stuff!
This is, again, stuff I'd searched for forever but never found the right collection of terms. This is what bash readline is REALLY all about.
if you make an inputrc file, you can get awesome readline capabilities (similar to ipython). This post at lifehacker showed me how to do it. This site has a more complete description of how to use readline/inputrc.
Even better: VI command line editing .
There's also something called magic space but I don't know what purpose it serves all the time... it will replace 'magic' characters like !* and !& and !! with whatever command they represent if you hit space immediately after typing them (explained better here).
Wow. Awesome.
This is, again, stuff I'd searched for forever but never found the right collection of terms. This is what bash readline is REALLY all about.
if you make an inputrc file, you can get awesome readline capabilities (similar to ipython). This post at lifehacker showed me how to do it. This site has a more complete description of how to use readline/inputrc.
Even better: VI command line editing .
There's also something called magic space but I don't know what purpose it serves all the time... it will replace 'magic' characters like !* and !& and !! with whatever command they represent if you hit space immediately after typing them (explained better here).
Wow. Awesome.
Thursday, November 20, 2008
MOVIES!
Making movies is surprisingly difficult. No matter what language you use, apparently true movie files can only be made by stringing together images, i.e. there is no native movie-producing feature. Gnuplot can do some simple animations but to do anything sophisticated you need to start delving into variables, and for that I switch to a real language.
So, I returned to python. As usual, it took no more than a few hours of coding and learning to come up with something. But it bothers me that it took that long: I still think python is most deficient in its failure to create a default column-text reader like 'readcol' in IDL. I can't complain that much, though: I wrote my own in about 5 minutes.
Anyway, the key is to use the
I use
ImageMagick's convert can be used to stitch any kind of image into a movie given that you've installed an mpeg2 encoder (fink gave me
I had to use a smaller image size because a series of 1000x12kb files somehow chomped ~6-8 GB of RAM and swap space.
So, I returned to python. As usual, it took no more than a few hours of coding and learning to come up with something. But it bothers me that it took that long: I still think python is most deficient in its failure to create a default column-text reader like 'readcol' in IDL. I can't complain that much, though: I wrote my own in about 5 minutes.
Anyway, the key is to use the
.set_xdata and .set_ydata functions of a plot to update a canvas. I still don't have nearly as high a plotting speed as I'd like, but it works alright if I don't display to screen. Probably a different backend would be more effective but I don't like to mess with backends.I use
savefig(filename,dpi=50) to reduce the image quality so that it's easier for the animator to handle.ImageMagick's convert can be used to stitch any kind of image into a movie given that you've installed an mpeg2 encoder (fink gave me
mpeg2vidcodec ). The command is very simple:
convert -size 300x300 *.png movie.mpg
I had to use a smaller image size because a series of 1000x12kb files somehow chomped ~6-8 GB of RAM and swap space.
Monday, November 17, 2008
Screen cont'd
Guide to screen
Particularly useful features:
Scrollback:
Switch to a numbered screen (doesn't work for me so far):
Particularly useful features:
Scrollback:
ctrl-a [
[scroll keys]
Switch to a numbered screen (doesn't work for me so far):
ctrl-a [number]
Sunday, November 16, 2008
Ultimate today
We went 2-1 in the tournament, losing the final in a windy game on universe point, 15-14. We won the first game 13-5 or 13-6, and the second by a similarly wide margin. I think I scored 15-20 points on the day and threw an additional few. 2-3 drops and 3-4 throwaways to balance that out. Pretty good on the whole for what was probably my last ultimate of the season.
Of course, that made me very exhausted for the review session this evening. Oh well, I think it worked out. I need to remember to come up with (and post!) a clear explanation of why you can use that silly little triangle to come up with the translation between angle on the sky and image size.
Of course, that made me very exhausted for the review session this evening. Oh well, I think it worked out. I need to remember to come up with (and post!) a clear explanation of why you can use that silly little triangle to come up with the translation between angle on the sky and image size.
Thursday, November 13, 2008
Eta Carinae
Post #4 and I'm diverging from the purported theme. But it's Eta Carinae, it must be discussed!
Eta Carinae is the most amazing object in the sky, in my opinion. It has been the most active nearby (around 2 kiloparsecs, or 6000 light-years, distant - close only relative to, say, the Galactic Center or other galaxies) object in the past few centuries. In 1840 it brightened by a few magnitudes, and it has been discovered since then that this event was a huge explosion.
But it wasn't a 'normal' explosion. Most exploding stars blow up and die in supernovae. Those are awesome events that sometimes completely annihilate a star and sometimes leave a black hole or neutron star behind.
This explosion was nearly as powerful as a supernova, but the star is still there. And it turns out that there are TWO stars there. There are a number of different theories going around, but this system is incredibly difficult to observe because there's just so much stuff there.
The prevailing theory is that there are two stars, one around 80 solar masses, the other around 20. Either one of these on its own would be pretty big, hot, and impressive. Two together is something pretty impressive - but actually not that uncommon.
Other pairs of hot, massive stars are observed frequently. One example is WR 40, in which two 80 solar mass stars are orbiting one another. The curious thing is that most of these massive stars have very circular orbits, but in Eta Carinae the small star has a very eccentric orbit around the large star.
One explanation for this - which is not widely accepted yet because it's very hard to confirm, but it's plausible - is that Eta Carinae used to be a nearly equal-mass binary system like WR40, but one of the stars blew off its outer envelope (the outer part of the star) in a huge explosion, and when it lost this mass its orbit changed. One bit of evidence in favor of this hypothesis is that the small star has a very fast wind - 3000 km/s, which is faster than any normal star's wind. Wolf Rayet stars have winds that fast, though, and they are stars with blown-off envelopes. It's going to take a lot of skilled observations to confirm this, though.
Another neat thing is that the smaller star disappears on its closest approach to the large star. Why? The large star is losing mass at a very fast rate - about a Jupiter mass every year! That's more than a million times as fast as the sun is losing its mass in the Solar Wind. All of this mass loss has an interesting result - you can't see the star, you only see the outside of its wind. But the small star passes inside the 'sphere' that is producing the wind, and so it's as though the small star was going inside the large one! You don't see it at all during closest approach.
I get excited about this star... I'd like to add a lot more on it, and so you might see Luminous Blue Variables as a recurring theme on this blog. I hope to post a lot of pictures and add some links to references too.
Eta Carinae is the most amazing object in the sky, in my opinion. It has been the most active nearby (around 2 kiloparsecs, or 6000 light-years, distant - close only relative to, say, the Galactic Center or other galaxies) object in the past few centuries. In 1840 it brightened by a few magnitudes, and it has been discovered since then that this event was a huge explosion.
But it wasn't a 'normal' explosion. Most exploding stars blow up and die in supernovae. Those are awesome events that sometimes completely annihilate a star and sometimes leave a black hole or neutron star behind.
This explosion was nearly as powerful as a supernova, but the star is still there. And it turns out that there are TWO stars there. There are a number of different theories going around, but this system is incredibly difficult to observe because there's just so much stuff there.
The prevailing theory is that there are two stars, one around 80 solar masses, the other around 20. Either one of these on its own would be pretty big, hot, and impressive. Two together is something pretty impressive - but actually not that uncommon.
Other pairs of hot, massive stars are observed frequently. One example is WR 40, in which two 80 solar mass stars are orbiting one another. The curious thing is that most of these massive stars have very circular orbits, but in Eta Carinae the small star has a very eccentric orbit around the large star.
One explanation for this - which is not widely accepted yet because it's very hard to confirm, but it's plausible - is that Eta Carinae used to be a nearly equal-mass binary system like WR40, but one of the stars blew off its outer envelope (the outer part of the star) in a huge explosion, and when it lost this mass its orbit changed. One bit of evidence in favor of this hypothesis is that the small star has a very fast wind - 3000 km/s, which is faster than any normal star's wind. Wolf Rayet stars have winds that fast, though, and they are stars with blown-off envelopes. It's going to take a lot of skilled observations to confirm this, though.
Another neat thing is that the smaller star disappears on its closest approach to the large star. Why? The large star is losing mass at a very fast rate - about a Jupiter mass every year! That's more than a million times as fast as the sun is losing its mass in the Solar Wind. All of this mass loss has an interesting result - you can't see the star, you only see the outside of its wind. But the small star passes inside the 'sphere' that is producing the wind, and so it's as though the small star was going inside the large one! You don't see it at all during closest approach.
I get excited about this star... I'd like to add a lot more on it, and so you might see Luminous Blue Variables as a recurring theme on this blog. I hope to post a lot of pictures and add some links to references too.
BASH discoveries, readline questions
1.
Maybe my hostname completion worked and then stopped working because the bash option hostcomplete was not set. Duh! Why? I don't know. Anyway,
2. it's really hard to search for readline stuff on google. Can anyone explain to me how BASH readline works? I would REALLY like to make bash readline work like ipython, in which you can start typing a command and hit the 'up' key to search through the history for anything beginning with the stuff you've typed up to that point. But I can't even find documentation for the ipython readline! Any hints, anyone?
3. my desktop at work blocks ssh connections. I can ssh into some computers and then into it, but not directly into it.
...as usual, I made a list where not-a-list would have sufficed, and I had to add the last thing because a 2-item list is dumb.
shopt -p Maybe my hostname completion worked and then stopped working because the bash option hostcomplete was not set. Duh! Why? I don't know. Anyway,
shopt -s hostcomplete solves the problem. nocaseglob is also pretty cool (case insensitive tab completion)2. it's really hard to search for readline stuff on google. Can anyone explain to me how BASH readline works? I would REALLY like to make bash readline work like ipython, in which you can start typing a command and hit the 'up' key to search through the history for anything beginning with the stuff you've typed up to that point. But I can't even find documentation for the ipython readline! Any hints, anyone?
3. my desktop at work blocks ssh connections. I can ssh into some computers and then into it, but not directly into it.
...as usual, I made a list where not-a-list would have sufficed, and I had to add the last thing because a 2-item list is dumb.
Wednesday, November 12, 2008
Why do astronomers have such a strong presence on the web?
I'm not making this a complete post, just a few examples of blogs and websites I'm aware of. But we do have a strong presence on the web - astronomers have an unusually high google ranking etc. Is it just because 'we' were here first ('we' excludes me, I'm just jumping on the bandwagon and getting a free ride)?
Examples:
Pamela Gay
Dr Lisa
Science writers who write on astronomy:
Dave Mosher
Examples:
Pamela Gay
Dr Lisa
Science writers who write on astronomy:
Dave Mosher
Sunday, November 9, 2008
Galactic Center Photo Contest
The NRAO Photo Contest
The Galactic Center and the surrounding Central Molecular Zone comprise the most active star formation region in the Milky Way. This 2 x 1 degree field was imaged at 20 cm (purple) with the NRAO Very Large Array, tracing H II regions that are illuminated by hot, massive stars, supernova remnants, and synchrotron emission. Emission at 1.1 mm (orange) was observed with the Caltech Submillimeter Observatory and highlights cold (20-30 K) dust associated with molecular gas. Some of this material will form stars within in the next few million years; the remainder will be blown away. The diffuse cyan and colored star images are from the Spitzer Space Observatory's Infrared Array Camera. The cyan is primarily emission from stars, the point sources, and from polycyclic aromatic hydrocarbons (PAHs), the diffuse component.
My view of the Galactic Center is the winner of this year's NRAO photo contest. It displays a multiwavelength view of the galactic center. The caption at the NRAO photo contest (reproduced above) describes it well, but I will go into further depth here.
The 'arches' are the large purple filamentary structures seen near the brightest point in the map, Sgr A*, the center of our galaxy. They are thought to be large scale magnetic fields possibly generated by winds from hot stars. Farhad Yusef-Zadeh, my collaborator on the photo, leads the study of radio emission from the Galactic Center.
The orange color is 1.1mm emission, which means that we're seeing something completely different. This light actually comes from dust. We usually see dust in absorption because in the optical, it blocks light. There is so much dust between us and the Galactic Center that we can't see it at all at optical wavelengths.
Sagittarius B2 is the bright blob just to the left of center. It is the most massive collection of gas and dust in our galaxy. It is full of different types of molecules including complex organic molecules. A new generation of massive stars is thought to be forming there.
Spitzer's 8 micron view is presented in cyan. Mostly this comes from PAH, Polycyclic Aromatic Hydrocarbon, emission. PAHs are big molecules with lots of Carbon and Hydrogen chained together. Stars are visible at these wavelengths because the dust isn't as effective at blocking out infrared light - but look closely at the orange regions on the right side of the image. You might notice that not only is there dust there, but there's also less PAH emission visible. This is because on the right side of the Galactic Center, the spiral arm is behind the dust, and on the left side it is in front. When dust blocks out infrared light, the object we see is called an 'infrared dark cloud'. These clouds are usually places where star formation has not yet started, but will soon.
This image was created using GIMP, the Gnu Image Manipulation Program, using a variety of layers, opacities, etc.
Wednesday, November 5, 2008
Tuesday, October 28, 2008
Python magic / advanced numpy indexing
Yeah, indexing python arrays should really be easy.
Stefan van der Walt's page
Stefan van der Walt's page
In [85]: bi = (f.bolo_indices[np.newaxis,:] + zeros([7751,1])).astype('int')
In [86]: whc = (whscan[:,np.newaxis] + zeros([1,107])).astype('int')
In [87]: array2d[whc,bi] = temp2d
Python: one-line arrays
Ahhh, refreshing:
Not like IDL, which takes at least 4 lines b/c of the variable declaration. There's probably a better way to do that too.
whscan = asarray([arange(scanlen)+i for i in scans_info[:,0]]).ravel()
Not like IDL, which takes at least 4 lines b/c of the variable declaration. There's probably a better way to do that too.
Port vs Fink
Port failed to install:
I don't want to deal with it so I'm just using Fink.
adam-macbook ~$ sudo port install py25-matplotlib
---> Building py25-hashlib with target build
---> Staging py25-hashlib into destroot
---> Installing py25-hashlib 2.5.2_0
---> Activating py25-hashlib 2.5.2_0
---> Cleaning py25-hashlib
---> Fetching py25-zlib
---> Verifying checksum(s) for py25-zlib
---> Extracting py25-zlib
---> Configuring py25-zlib
---> Building py25-zlib with target build
---> Staging py25-zlib into destroot
---> Installing py25-zlib 2.5.2_0
---> Activating py25-zlib 2.5.2_0
---> Cleaning py25-zlib
---> Fetching py25-setuptools
---> Attempting to fetch setuptools-0.6c8.tar.gz from http://cheeseshop.python.org/packages/source/s/setuptools/
---> Verifying checksum(s) for py25-setuptools
---> Extracting py25-setuptools
---> Applying patches to py25-setuptools
---> Configuring py25-setuptools
---> Building py25-setuptools with target build
---> Staging py25-setuptools into destroot
---> Installing py25-setuptools 0.6c8_2
---> Activating py25-setuptools 0.6c8_2
---> Cleaning py25-setuptools
---> Fetching py25-tz
---> Attempting to fetch pytz-2007k.tar.bz2 from http://cheeseshop.python.org/packages/source/p/pytz/
---> Verifying checksum(s) for py25-tz
---> Extracting py25-tz
---> Configuring py25-tz
---> Building py25-tz with target build
---> Staging py25-tz into destroot
---> Installing py25-tz 2007k_0
---> Activating py25-tz 2007k_0
---> Cleaning py25-tz
---> Fetching py25-dateutil
---> Attempting to fetch python-dateutil-1.4.1.tar.gz from http://labix.org/download/python-dateutil
---> Verifying checksum(s) for py25-dateutil
---> Extracting py25-dateutil
---> Configuring py25-dateutil
---> Building py25-dateutil with target build
---> Staging py25-dateutil into destroot
---> Installing py25-dateutil 1.4.1_0
---> Activating py25-dateutil 1.4.1_0
---> Cleaning py25-dateutil
---> Fetching fftw-3
---> Attempting to fetch fftw-3.1.3.tar.gz from http://www.fftw.org/
---> Verifying checksum(s) for fftw-3
---> Extracting fftw-3
---> Configuring fftw-3
---> Building fftw-3 with target all
---> Staging fftw-3 into destroot
---> Installing fftw-3 3.1.3_0
---> Activating fftw-3 3.1.3_0
---> Cleaning fftw-3
---> Fetching py25-nose
---> Attempting to fetch nose-0.10.4.tar.gz from http://somethingaboutorange.com/mrl/projects/nose/
---> Verifying checksum(s) for py25-nose
---> Extracting py25-nose
---> Configuring py25-nose
---> Building py25-nose with target build
---> Staging py25-nose into destroot
---> Installing py25-nose 0.10.4_0
---> Activating py25-nose 0.10.4_0
---> Cleaning py25-nose
---> Fetching py25-numpy
---> Attempting to fetch numpy-1.2.0.tar.gz from http://downloads.sourceforge.net/numpy
---> Verifying checksum(s) for py25-numpy
---> Extracting py25-numpy
---> Applying patches to py25-numpy
---> Configuring py25-numpy
---> Building py25-numpy with target build
Error: Target org.macports.build returned: shell command " cd "/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_python_py25-numpy/work/numpy-1.2.0" && /opt/local/bin/python2.5 setup.py build " returned error 1
Command output: setup_package()
File "setup.py", line 89, in setup_package
configuration=configuration )
File "/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_python_py25-numpy/work/numpy-1.2.0/numpy/distutils/core.py", line 184, in setup
return old_setup(**new_attr)
File "/opt/local/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/distutils/core.py", line 151, in setup
dist.run_commands()
File "/opt/local/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/distutils/dist.py", line 974, in run_commands
self.run_command(cmd)
File "/opt/local/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/distutils/dist.py", line 994, in run_command
cmd_obj.run()
File "/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_python_py25-numpy/work/numpy-1.2.0/numpy/distutils/command/build.py", line 37, in run
old_build.run(self)
File "/opt/local/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/distutils/command/build.py", line 112, in run
self.run_command(cmd_name)
File "/opt/local/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/distutils/cmd.py", line 333, in run_command
self.distribution.run_command(command)
File "/opt/local/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/distutils/dist.py", line 994, in run_command
cmd_obj.run()
File "/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_python_py25-numpy/work/numpy-1.2.0/numpy/distutils/command/build_src.py", line 130, in run
self.build_sources()
File "/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_python_py25-numpy/work/numpy-1.2.0/numpy/distutils/command/build_src.py", line 147, in build_sources
self.build_extension_sources(ext)
File "/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_python_py25-numpy/work/numpy-1.2.0/numpy/distutils/command/build_src.py", line 250, in build_extension_sources
sources = self.generate_sources(sources, ext)
File "/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_python_py25-numpy/work/numpy-1.2.0/numpy/distutils/command/build_src.py", line 307, in generate_sources
source = func(extension, build_dir)
File "numpy/core/setup.py", line 87, in generate_config_h
raise SystemError,"Failed to test configuration. "\
SystemError: Failed to test configuration. See previous error messages for more information.
Error: The following dependencies failed to build: py25-numpy py25-tkinter tcl tk
Error: Status 1 encountered during processing.
I don't want to deal with it so I'm just using Fink.
Sunday, October 26, 2008
Reinstalling a system
Python died tonight when I foolishly tried to update numpy or matplotlib. One was not compatible with the other in upgraded form, and the update for matplotlib wouldn't install because of GTK issues that were totally opaque.
So, I reinstalled EVERYTHING - fink AND macports - from scratch. If I don't get things running by morning I'm going to have to do a restore, which is ugly as death.
Things I need to do post-install:
test ipython
test matplotlib
test numpy
So, I reinstalled EVERYTHING - fink AND macports - from scratch. If I don't get things running by morning I'm going to have to do a restore, which is ugly as death.
Things I need to do post-install:
sudo port install python_select
sudo python_select python25
test ipython
test matplotlib
test numpy
Python: setting matplotlib defaults
Setting matplotlib defaults is a lot more difficult than it should be.
matplotlib.defaultParams['image.origin']='lower'
matplotlib.defaultParams['image.interpolation']='nearest'
Monday, October 20, 2008
Observing 10/20
I don't have a better place to post this one, so here it is:
http://sites.google.com/site/iras05358/figure-discussion/figures/iras05358_triple_17.29.png?attredirects=0
My automated fitter (Gaussfitting Cube Collapser) has come a long way. I now adaptively choose to fit 1, 2, or 3 Gaussian components to output to a data cube.
The purpose of that code is primarily to find a two-dimensional way to display information about the 3D structure, specifically about the presence/absence of outflows. Outflows will inevitably be confused with multiple velocity components, but they are also likely to be convolved with them.
http://sites.google.com/site/iras05358/figure-discussion/figures/iras05358_triple_17.29.png?attredirects=0
My automated fitter (Gaussfitting Cube Collapser) has come a long way. I now adaptively choose to fit 1, 2, or 3 Gaussian components to output to a data cube.
The purpose of that code is primarily to find a two-dimensional way to display information about the 3D structure, specifically about the presence/absence of outflows. Outflows will inevitably be confused with multiple velocity components, but they are also likely to be convolved with them.
Sunday, October 19, 2008
Screen, nohup, ssh, scp
I learned a lot about the above in the past day, but I didn't keep track of the links.
First, screen is very useful: it allows you to run any task, detach the screen, and let it run in the background. You can resume it later.
Example:
Second, it's a huge pain to type a password every time I use scp and ssh. The solution is to make a key on your computer and put it in the authorized_keys file.
Use ssh -v to figure out why it fails if it still asks you for a password. In one case, this happened because the computer I was using expected the id to be in ~/.ssh/identity instead of ~/.ssh/id_rsa. There may also be permissions issues (i.e. you want restrictive permissions on your ssh keys).
With the latter, you can still use nohup, which is helpful if you want to pipe your output to a log file.
First, screen is very useful: it allows you to run any task, detach the screen, and let it run in the background. You can resume it later.
Example:
screen
ipython run_fitter.py
<ctrl-a> d
screen -r
Second, it's a huge pain to type a password every time I use scp and ssh. The solution is to make a key on your computer and put it in the authorized_keys file.
ssh-keygen -t rsa
scp ~/.ssh/id_rsa ginsbura@milkyway.colorado.edu:.
ssh ginsbura@milkyway.colorado.edu
cat id_rsa >> ~/.ssh/authorized_keys
ssh -v ginsbura@milkyway.colorado.edu
Use ssh -v to figure out why it fails if it still asks you for a password. In one case, this happened because the computer I was using expected the id to be in ~/.ssh/identity instead of ~/.ssh/id_rsa. There may also be permissions issues (i.e. you want restrictive permissions on your ssh keys).
With the latter, you can still use nohup, which is helpful if you want to pipe your output to a log file.
Friday, October 17, 2008
Listing variables (e.g IDL help) in Python
Again, IDL has the simple 'help' command to tell you all variables in your namespace. Python has the same thing, but the namespace tends to be cluttered with imported functions. The commands who, who_ls, and whos are meant for interactive use. They are a hell of a lot more useful than var, locals, globals, and dir.
examples:
I'm afraid I don't know how to make the last two lines into a one-liner, as would be desirable.
examples:
whos float
whos ndarray
who module
float_vars = %who_ls float
grep('x',float_vars)
I'm afraid I don't know how to make the last two lines into a one-liner, as would be desirable.
SAVE / RESTORE in Python
Save/Restore is probably the single best feature of IDL that, sadly, is very poorly replicated in Python. For 1 or 2 dimensional variables, you can use Pylab's
e.g.:
save/load, but I never use such piddling tiny arrays. For higher dimensional objects, either using FITS files (a pain because of header definitions) or pickling ought to work.e.g.:
import numpy
import pickle
x=ones([10,10,10,10],dtype='float64')
pickle.dump(x,open('x.pysav','w'))
X = pickle.load(open('x.pysav','r'))
Wednesday, October 15, 2008
Last month's observations in L111
Mapped sources:
SH 2-156 [no high resolution observations; covered by 1.2m CO 1-0 survey]
SH 2-157 Yancy CS 5-4 , CO, 13CO 1-0, CS 2-1 , HCN & SCUBA too....
IRAS 23151 (two pointings)
Unmapped:
BFS 18 = IRAS 23042+6000 = IRCO 859
Wu et al mapped this in CO 2-1 and detected no outflows. They used the NRAO 12m with a heterodyne array with very high Tsys ~700-1100K.
SH 2-156 [no high resolution observations; covered by 1.2m CO 1-0 survey]
SH 2-157 Yancy CS 5-4 , CO, 13CO 1-0, CS 2-1 , HCN & SCUBA too....
IRAS 23151 (two pointings)
Unmapped:
BFS 18 = IRAS 23042+6000 = IRCO 859
Wu et al mapped this in CO 2-1 and detected no outflows. They used the NRAO 12m with a heterodyne array with very high Tsys ~700-1100K.
Sunday, October 12, 2008
How do I print a multi-page postscript?
I want to print a multi-page postscript duplexed and, if possible, two per page. pstopng, pstopdf, and the variants failed me horribly by making postscripts that don't fit the bounding box at all. Any ideas?
So far the best I've come up with is switching gnuplot's output to pdf, which is a workaround rather than a solution, but is useful nonetheless:
So far the best I've come up with is switching gnuplot's output to pdf, which is a workaround rather than a solution, but is useful nonetheless:
set terminal pdf enhanced dashed
Mathematica & Series expansions
Mathematica does something stupid and incomprehensible (to me) with their series expansions. For no clear reason, it makes the O[x]n uneditable, and adds a bunch of formatting things that prevent you from copying & pasting the expression and then solving for a variable. Argh. While this is repairable using the "show expression" tool (Shift-Apple-E), it's a pain and "show expression" is not meant to be human readable.
Saturday, October 11, 2008
Tuesday, October 7, 2008
gardening
I realized I don't have a gardening category and I should.
I have eaten 2 garden tomatoes. There have been a total of 3 that ripened. There is a fourth on the vine. There are ~40-50 big green ones being big and green.
I ate basil in my pasta sauce yesterday. I need to buy some pine nuts. But I'm not pining for nuts.
Next year, start growing in March!
I have eaten 2 garden tomatoes. There have been a total of 3 that ripened. There is a fourth on the vine. There are ~40-50 big green ones being big and green.
I ate basil in my pasta sauce yesterday. I need to buy some pine nuts. But I'm not pining for nuts.
Next year, start growing in March!
sage casjobs
annoying feature of sage casjobs: no spaces are allowed, no matter what, even if quoted. so when importing must remove all spaces. gurrh.
mathematica
//. and /. use 'rules' of the form {x->y} for substitutions
how do you make the layout of a page look nice? How do you change line spacing?
how do you make the layout of a page look nice? How do you change line spacing?
Saturday, September 27, 2008
Observing run summary
3 nights of worse than grade 4 weather (tau ~.4, i.e. the atmosphere is nearly opaque) have given me a few interesting results.
1. Non-detections of CS 5-4 in a few bright BGPS cores. My best guess is that the pointing was off, though, because I definitely got down to a noise level around 50mK.
2. Lots of CO 2-1 mapping, 13CO/C18O too. That means I can now use this plot: http://dustem.astro.umd.edu/pdrt/co3221web.html to try to measure the incident flux.
3. I can't do homework at the summit. Maybe it's just that I can't do stellar structures homework, but I really am not processing this stuff well. Damn.
4. Running at altitude uphill is really, really hard. All the running I've been doing the past N months has not even come close to preparing me. Need to correct that...
5. "I am inappropriate." I probably violate a bunch of those rules... the important ones...
6. gotta get back to the mountains. Really. Seriously. OMFG.
7. ARGH those papers need writing! Lots of it!
For the next observing run, need the following:
1. Pre-prepared A3 rasters
2. Larger selection of summit music
3. Better sleep schedule preparation
4. More caffeinated drinks per night
5. Cloud cannons.
1. Non-detections of CS 5-4 in a few bright BGPS cores. My best guess is that the pointing was off, though, because I definitely got down to a noise level around 50mK.
2. Lots of CO 2-1 mapping, 13CO/C18O too. That means I can now use this plot: http://dustem.astro.umd.edu/pdrt/co3221web.html to try to measure the incident flux.
3. I can't do homework at the summit. Maybe it's just that I can't do stellar structures homework, but I really am not processing this stuff well. Damn.
4. Running at altitude uphill is really, really hard. All the running I've been doing the past N months has not even come close to preparing me. Need to correct that...
5. "I am inappropriate." I probably violate a bunch of those rules... the important ones...
6. gotta get back to the mountains. Really. Seriously. OMFG.
7. ARGH those papers need writing! Lots of it!
For the next observing run, need the following:
1. Pre-prepared A3 rasters
2. Larger selection of summit music
3. Better sleep schedule preparation
4. More caffeinated drinks per night
5. Cloud cannons.
Thursday, September 25, 2008
running at altitude
wow dude. I suck at that.
I swear there were some other things I had to write down but I can't remember them now.
precipitable water vapor precipitates and that can be good
I swear there were some other things I had to write down but I can't remember them now.
precipitable water vapor precipitates and that can be good
pwv
precipitable water vapor, precipitable water vapor, go away, come again some other observer's day
little HARPy wants to receive
so clear up before I leave
stupid bumping-head old man.
little HARPy wants to receive
so clear up before I leave
stupid bumping-head old man.
Wednesday, September 24, 2008
Marc Pound's toolbox
This is neat:
Line Ratio as a function of density, incident flux from Marc Pound. Sadly, straight CO 3-2 intensity isn't included.
Line Ratio as a function of density, incident flux from Marc Pound. Sadly, straight CO 3-2 intensity isn't included.
python stuff
to get imshow defaults to be nearest neighbor, need to edit ~/.matplotlib/matplotlibrc. Still don't know how to change default command-line output format.
Do not install 10.5.5
Well, I didn't obey that excellent advice from Charles. I installed the 10.5.5 update. Now DS9 - in my opinion, one of the heartiest and most reliable programs I've dealt with in the astronomy world - won't work because TCL is broken. I doubt TUI or most of the other things I use will work either. That's a tragedy of significant proportions - I NEED ds9 to make my observing files! I can only pray that an update of TCL using port will work... if not I'm basically screwed.
Saturday, September 13, 2008
PHP and CGI working in the same directory
While it's not possible to get CGI and PHP to work within the same file, you can get .html and .htm files to use cgi/php scripts alternately in the same directory with the following .htaccess file:
If anyone can prove the previous statement false, please do! I spent hours googling and hacking to try to get php and cgi to work within the same file.
Also, the way the above file is parsed makes no sense to me. It is inverted from how I would have thought it should work.
RemoveHandler .html
AddType text/html .htm
AddHandler application/x-httpd-php .php .html
AddHandler server-parsed .htm .html
If anyone can prove the previous statement false, please do! I spent hours googling and hacking to try to get php and cgi to work within the same file.
Also, the way the above file is parsed makes no sense to me. It is inverted from how I would have thought it should work.
Thursday, September 11, 2008
Wednesday, September 10, 2008
no-coffee day make me useless
I had a few productive hours after my run, then things kinda went downhill.... in the morning, I wrote two guides and then taught someone else to use that stuff and then I tried to write up my lecture for tomorrow and totally failed. I didn't accomplish anything after about 1 pm. Lesson to me: don't ride no-coffee days at work.
Tuesday, September 9, 2008
Unhealthy outlier
162.7 lbs - only a 2-sigma outlier from my weight this past month, but definitely unhealthily low. Probably says something about whether I should have played disc tonight. Or something else.
CGI errors on Origins
First step:
Then, make sure .htaccess has the right lines.
tail -f /usr/local/apache/logs/casa/casa_error_log
Then, make sure .htaccess has the right lines.
AddType text/html .htm
AddHandler application/x-httpd-php .php .htm
AddHandler server-parsed .htm .html
AddHandler cgi-script .cgi .pl
Sunday, September 7, 2008
Baking again
Since I totally ran out of food last night (ate my last frozen pizza), I decided it's time to try my hand at cooking things again. That will start with making a bunch of dough, possibly shaping it, and freezing some, fridging some, etc. Homemade pizza is way better.
Need to get pine nuts. 'n stuff.
Need to get pine nuts. 'n stuff.
Friday, September 5, 2008
Hi Kyle!
How's google reader treating you today? I hope that's what you're using.
When you come back to observe at 11:30, can you bring extra beer?
When you come back to observe at 11:30, can you bring extra beer?
Thursday, September 4, 2008
Non-greedy vim
To make .* in VIM be non-greedy (i.e. match just <a href=x> out of "<a href=x> </a>"), use:
.\{-}
List
I have something of a list going of 'hottest girls who won't talk to me'. I think I'll award prizes.
1st prize - the hoppiest beer I can find
2nd prize - baking chocolate
3rd prize - extra dark, extra fine-ground coffee
4th prize - angosteen bitters
1st prize - the hoppiest beer I can find
2nd prize - baking chocolate
3rd prize - extra dark, extra fine-ground coffee
4th prize - angosteen bitters
Wednesday, September 3, 2008
Using AWK to fixed-format print in VIM
24,42!awk '{printf("\%30s\%20f\%15g\%15g\%15g\%15g\%15g\%15g\%15g\%5i\n" , $1,$2,$3,$4,$5,$6,$7,$8,$9,$10)}'It's important that the %'s must be escaped, otherwise they print the current filename.
IDL colors
To get IDL colors to 'work right' (i.e. a colortable has 255 colors etc), put this in a startup file:
device, true_color = 24, retain = 2, decompose = 0
red = [0,1,1,0,0,1]
green = [0,1,0,1,0,1]
blue = [0,1,0,0,1,0]
if not strcmp(getenv('DISPLAY'),'') then $
tvlct, 255*red, 255*green, 255*blue
Monday, September 1, 2008
1 year of data
Weight in pounds as a function of time, starting 8/12/07, with nearly uniform sampling immediately before bed each night.
Sunday, August 31, 2008
Saturday, August 30, 2008
Saturday, August 23, 2008
Merging postscripts
This is essential and really difficult to find answers to, but this guy gave it:
http://ludo.qix.it/archive/2005/08/merge-postscript-files.html
The keywords I would have liked to see:
"merge postscripts into multi-page document"
or "combine postscript multiple page"
e.g.
http://ludo.qix.it/archive/2005/08/merge-postscript-files.html
The keywords I would have liked to see:
"merge postscripts into multi-page document"
or "combine postscript multiple page"
e.g.
gs -sDEVICE=pswrite -sOutputFile=output.ps -dNOPAUSE -dBATCH file1.ps file2.ps file3.ps
or
gs -sDEVICE=pswrite -sOutputFile=05358spectra.ps -dNOPAUSE -dEPSFitPage -dBATCH `ls 05358_*.eps`
(the added option is to make sure the .eps isn't cropped)
NSF eligibility
Am I eligible for the NSF GRFP? Sam says yes, I thought no, but re-reading.... umm... isn't this yes?
No more than 12 months means 12 months is OK?
What was I thinking!?
Applicants must have completed no more than twelve (12) months of full-time graduate
study or its equivalent as of August 31, 2008. Full time graduate study is as defined by the
universities attended.
No more than 12 months means 12 months is OK?
What was I thinking!?
Friday, August 22, 2008
Back to site editing
I'm editing class web pages, php, cgi, etc. again. I don't think I ever really understood how to do this, but I'd like to know how to use PHP to fold text. That would be sweet.
I found it by googling this time.
This is really the right way to do things:
http://help.blogger.com/bin/answer.py?hl=en&answer=42214
Of course, it's pretty disappointing that it doesn't work with the non-classic template for blogger.
I found it by googling this time.
This is really the right way to do things:
http://help.blogger.com/bin/answer.py?hl=en&answer=42214
Of course, it's pretty disappointing that it doesn't work with the non-classic template for blogger.
Thursday, August 21, 2008
Extracting figures from papers
I want to automatically extract figures from papers, but of course most papers come in .pdf, not .ps form...
This site has a great tool for extracting figures from .ps papers:
http://www-hermes.desy.de/notes/extras/postscript.html
If anyone can tell me how to do the same thing for .pdfs (automatically extract figures from pdfs) I'd appreciate it.
Update 12/13/08: One solution is listed at http://ubuntuforums.org/showthread.php?t=884002: use pdfimages, part of the xpdf package.
This site has a great tool for extracting figures from .ps papers:
http://www-hermes.desy.de/notes/extras/postscript.html
If anyone can tell me how to do the same thing for .pdfs (automatically extract figures from pdfs) I'd appreciate it.
Update 12/13/08: One solution is listed at http://ubuntuforums.org/showthread.php?t=884002: use pdfimages, part of the xpdf package.
Sunday, August 10, 2008
Movie Review: Dark Knight
Definitely the best movie in the past... 3 years? IMAX didn't add much to the experience though.
A little bit more: most people I've talked to about it (ummm... 2 grad students and 1 blog?) thought the movie was overhyped and didn't live up to the hype. I liked it a lot, and I never saw a single preview. Lesson? Don't see previews.
A little bit more: most people I've talked to about it (ummm... 2 grad students and 1 blog?) thought the movie was overhyped and didn't live up to the hype. I liked it a lot, and I never saw a single preview. Lesson? Don't see previews.
Saturday, August 9, 2008
Blogroll
Since you're now subscribing (I'm assuming, at the moment and for the next... say... year, Kyle's the only one who will see this), here's a list of other blogs to check out:
http://sketchbrewing.blogspot.com/
http://thedrog.blogspot.com/
I'd send you to some other ones too, but they're not complete. I have vague ambitions to write a genuine science blog for the international year of astronomy and public outreach.... but I don't really know if that will go anywhere.
http://sketchbrewing.blogspot.com/
http://thedrog.blogspot.com/
I'd send you to some other ones too, but they're not complete. I have vague ambitions to write a genuine science blog for the international year of astronomy and public outreach.... but I don't really know if that will go anywhere.
Movie Review: Hancock
Saw Hancock on Tuesday. Basically, as hinted at by friends, didn't live up to expectations. The funny quit early, and they tried to insert some 'dark superhero movie' into it in the second half. From the ads, and the presence of Jason Bateman & his character, this was inappropriate. Not a great movie, stick to the Dark Knight.
More beer tastes
Stone Vertical Epic 8: described as an overhopped Belgian, and I think they hit the mark exactly. It was pleasant and drinkable, with a strong hop flavor but weak hop bitterness. 8.6% ABV
Fort Collins brewery double chocolate stout: a reasonable dessert beer. Not nearly as dry as Stone's Oatmeal Chocolate, but not as sweet as Young's double chocolate. Chocolate undertones stayed undertones, which is not what I expect from a double chocolate stout. 8.1% ABV
Avery / Liquor Mart 40th anniversary: Just a hoppy ale, nothing special. Had a stronger hop aftertaste than most, not up to Avery's standard. No ABV label.
Steamworks Kolsch: Very light kolsch. I think it may have had a crisp flavor to it, but it was the last in line and unremarkable enough that my comments on it are untrustworthy. 4.86% ABV.
Also, talked to Bobby again, he recommends St Bernard as a good Belgian. That's on the list.
Fort Collins brewery double chocolate stout: a reasonable dessert beer. Not nearly as dry as Stone's Oatmeal Chocolate, but not as sweet as Young's double chocolate. Chocolate undertones stayed undertones, which is not what I expect from a double chocolate stout. 8.1% ABV
Avery / Liquor Mart 40th anniversary: Just a hoppy ale, nothing special. Had a stronger hop aftertaste than most, not up to Avery's standard. No ABV label.
Steamworks Kolsch: Very light kolsch. I think it may have had a crisp flavor to it, but it was the last in line and unremarkable enough that my comments on it are untrustworthy. 4.86% ABV.
Also, talked to Bobby again, he recommends St Bernard as a good Belgian. That's on the list.
Wednesday, August 6, 2008
Applescripting out that VIM error
I have a frequent problem where my VIM window is too large for my macbook screen but it fits on my external monitor, so when I unplug the external the VIM file bar gets stuck behind the Mac OS menu bar. It's a huge pain to fix this normally, but I wrote/stole an applescript to fix the problem:
try
tell application "Vim"
activate
end tell
tell application "System Events"
tell process "Vim"
set size of the first window to {1000, 200}
set position of the first window to {50, 50}
end tell
end tell
end try
Tuesday, August 5, 2008
Defining a Molecular Cloud
I'm going to try to answer a few questions [this post is in-progress and will be updated]:
1. Molecular clouds are regions in space with densities high enough and temperatures low enough that molecules can form. The most common molecule is Hydrogen (H2), simply because hydrogen is the most common atom. The next most common is Carbon Monoxide, CO, that is ~104 times less common. It is the easiest to observe, though, because H2 can only emit light and energy in 'forbidden' transitions that don't happen very often.
2. Size, temperature, and composition all distinguish molecular clouds from other nebulae.
3. It depends on what wavelength you're looking in. In the optical, where our eyes can see, they look dark - they're seen because they absorb light from behind them. In the infrared, the hot ones glow, but the cold ones are still invisible. If you go all the way to the millimeter, all molecular clouds glow, no matter how cold, but in order to see them there has to be a lot of material. I'll add some pictures here later.
4. Lots of ways. Optical telescopes aren't the best choice, though. With near-infrared images, we can detect molecular clouds by counting stars and noticing when there aren't as many in some regions. This technique is called "NICE" and works because there are so many stars in the sky and dust in molecular clouds includes some, but not all, of their light. In the far infrared, we can see hot dust glowing, but this is very difficult because our atmosphere is opaque in the infrared - it's like trying to look through a brick wall. We need telescopes in space to be able to see anything at these wavelengths. In the millimeter, there are some 'holes' in the atmospheric absorption, sort of like color filters, that we can see through to detect the coldest emission.
- Why are molecular clouds called molecular clouds?
- What distinguishes a molecular cloud from similar objects (e.g. cores, clumps, HII regions)?
- What do molecular clouds look like?
- How do astronomers look at molecular clouds?
- Why are astronomers look at molecular clouds (why are they interesting)?
- What do molecular clouds have to do with star formation and planet formation?
- How are molecular clouds involved in galactic evolution?
1. Molecular clouds are regions in space with densities high enough and temperatures low enough that molecules can form. The most common molecule is Hydrogen (H2), simply because hydrogen is the most common atom. The next most common is Carbon Monoxide, CO, that is ~104 times less common. It is the easiest to observe, though, because H2 can only emit light and energy in 'forbidden' transitions that don't happen very often.
2. Size, temperature, and composition all distinguish molecular clouds from other nebulae.
- An H II region is very hot - in it, hydrogen atoms are ionized by radiation from a nearby O-type star. In regions hot enough to ionize hydrogen, all molecules are destroyed unless they have already collected into larger dust particles.
- A core is a much denser bunch of gas than a molecular cloud. Cores are dense enough that their outsides shield their insides from the radiation of the rest of the universe and their insides cool to very low temperatures. At low temperatures, the pressure supporting the cloud against collapse is lower, and stars can form via gravitational collapse.
- I don't really know what clumps are... they're supposed to be something intermediate between cores and clouds, but I don't know what their distinguishing features are.
- A core is a much denser bunch of gas than a molecular cloud. Cores are dense enough that their outsides shield their insides from the radiation of the rest of the universe and their insides cool to very low temperatures. At low temperatures, the pressure supporting the cloud against collapse is lower, and stars can form via gravitational collapse.
3. It depends on what wavelength you're looking in. In the optical, where our eyes can see, they look dark - they're seen because they absorb light from behind them. In the infrared, the hot ones glow, but the cold ones are still invisible. If you go all the way to the millimeter, all molecular clouds glow, no matter how cold, but in order to see them there has to be a lot of material. I'll add some pictures here later.
4. Lots of ways. Optical telescopes aren't the best choice, though. With near-infrared images, we can detect molecular clouds by counting stars and noticing when there aren't as many in some regions. This technique is called "NICE" and works because there are so many stars in the sky and dust in molecular clouds includes some, but not all, of their light. In the far infrared, we can see hot dust glowing, but this is very difficult because our atmosphere is opaque in the infrared - it's like trying to look through a brick wall. We need telescopes in space to be able to see anything at these wavelengths. In the millimeter, there are some 'holes' in the atmospheric absorption, sort of like color filters, that we can see through to detect the coldest emission.
Molecular Clouds
I'm starting a blog on Molecular Clouds partly for the International Year of Astronomy. Ideally, I'd like to present a tour of the W3/4/5 region, but at the moment that's more of a research blog than a public blog. This is partly in response to Pamela Gay's request for astronomy blog feeds.
I've generally spent more time maintaining my personal web page than my blog pages. The most interesting are:
W5
P Cygni
Some pictures
but feel free to browse around.
I've generally spent more time maintaining my personal web page than my blog pages. The most interesting are:
W5
P Cygni
Some pictures
but feel free to browse around.
Sunday, August 3, 2008
Review: Fruition restaurant
I ate at Fruition restaurant last night with the family. They serve haute cuisine. I had the beef coulotte. It was pretty tasty, but overall a bit too creamy for my taste. The wine had a hint of cherry flavor to it. The dinner was pretty good, but overall a little disappointing for the price. The French press coffee was a bit much as I would have preferred a small cup of cappuccino. I wouldn't object to going back when they switch to a new seasonal menu, but I'm not going out of my way.
Saturday, August 2, 2008
Photos featured in an article
My photos, article by Jeff Marlow:
http://www.atmagazine.co.uk/content/view/853/94/
http://www.atmagazine.co.uk/content/view/853/94/
Friday, August 1, 2008
Movie Review: The Namesake
The Namesake was pretty damned boring. The rich white girl might have been flaky and a touch rich-bitchy, but she was attractive, she loved him, and she did nothing wrong except perhaps out of ignorance. Coming of age stories like this one are dumb and take away my respect for ultraconservative, traditional cultures.
The one thing I can really say for the movie, though, is that it had awesome music. Federico Aubele's Postales showed up at one point, and there was another techno/ambient style song I couldn't pin down.
The one thing I can really say for the movie, though, is that it had awesome music. Federico Aubele's Postales showed up at one point, and there was another techno/ambient style song I couldn't pin down.
Thursday, July 31, 2008
Postscript to PDF conversion
ls *.ps | sed 's/\(.*\).ps/ps2pdf \1.ps \1.pdf/' | bashbecause it's impossible to view multiple postscripts in a single window.... how annoying.
Corollary:
ls *.pdf | sed 's/\(.*\).pdf/pdf2jpeg \1.pdf \1.jpg/' | bash
Wednesday, July 30, 2008
Boulder Obovoid
Obovoid is another imperial stout - "empirical stout" - and is milder than some of the others I've had lately, but still a touch bitter. Nice but not worth repeating much.
IDL string conversion
Just a reminder to me:
to convert an IDL number to a string, use string!
Also,
to convert an IDL number to a string, use string!
s = string(x,format='(F10.5)')Also,
Julian Date conversion
Tuesday, July 29, 2008
find and xargs
pretty useful info on find and xargs, e.g.
find /scratch/adam_work/ -name 'mmm*sav' | xargs rm
SETI
I read a paper on ArXiv today about SETI, and talked to Jeremy Darling about water masers. This got me thinking: the best way to signal your presence to the universe as an intelligent civilization is to take advantage of natural amplification. So, say, you go find a big cloud of water molecules in the right state to start masing, and blast them with some of the right frequency radio emission: you get a strongly enhanced isotropic emitter. Maybe SETI should look for those, then search for nearby coherent sources....
EDIT: Isotropic? No.
EDIT: Isotropic? No.
Monday, July 28, 2008
Hebrew Jewbilation 11
Wow... 11 malts and 11 hops, and it tastes like at least 3 different types of beer. There was an ale start and a porter finish to the first few sips, but later ones had a fuller, maltier flavor at the top of the mouth. It has some Irish Stout smoothness too. There's some residual porter-esque bitterness. A lot of depth and complexity that keeps going. Whew.
Sunday, July 27, 2008
Delirium Nocturne
Delirium Tremens is a tasty beer itself, I think. I'm drinking Delirium Nocturne right now. It has a cool (opaque!) bottle, plus a pink elephant on the label. At 8.5% it's pretty strong, but very drinkable with no taste of alcohol. It's definitely my favorite Belgian brown now. It still has that wheat-soy-sauce taste to it, but very weakly. It poured with ridiculous head - I could not control it. Overall, great beer, worth the ridiculous price I paid for it, whatever that was.
Saturday, July 26, 2008
More tastings
Deschutes Black Butte XX - good for a porter, reasonably drinkable, but still a porter
Rogue Chocolate Stout - a bit bitterer than I want in a chocolate stout
Victory 12 - pretty good, tasted a bit like a trippel
Sketch Wheat - much more drinkable than last time; "weirdness" gone
Rogue Chocolate Stout - a bit bitterer than I want in a chocolate stout
Victory 12 - pretty good, tasted a bit like a trippel
Sketch Wheat - much more drinkable than last time; "weirdness" gone
Thursday, July 24, 2008
Guinness and Dogfish Head
First, the Guinness: I'm drinking a Guinness Extra Stout right now. I'm a big fan of the normal Guinness Stout, but the extra stout is bitter - bridging on skunky. I don't like it. It has a sour edge.
Last night, tried the Dogfish Head Midas Touch. It's interesting - lots of flavors, topped by grape juice. It's a champagne-meade-beer. OK, but not spectacular. Pretty drinkable.
Dogfish Head Santo Palo tonight. It's a brown ale that looks more like a stout. I think it needs to mellow - it has an astringent alcohol flavor, but depth once you get further in. It was a lot better when it poured with some head than not. A little chocolate at the end, not much though.
Also tried the Avery "Ale to the Chief" Presidential Pale Ale. A little like an IPA, but really a hoppy pale ale. Not bad, not my favorite.
Those Dogfish beers are too expensive. I ought to stick to Avery and homebrews....
Last night, tried the Dogfish Head Midas Touch. It's interesting - lots of flavors, topped by grape juice. It's a champagne-meade-beer. OK, but not spectacular. Pretty drinkable.
Dogfish Head Santo Palo tonight. It's a brown ale that looks more like a stout. I think it needs to mellow - it has an astringent alcohol flavor, but depth once you get further in. It was a lot better when it poured with some head than not. A little chocolate at the end, not much though.
Also tried the Avery "Ale to the Chief" Presidential Pale Ale. A little like an IPA, but really a hoppy pale ale. Not bad, not my favorite.
Those Dogfish beers are too expensive. I ought to stick to Avery and homebrews....
Tuesday, July 22, 2008
Monday, July 21, 2008
Parallactic Angle
So it turns out parallactic angle and field rotation are essentially the same thing. While working on the data pipeline for the Bolocam Galactic Plane Survey, I've needed to deal with this. The parallactic angle is dependent on the alt/az, and therefore the time, of the observation. So using an RA/Dec in J2000 coordinates to get alt/az isn't going to work, because my observations weren't in 2000.0.... my best hope is that the PA in the NCDF files is actually the correct parallactic angle.
YARGH!
YARGH!
Stone 12th Anniversary Bitter Chocolate Stout
Tried this 9.2% bomber last night, and even drinking it over the course of an hour and a half I could feel my nose not feeling afterwards.
The pour was probably the most exciting part of this beer. It came out an extremely dark, rich, thick motor-oil looking thing with nearly no head. The flavor, though, is of a very bitter and very dry imperial stout. I didn't get much of the chocolate tone, because I only taste chocolate when it is sweetened - bakers chocolate doesn't taste like chocolate to me. I think this beer would have been better if kept a little bit sweeter, much unlike the other Imperial stouts I've tried lately.
The pour was probably the most exciting part of this beer. It came out an extremely dark, rich, thick motor-oil looking thing with nearly no head. The flavor, though, is of a very bitter and very dry imperial stout. I didn't get much of the chocolate tone, because I only taste chocolate when it is sweetened - bakers chocolate doesn't taste like chocolate to me. I think this beer would have been better if kept a little bit sweeter, much unlike the other Imperial stouts I've tried lately.
Sunday, July 20, 2008
Friday, July 18, 2008
Observing this weekend
whoops, never made this post. Turns out I got rained out - 3 inches of rain on my night. Sad, no DIS time for P Cyg.
Lunar Occultation Hi-Res measurements
The folks at the VLT have come up with a means of achieving extremely high resolution from the ground: Observe when the moon occults a source and how long the occultation takes. Sweet.
This was applied in the galactic center first, of course:
Observations of binaries in GC
Additional observations in GC
This was applied in the galactic center first, of course:
Observations of binaries in GC
Additional observations in GC
ds9 commands to test distortion mapping
trying to figure out whether I'm screwing something up in the distortion map phase... the -scale limits command here is new and useful
ds9 testfields_*pix5*_map0.fits -cmap sls -scale limits 0 .005 -zoom 4 -match frames wcs -match scales -match colorbars &
ds9 testfields_*pix15*_map0.fits -cmap sls -scale limits 0 .003 -zoom 4 -match frames wcs -match scales -match colorbars &
ds9 testfields_*pix10*_map0.fits -cmap sls -scale limits 0 .004 -zoom 4 -match frames wcs -match scales -match colorbars &
Russell Owen, DS9
As I might have expected, Russell Owen from the University of Washington - the man responsible for the Telescope User Interface, TUI, on the Apache Point Observatory telescope - has written up a help page for XPA access points on ds9. He knows his stuff... I think I need to aspire to that.
World Wide Stout
Last night, had ~4oz of Dogfish Head's World Wide Stout. It was the Nov 2007 batch and was WAY too sweet. Just like the Stone Russian Imperial Stout, it was almost undrinkably sweet, so I recapped it and am going to "cellar" it (at 78 degrees, yikes) for a few months. The Stone brew turned out a lot better on the second try, I hope the Dogfish Head does too.
I also had a bottle of Young's Double Chocolate Stout. Oddly, they don't have it on their website. It has a strong chocolate flavor at the start, but it fades quickly into a generic, strongly carbonated stout - it's nothing special after the first taste.
I also had a bottle of Young's Double Chocolate Stout. Oddly, they don't have it on their website. It has a strong chocolate flavor at the start, but it fades quickly into a generic, strongly carbonated stout - it's nothing special after the first taste.
Thursday, July 17, 2008
Mac terminal stuff
xterm-color has some annoying features when SSHing to other compys, including 'top' not functioning properly. The error I receive:
top: no termcap entry for a `xterm-color' terminal
The solution:
From MacOSXhints, in bash just add export TERM="vt100" to your .profile file.
top: no termcap entry for a `xterm-color' terminal
The solution:
From MacOSXhints, in bash just add export TERM="vt100" to your .profile file.
Awk sexagesimal to decimal conversion
In VIM I often need to convert columns of RA/Dec from Sexagesimal into Decimal format.
%!awk '{ra = ($2+$3/60+$4/3600)*15; dec = $6+$7/60+$8/3600; print $1,"ra=",ra,", dec= ",$5,dec}'
The far more irritating inverse operation:
%!awk '{h=($2/15); h=h-(h\%1); m=($2-h*15)/15*60; m=m-(m\%1); s=($2-h*15-m*15/60)/15*3600; d=-($3-$3\%1); am=(-$3-d)*60; am=am-(am\%1); as=(-$3-d-am/60)*3600; printf "\%s \%02i:\%02i:\%02.2f
, -\%02i:\%02i:\%02.2f \%s \%s \%s\n" , $1,h,m,s,d,am,as,$4,$5,$6}'
%!awk '{ra = ($2+$3/60+$4/3600)*15; dec = $6+$7/60+$8/3600; print $1,"ra=",ra,", dec= ",$5,dec}'
The far more irritating inverse operation:
%!awk '{h=($2/15); h=h-(h\%1); m=($2-h*15)/15*60; m=m-(m\%1); s=($2-h*15-m*15/60)/15*3600; d=-($3-$3\%1); am=(-$3-d)*60; am=am-(am\%1); as=(-$3-d-am/60)*3600; printf "\%s \%02i:\%02i:\%02.2f
, -\%02i:\%02i:\%02.2f \%s \%s \%s\n" , $1,h,m,s,d,am,as,$4,$5,$6}'
Wednesday, July 16, 2008
Outreach Observing
Some notes for this Friday...
the ISS will rise at about 10:55 on Friday, July 18th, and 9:43 on Saturday, July 19th.
Jupiter will be in the sky, but Venus and Mercury will not. Jupiter will rise in the Southeast
the ISS will rise at about 10:55 on Friday, July 18th, and 9:43 on Saturday, July 19th.
Jupiter will be in the sky, but Venus and Mercury will not. Jupiter will rise in the Southeast
Denver, CO 2008/07/16 10:15 pm Sunrise/set: 05:46 am / 08:26 pm
Predicted passes for: ISS daylight-saving time enabled
Local Duration
Date Time (min) Approach Max. elevation Departure
------------------------------------------------------------------------
2008/07/17 02:57 am <1 10° above N 10° above N 10° above NNE
2008/07/17 04:30 am 5 12° above NNW 31° above NE 10° above E
2008/07/17 10:30 pm <1 13° above WSW 15° above WSW 15° above WSW
2008/07/18 03:19 am 4 10° above NNW 16° above NNE 10° above ENE
2008/07/18 04:54 am 5 15° above NW 78° above W 10° above SE
2008/07/18 09:19 pm 5 13° above SSW 29° above SE 11° above ENE
2008/07/18 10:55 pm 4 17° above WNW 25° above NNW 10° above NNE
2008/07/19 03:42 am 5 13° above NNW 32° above NNE 12° above E
2008/07/19 05:18 am 4 16° above W 22° above SW 11° above S
2008/07/19 09:43 pm 5 19° above WSW 71° above N 10° above NE
google accident
I'm still trying to figure out why google doesn't link to my page for the term died in a google accident, but we'll see if this page ever changes that. I link to XKCD #369 from that page with that phrase.
meta
Why a blog? I've been strongly opposed to blogs for the longest time (after all, Desert Strike coined the word BLOG as an acronym for building), so why use one now?
1. link power - more links to my own pages if done correctly
2. google owns the blogger
2.a. google made it really easy
3. it's easier to add content here than to my real 'web page' which is already too content-heavy. Still, my web page is more useful for displaying code, my cv, etc.
4. needed a place to post beer reviews / brewing related content that doesn't fit on a normal web page unless I made a blog there
5. easiest way to take notes that I can't just accidentally delete...
and now that I've made this post,
6. because I like metahumour
1. link power - more links to my own pages if done correctly
2. google owns the blogger
2.a. google made it really easy
3. it's easier to add content here than to my real 'web page' which is already too content-heavy. Still, my web page is more useful for displaying code, my cv, etc.
4. needed a place to post beer reviews / brewing related content that doesn't fit on a normal web page unless I made a blog there
5. easiest way to take notes that I can't just accidentally delete...
and now that I've made this post,
6. because I like metahumour
DS9 - crosshairs on command line
ds9 pmm*_map0.fits -zscale -match scales -zoom 4 -match frames wcs -crosshair 17:33:02.7 -13:04:49.5 wcs fk5 -lock crosshairs wcs &
ds9 pmp*_map0.fits -zscale -match scales -zoom 4 -match frames wcs -crosshair 17:33:02.7 -13:04:49.5 wcs fk5 -lock crosshairs wcs &
ds9 mpp*_map0.fits -zscale -match scales -zoom 4 -match frames wcs -crosshair 17:33:02.7 -13:04:49.5 wcs fk5 -lock crosshairs wcs &
ds9 mpm*_map0.fits -zscale -match scales -zoom 4 -match frames wcs -crosshair 17:33:02.7 -13:04:49.5 wcs fk5 -lock crosshairs wcs &
DS9's crosshairs are extremely useful for checking on WCS coordinate matching, especially when they can be set precisely using the command line. I don't know how to set the coordinates exactly interactively... that may come later.
ds9 pmp*_map0.fits -zscale -match scales -zoom 4 -match frames wcs -crosshair 17:33:02.7 -13:04:49.5 wcs fk5 -lock crosshairs wcs &
ds9 mpp*_map0.fits -zscale -match scales -zoom 4 -match frames wcs -crosshair 17:33:02.7 -13:04:49.5 wcs fk5 -lock crosshairs wcs &
ds9 mpm*_map0.fits -zscale -match scales -zoom 4 -match frames wcs -crosshair 17:33:02.7 -13:04:49.5 wcs fk5 -lock crosshairs wcs &
DS9's crosshairs are extremely useful for checking on WCS coordinate matching, especially when they can be set precisely using the command line. I don't know how to set the coordinates exactly interactively... that may come later.
Tuesday, July 15, 2008
Mountain Sun
Went to the Mountain Sun, had the Thunderhead Stout and Stone's Levitation IPA.
The Thunderhead Stout is one of my favorites. It's rich, slightly sweet, full bodied. I can't identify any particular flavors, but it's delicious.
Stone's Levitation IPA is dark brown for an IPA. It's hoppy and pretty tasty, but nowhere near as good as the Dogfish Head IPAs.
The Thunderhead Stout is one of my favorites. It's rich, slightly sweet, full bodied. I can't identify any particular flavors, but it's delicious.
Stone's Levitation IPA is dark brown for an IPA. It's hoppy and pretty tasty, but nowhere near as good as the Dogfish Head IPAs.
SSH
Learned something knew about ssh... keys are stored in ~/.ssh/known_hosts. If a remote key changes, you may be denied access if your saved key is wrong.
The System
Cat and Girl is a webcomic probably best described as being subversive.
Today's comic reveals that my role is to perpetuate a system. It's always nice to analyze your workplace and world on a deeper level. I am perpetuating evil, but the alternatives are scary.
Also, just to make sure there's a tie back to beer, Goats has a litany for beer.
Today's comic reveals that my role is to perpetuate a system. It's always nice to analyze your workplace and world on a deeper level. I am perpetuating evil, but the alternatives are scary.
Also, just to make sure there's a tie back to beer, Goats has a litany for beer.
Explosion
Our "Workie Ticket", which was doing its primary in the carboy, popped its top. It made a mess but should still be healthy... the yeast was going extremely strong at the way-too-hot 78 degrees it's sitting at this morning. It already has a fruity aroma, though, so it's likely that all of our brews are going to be extremely fruity. Bleh.
Our Honey Brown , which was excellent at the start, has gotten a bit too dry and alcoholic. The first few sips have a strong flavor of alcohol to them. However, once you get into the beer, it gains a little bit of depth and isn't too bad... probably the alcohol talking. I think most of the bottles are probably overcarbonated at this point, but we only have a few left. Overall point: Honey beers need to be drunk within about a week of bottling for peak deliciousness; apparently honey dries out very strongly.
Our Honey Brown , which was excellent at the start, has gotten a bit too dry and alcoholic. The first few sips have a strong flavor of alcohol to them. However, once you get into the beer, it gains a little bit of depth and isn't too bad... probably the alcohol talking. I think most of the bottles are probably overcarbonated at this point, but we only have a few left. Overall point: Honey beers need to be drunk within about a week of bottling for peak deliciousness; apparently honey dries out very strongly.
Monday, July 14, 2008
Pumphouse Brewery Review
Went to the pumphouse brewery last night and 'tasted' 4 beers.
First, for a brewpub, it's pretty weird to serve tastes in single-sip servings, but that's what they did. They didn't charge, so I didn't object, but it was difficult to get a strong sense for what the beers were meant to be.
Beers:
Cherry Saison: light colored, mild cherry flavor, no sweetness or depth to match the cherry flavor
Ginger saison: strong ginger flavor, no depth or interesting beer flavor. Like a dry but strong ginger ale
Saison: just... average
Copper: Weak flavored reddish beer
First, for a brewpub, it's pretty weird to serve tastes in single-sip servings, but that's what they did. They didn't charge, so I didn't object, but it was difficult to get a strong sense for what the beers were meant to be.
Beers:
Cherry Saison: light colored, mild cherry flavor, no sweetness or depth to match the cherry flavor
Ginger saison: strong ginger flavor, no depth or interesting beer flavor. Like a dry but strong ginger ale
Saison: just... average
Copper: Weak flavored reddish beer
Tuesday, July 8, 2008
Buffalo the Destroyer!
Buffalo the Destroyer!
Adam Ginsburg
I'm not going to explain this nickname any more than I explained my other nickname, Keflavich
Adam Ginsburg
I'm not going to explain this nickname any more than I explained my other nickname, Keflavich
Subscribe to:
Posts (Atom)