BEGIN{ printf("EFLD 1.0\n") printf("# input catalog for sky-camera mapping\n") printf("LABEL %s\n", fieldname) printf("UTDATE 2008 08 14.000\n") printf("CENTRE %s %s %s %s %s %s\n", rah, ram, ras, deh, dem, des) printf("EQUINOX J2000.0\n") printf("AGMODE Point-&-Stare\n") printf("OBS_DURATION 10\n") printf("POS_ANGLE 0\n") printf("WAVELENGTH 900\n") printf("\n") } { if($1>0 && $2<30.0 && $14<14){ printf("%s ",$1) RA = $10 DE = $11 dRA = $17*7.9 /1000/3600 dDE = $18*7.9 /1000/3600 # with proper motion correction # RAnew = RA + dRA/cos(DE*3.1415926/180.0) # DEnew = DE + dDE RAnew = RA DEnew = DE if(DEnew<0){DEnew=-DEnew; DEsig=-1} else{DEnew=DEnew; DEsig=1} RAh = int( RAnew / 15 ) RAm = int( (RAnew/15 - RAh) * 60.0 ) RAs = ((RAnew/15-RAh)*60.0-RAm)*60.0 DEh = int( DEnew / 1 ) DEm = int( (DEnew - DEh) * 60.0 ) DEs = ((DEnew-DEh)*60.0-DEm)*60.0 if(RAh>=10){ printf("%2d ", RAh) } else{ printf("0%1d ", RAh) } if(RAm>=10){ printf("%2d ", RAm) } else{ printf("0%1d ", RAm) } if(RAs>=10){ printf("%6.3f ", RAs) } else{ printf("0%5.3f ", RAs) } if(DEsig=="-1"){ printf("-") } else{ printf(" ") } if(DEh>=10){ printf("%2d ", DEh) } else{ printf("0%1d ", DEh) } if(DEm>=10){ printf("%2d ", DEm) } else{ printf("0%1d ", DEm) } if(DEs>=10){ printf("%6.3f ", DEs) } else{ printf("0%5.3f ", DEs) } # calibration stars if($2<5 && 13.0<$14 && $14<14){ printf("C 1 ") } else if($2>15 && $14<14.0){ printf("F 1 ") } else if(10<$14 && $14<12){ printf("F 1 ") } else if(12<$14 && $14<14){ printf("F 3 ") } else{ printf("F 6 ") } printf("%5.2f ",$14) printf("* ") printf("\n") } }