diff --git a/README.md b/README.md index 3f253dc..210fbf9 100644 --- a/README.md +++ b/README.md @@ -1,12 +1,28 @@ Satellite Tracking Toolkit +========= -Required libraries: -qfits-5.2.0: ftp://ftp.eso.org/pub/qfits/qfits-5.2.0.tar.gz -pgplot-5.2.2: http://www.astro.caltech.edu/~tjp/pgplot/ -gsl-1.15: ftp://ftp.gnu.org/gnu/gsl/gsl-1.15.tar.gz -wcslib-2.9 +Sattools is a collection of tools to facilitate Photographic and Video satellite tracking. -Required environment variables: -ST_COSPAR: COSPAR number -ST_DATADIR: path to sattools directory -ST_TLEDIR: path to TLE directory +Build +------ +* Clone locally the code repository +* Install common dependencies + * gfortran + * gcc + * libpng-dev + * libx11-dev + * libjpeg-dev + * libexif-dev +* Build & install required libraries + * qfits-5.2.0: ftp://ftp.eso.org/pub/qfits/qfits-5.2.0.tar.gz + * pgplot-5.2.2: http://www.astro.caltech.edu/~tjp/pgplot/ + * gsl-1.15: ftp://ftp.gnu.org/gnu/gsl/gsl-1.15.tar.gz + * wcslib-2.9: http://www.epta.eu.org/~bassa/wcslib-2.9.tar +* Run `make` on the sattools folder + +Run +----- +You will need to set the following environment variables to run sattools: +* `ST_COSPAR` COSPAR number +* `ST_DATADIR` path to sattools directory +* `ST_TLEDIR` path to TLE directory diff --git a/calibrate.c b/calibrate.c index eede038..f5d08f1 100644 --- a/calibrate.c +++ b/calibrate.c @@ -86,63 +86,68 @@ int main(int argc,char *argv[]) sprintf(starfile,"%s/data/tycho2.dat",env); // Decode options - while ((arg=getopt(argc,argv,"f:R:D:s:q:m:r:M:"))!=-1) { - switch(arg) { - - case 'f': - fitsfile=optarg; - break; - - case 'R': - strcpy(sra,optarg); - if (strchr(sra,':')!=NULL) - ra=15.0*sex2dec(sra); - else - ra=atof(sra); - break; - - case 'D': - strcpy(sde,optarg); - if (strchr(sde,':')!=NULL) - de=sex2dec(sde); - else - de=atof(sde); - break; - - case 's': - if (strchr(optarg,',')!=NULL) { - sscanf(optarg,"%f,%f",&sx,&sy); - } else { - sx=-atof(optarg); - sy=-sx; + if (argc>1) { + while ((arg=getopt(argc,argv,"f:R:D:s:q:m:r:M:"))!=-1) { + switch(arg) { + + case 'f': + fitsfile=optarg; + break; + + case 'R': + strcpy(sra,optarg); + if (strchr(sra,':')!=NULL) + ra=15.0*sex2dec(sra); + else + ra=atof(sra); + break; + + case 'D': + strcpy(sde,optarg); + if (strchr(sde,':')!=NULL) + de=sex2dec(sde); + else + de=atof(sde); + break; + + case 's': + if (strchr(optarg,',')!=NULL) { + sscanf(optarg,"%f,%f",&sx,&sy); + } else { + sx=-atof(optarg); + sy=-sx; + } + break; + + case 'q': + q=atof(optarg); + break; + + case 'm': + mag=atof(optarg); + break; + + case 'M': + rmask=atof(optarg); + break; + + case 'r': + rmax=atof(optarg); + break; + + case 'h': + usage(); + return 0; + break; + + default: + usage(); + return 0; } - break; - - case 'q': - q=atof(optarg); - break; - - case 'm': - mag=atof(optarg); - break; - - case 'M': - rmask=atof(optarg); - break; - - case 'r': - rmax=atof(optarg); - break; - - case 'h': - usage(); - return 0; - break; - - default: - usage(); - return 0; - } + } + } else { + usage(); + return 0; } // Read image @@ -224,6 +229,22 @@ int main(int argc,char *argv[]) if (c=='q') break; + // Help + if (c=='h') { + printf("Calibrates astrometry. Initially requires manual matching of at least three stars. Use 'a' to select star on the image, then 'b' to select star from the catalog (blue circles). If at least three stars are matched, use 'f' and 'm' to converge the calibration.\n\n"); + printf("q Quit\n"); + printf("a Select star on image\n"); + printf("b Select star from catalog\n"); + printf("c Center image on pixel\n"); + printf("f Fit calibration\n"); + printf("m Match stars using current calibration\n"); + printf("z/+ Zoom in on cursor\n"); + printf("x/- Zoom out on cursor\n"); + printf("p Plot sextractor catalog\n"); + printf("r Reset zoom\n"); + + } + // Select pixel catalog if (c=='a' && click==0) { i=select_nearest(cat,x,y); @@ -259,7 +280,7 @@ int main(int argc,char *argv[]) } // Zoom - if (c=='z' || c=='+') { + if (c=='z' || c=='+' || c=='=') { width/=1.25; xmin=x-0.5*width; xmax=x+0.5*width; @@ -276,7 +297,7 @@ int main(int argc,char *argv[]) } // Unzoom - if (c=='x' || c=='+' || c=='=') { + if (c=='x' || c=='-') { width*=1.25; xmin=x-0.5*width; xmax=x+0.5*width; diff --git a/data/desig.txt b/data/desig.txt index 8c5645a..63fc867 100644 --- a/data/desig.txt +++ b/data/desig.txt @@ -240,7 +240,6 @@ 00829 64038A 00831 64038C 00836 64040A -00858 64047A 00862 64047B 00870 64051A 00871 64051B @@ -415,7 +414,6 @@ 02150 66034A 02153 65098J 02154 66008E -02161 66013G 02173 66040A 02174 66040B 02176 66041A @@ -480,7 +478,6 @@ 02520 66089B 02529 65034D 02530 63014AJ -02531 63014AK 02580 66076C 02608 66110A 02610 66111A @@ -581,6 +578,7 @@ 03081 67127A 03093 68002A 03094 68002B +03121 63030F 03122 67092D 03123 67114D 03124 61015HX @@ -589,6 +587,7 @@ 03127 68002D 03129 68011A 03131 68011B +03132 63030G 03133 68012A 03137 68012B 03158 68023A @@ -714,6 +713,7 @@ 03841 69010D 03842 61015JF 03843 68097BG +03848 68055C 03874 68091BH 03889 69036A 03890 69037A @@ -893,7 +893,6 @@ 04643 70025AB 04644 70025AC 04645 70025AD -04648 70025AG 04649 70025AH 04650 70025AJ 04651 70025AK @@ -909,7 +908,6 @@ 04667 70089K 04668 70089L 04672 70025AV -04673 70025AW 04675 70025AY 04676 70025AZ 04677 70025BA @@ -1031,7 +1029,6 @@ 04934 70025GA 04935 70025GB 04938 69082HJ -04941 70091AJ 04942 70091AK 04943 70025GD 04945 70025GF @@ -1271,7 +1268,6 @@ 05986 65108Q 05987 67001S 05991 65108R -05992 65108S 05998 65108T 06014 70025MR 06019 72035A @@ -1405,7 +1401,6 @@ 06920 73086A 06921 73086B 06938 73088D -06939 73084D 06965 73098A 06966 73098B 06973 73100A @@ -1700,7 +1695,6 @@ 07902 75042B 07907 72058BS 07924 75052A -07930 72058BY 07933 72058CB 07937 75051C 07938 75051D @@ -2161,7 +2155,6 @@ 10234 77065BM 10235 77065BN 10240 77065BR -10246 77065BX 10250 77065CB 10251 77065CC 10252 77065CD @@ -2479,7 +2472,6 @@ 10991 78074A 10992 78074B 10998 78066D -11007 78080A 11027 78087A 11028 78087B 11042 78091A @@ -2645,7 +2637,6 @@ 11822 80044B 11841 80049A 11847 80050B -11849 80051B 11852 80052C 11862 80049F 11869 80056A @@ -2679,7 +2670,6 @@ 12032 80085A 12046 80087A 12065 80091A -12072 80093B 12078 80095A 12086 80095E 12087 80097A @@ -2724,7 +2714,6 @@ 12179 78026U 12181 78026W 12183 78026Y -12184 78026Z 12185 78026AA 12189 78026AE 12190 78026AF @@ -2806,7 +2795,6 @@ 12376 81031A 12384 81031D 12388 81033A -12389 81033B 12404 80022C 12406 78014F 12409 81037A @@ -2998,7 +2986,6 @@ 12849 81094B 12850 81069F 12851 81061F -12854 70103G 12855 81096A 12857 77012H 12858 77106C @@ -3020,11 +3007,9 @@ 12893 81031F 12894 81031G 12897 81102A -12903 81103A 12904 81103B 12906 77068F 12907 78083E -12909 79058F 12930 81107A 12932 81107C 12933 81108A @@ -3213,7 +3198,6 @@ 13518 70025NR 13524 73086GM 13525 73086GN -13526 74089EM 13527 74089EN 13528 74089EP 13529 76077EB @@ -3340,7 +3324,6 @@ 13908 81088H 13909 81088J 13910 81088K -13911 81088L 13912 65108V 13916 83021A 13917 83021B @@ -3497,7 +3480,6 @@ 14478 70025NW 14480 73086GR 14483 83111A -14484 83111B 14493 77065FJ 14496 77065FM 14497 77065FN @@ -3670,7 +3652,6 @@ 15139 84028F 15141 83066F 15144 84078A -15147 84079A 15152 84080A 15158 84081A 15159 84081B @@ -4032,7 +4013,6 @@ 16860 86052A 16864 86052B 16865 86052C -16866 86052D 16867 86052E 16870 86007F 16881 86055A @@ -4058,9 +4038,7 @@ 17024 77065FT 17025 81053KJ 17026 81053KK -17027 78026GU 17035 86062C -17037 86078D 17046 86082A 17065 86082F 17066 86086A @@ -4118,6 +4096,7 @@ 17240 86100B 17241 86101A 17242 86101B +17253 86026E 17254 86026F 17255 85056C 17268 86101C @@ -4512,7 +4491,6 @@ 18959 88020B 18960 88021A 18961 88021B -18962 63014DN 18964 65027AM 18965 65027AN 18967 65027AP @@ -4539,6 +4517,7 @@ 19048 88029C 19051 63014EA 19054 77065FX +19059 63014EG 19070 88033A 19071 88033B 19072 88033C @@ -4627,6 +4606,7 @@ 19306 76077FJ 19309 68097ES 19311 75004HP +19314 77012L 19316 76077FK 19317 68091DH 19324 88062A @@ -4801,12 +4781,12 @@ 20012 76077FM 20014 76067BT 20016 78100AQ -20023 63014EZ 20024 89039A 20025 89039B 20026 89039C 20028 89039E 20029 78026HS +20030 76067BU 20031 65027BB 20032 77121BW 20040 89041A @@ -4889,7 +4869,6 @@ 20330 89091A 20333 89091D 20344 89061D -20347 89080J 20348 72049C 20349 76077FR 20355 89090B @@ -5039,7 +5018,6 @@ 20801 78062D 20804 90083A 20805 90083B -20811 90065C 20826 90086A 20827 90086B 20830 90088A @@ -5057,7 +5035,6 @@ 20853 90081R 20855 90081T 20858 90081W -20859 90081X 20861 90081Z 20862 90081AA 20864 90081AC @@ -5112,7 +5089,6 @@ 20962 75100F 20963 90097B 20969 90081CB -20971 90081CD 20972 90081CE 20973 90081CF 20974 90081CG @@ -5390,6 +5366,7 @@ 21473 75052FH 21475 75052FK 21492 75052FM +21493 75052FN 21496 75052FR 21497 75052FS 21498 75052FT @@ -5441,7 +5418,6 @@ 21601 75052HV 21602 75052HW 21610 91050F -21622 87036K 21623 87036L 21634 90065L 21635 90065M @@ -5523,7 +5499,6 @@ 21820 91086B 21821 91087A 21824 91087D -21825 91082D 21826 91086C 21827 91086D 21833 91088A @@ -5570,7 +5545,6 @@ 21941 92021C 21942 92021D 21943 91009CE -21944 80032C 21948 90065S 21949 92023A 21955 91009CF @@ -5600,9 +5574,9 @@ 22007 92036B 22009 92037A 22011 92037C +22014 92039A 22015 92039B 22017 92040A -22020 92040D 22025 91009CJ 22027 92041A 22028 92041B @@ -5678,7 +5652,6 @@ 22245 92082A 22248 92082D 22253 92084A -22254 92084B 22266 84080E 22269 92088A 22272 92088D @@ -5770,7 +5743,6 @@ 22408 92093CZ 22409 92093DA 22410 92093DB -22411 92093DC 22415 92093DG 22416 92093DH 22417 92093DJ @@ -5814,7 +5786,6 @@ 22478 92093FT 22480 92093FV 22481 92093FW -22485 92093GA 22487 93008A 22488 93008B 22489 93009A @@ -5824,13 +5795,11 @@ 22494 92093GD 22495 92093GE 22496 92093GF -22498 92093GH 22499 92093GJ 22500 92093GK 22502 92093GM 22503 92093GN 22504 92093GP -22506 92093GR 22509 92093GU 22512 93010A 22513 93010B @@ -5848,7 +5817,6 @@ 22549 92093HK 22551 92093HM 22552 92093HN -22553 92093HP 22554 92093HQ 22557 93013A 22561 93014A @@ -5939,7 +5907,6 @@ 22783 93055B 22784 93055C 22787 93056A -22788 93056B 22796 93058B 22797 93058D 22801 93050C @@ -5965,7 +5932,6 @@ 22841 91009CL 22842 91009CM 22843 91009CN -22845 91086G 22849 92089D 22850 92021E 22871 93066A @@ -6028,7 +5994,6 @@ 23093 94024B 23094 94023F 23097 94026A -23104 92094C 23105 94029A 23106 94029B 23108 94030A @@ -6067,7 +6032,6 @@ 23230 94056A 23232 91015H 23233 94057A -23234 94057B 23237 92021G 23238 92021H 23239 92021J @@ -6086,6 +6050,7 @@ 23260 90081CR 23261 90081CS 23263 91009CU +23265 91009CW 23267 94060A 23270 94060D 23275 92093JF @@ -6095,6 +6060,7 @@ 23282 92093JG 23292 85118P 23293 85118Q +23297 85118U 23301 92072D 23302 92093JH 23303 93059G @@ -6194,7 +6160,6 @@ 23527 95012B 23528 95013A 23533 95015A -23534 95015B 23535 95015C 23536 95016A 23537 95016B @@ -6306,14 +6271,13 @@ 23814 96014A 23815 96014B 23816 96015A +23820 66077D 23827 96017A 23828 96017B -23830 66077G 23833 96019A 23836 78007D 23839 96020A 23840 96020B -23841 66077L 23842 96021A 23845 96021D 23846 96022A @@ -6352,12 +6316,8 @@ 24051 94029CY 24082 94029ED 24091 94029EN -24092 94029EP -24094 94029ER 24097 94029EU 24102 94029EZ -24114 94029FM -24115 94029FN 24126 94029FZ 24129 94029GC 24138 94029GM @@ -6371,7 +6331,6 @@ 24227 94029LA 24236 94029LK 24239 94029LN -24243 94029LS 24277 96046A 24278 96046B 24279 96046C @@ -6412,6 +6371,7 @@ 24665 96067A 24674 96070A 24675 96070B +24676 66077Q 24677 96071A 24678 96071B 24679 94085Y @@ -6518,7 +6478,6 @@ 24932 97049B 24936 97050A 24937 97050B -24938 92047K 24944 97051A 24945 97051B 24946 97051C @@ -6571,8 +6530,6 @@ 25053 97070E 25054 97070F 25057 90036C -25063 97074A -25064 97074B 25067 97075A 25068 97075B 25071 97076A @@ -6657,7 +6614,6 @@ 25289 98021E 25290 98021F 25291 98021G -25300 87062C 25306 98023A 25307 98023B 25308 98023C @@ -6830,7 +6786,6 @@ 25662 99014B 25666 99016A 25667 99016B -25668 92060E 25669 99017A 25673 99018A 25676 99019A @@ -6877,23 +6832,18 @@ 25793 61018F 25795 80089CJ 25798 89001K -25803 87079N 25804 87079P +25807 87079S 25808 87079T 25809 87079U -25818 89101Q 25822 84023D 25823 84023E 25824 84023F 25825 66040H -25826 66077T -25833 78100AV 25834 78100AW -25837 79104L 25838 81071G 25839 81071H 25840 81071J -25843 89062G 25844 89062H 25847 99036A 25850 99036D @@ -7018,7 +6968,6 @@ 26045 99073D 26048 65108AG 26049 65108AH -26050 65108AJ 26052 00001A 26054 00001C 26056 00002A @@ -7117,13 +7066,11 @@ 26222 99057DL 26225 99057DP 26227 99057DR -26228 99057DS 26231 99057DV 26233 99057DX 26234 99057DY 26235 99057DZ 26236 99057EA -26238 99057EC 26239 99057ED 26241 99057EF 26243 00019A @@ -7139,7 +7086,6 @@ 26272 99057FJ 26273 99057FK 26275 99057FM -26277 99057FP 26280 99057FS 26283 99057FV 26284 99057FW @@ -7150,7 +7096,6 @@ 26300 99057GL 26306 99057GP 26310 99057GT -26311 99057GU 26313 99057GW 26321 99057HE 26323 99057HG @@ -7218,7 +7163,6 @@ 26476 00048A 26477 00049A 26480 00049D -26481 00050A 26483 00051A 26486 00051D 26487 00052A @@ -7240,6 +7184,7 @@ 26512 81021P 26513 84069E 26514 84069F +26516 84069H 26517 84069J 26519 84069L 26520 84069M @@ -7248,12 +7193,10 @@ 26523 84069Q 26524 65108AK 26525 65108AL -26526 65108AM 26527 65108AN 26528 65108AP 26529 65108AQ 26530 65108AR -26532 79104N 26534 79104Q 26535 79104R 26536 00055A @@ -7303,9 +7246,7 @@ 26638 00081A 26639 00081B 26643 00082A -26645 84069S 26646 84069T -26647 84069U 26650 99057LG 26654 99057LL 26658 99057LQ @@ -7313,8 +7254,6 @@ 26666 01002A 26669 66040J 26670 66040K -26672 68055F -26675 80050E 26676 80050F 26677 80050G 26678 80085E @@ -7358,7 +7297,6 @@ 26761 01018A 26762 01018B 26763 65108AS -26764 65108AT 26765 65108AU 26766 01019A 26769 01019D @@ -7375,7 +7313,6 @@ 26787 85056M 26788 86010E 26789 99055D -26792 92009D 26794 91015L 26795 90102G 26799 89062J @@ -7399,6 +7336,7 @@ 26828 75052KG 26829 75052KH 26830 75052KJ +26831 75052KK 26832 75052KL 26833 75052KM 26834 75052KN @@ -7440,7 +7378,6 @@ 26906 01040B 26907 01040C 26911 81021Q -26916 83072D 26918 80085G 26919 80085H 26921 85093C @@ -7451,7 +7388,6 @@ 26930 01043B 26931 01043C 26932 01043D -26934 01044A 26936 01045A 26939 01045D 26948 01046A @@ -7461,16 +7397,12 @@ 26960 01049D 26962 66040M 26964 77047E -26967 93017E 26968 93016AJ 26975 78066F 26977 80057H -26978 80057J 26980 92093JK -26981 92072L 26982 87020DU 26985 01052A -26986 01052B 26987 01053A 26988 01053B 26989 01053C @@ -7494,6 +7426,7 @@ 27032 89006M 27033 89006N 27034 89006P +27035 89006Q 27036 89006R 27037 89006S 27038 89006T @@ -7517,7 +7450,6 @@ 27059 01058E 27060 01058F 27061 01058G -27063 01049F 27089 01049AH 27097 01049AR 27105 01049AZ @@ -7555,7 +7487,6 @@ 27298 02002A 27299 02002B 27305 01049KB -27317 01049KP 27353 01049MA 27367 02003A 27368 02003B @@ -7618,6 +7549,7 @@ 27476 02037G 27478 61028H 27479 63038P +27481 68055J 27482 90081CU 27489 97005G 27491 98043L @@ -7651,11 +7583,10 @@ 27561 02054C 27566 02055A 27567 02055B -27568 80034E 27569 81021R 27570 81021S 27571 81021T -27574 75052LC +27573 84069X 27575 77065GE 27576 75052LD 27577 75052LE @@ -7665,7 +7596,6 @@ 27583 01049NB 27585 01049ND 27586 01049NE -27590 01049NH 27593 88040C 27594 88040D 27595 88040E @@ -7691,14 +7621,12 @@ 27618 02060B 27619 02060C 27620 02060D -27624 68055K 27625 88098H 27626 88098J 27628 87068CK 27629 87068CL 27632 02062A 27633 02062B -27635 66077AA 27638 84069Y 27639 87068CM 27640 03001A @@ -7756,6 +7684,7 @@ 27727 76077FS 27728 76077FT 27729 77065GF +27731 78087C 27732 79104S 27733 81053MU 27734 81053MV @@ -7763,7 +7692,6 @@ 27736 81120H 27737 81120J 27738 81120K -27739 84069Z 27740 84114E 27741 85090C 27742 87020DV @@ -7778,6 +7706,7 @@ 27751 92093JN 27752 93014G 27760 94077M +27761 96056D 27762 99057MV 27763 00015B 27765 00040D @@ -7791,16 +7720,13 @@ 27785 68063B 27786 69036B 27787 70069B -27788 65089E 27789 65089F -27790 65108BE -27791 65108BF 27792 65108BG +27793 65108BH 27794 65108BJ 27795 79104T 27796 81031H 27797 81053MX -27798 81071M 27800 81071P 27801 81071Q 27802 81071R @@ -7836,8 +7762,6 @@ 27858 03036A 27859 03036B 27864 65108BK -27865 65108BL -27866 65108BM 27867 65108BN 27868 03037A 27869 03037B @@ -7887,6 +7811,7 @@ 27954 03044A 27958 87079AX 27976 88085W +27998 88085AU 28039 91003M 28040 91015P 28041 91075L @@ -7894,12 +7819,13 @@ 28050 03046B 28051 03046A 28054 03048A +28055 03048B +28056 03048C 28057 03049A 28058 03049B 28059 03049C 28060 03050A 28061 03050B -28064 99057ST 28067 99057MX 28068 99057MY 28069 99057MZ @@ -7917,7 +7843,6 @@ 28095 03054A 28096 03054B 28097 03054C -28098 03055A 28112 03056A 28113 03056B 28114 03056C @@ -7939,7 +7864,6 @@ 28156 00024E 28158 04004A 28161 04004D -28163 04005A 28164 04005B 28184 04007A 28187 04008A @@ -7963,7 +7887,6 @@ 28221 04012B 28222 04012C 28225 04012F -28226 04012G 28230 04014A 28232 04012J 28234 04015A @@ -7971,16 +7894,13 @@ 28238 04016A 28240 04015D 28243 02003E -28248 64083M 28252 04017A 28254 04018A 28255 04018B 28256 04010F 28265 66052E -28266 68055L 28270 81016L 28271 84069AG -28273 84069AJ 28274 84069AK 28275 88019F 28277 91068AS @@ -8011,7 +7931,6 @@ 28321 91082AJ 28322 91082AK 28323 91082AL -28324 91082AM 28325 91082AN 28326 91082AP 28327 91082AQ @@ -8132,21 +8051,14 @@ 28545 05005D 28546 04050B 28551 65108BS -28552 65108BT 28553 65108BU 28554 79095C 28555 79095D 28556 79104V -28557 79104W -28558 81071S 28562 88040F 28563 88098L -28564 89006AF 28565 89006AG 28566 89006AH -28567 89006AJ -28568 89006AK -28569 89006AL 28570 89006AM 28573 91003N 28574 91003P @@ -8166,12 +8078,10 @@ 28598 82116E 28599 81043H 28600 81043J -28601 81043K 28603 85013C 28608 68055M 28609 66040Y 28610 66040Z -28611 66040AA 28612 66040AB 28613 66040AC 28615 66040AE @@ -8189,7 +8099,6 @@ 28634 05010F 28637 05011B 28638 05012A -28642 05014A 28644 05015A 28645 05015B 28646 05016A @@ -8203,10 +8112,11 @@ 28660 05019B 28661 05019C 28664 66077AD +28671 68055U 28672 88074H +28673 88074J +28675 88074L 28676 88074M -28678 92047AB -28685 92047AJ 28702 05022A 28704 05023H 28706 90017C @@ -8221,27 +8131,19 @@ 28786 05028A 28787 05028B 28790 05030A -28793 65063C -28795 65063E 28796 65063F -28798 65063H 28799 65063J -28804 65108BW 28805 65108BX -28807 65108BZ -28808 65108CA +28806 65108BY 28809 05031A 28810 05031B 28811 05031C 28812 05031D 28813 05031E +28814 65108CB 28815 65108CC 28817 65108CE -28820 65108CH -28821 65108CJ 28822 05032A -28834 65108CL -28835 65108CM 28837 65108CP 28838 65108CQ 28839 65108CR @@ -8313,7 +8215,6 @@ 28956 93014AH 28958 93014AK 28961 93014AN -28965 93014AS 28972 79095H 28973 79095J 28974 79095K @@ -8332,8 +8233,6 @@ 28989 88005D 28990 88005E 28991 88005F -28992 88005G -28993 88005H 28994 88005J 28995 88005K 28998 03024C @@ -8344,18 +8243,14 @@ 29003 79101F 29004 79101G 29006 62047F -29011 65108CT -29012 65108CU -29013 65108CV +29008 64038G +29009 65063P 29014 77092K 29019 80034F -29020 81021V -29021 84069AL 29023 89006AP -29024 91003Q 29026 92021S 29029 93014AW -29032 93014AZ +29030 93014AX 29033 94029AEQ 29036 79095N 29039 99057NJ @@ -8379,7 +8274,6 @@ 29061 65027BL 29064 73088H 29065 78100AX -29066 88005L 29067 88005M 29068 91082BP 29069 91082BQ @@ -8389,7 +8283,6 @@ 29073 91082BU 29074 91082BV 29075 91082BW -29076 91082BX 29077 91082BY 29078 91082BZ 29079 06014A @@ -8409,7 +8302,6 @@ 29096 93016AK 29097 93016AL 29098 93016AM -29099 93016AN 29100 93016AP 29101 93016AQ 29102 93016AR @@ -8421,6 +8313,7 @@ 29110 06016D 29126 85108Q 29136 85100C +29137 85100D 29138 66076J 29139 74015F 29155 06018A @@ -8438,20 +8331,16 @@ 29185 68002J 29186 68002K 29187 68002L -29189 68002N 29194 76103G 29195 76103H 29196 76103J 29197 78047C +29198 79095P 29199 80073E 29200 80073F 29201 80073G 29203 81037E -29204 82099F -29208 87080L 29209 87080M -29211 87080P -29212 88033H 29213 99057NL 29215 02024F 29216 02024G @@ -8504,7 +8393,6 @@ 29317 89039Y 29325 89039AG 29328 89039AK -29347 89039BE 29349 06034A 29350 06034B 29355 87020DW @@ -8514,18 +8402,22 @@ 29359 87020EA 29360 87020EB 29362 88033J +29363 00075F 29393 06037A 29398 06038A 29414 89039BS +29416 89039BU +29422 04023D 29441 89039CQ 29444 81021W 29445 84035D 29446 84035E 29447 86010F +29448 86010G 29449 89039CT +29450 89039CU 29479 06041A 29486 06042A -29489 89039DZ 29494 06043A 29495 06043B 29498 06043E @@ -8600,7 +8492,6 @@ 29754 99025AU 29756 99025AW 29758 99025AY -29759 99025AZ 29760 99025BA 29762 99025BC 29763 99025BD @@ -8665,7 +8556,6 @@ 29835 99025ED 29836 99025EE 29837 99025EF -29838 99025EG 29840 99025EJ 29841 99025EK 29842 99025EL @@ -8689,7 +8579,6 @@ 29860 99025FE 29865 99025FK 29866 99025FL -29867 99025FM 29868 99025FN 29869 99025FP 29870 99025FQ @@ -8855,7 +8744,6 @@ 30052 99025PD 30053 99025PE 30054 99025PF -30055 99025PG 30056 99025PH 30057 99025PJ 30058 99025PK @@ -8973,13 +8861,11 @@ 30181 99025UN 30182 99025UP 30183 99025UQ -30189 99025UW 30190 99025UX 30191 99025UY 30192 99025UZ 30193 99025VA 30194 99025VB -30195 99025VC 30196 99025VD 30197 99025VE 30198 99025VF @@ -9003,7 +8889,6 @@ 30222 99025WF 30224 99025WH 30225 99025WJ -30226 99025WK 30227 99025WL 30228 99025WM 30229 99025WN @@ -9058,7 +8943,6 @@ 30289 99025ZA 30290 99025ZB 30291 99025ZC -30292 99025ZD 30293 99025ZE 30294 99025ZF 30295 99025ZG @@ -9170,7 +9054,6 @@ 30418 99025AEH 30419 99025AEJ 30420 99025AEK -30421 99025AEL 30423 99025AEN 30429 99025AEQ 30431 99025AES @@ -9377,7 +9260,6 @@ 30689 99025AQX 30690 99025AQY 30691 99025AQZ -30692 99025ARA 30693 99025ARB 30694 99025ARC 30695 99025ARD @@ -9395,7 +9277,6 @@ 30708 99025ARS 30709 99025ART 30710 99025ARU -30712 99025ARW 30714 99025ARY 30715 99025ARZ 30716 99025ASA @@ -9427,7 +9308,6 @@ 30744 99025ATE 30745 99025ATF 30746 99025ATG -30748 99025ATJ 30749 99025ATK 30750 99025ATL 30751 99025ATM @@ -9475,7 +9355,6 @@ 30815 99025AUZ 30816 99025AVA 30817 99025AVB -30819 99025AVD 30820 99025AVE 30821 99025AVF 30822 99025AVG @@ -9486,7 +9365,6 @@ 30829 99025AVP 30830 99025AVQ 30831 99025AVR -30832 99025AVS 30833 99025AVT 30835 99025AVV 30836 99025AVW @@ -9544,8 +9422,6 @@ 30905 99025AYT 30906 99025AYU 30907 99025AYV -30909 99025AYX -30910 99025AYY 30912 99025AZA 30913 99025AZB 30914 99025AZC @@ -9590,7 +9466,6 @@ 30959 99025BAZ 30960 99025BBA 30961 99025BBB -30963 99025BBD 30964 99025BBE 30966 99025BBG 30967 99025BBH @@ -9606,7 +9481,6 @@ 30980 99025BBW 30981 99025BBX 30982 99025BBY -30983 99025BBZ 30984 99025BCA 30985 99025BCB 30986 99025BCC @@ -9614,7 +9488,6 @@ 30989 99025BCF 30993 99025BCK 30994 99025BCL -30995 99025BCM 30996 99025BCN 30997 99025BCP 30998 99025BCQ @@ -9682,7 +9555,7 @@ 31073 99025BFT 31074 99025BFU 31077 99025BFX -31079 99025BFZ +31078 99025BFY 31080 99025BGA 31081 99025BGB 31082 99025BGC @@ -9730,7 +9603,6 @@ 31137 95021D 31138 79101H 31139 03049E -31140 07014A 31142 66040AQ 31146 99025BGY 31147 99025BGZ @@ -9752,7 +9624,6 @@ 31167 99025BHV 31168 99025BHW 31169 99025BHX -31170 99025BHY 31171 99025BHZ 31172 99025BJA 31173 99025BJB @@ -9817,7 +9688,6 @@ 31241 99025BLX 31242 99025BLY 31243 99025BLZ -31244 99025BMA 31246 99025BMC 31247 99025BMD 31248 99025BME @@ -9839,15 +9709,14 @@ 31265 99025BMX 31266 99025BMY 31267 99025BMZ -31268 99025BNA 31269 99025BNB 31270 99025BNC 31272 99025BNE 31273 99025BNF 31274 99025BNG 31275 99025BNH -31276 99025BNJ 31279 99025BNM +31280 99025BNN 31281 99025BNP 31282 99025BNQ 31284 99025BNS @@ -9896,6 +9765,7 @@ 31331 99025BQK 31332 99025BQL 31333 99025BQM +31334 99025BQN 31335 99025BQP 31336 99025BQQ 31337 99025BQR @@ -9928,6 +9798,7 @@ 31370 99025BSA 31371 99025BSB 31372 99025BSC +31374 99025BSE 31376 99025BSG 31378 99025BSJ 31379 99025BSK @@ -9957,7 +9828,6 @@ 31410 99025BSW 31411 99025BSX 31412 99025BSY -31413 99025BSZ 31414 99025BTA 31415 99025BTB 31416 99025BTC @@ -9965,7 +9835,6 @@ 31418 99025BTE 31419 99025BTF 31420 99025BTG -31421 99025BTH 31422 99025BTJ 31423 99025BTK 31424 99025BTL @@ -10050,7 +9919,6 @@ 31516 99025BXF 31519 99025BXJ 31520 99025BXK -31521 99025BXL 31522 99025BXM 31523 99025BXN 31524 99025BXP @@ -10073,7 +9941,6 @@ 31546 99025BYM 31547 99025BYN 31548 99025BYP -31549 99025BYQ 31550 99025BYR 31551 99025BYS 31552 99025BYT @@ -10101,27 +9968,22 @@ 31577 07021A 31579 66040AR 31581 99057PC -31582 99057PD 31583 99057PE -31586 99057PH 31587 99057PJ 31591 99057PN 31592 99057PP 31594 99057PR 31598 07023A +31601 07025A 31608 99025BZN 31609 99025BZP 31611 99025BZR 31612 99025BZS 31613 99025BZT -31614 99025BZU 31615 99025BZV -31616 99025BZW -31617 99025BZX 31618 99025BZY 31619 99025BZZ 31620 99025CAA -31621 99025CAB 31622 99025CAC 31623 99025CAD 31624 99025CAE @@ -10134,6 +9996,7 @@ 31631 99025CAM 31632 99025CAN 31633 99025CAP +31634 99025CAQ 31635 99025CAR 31636 99025CAS 31637 99025CAT @@ -10148,7 +10011,6 @@ 31649 99025CBF 31650 99025CBG 31652 99025CBJ -31653 99025CBK 31655 99025CBM 31656 99025CBN 31657 99025CBP @@ -10170,6 +10032,7 @@ 31678 99025CCL 31679 99025CCM 31680 99025CCN +31682 99025CCQ 31683 99025CCR 31685 99025CCT 31686 99025CCU @@ -10178,7 +10041,6 @@ 31690 99025CCY 31691 99025CCZ 31692 99025CDA -31694 99025CDC 31695 99025CDD 31696 99025CDE 31698 07026A @@ -10189,15 +10051,14 @@ 31709 99025CDG 31710 99025CDH 31711 99025CDJ -31712 99025CDK 31714 99025CDM 31715 99025CDN 31716 99025CDP 31717 99025CDQ +31719 99025CDS 31720 99025CDT 31721 99025CDU 31722 99025CDV -31723 99025CDW 31724 99025CDX 31725 99025CDY 31726 99025CDZ @@ -10219,8 +10080,6 @@ 31744 99025CET 31745 99025CEU 31746 99025CEV -31750 99025CEZ -31751 99025CFA 31752 99025CFB 31753 99025CFC 31754 99025CFD @@ -10241,7 +10100,6 @@ 31771 99025CFW 31772 99025CFX 31773 99025CFY -31774 99025CFZ 31776 99025CGB 31779 99025CGE 31780 99025CGF @@ -10249,11 +10107,9 @@ 31782 99025CGH 31783 99025CGJ 31784 99025CGK -31786 99025CGM 31787 99025CGN 31788 99025CGP 31789 07028A -31790 07028B 31791 07028C 31792 07029A 31793 07029B @@ -10264,18 +10120,15 @@ 31799 07029F 31800 07031A 31802 99025CGQ -31803 99025CGR 31804 99025CGS 31805 99025CGT 31806 99025CGU 31807 99025CGV 31808 99025CGW -31809 99025CGX 31810 99025CGY 31811 99025CGZ 31812 99025CHA 31813 99025CHB -31814 99025CHC 31815 99025CHD 31816 99025CHE 31817 99025CHF @@ -10290,13 +10143,12 @@ 31829 99025CHT 31830 99025CHU 31831 99025CHV -31832 99025CHW 31833 99025CHX 31834 99025CHY 31836 99025CJA 31837 99025CJB 31838 99025CJC -31839 99025CJD +31840 99025CJE 31841 99025CJF 31842 99025CJG 31843 99025CJH @@ -10308,19 +10160,16 @@ 31851 99025CJR 31852 99025CJS 31854 99025CJU -31855 99025CJV +31856 99025CJW 31857 99025CJX -31858 99025CJY 31861 99025CKB 31862 07032A 31863 07032B 31864 07032C -31867 99057PS 31870 99057PV 31871 99057PW 31872 99057PX 31873 99057PY -31875 99057QA 31876 99057QB 31877 99025CKC 31878 99025CKD @@ -10345,7 +10194,6 @@ 31904 99025CLF 31905 99025CLG 31906 99025CLH -31908 99025CLK 31910 99025CLM 31911 99025CLN 31912 99025CLP @@ -10366,7 +10214,6 @@ 31934 99025CMK 31935 99025CML 31938 99025CMP -31939 99025CMQ 31942 99025CMT 31944 99025CMV 31945 99025CMW @@ -10374,7 +10221,6 @@ 31947 99025CMY 31948 99025CMZ 31949 99025CNA -31950 99025CNB 31952 99025CND 31953 99025CNE 31954 99025CNF @@ -10423,7 +10269,6 @@ 32021 07036D 32025 06006C 32026 06006D -32027 06006E 32028 06006F 32029 06006G 32030 06006H @@ -10435,7 +10280,6 @@ 32036 06006P 32037 06006Q 32038 06006R -32040 06006T 32041 06006U 32042 06006V 32043 06006W @@ -10450,15 +10294,11 @@ 32066 99057QE 32071 99057QK 32072 99057QL -32074 99057QN 32075 99057QP -32076 99057QQ -32078 99057QS 32079 99057QT 32080 99057QU 32081 99057QV 32082 99057QW -32083 99057QX 32084 99057QY 32085 99057QZ 32086 99025CQC @@ -10471,12 +10311,12 @@ 32094 99025CQL 32095 99025CQM 32096 99025CQN -32097 99025CQP 32098 99025CQQ 32099 99025CQR 32101 99025CQT 32102 99025CQU 32103 99025CQV +32104 99025CQW 32105 99025CQX 32106 99025CQY 32109 99025CRB @@ -10486,7 +10326,6 @@ 32113 99025CRF 32114 99025CRG 32115 99025CRH -32116 99025CRJ 32117 99025CRK 32118 99025CRL 32119 99025CRM @@ -10505,8 +10344,8 @@ 32136 99025CSE 32137 99025CSF 32138 99025CSG -32139 99025CSH 32141 99025CSK +32142 99025CSL 32143 99025CSM 32146 99025CSQ 32148 99025CSS @@ -10533,6 +10372,7 @@ 32172 99025CTS 32173 99025CTT 32174 99025CTU +32175 99025CTV 32176 99025CTW 32177 99025CTX 32178 99025CTY @@ -10550,8 +10390,10 @@ 32190 99025CUL 32191 99025CUM 32192 99025CUN +32193 99025CUP +32194 99025CUQ +32195 99025CUR 32197 99025CUT -32198 99025CUU 32199 99025CUV 32200 99025CUW 32202 99025CUY @@ -10559,7 +10401,6 @@ 32204 99025CVA 32205 99025CVB 32207 99025CVD -32208 99025CVE 32209 99025CVF 32210 99025CVG 32212 99025CVJ @@ -10580,7 +10421,6 @@ 32233 99025CWF 32234 99025CWG 32235 99025CWH -32237 99025CWK 32238 99025CWL 32241 99025CWP 32242 99025CWQ @@ -10603,6 +10443,7 @@ 32276 07052B 32277 07052C 32278 07052D +32280 07052F 32281 07052G 32283 07053A 32287 07054A @@ -10619,12 +10460,12 @@ 32305 99025CWW 32308 99025CWZ 32311 99025CXC -32312 99025CXD 32313 99025CXE 32314 99025CXF 32315 99025CXG 32316 99025CXH 32317 99025CXJ +32318 99025CXK 32319 99025CXL 32320 99025CXM 32322 99025CXP @@ -10661,6 +10502,7 @@ 32364 99025CZH 32366 99025CZK 32368 99025CZM +32370 99025CZP 32371 99025CZQ 32372 99025CZR 32373 07058A @@ -10686,7 +10528,6 @@ 32406 99025CZS 32407 99025CZT 32409 99025CZV -32410 99025CZW 32411 99025CZX 32414 99025DAA 32416 99025DAC @@ -10696,20 +10537,20 @@ 32425 99025DAM 32427 99025DAP 32429 99025DAR +32430 99025DAS 32431 99025DAT 32432 99025DAU -32433 99025DAV 32435 99025DAX -32440 99025DBC +32438 99025DBA 32443 99025DBF 32444 99025DBG 32445 99025DBH 32448 99025DBL -32450 99025DBN 32451 99025DBP -32452 99025DBQ 32453 99025DBR +32454 99025DBS 32455 99025DBT +32456 99025DBU 32457 99025DBV 32458 99025DBW 32459 99025DBX @@ -10718,6 +10559,7 @@ 32464 99025DCC 32467 99057RB 32468 99057RC +32470 99057RE 32473 07042C 32474 07042D 32475 66040AS @@ -10729,17 +10571,12 @@ 32483 79009D 32487 08006A 32488 08006B -32492 06006AF 32493 06006AG 32494 06006AH -32495 06006AJ -32496 06006AK -32499 06006AN 32500 08007A 32501 08007B 32705 77093E 32706 08010A -32707 08010B 32708 08011A 32709 08011B 32711 08012A @@ -10748,7 +10585,6 @@ 32730 08013B 32750 08014A 32763 08016A -32765 08017A 32766 08017B 32767 08018A 32768 08018B @@ -10777,7 +10613,6 @@ 32957 08025E 32958 08026A 32959 08026B -32974 78083S 32977 79058S 32982 79077H 32984 79101K @@ -10791,7 +10626,6 @@ 32993 79101U 32998 91075P 32999 92021T -33000 78012D 33051 08028A 33053 08029A 33055 08030A @@ -10904,11 +10738,11 @@ 33505 09003B 33508 67048G 33509 68081K -33510 68081L 33511 68081M 33512 68081N 33513 68081P 33514 71024E +33519 77092L 33520 81031M 33521 81049D 33524 81071Z @@ -10926,16 +10760,12 @@ 33544 88109T 33545 88109U 33546 88109V -33547 88109W 33548 89006AQ 33549 89006AR 33551 89006AT -33554 89006AW -33555 89006AX 33556 89006AY 33557 89006AZ 33558 89006BA -33559 89006BB 33560 89006BC 33561 89006BD 33562 89006BE @@ -10944,12 +10774,11 @@ 33567 89006BK 33568 89006BL 33569 89006BM -33570 89006BN -33572 89006BQ 33574 89006BS 33578 89006BW 33579 89006BX 33580 89006BY +33581 89006BZ 33582 89006CA 33584 89027E 33588 91003R @@ -10964,8 +10793,6 @@ 33601 99025DCJ 33603 99025DCL 33605 99025DCN -33606 99025DCP -33612 99025DCV 33613 99025DCW 33614 99025DCX 33617 99025DDA @@ -10977,6 +10804,7 @@ 33625 99025DDJ 33628 99025DDM 33629 99025DDN +33630 99025DDP 33632 99025DDR 33633 99025DDS 33636 99025DDV @@ -11001,7 +10829,6 @@ 33663 99025DEY 33664 99025DEZ 33665 99025DFA -33666 99025DFB 33668 99025DFD 33673 99025DFJ 33674 99025DFK @@ -11044,7 +10871,6 @@ 33723 99025DHL 33724 99025DHM 33726 99025DHP -33727 99025DHQ 33728 99025DHR 33730 99025DHT 33731 99025DHU @@ -11058,8 +10884,6 @@ 33740 99025DJD 33741 99025DJE 33745 99025DJJ -33746 99025DJK -33747 99025DJL 33748 99025DJM 33749 09008A 33750 09008B @@ -11089,7 +10913,6 @@ 33784 93036Z 33785 93036AA 33789 93036AE -33790 93036AF 33791 93036AG 33792 93036AH 33793 93036AJ @@ -11135,7 +10958,6 @@ 33853 97051W 33854 97051X 33855 97051Y -33856 97051Z 33858 97051AB 33859 97051AC 33860 97051AD @@ -11155,8 +10977,6 @@ 33878 97051AX 33879 97051AY 33881 97051BA -33882 97051BB -33883 97051BC 33884 97051BD 33886 97051BF 33887 97051BG @@ -11176,7 +10996,6 @@ 33907 93036DL 33908 93036DM 33909 93036DN -33910 93036DP 33911 93036DQ 33912 93036DR 33913 93036DS @@ -11206,7 +11025,6 @@ 33942 93036EX 33944 93036EZ 33945 93036FA -33947 93036FC 33948 93036FD 33950 97051BK 33951 97051BL @@ -11226,7 +11044,6 @@ 33971 93036FE 33972 93036FF 33973 93036FG -33974 93036FH 33975 93036FJ 33980 93036FP 33981 93036FQ @@ -11297,7 +11114,6 @@ 34069 93036KG 34070 93036KH 34071 97051CE -34073 97051CG 34075 97051CJ 34076 97051CK 34077 97051CL @@ -11313,7 +11129,6 @@ 34097 97051DG 34098 97051DH 34099 97051DJ -34100 97051DK 34101 97051DL 34102 97051DM 34103 97051DN @@ -11352,15 +11167,11 @@ 34148 97051EB 34150 97051ED 34155 97051EJ -34156 97051EK -34157 97051EL 34159 97051EN 34160 97051EP -34164 65108CX -34167 65108DA +34165 65108CY 34169 65108DC 34170 65108DD -34173 65108DG 34176 87011G 34178 92021U 34179 92021V @@ -11368,11 +11179,11 @@ 34181 92041K 34184 92072Q 34186 92072S -34188 92072U 34189 92072V 34190 92072W 34191 92072X 34192 92072Y +34193 92072Z 34194 92072AA 34195 92072AB 34196 92072AC @@ -11387,31 +11198,27 @@ 34213 99057RQ 34214 99057RR 34215 99057RS -34217 99057RU 34221 01049PR 34223 02024H 34224 06006AP -34225 06006AQ 34226 06006AR 34229 06006AU 34230 06006AV -34231 06006AW -34234 06006AZ -34235 06006BA -34236 06006BB +34233 06006AY +34238 06006BD 34239 06006BE 34241 06006BG -34242 06006BH -34243 06006BJ 34246 06006BM +34247 06006BN 34248 06006BP 34250 06006BR 34252 06006BT 34254 06006BV 34256 06006BX 34257 06006BY +34258 06006BZ 34259 06006CA -34261 06006CC +34260 06006CB 34263 06006CE 34264 09010A 34265 09010B @@ -11451,9 +11258,7 @@ 34331 93036PD 34333 93036PF 34337 93036PK -34339 93036PM 34340 93036PN -34341 93036PP 34342 93036PQ 34343 93036PR 34344 93036PS @@ -11461,18 +11266,13 @@ 34350 97051ES 34351 97051ET 34354 97051EW -34356 97051EY 34358 97051FA 34359 97051FB -34360 97051FC 34361 97051FD 34363 97051FF -34364 97051FG 34366 97051FJ 34367 97051FK 34368 97051FL -34373 97051FR -34374 97051FS 34375 97051FT 34376 97051FU 34378 97051FW @@ -11498,7 +11298,6 @@ 34420 93036RM 34421 93036RN 34422 93036RP -34423 93036RQ 34424 93036RR 34425 93036RS 34426 93036RT @@ -11510,8 +11309,6 @@ 34432 93036RZ 34433 93036SA 34434 93036SB -34435 93036SC -34436 93036SD 34439 93036SG 34442 93036SK 34445 93036SN @@ -11544,7 +11341,6 @@ 34488 97051GD 34489 97051GE 34490 97051GF -34491 97051GG 34492 97051GH 34493 97051GJ 34496 97051GM @@ -11552,7 +11348,6 @@ 34503 97051GU 34508 97051GZ 34511 97051HC -34515 97051HG 34517 97051HJ 34518 97051HK 34520 97051HM @@ -11586,7 +11381,6 @@ 34570 93036VG 34571 93036VH 34572 93036VJ -34573 93036VK 34574 93036VL 34576 93036VN 34577 93036VP @@ -11624,13 +11418,11 @@ 34652 97051KD 34657 97051KJ 34661 09014A -34662 09014B 34667 91025L 34671 93036XR 34672 93036XS 34674 93036XU 34675 93036XV -34676 93036XW 34678 93036XY 34679 93036XZ 34680 93036YA @@ -11659,12 +11451,10 @@ 34724 93036YS 34725 93036YT 34726 93036YU -34727 93036YV 34728 93036YW 34729 93036YX 34730 93036YY 34732 93036ZA -34733 93036ZB 34735 93036ZD 34736 93036ZE 34737 93036ZF @@ -11704,22 +11494,17 @@ 34791 93036AAU 34792 93036AAV 34793 93036AAW -34794 93036AAX -34795 93036AAY 34796 93036AAZ 34797 93036ABA 34807 09019A -34809 09019C 34810 09020A 34811 09020B -34812 93036ABB 34813 93036ABC 34814 93036ABD 34815 93036ABE 34816 93036ABF 34818 93036ABH 34819 93036ABJ -34820 93036ABK 34822 93036ABM 34823 93036ABN 34825 97051LZ @@ -11733,14 +11518,12 @@ 34852 93036ABW 34853 93036ABX 34854 93036ABY -34855 93036ABZ 34856 93036ACA 34857 93036ACB 34858 93036ACC 34859 93036ACD 34860 93036ACE 34861 93036ACF -34863 93036ACH 34868 97051MN 34869 97051MP 34870 97051MQ @@ -11775,7 +11558,6 @@ 34916 93036ADM 34918 93036ADP 34920 93036ADR -34921 93036ADS 34922 93036ADT 34924 93036ADV 34926 97051NF @@ -11783,9 +11565,7 @@ 34931 97051NL 34941 09027A 34942 09027B -34943 09027C 34944 93036ADX -34945 93036ADY 34946 93036ADZ 34948 93036AEB 34949 93036AEC @@ -11809,7 +11589,6 @@ 34977 93036AFG 34978 93036AFH 34979 93036AFJ -34982 97051NN 34983 97051NP 34985 97051NR 34986 97051NS @@ -11817,7 +11596,6 @@ 34989 93036AFM 34990 93036AFN 34991 93036AFP -34992 93036AFQ 34993 93036AFR 34995 93036AFT 34997 93036AFV @@ -11833,7 +11611,6 @@ 35018 93036AGF 35019 93036AGG 35020 93036AGH -35021 93036AGJ 35022 93036AGK 35023 93036AGL 35025 93036AGN @@ -11842,7 +11619,6 @@ 35031 93036AGU 35032 93036AGV 35033 93036AGW -35034 93036AGX 35037 93036AHA 35038 93036AHB 35039 93036AHC @@ -11850,7 +11626,6 @@ 35043 93036AHG 35045 93036AHJ 35047 93036AHL -35048 97051NU 35049 97051NV 35050 97051NW 35051 97051NX @@ -11878,10 +11653,11 @@ 35093 99025DJV 35094 99025DJW 35095 99025DJX -35097 99025DJZ +35099 99025DKB 35102 99025DKE 35104 99025DKG 35105 99025DKH +35108 99025DKL 35109 99025DKM 35111 99025DKP 35114 99025DKS @@ -11889,12 +11665,12 @@ 35117 99025DKV 35118 99025DKW 35119 99025DKX +35120 99025DKY 35121 99025DKZ 35122 99025DLA 35123 99025DLB 35125 99025DLD 35128 99025DLG -35129 99025DLH 35130 99025DLJ 35134 99025DLN 35135 99025DLP @@ -11917,7 +11693,7 @@ 35159 99025DMP 35160 99025DMQ 35161 99025DMR -35162 99025DMS +35163 99025DMT 35164 99025DMU 35165 99025DMV 35167 99025DMX @@ -11927,6 +11703,7 @@ 35171 99025DNB 35173 99025DND 35174 99025DNE +35175 99025DNF 35176 99025DNG 35177 99025DNH 35180 99025DNL @@ -11959,20 +11736,18 @@ 35221 99025DQD 35222 99025DQE 35223 99025DQF -35224 99025DQG 35225 99025DQH 35226 99025DQJ 35227 99025DQK -35228 99025DQL 35229 99025DQM 35230 99025DQN 35231 99025DQP 35233 99025DQR 35234 99025DQS -35235 99025DQT 35236 99025DQU 35237 99025DQV 35239 99025DQX +35240 99025DQY 35241 99025DQZ 35242 99025DRA 35243 99025DRB @@ -11981,7 +11756,6 @@ 35248 07055C 35249 07055D 35251 07055F -35252 07055G 35253 79095U 35254 79095V 35255 79095W @@ -12031,7 +11805,6 @@ 35331 99057SK 35332 99057SL 35334 99057SN -35336 01049PU 35339 01049PX 35341 01049PZ 35342 01049QA @@ -12053,7 +11826,6 @@ 35364 88005Y 35365 88005Z 35366 88005AA -35367 88005AB 35368 88005AC 35369 88005AD 35370 88005AE @@ -12070,15 +11842,14 @@ 35382 03049J 35383 03049K 35384 99057SP -35386 99057SR 35387 99057SS -35388 08064C 35391 95021E 35392 95021F 35393 95021G 35397 07042E 35398 07042F 35401 93041E +35403 93041G 35404 02024J 35405 02024K 35406 82116F @@ -12094,7 +11865,8 @@ 35420 01056M 35422 75124D 35423 77002F -35428 04045D +35424 81043M +35426 80073J 35429 93036AJX 35430 93036AJY 35431 93036AJZ @@ -12115,10 +11887,9 @@ 35453 93036AKX 35455 93036AKZ 35456 93036ALA -35459 93036ALD +35460 93036ALE 35461 93036ALF 35462 93036ALG -35463 93036ALH 35464 93036ALJ 35465 93036ALK 35467 93036ALM @@ -12132,7 +11903,6 @@ 35479 97051PZ 35480 97051QA 35483 97051QD -35484 97051QE 35487 97051QH 35488 97051QJ 35491 09033A @@ -12159,7 +11929,6 @@ 35523 87068CS 35524 87068CT 35525 87068CU -35526 88074P 35528 88074R 35530 91068AX 35532 91082CP @@ -12168,20 +11937,23 @@ 35535 75052LH 35536 01049QL 35538 06006CF -35543 06006CL +35540 06006CH +35542 06006CK +35544 06006CM 35545 06006CN +35547 07003F 35548 07003G -35549 07003H 35550 07003J 35551 07003K 35552 07003L -35553 07003M 35554 07003N 35555 07003P 35556 07003Q 35557 07003R 35558 07003S +35560 07003U 35563 07003X +35564 07003Y 35567 07042G 35568 07047D 35571 08053G @@ -12197,13 +11969,14 @@ 35586 93036AMD 35587 93036AME 35589 93036AMG +35591 93036AMJ 35592 93036AMK 35593 93036AML 35594 93036AMM -35596 93036AMP 35598 93036AMR 35600 93036AMT 35602 93036AMV +35603 93036AMW 35604 93036AMX 35605 93036AMY 35606 93036AMZ @@ -12217,8 +11990,8 @@ 35616 97051QN 35617 97051QP 35618 97051QQ +35620 97051QS 35622 97051QU -35623 97051QV 35624 97051QW 35625 97051QX 35627 97051QZ @@ -12235,10 +12008,10 @@ 35646 93036ANM 35647 93036ANN 35650 93036ANR -35652 93036ANT 35653 93036ANU 35655 93036ANW 35657 93036ANY +35658 93036ANZ 35660 93036APB 35661 93036APC 35663 93036APE @@ -12275,8 +12048,8 @@ 35706 93036AQB 35707 93036AQC 35709 93036AQE -35710 93036AQF 35711 93036AQG +35713 93036AQJ 35716 93036AQM 35717 93036AQN 35718 93036AQP @@ -12288,18 +12061,15 @@ 35729 93036ARA 35730 93036ARB 35732 97051RK -35733 97051RL 35735 97051RN 35737 97051RQ 35739 97051RS -35741 97051RU 35742 97051RV 35744 97051RX 35745 97051RY 35747 97051SA 35748 97051SB 35749 97051SC -35750 97051SD 35752 09043A 35754 09043C 35755 09044A @@ -12334,13 +12104,11 @@ 35799 97051SH 35800 97051SJ 35802 97051SL -35803 97051SM 35805 97051SP 35806 97051SQ 35809 97051ST 35812 09046A 35815 09047A -35821 93036AST 35822 93036ASU 35824 93036ASW 35825 93036ASX @@ -12382,7 +12150,6 @@ 35885 93036ATY 35887 93036AUA 35888 93036AUB -35889 93036AUC 35890 93036AUD 35891 93036AUE 35893 93036AUG @@ -12397,12 +12164,9 @@ 35905 93036AUU 35906 93036AUV 35911 97051TU -35914 97051TX 35915 97051TY -35916 97051TZ 35917 97051UA 35918 97051UB -35920 97051UD 35921 97051UE 35922 97051UF 35925 97051UJ @@ -12426,7 +12190,6 @@ 35954 09057D 35955 93061K 35956 93036AUW -35957 93036AUX 35958 93036AUY 35959 93036AUZ 35960 93036AVA @@ -12438,6 +12201,7 @@ 35970 93036AVL 35974 93036AVQ 35975 93036AVR +35976 93036AVS 35979 93036AVV 35980 93036AVW 35982 93036AVY @@ -12447,10 +12211,9 @@ 35987 93036AWD 35989 93036AWF 35991 93036AWH -35992 93036AWJ 35993 93036AWK -35996 93036AWN 35997 93036AWP +36000 93036AWS 36002 93036AWU 36004 93036AWW 36005 93036AWX @@ -12462,8 +12225,6 @@ 36021 97051VC 36023 97051VE 36025 97051VG -36026 97051VH -36027 97051VJ 36028 97051VK 36030 96051M 36031 96051N @@ -12481,7 +12242,6 @@ 36048 93036AXK 36049 93036AXL 36050 93036AXM -36051 93036AXN 36052 93036AXP 36053 93036AXQ 36055 93036AXS @@ -12545,6 +12305,7 @@ 36140 99025DRM 36141 99025DRN 36142 99025DRP +36147 99025DRU 36148 99025DRV 36149 99025DRW 36152 99025DRZ @@ -12558,9 +12319,8 @@ 36163 99025DSL 36166 99025DSP 36167 99025DSQ -36169 99025DSS +36170 99025DST 36172 99025DSV -36173 99025DSW 36177 99025DTA 36180 99025DTD 36181 99025DTE @@ -12575,15 +12335,16 @@ 36198 99025DTX 36199 99025DTY 36200 99025DTZ +36201 99025DUA 36203 99025DUC 36204 99025DUD 36206 99025DUF 36207 99025DUG 36211 99025DUL -36214 99025DUP 36216 99025DUR 36218 99025DUT 36222 99025DUX +36223 99025DUY 36224 99025DUZ 36225 99025DVA 36226 99025DVB @@ -12616,6 +12377,7 @@ 36269 99025DWW 36270 99025DWX 36271 99025DWY +36272 99025DWZ 36273 99025DXA 36274 99025DXB 36278 99025DXF @@ -12623,7 +12385,6 @@ 36280 99025DXH 36281 99025DXJ 36283 99025DXL -36284 99025DXM 36287 10001A 36288 10001B 36289 08053K @@ -12636,6 +12397,7 @@ 36301 69009D 36302 75087C 36304 99057SV +36307 80073K 36308 80073L 36309 80073M 36318 82025C @@ -12643,16 +12405,17 @@ 36320 82025E 36321 82025F 36322 82025G -36323 82025H -36324 82025J 36325 82025K 36326 82025L +36327 82025M 36328 82025N 36329 82025P 36330 82025Q 36331 82025R +36333 82025T 36334 82025U 36335 82025V +36336 82025W 36337 82025X 36338 82025Y 36339 82025Z @@ -12661,7 +12424,6 @@ 36343 82025AD 36344 82025AE 36345 82025AF -36346 82025AG 36347 82025AH 36348 82025AJ 36349 82025AK @@ -12669,6 +12431,7 @@ 36351 82025AM 36354 82025AQ 36355 82025AR +36357 82025AT 36358 10002A 36359 10002B 36360 10002C @@ -12676,12 +12439,10 @@ 36365 06015J 36366 93036AYT 36367 93036AYU -36368 93036AYV 36369 93036AYW 36371 93036AYY 36372 93036AYZ 36375 93036AZC -36376 93036AZD 36378 93036AZF 36379 93036AZG 36380 93036AZH @@ -12699,6 +12460,7 @@ 36401 10007B 36402 10007C 36403 10007D +36406 10007G 36407 10007H 36411 10008A 36412 10008B @@ -12724,21 +12486,19 @@ 36451 93036BBA 36453 93036BBC 36454 93036BBD +36459 93036BBJ 36461 93036BBL 36462 93036BBM 36467 93036BBS 36470 93036BBV -36471 93036BBW -36474 93036BBZ 36476 93036BCB 36479 93036BCE -36480 93036BCF -36482 97051VY 36483 97051VZ 36485 97051WB 36486 97051WC 36487 97051WD 36488 97051WE +36489 97051WF 36490 97051WG 36491 97051WH 36492 97051WJ @@ -12759,6 +12519,7 @@ 36528 93036BCN 36529 93036BCP 36530 93036BCQ +36531 93036BCR 36532 93036BCS 36533 93036BCT 36539 93036BCZ @@ -12766,16 +12527,15 @@ 36541 93036BDB 36542 93036BDC 36543 93036BDD -36544 93036BDE 36548 93036BDJ 36550 93036BDL 36551 93036BDM 36553 93036BDP -36556 93036BDS +36554 93036BDQ 36557 93036BDT 36558 93036BDU 36559 93036BDV -36561 93036BDX +36563 97051WS 36566 97051WV 36568 97051WX 36571 06002AA @@ -12803,7 +12563,6 @@ 36605 10030A 36606 10030B 36608 10031A -36613 93036BEB 36614 93036BEC 36615 93036BED 36617 93036BEF @@ -12814,7 +12573,6 @@ 36632 93036BEW 36642 97051XC 36644 99025DXN -36645 99025DXP 36652 99025DXW 36654 99025DXY 36656 99025DYA @@ -12840,6 +12598,7 @@ 36686 99025DZG 36687 99025DZH 36689 99025DZK +36691 99025DZM 36692 99025DZN 36694 99025DZQ 36695 99025DZR @@ -12853,6 +12612,7 @@ 36706 99025EAC 36707 99025EAD 36710 99025EAG +36711 99025EAH 36713 99025EAK 36714 99025EAL 36716 99025EAN @@ -12868,7 +12628,6 @@ 36735 99025EBH 36736 99025EBJ 36738 99025EBL -36740 99025EBN 36741 99025EBP 36743 99025EBR 36744 10032A @@ -12892,68 +12651,49 @@ 36832 10037C 36833 10037D 36834 10038A -36835 10038B 36858 82025AV -36859 82025AW +36860 82025AX 36861 82025AY 36862 82025AZ 36863 82025BA 36864 82025BB 36865 82025BC -36867 82025BE +36866 82025BD 36868 10039A 36875 81043N -36876 81043P -36877 81043Q 36878 81043R 36879 81043S 36881 05011C -36882 89001V 36886 79095AB -36888 79095AD -36890 79095AF 36891 79095AG 36894 79095AK -36897 79095AN 36898 79095AP 36899 79095AQ 36901 79095AS 36902 79095AT 36903 79095AU -36904 79095AV -36905 79095AW 36906 79095AX -36907 79095AY 36908 79095AZ 36911 79095BC -36912 79095BD 36913 79095BE 36914 79095BF 36915 79095BG 36916 79095BH 36918 79095BK -36919 79095BL -36920 79095BM 36921 79095BN 36922 79095BP 36923 79095BQ -36930 71120T 36935 87011P 36936 87011Q 36937 87011R 36939 87011T 36940 87011U -36941 87011V 36942 87011W 36943 87011X 36944 87011Y 36946 87011AA 36947 87011AB -36948 87011AC -36949 87011AD -36951 87011AF 36959 06006CR -36963 06006CV 36964 06006CW 36965 07003AB 36966 07003AC @@ -12961,48 +12701,40 @@ 36968 07003AE 36970 07003AG 36971 07003AH -36972 81043T -36981 65108DM -36982 65108DN +36979 65108DK 36983 65108DP 36984 65108DQ 36985 10040A 36986 99025EBS +36987 99025EBT 36989 99025EBV 36990 99025EBW 36991 99025EBX 36994 99025ECA 36995 99025ECB 36996 99025ECC -36997 99025ECD +36999 99025ECF 37001 99025ECH 37003 99025ECK 37006 99025ECN 37007 99025ECP -37008 99025ECQ -37009 99025ECR 37011 99025ECT 37012 99025ECU -37013 99025ECV 37014 99025ECW -37015 99025ECX -37016 99025ECY 37017 99025ECZ 37018 99025EDA 37019 99025EDB 37021 99025EDD 37022 99025EDE -37024 99025EDG 37026 99025EDJ -37030 99025EDN +37028 99025EDL +37032 99025EDQ 37035 99025EDT 37036 99025EDU 37040 99025EDY 37041 99025EDZ 37042 99025EEA -37043 99025EEB 37044 99025EEC -37045 99025EED 37046 99025EEE 37047 99025EEF 37048 99025EEG @@ -13014,15 +12746,16 @@ 37056 99025EEQ 37057 99025EER 37058 99025EES -37059 99025EET 37060 99025EEU -37061 99025EEV 37064 99025EEY 37065 99025EEZ 37067 99025EFB +37068 99025EFC 37069 99025EFD 37077 99025EFM 37079 99025EFP +37080 99025EFQ +37081 99025EFR 37084 93036BFE 37085 93036BFF 37086 93036BFG @@ -13034,18 +12767,16 @@ 37094 93036BFQ 37095 93036BFR 37101 93036BFX -37102 93036BFY 37103 93036BFZ 37105 93036BGB 37108 93036BGE 37110 93036BGG 37111 93036BGH +37114 93036BGL 37118 93036BGQ 37120 93036BGS -37122 93036BGU 37123 93036BGV 37124 93036BGW -37125 93036BGX 37126 93036BGY 37127 93036BGZ 37131 93036BHD @@ -13060,7 +12791,6 @@ 37144 10041H 37148 74089FH 37150 10042A -37151 10042B 37152 10043A 37153 10043B 37154 10043C @@ -13099,7 +12829,6 @@ 37206 10056A 37207 10056B 37208 10056C -37209 10056D 37210 10057A 37212 10058A 37213 10058B @@ -13136,7 +12865,6 @@ 37275 93036BHP 37277 93036BHR 37278 93036BHS -37280 93036BHU 37282 93036BHW 37284 93036BHY 37287 93036BJB @@ -13148,12 +12876,14 @@ 37296 93036BJL 37297 93036BJM 37299 93036BJP -37302 93036BJS 37303 93036BJT +37305 93036BJV 37306 93036BJW +37307 93036BJX 37309 93036BJZ 37312 93036BKC 37314 93036BKE +37315 93036BKF 37316 93036BKG 37320 93036BKL 37321 93036BKM @@ -13169,7 +12899,6 @@ 37335 93036BLB 37336 93036BLC 37337 93036BLD -37338 93036BLE 37340 93036BLG 37341 93036BLH 37342 93036BLJ @@ -13201,16 +12930,17 @@ 37398 11018A 37399 11018B 37400 99025EFS +37401 99025EFT 37402 99025EFU 37403 99025EFV 37404 99025EFW 37405 99025EFX 37406 99025EFY 37408 99025EGA +37409 99025EGB 37410 99025EGC -37412 99025EGE +37411 99025EGD 37413 99025EGF -37414 99025EGG 37415 99025EGH 37416 99025EGJ 37417 99025EGK @@ -13220,9 +12950,7 @@ 37422 99025EGQ 37423 99025EGR 37424 99025EGS -37425 99025EGT 37427 99025EGV -37428 99025EGW 37429 99025EGX 37430 99025EGY 37432 99025EHA @@ -13232,34 +12960,25 @@ 37438 99025EHG 37439 99025EHH 37440 99025EHJ -37441 99025EHK 37442 99025EHL -37443 99025EHM 37444 99025EHN 37445 99025EHP 37447 99025EHR 37448 99025EHS 37449 99025EHT -37450 99025EHU -37452 99025EHW 37453 99025EHX 37454 99025EHY 37455 99025EHZ 37456 99025EJA 37457 99025EJB 37460 99025EJE -37461 99025EJF 37462 99025EJG 37463 99025EJH 37464 99025EJJ 37465 99025EJK -37466 99025EJL -37467 99025EJM -37468 99025EJN 37469 99025EJP 37470 99025EJQ 37471 99025EJR -37472 99025EJS 37473 99025EJT 37474 99025EJU 37475 99025EJV @@ -13287,7 +13006,6 @@ 37506 93036BMK 37508 93036BMM 37510 93036BMP -37511 93036BMQ 37515 93036BMU 37518 93036BMX 37520 93036BMZ @@ -13334,14 +13052,13 @@ 37585 91082CR 37586 04025L 37588 99025EKG -37590 99025EKJ +37589 99025EKH 37591 99025EKK 37593 99025EKM 37594 99025EKN 37596 99025EKQ 37597 99025EKR 37598 99025EKS -37600 99025EKU 37602 11021A 37603 11021B 37604 11021C @@ -13353,6 +13070,7 @@ 37612 90081CY 37613 90081CZ 37614 90081DA +37615 90081DB 37616 90081DC 37617 90081DD 37619 90081DF @@ -13390,7 +13108,6 @@ 37661 78026JY 37662 78026JZ 37663 78026KA -37664 78026KB 37665 78026KC 37666 78026KD 37667 78026KE @@ -13409,30 +13126,23 @@ 37688 99025ELD 37689 99025ELE 37690 99025ELF -37691 99025ELG -37692 99025ELH 37693 99025ELJ 37694 99025ELK -37697 99025ELN +37695 99025ELL 37698 99025ELP 37699 99025ELQ 37701 99025ELS 37703 99025ELU 37704 99025ELV -37705 99025ELW 37706 99025ELX 37707 99025ELY 37708 99025ELZ -37710 99025EMB 37712 99025EMD 37713 99025EME 37714 99025EMF -37715 99025EMG 37716 99025EMH 37717 99025EMJ -37718 99025EMK 37719 99025EML -37722 99025EMP 37728 11029A 37730 11030A 37731 11030B @@ -13483,7 +13193,6 @@ 37794 11044G 37795 11044H 37796 11044J -37797 11044K 37799 11045B 37804 11047A 37806 11048A @@ -13492,7 +13201,6 @@ 37809 11049A 37810 11049B 37811 11049C -37812 11049D 37813 11050A 37814 11050B 37816 11051A @@ -13535,7 +13243,6 @@ 37869 11064C 37870 11064D 37871 11064E -37874 11066A 37875 11066B 37884 87079BP 37885 90105AF @@ -13548,11 +13255,12 @@ 37892 66040AX 37893 66040AY 37894 99025EMT -37896 99025EMV +37895 99025EMU 37899 99025EMY 37900 99025EMZ 37901 99025ENA 37902 99025ENB +37903 99025ENC 37904 99025END 37905 99025ENE 37906 99025ENF @@ -13561,8 +13269,8 @@ 37909 99025ENJ 37911 99025ENL 37912 99025ENM +37913 99025ENN 37914 99025ENP -37915 99025ENQ 37917 99025ENS 37918 99025ENT 37919 99025ENU @@ -13574,7 +13282,6 @@ 37926 99025EPB 37927 99025EPC 37928 99025EPD -37929 99025EPE 37930 11068A 37931 11068B 37932 11068C @@ -13589,7 +13296,6 @@ 37951 11074B 37952 11074C 37954 11075A -37955 11075B 37956 93036BPF 37957 93036BPG 37958 93036BPH @@ -13600,8 +13306,6 @@ 37966 93036BPR 37967 93036BPS 37968 93036BPT -37969 93036BPU -37970 93036BPV 37971 93036BPW 37973 93036BPY 37974 93036BPZ @@ -13610,6 +13314,7 @@ 37978 93036BQD 37979 93036BQE 37983 93036BQJ +37984 93036BQK 37985 93036BQL 37986 93036BQM 37988 93036BQP @@ -13618,7 +13323,6 @@ 37994 93036BQV 37995 93036BQW 37996 93036BQX -37997 93036BQY 37998 93036BQZ 37999 93036BRA 38000 93036BRB @@ -13636,13 +13340,12 @@ 38015 11077B 38016 97051YC 38017 97051YD -38018 97051YE 38019 97051YF 38020 97051YG +38022 97051YJ 38023 97051YK 38024 97051YL 38025 97051YM -38026 97051YN 38028 97051YQ 38030 97051YS 38033 97051YV @@ -13659,7 +13362,6 @@ 38047 12001B 38049 12002A 38050 12002B -38053 93036BRK 38054 93036BRL 38055 93036BRM 38056 93036BRN @@ -13697,7 +13399,6 @@ 38107 12013A 38109 12014A 38110 99025EPF -38112 99025EPH 38113 99025EPJ 38114 99025EPK 38116 99025EPM @@ -13708,22 +13409,20 @@ 38124 99025EPV 38125 99025EPW 38127 99025EPY +38129 99025EQA 38130 99025EQB 38131 99025EQC 38132 99025EQD 38133 99025EQE 38134 99025EQF -38136 99025EQH +38135 99025EQG 38137 99025EQJ 38138 99025EQK -38139 99025EQL 38140 99025EQM 38142 99025EQP 38143 99025EQQ 38144 99025EQR -38145 99025EQS 38147 99025EQU -38148 99025EQV 38149 99025EQW 38150 99025EQX 38154 99025ERB @@ -13738,7 +13437,6 @@ 38164 99025ERM 38165 99025ERN 38166 99025ERP -38167 99025ERQ 38168 99025ERR 38169 99025ERS 38170 99025ERT @@ -13769,7 +13467,6 @@ 38205 93036BSV 38208 93036BSY 38210 93036BTA -38211 93036BTB 38212 93036BTC 38213 93036BTD 38214 93036BTE @@ -13793,6 +13490,7 @@ 38237 97051ZM 38241 97051ZR 38243 97051ZT +38244 97051ZU 38245 12016A 38246 12016B 38247 12016C @@ -13811,7 +13509,6 @@ 38265 09069B 38266 09069C 38267 10040B -38268 10040C 38269 12020B 38270 12020C 38271 79095BR @@ -13824,10 +13521,11 @@ 38279 79095BZ 38280 79095CA 38282 88005AK -38283 88005AL -38285 88005AN +38284 88005AM 38286 88005AP +38287 88005AQ 38288 81043W +38290 82116J 38293 10059C 38294 10059D 38295 10059E @@ -13851,14 +13549,12 @@ 38314 11043D 38315 11043E 38316 11043F -38317 12001D 38318 12001E 38323 08056J 38324 08056K 38325 09072D 38326 10051H 38327 11068J -38328 01049QP 38329 91082CS 38330 88005AR 38331 12023A @@ -13885,14 +13581,12 @@ 38358 12031A 38360 65027BM 38361 65027BN -38362 65027BP 38363 65027BQ 38364 65027BR 38365 65027BS 38368 65027BV 38369 65027BW 38370 65027BX -38371 65027BY 38372 65027BZ 38373 65027CA 38374 65027CB @@ -13909,10 +13603,12 @@ 38386 65027CP 38387 65027CQ 38388 65027CR +38389 65027CS 38390 65027CT 38391 65027CU 38392 65027CV 38393 65027CW +38394 65027CX 38395 65027CY 38396 65027CZ 38397 65027DA @@ -13926,13 +13622,11 @@ 38405 65027DJ 38406 65027DK 38407 65027DL -38408 65027DM 38409 65027DN 38410 65027DP 38411 65027DQ 38412 65027DR 38413 65027DS -38414 65027DT 38415 65027DU 38416 65027DV 38417 65027DW @@ -13950,6 +13644,7 @@ 38429 65027EJ 38430 65027EK 38431 65027EL +38432 65027EM 38433 65027EN 38434 65027EP 38436 65027ER @@ -13966,12 +13661,10 @@ 38448 65027FD 38449 65027FE 38450 65027FF -38451 65027FG 38452 65027FH 38453 65027FJ 38454 65027FK 38455 65027FL -38456 65027FM 38457 65027FN 38458 65027FP 38459 65027FQ @@ -13986,6 +13679,7 @@ 38477 97051AAE 38478 93036BTN 38479 93036BTP +38480 93036BTQ 38481 93036BTR 38483 93036BTT 38485 93036BTV @@ -13994,16 +13688,20 @@ 38489 99025ERY 38491 99025ESA 38492 99025ESB +38493 99025ESC +38494 99025ESD 38495 99025ESE 38496 99025ESF +38497 99025ESG +38498 99025ESH 38500 99025ESK -38501 99025ESL +38502 99025ESM 38503 99025ESN 38504 99025ESP -38510 99025ESV +38506 99025ESR +38508 99025EST 38512 99025ESX 38513 99025ESY -38515 99025ETA 38516 99025ETB 38517 99025ETC 38518 99025ETD @@ -14020,7 +13718,6 @@ 38530 67048H 38531 71046C 38532 75064E -38534 79095CB 38535 79095CC 38536 83109D 38537 83109E @@ -14029,16 +13726,14 @@ 38541 88005AU 38542 88005AV 38543 88005AW -38544 88005AX 38546 91032F 38547 91032G +38549 68014C 38551 12035A 38552 12035B 38553 12035C 38554 12035D -38557 06006CX 38559 06006CZ -38564 06006DE 38565 08011L 38566 08011M 38567 08011N @@ -14051,6 +13746,7 @@ 38575 08011W 38576 08011X 38577 08011Y +38578 08011Z 38579 08011AA 38580 08011AB 38581 08011AC @@ -14064,21 +13760,19 @@ 38592 08011AP 38593 08011AQ 38594 08011AR -38597 08011AU +38598 08011AV 38599 08011AW 38600 08011AX 38601 08011AY 38602 08011AZ 38603 08011BA 38604 08011BB -38605 08011BC 38606 08011BD 38607 08011BE 38608 08011BF 38609 08011BG 38610 08011BH 38612 08011BK -38613 08011BL 38614 08011BM 38615 08011BN 38617 08011BQ @@ -14086,13 +13780,10 @@ 38622 08011BV 38624 08011BX 38625 08011BY -38629 08011CC 38631 08011CE 38632 08011CF 38636 08011CK 38640 08011CP -38641 08011CQ -38642 08011CR 38643 08011CS 38645 08011CU 38646 08011CV @@ -14117,10 +13808,12 @@ 38665 71035E 38666 71035F 38667 74071E -38670 07003AJ 38673 86010H +38674 06042D 38682 79009E +38683 79009F 38684 80018E +38685 80018F 38687 79101W 38688 79101X 38689 79101Y @@ -14129,11 +13822,12 @@ 38692 68081S 38693 68081T 38694 68081U +38695 68081V +38696 68081W 38697 68081X 38698 68081Y +38699 68081Z 38700 68081AA -38701 68081AB -38702 68081AC 38703 68081AD 38704 68081AE 38705 68081AF @@ -14144,7 +13838,9 @@ 38711 12039E 38714 99025ETF 38715 99025ETG +38716 99025ETH 38717 99025ETJ +38718 99025ETK 38720 99025ETM 38721 99025ETN 38722 99025ETP @@ -14154,7 +13850,6 @@ 38726 99025ETT 38727 99025ETU 38728 99025ETV -38729 99025ETW 38730 12040A 38732 94089D 38733 12041A @@ -14188,35 +13883,39 @@ 38780 12051C 38781 12051D 38782 12052A -38783 99025ETX 38784 99025ETY +38785 99025ETZ 38786 99025EUA 38789 99025EUD 38790 99025EUE 38792 99025EUG 38793 99025EUH +38794 99025EUJ 38795 99025EUK 38797 99025EUM 38799 99025EUP 38801 99025EUR -38802 99025EUS 38803 99025EUT -38806 99025EUW +38804 99025EUU +38807 99025EUX 38808 99025EUY 38809 99025EUZ 38810 99025EVA 38811 99025EVB 38813 99025EVD +38814 99025EVE 38816 99025EVG 38817 99025EVH 38818 99025EVJ 38819 99025EVK 38820 99025EVL +38821 99025EVM +38822 99025EVN +38823 99025EVP 38824 99025EVQ 38825 99025EVR 38826 99025EVS 38827 99025EVT -38828 99025EVU 38829 99025EVV 38831 99025EVX 38832 99025EVY @@ -14228,7 +13927,6 @@ 38840 83022M 38841 83022N 38842 05017F -38843 74029D 38845 88019G 38857 12055A 38858 12055B @@ -14244,19 +13942,10 @@ 38884 12044R 38886 12044T 38888 12044V -38890 12044X -38894 12044AB -38895 12044AC -38898 12044AF -38902 12044AK -38905 12044AN 38906 12044AP -38914 12044AX -38918 12044BB 38921 12044BE 38922 12044BF 38926 12044BK -38941 12044CA 38948 12044CH 38953 12059A 38954 12059B @@ -14264,7 +13953,6 @@ 38960 12044CT 38961 12044CU 38971 12044DE -38972 12044DF 38977 12061A 38978 12061B 38979 12061C @@ -14278,7 +13966,6 @@ 38997 12064A 38998 12064B 38999 12064C -39000 12064D 39008 12065A 39009 12065B 39010 12065C @@ -14307,8 +13994,6 @@ 39043 88006J 39046 88006M 39052 99057SZ -39054 92047AV -39055 90081DP 39057 13001A 39058 13001B 39059 13001C @@ -14356,7 +14041,7 @@ 39110 93036BTZ 39111 93036BUA 39112 99025EVZ -39113 99025EWA +39114 99025EWB 39119 09049J 39120 13011A 39122 13012A @@ -14427,10 +14112,8 @@ 39223 13041B 39224 13035E 39225 13035F -39226 13030C 39227 13042A 39228 13042B -39229 13042C 39231 13037H 39232 13043A 39233 13044A @@ -14455,7 +14138,6 @@ 39256 13050A 39260 13052A 39261 13052B -39262 13053A 39265 13055A 39266 13055B 39267 13055C @@ -14466,14 +14148,9 @@ 39285 13056A 39286 13056B 39287 13056C -39290 13055Y -39292 13055AA -39293 13055AB 39294 66040BC 39295 66040BD 39296 68081AG -39297 68081AH -39298 68081AJ 39300 92093JR 39301 92093JS 39302 92093JT @@ -14502,19 +14179,19 @@ 39325 08011DC 39326 08011DD 39327 08011DE +39328 08011DF 39330 08011DH 39332 10042D 39333 10042E 39334 10042F 39335 10042G -39336 10042H 39337 10042J 39338 10057C 39339 10057D 39340 10057E 39341 10057F 39342 10057G -39343 10057H +39344 10057J 39345 10057K 39346 10057L 39347 10057M @@ -14537,38 +14214,16 @@ 39364 13059B 39369 13055AC 39371 13060B -39372 11066G 39375 13062A 39376 13062B 39377 13062C 39380 13064A 39381 13064B 39383 13064D -39385 13064F -39386 13064G -39387 13064H -39388 13064J -39389 13064K -39390 13064L -39391 13064M 39392 13064N -39393 13064P -39394 13064Q -39395 13064R -39396 13064S 39397 13064T -39398 13064U -39399 13064V -39400 13064W -39401 13064X -39402 13064Y 39403 13064Z 39404 13064AA -39405 13064AB -39406 13064AC -39407 13064AD -39408 13064AE -39409 13064AF 39410 13065A 39411 13065B 39416 13066A @@ -14603,13 +14258,13 @@ 39446 13066AG 39447 78083T 39448 13066AJ -39449 13066AK 39450 13066AL 39451 13067A 39452 13067B 39453 13067C 39454 13067D 39455 13068A +39459 13070B 39460 13071A 39462 13072A 39475 13072P @@ -14641,7 +14296,6 @@ 39524 14007C 39533 14008A 39534 14008B -39535 99025EWC 39536 99025EWD 39537 99025EWE 39539 99025EWG @@ -14656,9 +14310,9 @@ 39549 93036BUF 39550 93036BUG 39552 93036BUJ -39553 93036BUK 39554 93036BUL 39574 14009C +39582 93036BUM 39583 93036BUN 39584 93036BUP 39585 93036BUQ @@ -14672,7 +14326,6 @@ 39594 97051AAG 39595 97051AAH 39596 99057TA -39597 01049QQ 39598 07003AK 39599 08011DK 39600 10057X @@ -14729,6 +14382,7 @@ 39662 87060Z 39663 87060AA 39664 87060AB +39665 87060AC 39666 87060AD 39667 87060AE 39668 87060AF @@ -14762,8 +14416,8 @@ 39707 92093JX 39708 92093JY 39709 92093JZ -39710 92093KA 39711 92093KB +39712 92093KC 39713 92093KD 39714 92093KE 39715 92093KF @@ -14802,15 +14456,15 @@ 39779 97051AAL 39782 97051AAP 39783 97051AAQ -39785 97051AAS 39786 97051AAT 39787 97051AAU 39788 97051AAV 39789 97051AAW +39790 97051AAX 39791 97051AAY 39792 97051AAZ 39793 97051ABA -39795 97051ABC +39796 97051ABD 39797 08046J 39798 08046K 39799 08046L @@ -14828,22 +14482,20 @@ 39812 87020EL 39813 87020EM 39814 87020EN +39815 87020EP 39816 87020EQ 39817 87020ER 39818 87020ES 39819 87020ET 39820 87020EU 39821 87020EV -39822 87020EW 39823 87020EX -39826 87020FA 39827 87020FB -39829 87020FD +39828 87020FC 39830 87020FE 39832 87020FG 39833 87020FH 39835 91009DT -39836 91009DU 39837 91009DV 39839 91009DX 39840 91009DY @@ -14852,26 +14504,25 @@ 39843 92093KV 39844 92093KW 39845 92093KX +39846 92093KY 39847 92093KZ -39848 92093LA 39849 92093LB 39851 92093LD 39852 92093LE 39853 92093LF -39854 92093LG 39855 92093LH 39856 92093LJ 39857 92093LK 39858 92093LL 39860 92093LN -39861 92093LP +39862 92093LQ 39863 92093LR 39864 92093LS 39865 92093LT 39866 92093LU 39867 92093LV 39868 93016AS -39869 93016AT +39870 93016AU 39871 93016AV 39872 93016AW 39873 93016AX @@ -14880,23 +14531,21 @@ 39877 93016BB 39878 93016BC 39879 93016BD +39881 93016BF 39882 93016BG -39883 93016BH +39884 93016BJ 39885 93016BK 39886 93016BL 39887 93016BM 39888 93016BN 39889 93016BP 39890 93016BQ -39891 93016BR 39892 93016BS 39893 93016BT +39895 93016BV 39896 93016BW -39897 93016BX -39898 93016BY 39899 93016BZ 39900 93016CA -39901 93016CB 39902 93016CC 39903 93016CD 39905 93016CF @@ -14910,8 +14559,6 @@ 39914 94029AFM 39915 94029AFN 39916 94029AFP -39917 94029AFQ -39918 94029AFR 39919 94029AFS 39921 94029AFU 39922 94029AFV @@ -14924,8 +14571,9 @@ 39929 94029AGC 39930 94029AGD 39932 94029AGF +39933 06006DF 39935 81053NB -39937 81053ND +39936 81053NC 39938 81053NE 39939 81053NF 39941 81053NH @@ -14933,6 +14581,7 @@ 39943 81053NK 39944 81053NL 39945 81053NM +39946 81053NN 39947 81053NP 39948 81053NQ 39949 81053NR @@ -14941,10 +14590,8 @@ 39952 81053NU 39953 81053NV 39954 81053NW -39955 81053NX 39956 81053NY 39957 81053NZ -39958 81053PA 39959 81053PB 39960 81053PC 39961 81053PD @@ -14954,6 +14601,7 @@ 39965 81053PH 39966 81053PJ 39967 81053PK +39968 81053PL 39969 87068CV 39970 87068CW 39971 87068CX @@ -14963,8 +14611,8 @@ 39975 87068DB 39976 87068DC 39977 87068DD +39978 87068DE 39979 87068DF -39980 87068DG 39981 87068DH 39982 87068DJ 39983 87068DK @@ -14973,6 +14621,8 @@ 39988 08011DM 39989 08011DN 39993 08011DS +39995 08011DU +39998 02040E 39999 02040F 40000 04042D 40001 14032A @@ -15073,7 +14723,6 @@ 40130 14050C 40136 14051A 40137 14051B -40138 14051C 40141 14052A 40143 14053A 40144 14053B @@ -15081,8 +14730,11 @@ 40147 14054B 40148 14054C 40149 14054D +40152 07003AL 40156 07003AQ +40157 07003AR 40158 10042K +40159 10042L 40160 10042M 40161 10042N 40162 10042P @@ -15096,7 +14748,7 @@ 40171 10057Z 40172 10057AA 40174 10057AC -40175 10057AD +40176 10057AE 40177 10057AF 40178 10057AG 40179 11077C @@ -15109,8 +14761,8 @@ 40187 11077L 40188 11077M 40189 11077N -40190 11077P 40191 11077Q +40193 11077S 40194 12008D 40195 12008E 40196 12008F @@ -15124,22 +14776,23 @@ 40204 12008P 40206 12008R 40207 12008S +40208 14055A 40209 14055B 40213 79095CL 40214 79095CM 40215 79095CN -40216 86088E 40217 88056C 40218 93036BUX -40219 93036BUY 40220 93036BUZ 40221 93036BVA 40223 93036BVC +40224 93036BVD 40225 93036BVE 40226 93036BVF 40227 93036BVG 40228 93036BVH 40230 93036BVK +40231 93036BVL 40232 93036BVM 40233 93036BVN 40234 93036BVP @@ -15148,19 +14801,16 @@ 40237 93036BVS 40238 93036BVT 40239 93036BVU -40241 94029AGH 40242 94029AGJ -40244 94029AGL +40243 94029AGK 40245 08046T 40248 97082H 40249 97082J 40251 97082L 40252 97082M -40253 97082N 40254 97082P 40255 97082Q 40256 97082R -40257 97082S 40258 14058A 40259 14058B 40260 14058C @@ -15182,7 +14832,6 @@ 40276 14063B 40277 14064A 40278 14064B -40279 14064C 40286 14066A 40287 14066B 40288 14066C @@ -15202,7 +14851,6 @@ 40304 14070G 40305 14071A 40306 14071B -40309 14071E 40310 14072A 40311 14073A 40314 98067FL @@ -15210,8 +14858,6 @@ 40316 14075B 40317 84072C 40318 84072D -40324 02005G -40325 02005H 40326 02005J 40327 02005K 40328 14000A @@ -15241,7 +14887,6 @@ 40355 14085A 40356 14085B 40358 14086A -40359 14086B 40360 14087A 40362 14088A 40364 14089A @@ -15250,6 +14895,7 @@ 40367 14090A 40368 14090B 40369 14090C +40374 15002A 40375 15002B 40376 15003A 40377 15003B @@ -15261,9 +14907,6 @@ 40384 15005A 40385 15005B 40386 15005C -40389 98067FM -40391 15007B -40394 95015F 40395 95015G 40396 95015H 40397 95015J @@ -15286,50 +14929,26 @@ 40414 95015AB 40415 95015AC 40416 95015AD -40417 95015AE 40418 95015AF -40419 95015AG 40420 15009A -40421 15009B -40422 98067FN -40423 98067FP 40424 15010A -40425 15010C +40425 15010B 40426 15010C -40427 98067FQ -40428 98067FR -40429 98067FS -40430 98067FT 40431 95015AH -40432 95015AJ -40433 95015AK 40434 95015AL 40435 95015AM 40436 95015AN 40437 95015AP 40438 95015AQ 40439 95015AR -40440 95015AS 40441 95015AT 40442 95015AU 40443 95015AV 40444 95015AW -40445 95015AX -40446 95015AY 40447 95015AZ 40448 95015BA 40449 95015BB 40450 95015BC -40451 98067FU -40452 98067FV -40453 98067FW -40454 98067FX -40455 98067FY -40456 98067FZ -40457 98067GA -40458 98067GB -40459 98067GC -40460 98067GD 40461 95015BD 40462 95015BE 40463 95015BF @@ -15357,7 +14976,6 @@ 40485 15011D 40486 15011E 40487 95015CA -40488 95015CB 40489 95015CC 40490 95015CD 40491 95015CE @@ -15372,7 +14990,6 @@ 40500 95015CP 40501 95015CQ 40502 95015CR -40503 95015CS 40504 95015CT 40505 15012A 40506 15012B @@ -15381,7 +14998,6 @@ 40509 95015CV 40510 95015CW 40511 95015CX -40512 95015CY 40513 95015CZ 40514 95015DA 40515 95015DB @@ -15410,21 +15026,17 @@ 40538 15015A 40539 15015B 40541 15014C -40542 15016A -40543 15016B 40544 15017A 40545 15017B 40546 15017C 40547 15018A 40548 15018B 40549 15019A -40550 15019B 40552 15020A 40553 15020B 40554 15020C 40555 15020D 40556 15020E -40557 95015DQ 40558 95015DR 40559 95015DS 40560 95015DT @@ -15451,14 +15063,966 @@ 40581 95015EQ 40582 95015ER 40583 95015ES -40584 95015ET 40585 95015EU 40586 95015EV -40587 95015EW -40588 15021A -40589 15021B -40590 15021C -40591 15021D +40592 64051F +40594 14088D +40600 95015EX +40601 95015EY +40602 95015EZ +40603 95015FA +40604 95015FB +40605 95015FC +40606 95015FD +40607 95015FE +40608 95015FF +40609 95015FG +40610 95015FH +40611 95015FJ +40612 95015FK +40613 15022A +40614 15022B +40615 15022C +40616 15022D +40617 15023A +40618 15023B +40641 66040BE +40642 66040BF +40644 66040BH +40645 66040BJ +40646 66040BK +40647 66040BL +40648 66040BM +40650 94074AF +40651 15025A +40663 15026A +40664 15026B +40665 15026C +40666 15026D +40672 95015FL +40673 95015FM +40674 95015FN +40675 95015FP +40676 95015FQ +40677 95015FR +40678 95015FS +40679 95015FT +40680 95015FU +40681 95015FV +40684 99025EWQ +40685 99025EWR +40686 99025EWS +40687 99025EWT +40689 99025EWV +40690 99025EWW +40691 99025EWX +40693 99025EWZ +40694 99025EXA +40696 99025EXC +40697 15028A +40698 15028B +40699 15029A +40701 15030A +40703 84072E +40704 84072F +40705 84072G +40706 84072H +40707 84072J +40708 84072K +40709 84072L +40710 84072M +40711 93055D +40715 15032A +40716 15032B +40717 15032C +40718 15032D +40719 15032E +40720 15032F +40721 95015FW +40726 98067GJ +40729 98067GM +40730 15033A +40731 15033B +40732 15034A +40733 15034B +40734 15034C +40735 15034D +40737 98067GP +40738 98067GQ +40739 98067GR +40740 98067GS +40741 98067GT +40746 15036A +40747 15036B +40748 15037A +40749 15037B +40750 15037C +40751 15037D +40752 76077FV +40754 76077FX +40755 76077FY +40756 76077FZ +40757 76077GA +40758 76077GB +40759 76077GC +40761 76077GE +40762 76077GF +40763 76077GG +40764 76077GH +40765 76077GJ +40766 76077GK +40768 76077GM +40769 76077GN +40770 76077GP +40771 76077GQ +40773 76077GS +40774 79104AB +40775 79104AC +40778 79104AF +40780 79104AH +40781 79104AJ +40783 89006CC +40784 89006CD +40785 89006CE +40786 89006CF +40787 89006CG +40788 89006CH +40792 93036BVV +40793 93036BVW +40794 93036BVX +40795 93036BVY +40796 93036BVZ +40797 93036BWA +40798 93036BWB +40799 93036BWC +40800 93036BWD +40802 93036BWF +40803 93036BWG +40804 93036BWH +40805 93036BWJ +40806 93036BWK +40807 93036BWL +40808 93036BWM +40810 93036BWP +40811 93036BWQ +40813 07003AS +40814 10042X +40815 10042Y +40816 10057AH +40818 10057AK +40819 10057AL +40822 10057AP +40823 10057AQ +40824 10057AR +40825 10057AS +40827 10057AU +40829 10057AW +40831 10057AY +40834 10057BB +40835 10057BC +40838 11077T +40839 11077U +40840 11077V +40841 11077W +40842 11077X +40843 11077Y +40845 11077AA +40848 11077AD +40849 11077AE +40850 11077AF +40851 11077AG +40852 11077AH +40853 12008T +40855 12008V +40856 12008W +40857 12008X +40858 12008Y +40860 12008AA +40861 12008AB +40863 12008AD +40864 12008AE +40867 12008AH +40868 12008AJ +40871 12035E +40874 15039A +40875 15039B +40876 15039C +40877 15039D +40878 15040A +40879 15040B +40880 15041A +40881 15041B +40882 15042A +40883 15042B +40884 15042C +40887 15044A +40888 15044B +40889 15045A +40890 15045B +40891 15045C +40892 15046A +40894 15047A +40895 15048A +40896 15048B +40898 98067GY +40899 15049A +40900 15049B +40901 15049C +40902 15049D +40903 15049E +40904 15049F +40905 15049G +40906 15049H +40907 15049J +40908 15049K +40909 15049L +40910 15049M +40911 15049N +40912 15049P +40913 15049Q +40914 82118G +40915 82118H +40916 82118J +40917 82118K +40918 82118L +40919 82118M +40920 15050A +40921 15050B +40922 15050C +40923 15050D +40924 85102G +40925 15051A +40926 15051B +40927 15051C +40928 15051D +40930 15052A +40931 15052B +40932 15052C +40933 15052D +40934 15052E +40935 15052F +40936 15052G +40937 15052H +40938 15053A +40940 15054A +40941 15054B +40942 15054C +40943 15054D +40946 15056A +40947 15056B +40949 98067HA +40950 98067HB +40951 98067HC +40952 98067HD +40953 98067HE +40954 98067HF +40955 98067HG +40956 98067HH +40958 15057A +40959 15057B +40960 15057C +40961 15057D +40962 98067HK +40963 98067HL +40964 15058A +40978 15058Q +40979 98067HM +40980 98067HN +40981 15058R +40982 15059A +40984 15060A +40985 15060B +40986 15060C +40987 08066D +40988 15061A +40989 15034E +40990 15034F +40991 15044C +40992 97051ABE +40993 97051ABF +40994 97051ABG +40995 97051ABH +40996 97051ABJ +40997 97051ABK +40998 97051ABL +40999 99025EXD +41000 99025EXE +41002 99025EXG +41005 99025EXK +41006 99025EXL +41007 99025EXM +41008 99025EXN +41009 99025EXP +41010 99025EXQ +41011 99025EXR +41012 99025EXS +41013 99025EXT +41014 99025EXU +41015 99025EXV +41016 99025EXW +41019 15062A +41020 15062B +41021 15063A +41022 15063B +41023 15061B +41024 15061C +41025 15061D +41026 15064A +41028 15065A +41029 15065B +41030 15065C +41031 15065D +41032 15066A +41033 15066B +41034 15067A +41035 15067B +41036 15068A +41037 15068B +41038 15069A +41039 15069B +41040 15064C +41041 15064D +41042 15064E +41044 15070B +41045 00055C +41046 00055D +41047 00055E +41048 00055F +41049 00055G +41050 00055H +41051 00055J +41052 00055K +41053 00055L +41054 00055M +41055 00055N +41056 00055P +41057 00055Q +41058 00055R +41059 00055S +41060 00055T +41061 00055U +41062 00055V +41063 00055W +41064 00055X +41065 00055Y +41066 00055Z +41067 00055AA +41068 00055AB +41069 00055AC +41070 00055AD +41071 00055AE +41072 00055AF +41073 00055AG +41074 00055AH +41075 00055AJ +41076 00055AK +41077 00055AL +41078 00055AM +41079 00055AN +41080 00055AP +41081 00055AQ +41082 00055AR +41083 00055AS +41084 00055AT +41085 00055AU +41086 00055AV +41087 00055AW +41088 00055AX +41089 00055AY +41090 00055AZ +41091 00055BA +41092 00055BB +41093 00055BC +41094 00055BD +41095 00055BE +41096 00055BF +41097 00055BG +41099 15071B +41103 15073A +41105 15074A +41106 15074B +41108 15074D +41113 15074J +41117 80085N +41118 91015Q +41119 11077AJ +41120 11077AK +41121 15075A +41122 15075B +41123 15075C +41126 00055BH +41127 00055BJ +41128 00055BK +41129 00055BL +41130 00055BM +41131 00055BN +41132 00055BP +41133 00055BQ +41134 00055BR +41135 00055BS +41136 00055BT +41137 00055BU +41138 00055BV +41139 00055BW +41140 00055BX +41141 00055BY +41142 00055BZ +41143 00055CA +41144 00055CB +41145 00055CC +41146 00055CD +41147 00055CE +41148 00055CF +41149 00055CG +41150 00055CH +41151 00055CJ +41152 00055CK +41153 00055CL +41154 00055CM +41155 00055CN +41156 00055CP +41157 00055CQ +41158 00055CR +41159 00055CS +41160 00055CT +41161 00055CU +41162 00055CV +41163 00055CW +41164 00055CX +41165 00055CY +41166 15077A +41167 15077B +41168 15077C +41169 15077D +41170 15077E +41171 15077F +41172 15077G +41173 15078A +41174 15079A +41175 15079B +41176 15079C +41177 15080A +41179 15081A +41180 15081B +41181 15081C +41182 15081D +41183 15081E +41184 15081F +41185 15081G +41186 15081H +41187 15081J +41188 15081K +41189 15081L +41191 15082A +41192 15082B +41193 15082C +41194 15083A +41195 15083B +41196 00055CZ +41197 00055DA +41198 00055DB +41199 00055DC +41200 00055DD +41201 00055DE +41202 00055DF +41203 00055DG +41204 00055DH +41205 00055DJ +41206 00055DK +41207 00055DL +41208 00055DM +41209 00055DN +41210 00055DP +41211 00055DQ +41212 00055DR +41213 00055DS +41214 00055DT +41215 00055DU +41216 00055DV +41217 00055DW +41218 00055DX +41219 00055DY +41220 00055DZ +41221 00055EA +41222 00055EB +41223 00055EC +41224 00055ED +41225 00055EE +41226 00055EF +41227 00055EG +41228 00055EH +41229 00055EJ +41230 00055EK +41231 00055EL +41232 00055EM +41233 00055EN +41234 00055EP +41235 00055EQ +41236 00055ER +41237 00055ES +41238 16001A +41239 16001B +41240 16002A +41241 16003A +41242 16003B +41243 00055ET +41244 00055EU +41245 00055EV +41246 00055EW +41247 00055EX +41248 00055EY +41249 00055EZ +41250 00055FA +41251 00055FB +41252 00055FC +41253 00055FD +41254 00055FE +41255 00055FF +41256 00055FG +41257 00055FH +41258 00055FJ +41259 00055FK +41260 00055FL +41261 00055FM +41262 00055FN +41263 00055FP +41264 00055FQ +41265 00055FR +41266 00055FS +41267 00055FT +41268 00055FU +41269 00055FV +41270 00055FW +41271 00055FX +41272 00055FY +41273 00055FZ +41274 00055GA +41275 00055GB +41276 00055GC +41277 00055GD +41278 00055GE +41279 00055GF +41280 00055GG +41281 00055GH +41282 00055GJ +41283 00055GK +41284 00055GL +41285 00055GM +41286 00055GN +41287 00055GP +41288 00055GQ +41289 00055GR +41290 00055GS +41291 00055GT +41292 00055GU +41293 00055GV +41294 00055GW +41295 00055GX +41296 00055GY +41297 00055GZ +41298 00055HA +41299 00055HB +41300 00055HC +41301 00055HD +41302 00055HE +41303 00055HF +41304 00055HG +41305 00055HH +41306 00055HJ +41307 00055HK +41308 16004A +41309 16004B +41310 16005A +41311 16005B +41312 16005C +41313 98067HP +41314 98067HQ +41315 16006A +41316 16006B +41317 12026D +41318 12026E +41319 12026F +41320 12026G +41321 12026H +41322 12026J +41323 12026K +41324 12026L +41325 16006C +41328 16007A +41329 16007B +41330 16008A +41331 16008B +41332 16009A +41333 16009B +41334 16010A +41335 16011A +41336 16011B +41337 16012A +41338 16012B +41339 16012C +41340 16012D +41341 16012E +41342 16012F +41343 00055HL +41344 00055HM +41345 00055HN +41346 00055HP +41347 00055HQ +41348 00055HR +41349 00055HS +41350 00055HT +41351 00055HU +41352 00055HV +41353 00055HW +41354 00055HX +41355 00055HY +41356 00055HZ +41357 00055JA +41358 00055JB +41359 00055JC +41360 00055JD +41361 00055JE +41362 00055JF +41363 00055JG +41364 00055JH +41365 00055JJ +41366 00055JK +41367 00055JL +41368 00055JM +41369 00055JN +41370 00055JP +41371 00055JQ +41372 00055JR +41373 00055JS +41374 00055JT +41375 00055JU +41376 00055JV +41377 00055JW +41378 00055JX +41379 00055JY +41380 16013A +41381 16013B +41382 16014A +41383 16014B +41384 16015A +41385 16015B +41386 16016A +41387 16016B +41390 16017C +41391 16018A +41394 16020A +41395 16020B +41396 00055JZ +41397 00055KA +41398 00055KB +41399 00055KC +41400 00055KD +41401 00055KE +41402 00055KF +41403 00055KG +41404 00055KH +41405 00055KJ +41406 00055KK +41407 00055KL +41408 00055KM +41409 00055KN +41410 00055KP +41411 00055KQ +41412 00055KR +41413 00055KS +41414 00055KT +41415 00055KU +41416 00055KV +41417 00055KW +41418 00055KX +41419 00055KY +41420 00055KZ +41421 00055LA +41422 00055LB +41423 00055LC +41424 00055LD +41425 00055LE +41426 00055LF +41427 00055LG +41428 00055LH +41429 00055LJ +41430 00055LK +41431 00055LL +41432 00055LM +41433 00055LN +41434 16021A +41435 16021B +41436 16022A +41439 16012H +41440 16012J +41441 16012K +41442 16012L +41444 16012N +41445 16012P +41446 16012Q +41447 16012R +41455 16023E +41456 16025A +41457 16025B +41458 16025C +41459 16025D +41460 16025E +41461 16025F +41463 98067HT +41464 16026A +41465 16026B +41466 16026C +41469 16027A +41470 16027B +41471 16028A +41472 16028B +41473 16029A +41474 98067HU +41475 98067HV +41476 98067HW +41477 98067HX +41478 98067HY +41479 98067HZ +41480 98067JA +41481 98067JB +41482 98067JC +41483 98067JD +41484 98067JE +41485 98067JF +41486 98067JG +41487 98067JH +41488 98067JJ +41489 98067JK +41490 98067JL +41491 00055LP +41492 00055LQ +41493 00055LR +41494 00055LS +41495 00055LT +41496 00055LU +41497 00055LV +41498 00055LW +41499 00055LX +41500 00055LY +41501 00055LZ +41502 00055MA +41503 00055MB +41504 00055MC +41505 00055MD +41506 00055ME +41507 00055MF +41508 00055MG +41509 00055MH +41510 00055MJ +41511 00055MK +41512 00055ML +41513 00055MM +41514 00055MN +41515 00055MP +41516 00055MQ +41517 00055MR +41518 00055MS +41519 00055MT +41520 00055MU +41521 00055MV +41522 00055MW +41523 00055MX +41524 00055MY +41525 00055MZ +41526 00055NA +41527 00055NB +41528 00055NC +41529 00055ND +41530 00055NE +41531 00055NF +41532 00055NG +41533 00055NH +41534 00055NJ +41535 00055NK +41536 00055NL +41537 00055NM +41538 00055NN +41539 00055NP +41540 00055NQ +41541 14087C +41543 15075D +41544 15075E +41546 15075G +41547 15075H +41548 15075J +41549 16030A +41550 16030B +41551 16030C +41552 16031A +41553 16031B +41554 16032A +41555 16032B +41556 16033A +41557 16033B +41558 16033C +41559 16033D +41560 08067J +41561 08067K +41562 08067L +41563 98067JM +41564 98067JN +41565 98067JP +41566 98067JQ +41567 98067JR +41568 98067JS +41569 98067JT +41570 98067JU +41571 98067JV +41572 98067JW +41573 98067JX +41574 98067JY +41575 98067JZ +41576 98067KA +41577 98067KB +41578 98067KC +41579 16034A +41580 16034B +41581 16035A +41582 16035B +41583 16035C +41584 16036A +41585 16036B +41586 16037A +41587 16037B +41588 16038A +41589 16038B +41590 16038C +41591 16039A +41592 16039B +41593 16039C +41594 16039D +41595 16019B +41596 16019C +41597 16019D +41598 16019E +41599 16040A +41600 16040B +41601 16040C +41602 16040D +41603 16040E +41604 16040F +41605 16040G +41606 16040H +41607 16040J +41608 16040K +41609 16040L +41610 16040M +41611 16040N +41612 16040P +41613 16040Q +41614 16040R +41615 16040S +41616 16040T +41617 16040U +41618 16040V +41620 16040X +41621 16040Y +41622 16041A +41623 16041B +41625 16042B +41626 16042C +41627 16042D +41628 16042E +41629 16042F +41630 16042G +41631 16042H +41632 16042J +41633 16042K +41634 16043A +41635 16043B +41636 16042L +41637 16042M +41639 16044A +41640 16044B +41641 08067M +41642 08067N +41643 08067P +41644 08067Q +41645 08067R +41646 08067S +41647 08067T +41648 08067U +41649 08067V +41650 08067W +41651 08067X +41652 00055NR +41653 00055NS +41654 00055NT +41655 00055NU +41656 00055NV +41657 00055NW +41658 00055NX +41659 00055NY +41660 00055NZ +41661 00055PA +41662 00055PB +41663 00055PC +41664 00055PD +41665 00055PE +41666 00055PF +41667 00055PG +41668 00055PH +41669 00055PJ +41670 16045A +41671 16045B +41672 16046A +41673 16046B +41674 16046C +41675 94076J +41676 07052H +41677 07052J +41678 07052K +41679 07052L +41680 07052M +41681 07052N +41682 07052P +41683 07052Q +41684 07052R +41685 07052S +41686 07052T +41687 07052U +41688 07052V +41689 07052W +41690 07052X +41691 07052Y +41692 07052Z +41693 07052AA +41694 07052AB +41695 10007M +41696 10007N +41697 10007P +41698 10007Q +41699 10007R +41700 10007S +41701 10007T +41702 10007U +41703 10007V +41704 10007W +41705 00055PK +41706 00055PL +41707 00055PM +41708 00055PN +41709 00055PP +41710 00055PQ +41711 00055PR +41712 00055PS +41713 00055PT +41714 00055PU +41715 00055PV +41716 00055PW +41717 00055PX +41718 15020F +41719 15020G +41720 15020H +41721 15020J +41722 15020K +41723 15020L +41724 16047A 90004 99750A 90006 00596A 90007 90097E @@ -15486,7 +16050,6 @@ 90099 14792C 90100 14798A 90101 14817A -90102 14819A 90103 14820A 90105 14829A 90106 14844A @@ -15494,6 +16057,13 @@ 90108 14845A 90109 15523A 90110 15545A +90111 15607A +90112 15614A +90113 16535A +90114 16591A +90115 16612A +90116 16697A +90117 16707A 96031 96531A 96041 96541A 96044 96544A @@ -15501,13 +16071,10 @@ 96137 80060G 96168 96668A 99079 13825B -99091 14547A 99096 14552E 99208 13594A 99210 13685A 99212 13706A -99221 14544A 99224 14658A 99406 14602A -40651 15025A 99999 99000A diff --git a/data/sites.txt b/data/sites.txt index 9a45b8b..77e6409 100644 --- a/data/sites.txt +++ b/data/sites.txt @@ -15,6 +15,7 @@ 0434 IR -26.1030 27.9288 1646 Ian Roberts 0710 LS 52.3261 10.6756 85 Lutz Schindler 1056 MK 57.0122 23.9833 4 Martins Keruss +1086 NK 46.4778 30.7556 56 Nikolay Koshkin 1244 AM 44.3932 33.9701 69 Andriy Makeyev 1747 DD 45.7275 -72.3526 191 Daniel Deak 1775 KF 44.6062 -75.6910 200 Kevin Fetter @@ -33,6 +34,7 @@ 4541 AR 41.9639 12.4531 80 Alberto Rango 4542 AR 41.9683 12.4545 80 Alberto Rango 4641 AR 41.1060 16.9010 70 Alberto Rango +5555 SG 56.1019 94.5533 154 Sergey Guryanov 5918 BG 59.2985 18.1045 40 Bjorn Gimle 5919 BG 59.2615 18.6206 30 Bjorn Gimle 6226 SC 28.4861 -97.8194 110 Scott Campbell @@ -45,14 +47,17 @@ 8536 TL 36.8479 -76.5010 4 Tim Luton 8539 SN 39.4707 -79.3388 839 Steve Newcomb 8597 TB -34.9638 138.6333 100 Tony Beresford +8600 PC -32.9770 151.6477 18 Paul Camilleri 8730 EC 30.3086 -97.7279 150 Ed Cannon 8739 DB 37.1133 -121.7028 282 Derek Breit +9461 KO 47.9175 19.89365 938 Konkoly Obs 9633 JN 33.0206 -96.9982 153 Jim Nix 9730 MM 30.3150 -97.8660 280 Mike McCants 6242 JM 42.9453 -2.8284 623 Jon Mikel 6241 JM 42.9565 -2.8203 619 Jon Mikel 4160 BD 51.2793 5.4768 35 Bram Dorreman -9001 MS -32.3793 20.810682 1760 Master K95 +9000 PP 40.0848 21.1581 1433 Pierros Papadeas +9001 AF 44.5700 6.68 1850 Alain Figer +9002 AO 62.2544 26.5962 150 Arto Oksanen 9999 GR 47.348 5.5151 100 Graves -9000 NJ 51.4198 5.4086 25 Nico Janssen -9905 DR -26.8724 26.6481 1200 Deon van Rooyen + diff --git a/detect.c b/detect.c index 314dc2a..67d4e3c 100644 --- a/detect.c +++ b/detect.c @@ -684,13 +684,14 @@ int main(int argc,char *argv[]) printf("%d points above %g sigma\n",n,sigma); // Exit if too many + /* if (n>2500) { file=fopen("skipped.dat","a"); fprintf(file,"%s : %d points above %g sigma\n",fitsfile,n,sigma); fclose(file); return 0; } - + */ // Fill points p=(struct point *) malloc(sizeof(struct point)*n); for (i=0,l=0;i +#include +#include +#include +#include +#include +#include "sgdp4h.h" + +#define LIM 128 +#define D2R M_PI/180.0 +#define R2D 180.0/M_PI +#define XKMPER 6378.135 // Earth radius in km +#define XKMPAU 149597879.691 // AU in km +#define FLAT (1.0/298.257) +#define STDMAG 6.0 +#define MMAX 1024 + + +long Isat=0; +long Isatsel=0; +extern double SGDP4_jd0; +struct map { + double lat,lng; + float alt; + char observer[32]; + int site_id; +} m; +void get_site(int site_id); +double nfd2mjd(char *date); +double date2mjd(int year,int month,double day); +void nfd_now(char *s); +void obspos_xyz(double mjd,xyz_t *pos,xyz_t *vel); +void sunpos_xyz(double mjd,xyz_t *pos,double *ra,double *de); +double gmst(double mjd); +double dgmst(double mjd); +double modulo(double x,double y); +void equatorial2horizontal(double mjd,double ra,double de,double *azi,double *alt); +void mjd2date(double mjd,char *date); +int properties(kep_t K,xyz_t obspos,xyz_t sunpos,float radius,float t,double *ra,double *de,double *r,float *mag); +void dec2sex(double x,char *s,int f,int len); + +void usage(void) +{ + return; +} + +int main(int argc,char *argv[]) +{ + int arg=0,imode,i; + int satno=-1; + orbit_t orb; + FILE *fp; + char tlefile[256]="classfd.tle"; + char nfd[32]; + double mjd0,mjd,jd,tsince; + int rv,withvel; + kep_t K; + double radius=1080.0; + xyz_t obspos,obsvel,sunpos; + double p,pmin,p0,p1,r,ra,de,azi,alt,sazi,salt,altmin=10.0,saltmin=-6.0; + float mag,mmin; + int state,pstate,nstate; + float t,length=86400.0,dt=60.0; + char sra[16],sde[16],type[32]; + + // Initialize + nfd_now(nfd); + mjd0=nfd2mjd(nfd); + get_site(4171); + + // Decode options + while ((arg=getopt(argc,argv,"t:c:i:s:l:hS:A:r:"))!=-1) { + switch (arg) { + + case 't': + strcpy(nfd,optarg); + mjd0=nfd2mjd(nfd); + break; + + case 'c': + strcpy(tlefile,optarg); + break; + + case 's': + get_site(atoi(optarg)); + break; + + case 'i': + satno=atoi(optarg); + break; + + case 'r': + radius=atof(optarg); + break; + + case 'l': + length=atoi(optarg); + if (strchr(optarg,'h')!=NULL) + length*=3600; + else if (strchr(optarg,'m')!=NULL) + length*=60; + else if (strchr(optarg,'d')!=NULL) + length*=86400; + break; + + case 'S': + saltmin=atof(optarg); + break; + + case 'A': + altmin=atof(optarg); + break; + + case 'h': + usage(); + return 0; + break; + + default: + usage(); + return 0; + } + } + + // Error checking + if (satno<=0) { + fprintf(stderr,"ERROR: NORAD satellite number not provided!\n"); + return -1; + } + + // Get TLE + fp=fopen(tlefile,"rb"); + if (fp==NULL) { + fprintf(stderr,"ERROR: Failed to open file with TLEs (%s)!\n",tlefile); + return -1; + } + + // Read TLE + while (read_twoline(fp,satno,&orb)==0) { + Isat=orb.satno; + imode=init_sgdp4(&orb); + + if (imode==SGDP4_ERROR) + continue; + + } + fclose(fp); + + // Object found? + if (orb.satno!=satno) { + fprintf(stderr,"ERROR: Object %d not found in %s!\n",satno,tlefile); + return -1; + } + + // Object found? + if (orb.rev<10.0) { + fprintf(stderr,"ERROR: Object %d is not in a LEO orbit.\n",satno); + return -1; + } + + // Print header + printf("Observer: %s (%04d) [%+.4f, %+.4f, %.0fm]\n",m.observer,m.site_id,m.lat,m.lng,m.alt*1000.0); + printf("Elements: %s\n",tlefile); + printf("Object: %d\n",satno); + printf("Start UT Date/Time: %s for %g h \n\n",nfd,length/3600.0); + printf("UT Date/Time R.A. Decl. Azi. Alt. Range Mag Sun Alt. Type\n"); + printf(" (deg) (deg) (km) (deg)\n"); + printf("=====================================================================================\n"); + + + + for (t=0.0;tnstate) + strcpy(type,"Ingress"); + else + strcpy(type,"Maximum"); + ra=modulo(ra,360.0); + equatorial2horizontal(mjd,ra,de,&azi,&alt); + azi=modulo(azi+180.0,360.0); + mjd2date(mjd,nfd); + dec2sex(ra/15.0,sra,0,2); + dec2sex(de,sde,0,2); + if (alt>altmin && salttm_year+1900,ptm->tm_mon+1,ptm->tm_mday,ptm->tm_hour,ptm->tm_min,ptm->tm_sec); + + return; +} + +// Observer position +void obspos_xyz(double mjd,xyz_t *pos,xyz_t *vel) +{ + double ff,gc,gs,theta,s,dtheta; + + s=sin(m.lat*D2R); + ff=sqrt(1.0-FLAT*(2.0-FLAT)*s*s); + gc=1.0/ff+m.alt/XKMPER; + gs=(1.0-FLAT)*(1.0-FLAT)/ff+m.alt/XKMPER; + + theta=gmst(mjd)+m.lng; + dtheta=dgmst(mjd)*D2R/86400; + + pos->x=gc*cos(m.lat*D2R)*cos(theta*D2R)*XKMPER; + pos->y=gc*cos(m.lat*D2R)*sin(theta*D2R)*XKMPER; + pos->z=gs*sin(m.lat*D2R)*XKMPER; + vel->x=-gc*cos(m.lat*D2R)*sin(theta*D2R)*XKMPER*dtheta; + vel->y=gc*cos(m.lat*D2R)*cos(theta*D2R)*XKMPER*dtheta; + vel->z=0.0; + + return; +} + +// Solar position +void sunpos_xyz(double mjd,xyz_t *pos,double *ra,double *de) +{ + double jd,t,l0,m,e,c,r; + double n,s,ecl; + + jd=mjd+2400000.5; + t=(jd-2451545.0)/36525.0; + l0=modulo(280.46646+t*(36000.76983+t*0.0003032),360.0)*D2R; + m=modulo(357.52911+t*(35999.05029-t*0.0001537),360.0)*D2R; + e=0.016708634+t*(-0.000042037-t*0.0000001267); + c=(1.914602+t*(-0.004817-t*0.000014))*sin(m)*D2R; + c+=(0.019993-0.000101*t)*sin(2.0*m)*D2R; + c+=0.000289*sin(3.0*m)*D2R; + + r=1.000001018*(1.0-e*e)/(1.0+e*cos(m+c)); + n=modulo(125.04-1934.136*t,360.0)*D2R; + s=l0+c+(-0.00569-0.00478*sin(n))*D2R; + ecl=(23.43929111+(-46.8150*t-0.00059*t*t+0.001813*t*t*t)/3600.0+0.00256*cos(n))*D2R; + + *ra=atan2(cos(ecl)*sin(s),cos(s))*R2D; + *de=asin(sin(ecl)*sin(s))*R2D; + + pos->x=r*cos(*de*D2R)*cos(*ra*D2R)*XKMPAU; + pos->y=r*cos(*de*D2R)*sin(*ra*D2R)*XKMPAU; + pos->z=r*sin(*de*D2R)*XKMPAU; + + return; +} + +// Greenwich Mean Sidereal Time +double dgmst(double mjd) +{ + double t,dgmst; + + t=(mjd-51544.5)/36525.0; + + dgmst=360.98564736629+t*(0.000387933-t/38710000); + + return dgmst; +} + +// Greenwich Mean Sidereal Time +double gmst(double mjd) +{ + double t,gmst; + + t=(mjd-51544.5)/36525.0; + + gmst=modulo(280.46061837+360.98564736629*(mjd-51544.5)+t*t*(0.000387933-t/38710000),360.0); + + return gmst; +} + +// Return x modulo y [0,y) +double modulo(double x,double y) +{ + x=fmod(x,y); + if (x<0.0) x+=y; + + return x; +} + +// Convert equatorial into horizontal coordinates +void equatorial2horizontal(double mjd,double ra,double de,double *azi,double *alt) +{ + double h; + + h=gmst(mjd)+m.lng-ra; + + *azi=modulo(atan2(sin(h*D2R),cos(h*D2R)*sin(m.lat*D2R)-tan(de*D2R)*cos(m.lat*D2R))*R2D,360.0); + *alt=asin(sin(m.lat*D2R)*sin(de*D2R)+cos(m.lat*D2R)*cos(de*D2R)*cos(h*D2R))*R2D; + + return; +} + +// Compute Date from Julian Day +void mjd2date(double mjd,char *date) +{ + double f,jd,dday; + int z,alpha,a,b,c,d,e; + int year,month,day,hour,min; + float sec,x; + + jd=mjd+2400000.5; + jd+=0.5; + + z=floor(jd); + f=fmod(jd,1.); + + if (z<2299161) + a=z; + else { + alpha=floor((z-1867216.25)/36524.25); + a=z+1+alpha-floor(alpha/4.); + } + b=a+1524; + c=floor((b-122.1)/365.25); + d=floor(365.25*c); + e=floor((b-d)/30.6001); + + dday=b-d-floor(30.6001*e)+f; + if (e<14) + month=e-1; + else + month=e-13; + + if (month>2) + year=c-4716; + else + year=c-4715; + + day=(int) floor(dday); + x=24.0*(dday-day); + x=3600.*fabs(x)+0.0001; + sec=fmod(x,60.); + x=(x-sec)/60.; + min=fmod(x,60.); + x=(x-min)/60.; + hour=x; + sec=floor(1000.0*sec)/1000.0; + + sprintf(date,"%04d-%02d-%02dT%02d:%02d:%02.0f",year,month,day,hour,min,sec); + + return; +} + +int properties(kep_t K,xyz_t obspos,xyz_t sunpos,float radius,float t,double *ra,double *de,double *r,float *mag) +{ + double st,ct,sn,cn,si,ci; + xyz_t satpos; + double dx,dy,dz,rsun,rearth,psun,pearth,p; + float phase; + int state; + + // Angles + sn=sin(K.ascn); + cn=cos(K.ascn); + si=sin(K.eqinc); + ci=cos(K.eqinc); + st=sin(K.theta+t); + ct=cos(K.theta+t); + + satpos.x=-sn*ci*st+cn*ct; + satpos.y=cn*ci*st+sn*ct; + satpos.z=si*st; + satpos.x*=(radius+XKMPER); + satpos.y*=(radius+XKMPER); + satpos.z*=(radius+XKMPER); + + // Sun position from satellite + dx=-satpos.x+sunpos.x; + dy=-satpos.y+sunpos.y; + dz=-satpos.z+sunpos.z; + + // Distances + rsun=sqrt(dx*dx+dy*dy+dz*dz); + rearth=sqrt(satpos.x*satpos.x+satpos.y*satpos.y+satpos.z*satpos.z); + + // Angles + psun=asin(696.0e3/rsun)*R2D; + pearth=asin(6378.135/rearth)*R2D; + + pearth=asin(6378.135/rearth)*R2D; + p=acos((-dx*satpos.x-dy*satpos.y-dz*satpos.z)/(rsun*rearth))*R2D; + p-=pearth; + + // Position differences + dx=satpos.x-obspos.x; + dy=satpos.y-obspos.y; + dz=satpos.z-obspos.z; + + // Celestial position + *r=sqrt(dx*dx+dy*dy+dz*dz); + *ra=atan2(dy,dx)*R2D; + *de=asin(dz/ *r)*R2D; + + // Visibility + if (p<-psun) { + // strcpy(state,"eclipsed"); + state=0; + } else if (p>-psun && ppsun) { + // strcpy(state,"sunlit"); + state=2; + } + + // Phase + phase=acos(((obspos.x-satpos.x)*(sunpos.x-satpos.x)+(obspos.y-satpos.y)*(sunpos.y-satpos.y)+(obspos.z-satpos.z)*(sunpos.z-satpos.z))/(rsun* *r))*R2D; + + // Magnitude + if (state==2) + *mag=STDMAG-15.0+5*log10( *r)-2.5*log10(sin(phase*D2R)+(M_PI-phase*D2R)*cos(phase*D2R)); + else + *mag=15; + + return state; +} + +// Convert Decimal into Sexagesimal +void dec2sex(double x,char *s,int f,int len) +{ + int i; + double sec,deg,min; + char sign; + char *form[]={"::",",,","hms"," "}; + + sign=(x<0 ? '-' : ' '); + x=3600.*fabs(x); + + sec=fmod(x,60.); + x=(x-sec)/60.; + min=fmod(x,60.); + x=(x-min)/60.; + // deg=fmod(x,60.); + deg=x; + + if (len==7) sprintf(s,"%c%02i%c%02i%c%07.4f%c",sign,(int) deg,form[f][0],(int) min,form[f][1],sec,form[f][2]); + if (len==6) sprintf(s,"%c%02i%c%02i%c%06.3f%c",sign,(int) deg,form[f][0],(int) min,form[f][1],sec,form[f][2]); + if (len==5) sprintf(s,"%c%02i%c%02i%c%05.2f%c",sign,(int) deg,form[f][0],(int) min,form[f][1],sec,form[f][2]); + if (len==4) sprintf(s,"%c%02i%c%02i%c%04.1f%c",sign,(int) deg,form[f][0],(int) min,form[f][1],sec,form[f][2]); + if (len==2) sprintf(s,"%c%02i%c%02i%c%02i%c",sign,(int) deg,form[f][0],(int) min,form[f][1],(int) floor(sec),form[f][2]); + + return; +} diff --git a/residuals.c b/residuals.c index fcbdc3e..2117bb9 100644 --- a/residuals.c +++ b/residuals.c @@ -241,6 +241,12 @@ struct point decode_iod_observation(char *iod_line) sscanf(iod_line+16,"%4d",&site_id); s=get_site(site_id); + // Skip if site not found + if (s.id<0) { + fprintf(stderr,"Site %d not found!\n",site_id); + p.flag=0; + } + // Decode date/time sscanf(iod_line+23,"%4d%2d%2d%2d%2d%5s",&year,&month,&iday,&hour,&min,secbuf); sec=atof(secbuf); @@ -298,7 +304,7 @@ struct point decode_iod_observation(char *iod_line) de=sign*(de+dd/100+ds/10000); break; default : - printf("IOD Format not implemented\n"); + fprintf(stderr,"IOD Format not implemented\n"); p.flag=0; break; } @@ -315,7 +321,7 @@ struct point decode_iod_observation(char *iod_line) } else if (epoch==5) { mjd0=51544.5; } else { - printf("Observing epoch not implemented\n"); + fprintf(stderr,"Observing epoch not implemented\n"); p.flag=0; } @@ -373,6 +379,9 @@ struct site get_site(int site_id) } fclose(file); + if (id!=site_id) + s.id==-1; + return s; } // Return x modulo y [0,y) diff --git a/rv2tle.c b/rv2tle.c index 5ba1e36..bdd728a 100644 --- a/rv2tle.c +++ b/rv2tle.c @@ -7,7 +7,7 @@ #include "satutl.h" #include -#define LIM 128 +#define LIM 1024 #define XKE 0.07436680 // Guassian Gravitational Constant #define XKMPER 6378.135 #define AE 1.0 @@ -342,7 +342,7 @@ int fgetline(FILE *file,char *s,int lim) int main(int argc,char *argv[]) { - int imode,satno=99000,arg; + int imode,satno=99000,arg,gmat=0; FILE *file; orbit_t orb; xyz_t r,v; @@ -352,7 +352,7 @@ int main(int argc,char *argv[]) char *env; // Decode options - while ((arg=getopt(argc,argv,"p:i:d:"))!=-1) { + while ((arg=getopt(argc,argv,"p:i:d:g"))!=-1) { switch (arg) { case 'p': @@ -372,6 +372,10 @@ int main(int argc,char *argv[]) return 0; break; + case 'g': + gmat=1; + break; + default: usage(); return 0; @@ -384,6 +388,11 @@ int main(int argc,char *argv[]) while (fgetline(file,line,LIM)>0) { sscanf(line,"%lf %lf %lf %lf %lf %lf %lf",&mjd,&r.x,&r.y,&r.z,&v.x,&v.y,&v.z); + // Convert to MJD + if (gmat==1) + mjd+=29999.5; + + // Convert orb=rv2el(orb.satno,mjd,r,v); orb.satno=satno; diff --git a/satfit.c b/satfit.c index 16af9e0..cfde9ed 100644 --- a/satfit.c +++ b/satfit.c @@ -654,9 +654,12 @@ int main(int argc,char *argv[]) // Adjust if (adjust==1) { // Count observations - for (i=0,nobs=0;i or , in mag.arcsec-2 -ANALYSIS_THRESH 1.5 # or , in mag.arcsec-2 +DETECT_THRESH 2.5 # or , in mag.arcsec-2 +ANALYSIS_THRESH 2.5 # or , in mag.arcsec-2 FILTER Y # apply filter for detection ("Y" or "N")? FILTER_NAME $ST_DATADIR/sextractor/default.conv # name of the file containing the filter diff --git a/skymap.c b/skymap.c index 687192a..d3ecffc 100644 --- a/skymap.c +++ b/skymap.c @@ -1078,6 +1078,7 @@ double modulo(double x,double y) return x; } + // Read a line of maximum length int lim from file FILE into string s int fgetline(FILE *file,char *s,int lim) { diff --git a/tleinfo.c b/tleinfo.c index 47b13f2..e1d7a5f 100644 --- a/tleinfo.c +++ b/tleinfo.c @@ -12,12 +12,16 @@ #define AE 1.0 /* Earth radius in "chosen units". */ #define XKE 0.743669161e-1 #define CK2 5.413080e-4 /* (0.5 * XJ2 * AE * AE) */ +#define D2R M_PI/180.0 +#define R2D 180.0/M_PI + extern double SGDP4_jd0; void usage(void) { return; } +double modulo(double x,double y); // Compute Julian Day from Date double date2mjd(int year,int month,double day) @@ -151,6 +155,23 @@ void mjd2nfd(double mjd,char *nfd) return; } +float orbital_longitude_at_midnight(orbit_t orb,double mjd0) +{ + int rv,imode; + double jd,tsince,mjd; + kep_t K; + + imode=init_sgdp4(&orb); + + mjd=floor(mjd0); + + jd=mjd+2400000.5; + tsince=1440.0*(jd-SGDP4_jd0); + rv=sgdp4(tsince,1,&K); + + return modulo(K.theta*R2D,360.0); +} + int main(int argc,char *argv[]) { int arg=0,satno=0,header=0,oneline=0,no,name=0,desig=0; @@ -158,7 +179,7 @@ int main(int argc,char *argv[]) char line0[LIM],line1[LIM],line2[LIM],nfd[32]; FILE *file; orbit_t orb; - float aodp,perigee,apogee,period; + float aodp,perigee,apogee,period,lng; int info=0; double mjd; char *env; @@ -167,7 +188,7 @@ int main(int argc,char *argv[]) sprintf(tlefile,"%s/bulk.tle",env); // Decode options - while ((arg=getopt(argc,argv,"c:i:aH1ftnd"))!=-1) { + while ((arg=getopt(argc,argv,"c:i:aH1ftndb"))!=-1) { switch (arg) { case 'c': @@ -198,6 +219,10 @@ int main(int argc,char *argv[]) info=1; break; + case 'b': + info=2; + break; + case 'H': header=1; break; @@ -281,6 +306,10 @@ int main(int argc,char *argv[]) mjd2nfd(mjd,nfd); if (info==0) printf("%05d %10.4lf %8.4f %8.4f %8.4f %8.4f %8.6f %8.5f\n",orb.satno,mjd,DEG(orb.eqinc),DEG(orb.ascn),DEG(orb.argp),DEG(orb.mnan),orb.ecc,orb.rev); if (info==1) printf("%05d %6.0f x %6.0f x %6.2f %8.2f %8.6f %14.8lf\n",orb.satno,perigee,apogee,DEG(orb.eqinc),period,orb.ecc,mjd); + if (info==2) { + lng=orbital_longitude_at_midnight(orb,mjd); + printf("%05d %10.4lf %8.4f %8.4f %8.4f %8.4f %8.6f %8.5f %10.4lf %8.4f\n",orb.satno,mjd,DEG(orb.eqinc),DEG(orb.ascn),DEG(orb.argp),DEG(orb.mnan),orb.ecc,orb.rev,floor(mjd),lng); + } } fclose(file); } else if (oneline==2) { @@ -302,3 +331,12 @@ int main(int argc,char *argv[]) return 0; } + +// Return x modulo y [0,y) +double modulo(double x,double y) +{ + x=fmod(x,y); + if (x<0.0) x+=y; + + return x; +}