actually, i had forgotten what they are doing. they are calling DCRAW with -D, but then they are scaling the data somehow. perhaps it's only a right shift, the authors never elucidated what exactly they are doing. if its only a right shift, then its OK.
rob
edit: i did go back into the thread and the thing that got me started on this is that i noticed that saturated stars did not have the value 16383, which is what i expected from a 14-bit sensor. but if they are only shifting left 2 bits, i would expect to find 2'b1111_1111_1111_1100, or 65532. instead i was seeing values in the high 50k, which did not make sense to me. i started looking at DCRAW options and wondered if they were calling with -d, but someone set me straight - they are using -D but subsequently doing some kind of 14b->16b scaling in the background. still if it were as simple as a shift left, the 50k values do not make sense.