double2short
double2short(double *&F, short *&T) { *T++ = (short)dbl2long(*F++, 32767); }
double2short(double in) // limit double to short