Generating a schedule for NGC 2903 observations last update: 7Feb2005 by mph comments inserted by jai: March 1, 2005 (see *** comments below) What you need: IDL (we run version 6.0, which at AO is called "idl60"; idl is aliased to that for the ngc2903 account) Phil's routines (2005Feb version: wasinit2n)**** DIFFERENT**** The phil/gen routines are invoked in the "ngc2903" account by ~ngc2903/.idl_startup eggidl/gen routines (see: /share/ngc2903/procs/eggidl) These can be compiled by calling "egggeninit"; the path is included in ~ngc2903/.idl_startup First, identify the blocks of observing time scheduled for the project. Go to the NAIC page, telescope schedule, searchable schedules. Enter the project number (A1963) for the relevant time period. Note that since the observing time sometimes crosses 0H AST, a block may actually have two entries. We designate blocks by the date of the START time. Also note that this utility is a little funky; you need to double check with the posted web schedule and also the one in the control room, especially for blocks that may start just before midnight. In fact, for the block this week, I am not 100% of the start times; see what I assume below. The schedule is slightly ambiguous; there appears to be some dead time just before the run, and in fact, the "searchable schedules" gives two, inconsistent entries for Feb 12 (2300-2400 or 2330-2400). Someone needs to check. Make up a simple list of the blocks that A1963 is assigned, e.g.: Block Day AST range 05.02.10 R 10Feb 2315-0215 05.02.11 F 11Feb 2315-0215 05.02.12 S 12Feb 2330-0215 05.02.13 S 13Feb 2315-0215 ***COMMENT: We are trying to start the N2903 observations at the same LST time each night. This LST is 083200 (see previous examples of schedule). So run lsttoazrt iteratively until there is a scan that starts at this time. Then just pick up the schedule when the allocated block time allows. ***COMMENT: Occasionally, there is an error message in running lsttoazrt indicating that there is not enough time to rotate. If this occurs, adjust the relevant delta T in the procedure upwards a bit. You are now ready to use the EGGIDL routines to help with the scheduling. -- On a linux machine as user ngc2903, cd to /share/ngc2903/procs/eggidl. -- Start idl > idl > @wasinit2n > @egggeninit -- The routine "blocksched" will help you determine the corresponding LST time. It works in AST, so the date changes (sure, we could improve this...) and thus you have to run it twice to get the right times. The syntax of blocksched is: SYNTAX blocksched,yymmdd,ast1HHMM,ast2HHMM EXAMPLE blocksched,050124L,0200,0615 It then prints the output to the screen. > blocksched,050210L,2315,2400 % Compiled module: YYMMDDTOJULDAY. % Compiled module: JULDAY. almanac for day 2005 02 10 JD(AST=0): 2453411.67 Computation via blocksched (eggidl:mph) at: Mon Feb 7 10:09:53 2005 % Compiled module: JULDAYTOLMST. % Compiled module: MEAN. % Compiled module: UTCINFOINP. % Compiled module: CALDAT. % Compiled module: DMTODAYNO. % Compiled module: ISLEAPYEAR. % Compiled module: FILE_EXISTS. % Compiled module: DAYNOTODM. % Compiled module: UTCTOUT1. % Compiled module: FISECMIDHMS3. JD at 0h 2453411.7 LMST at 0h: 08:54:19 23:15:00 24:00:00 32:13:08 32:58:16 08h54 23h15-24h00 08h13-08h58 And, for the other part of the block, after midnight AST, > blocksched,050211L,0000,0215 almanac for day 2005 02 11 JD(AST=0): 2453412.67 Computation via blocksched (eggidl:mph) at: Mon Feb 7 10:11:07 2005 JD at 0h 2453412.7 LMST at 0h: 08:58:16 00:00:00 02:15:00 08:58:16 11:13:38 08h58 00h00-02h15 08h58-11h13 So the run will extend from 08h13m to 11h13m LST. -- Continue running blocksched for the relevant dates so that you can fill out the schedule. Block Day AST range 05.02.10 R 10Feb 2315-0215 08h13-11h13 05.02.11 F 11Feb 2315-0215 08h17-11h17 05.02.12 S 12Feb 2330-0215 08h36-11h21 05.02.13 S 13Feb 2315-0215 08h24-11h25 -- Now you are ready to calculate the drift settings. This uses the routine lsttoazRT.pro. ***THIS ROUTINE IS NOT GUARANTEED TO WORK** and it is best if you try to understand what it does. The syntax is SYNTAX: lsttoazRT,yymmdd,lmstHHMMSS,raHHMMSS,decDDMMSS,480. EXAMPLE: lsttoazRT,041202L,083200L,092810L,211326L,480. lsttoazRT,041202L,074000L,085248L,200631L,240. lsttoazRT,041202L,091000L,101955L,215931L,240. Lyle's note indicates we need to complete the following drift strips: Declination settings for Continuum beam maps: 0805+21 1019+21 +9'24" +211615 +220855 -9'24" +205727 +215007 +11'17" +211808 +221048 -11'17" +205534 +214814 Declination settings for NGC 2903 line maps: -14'15 +211811 +19'00 +215126 +23'45 +215611 -23'45 +210841 +28'30 +220056 The calculations should begin as follows: NGC2903 083200 0805+21 065000 1019+21 090800 1. Decide which date you want to observe the target dec. strip. 2. Make sure the JD at midnight is correct for this date; for Feb 10th start use Feb 11th as the input date for the calculation. 3. Run lsttoazRT for the declination for that date. 4. Pick out the lst times that work for our observing run. 5. Write out the schedule as I have done before. Continuum scans at +000924 First: 0805+21 08:05:38.5 +21:06:51 +210651 on-source ==>> 211615 IDL> lsttoazRT,050211L,065000L,080338L,211615L,240. almanac for day 2005 02 11 JD(AST=0): 2453412.67 Computation via lsttoazRT (eggidl:mph) at: Mon Feb 7 10:33:23 2005 Source ra,dec: 80338 211615 8.0605553 20.989166 16 2453412.638551 23h19m31s 08h17m40s 129.4340 4.161 66.4 1.91 59.9 19.7 120.0 17 2453412.642360 23h25m00s 08h23m10s 120.7090 5.229 76.1 1.87 46.7 15.7 90.0 18 2453412.646169 23h30m29s 08h28m40s 115.1103 6.380 82.5 1.85 48.8 13.8 90.0 19 2453412.649978 23h35m58s 08h34m10s 111.3044 7.577 87.3 1.84 49.9 12.8 90.0 20 2453412.653787 23h41m27s 08h39m40s 108.6103 8.803 90.9 1.82 50.7 12.1 90.0 21 2453412.657596 23h46m56s 08h45m10s 106.6341 10.050 93.7 1.82 51.2 11.6 90.0 22 2453412.661405 23h52m25s 08h50m40s 105.1348 11.308 95.7 1.81 51.4 11.2 90.0 23 2453412.665214 23h57m54s 08h56m10s 103.9723 12.571 97.6 1.81 51.6 11.1 90.0 24 2453412.669023 00h03m24s 09h01m40s 103.0688 13.840 99.2 1.80 51.7 11.0 90.0 Move to second source: 1021+21 10:21:54.6 +21:59:31 +215931 on-source ==>> 220855 You may not make it in time to do the first drift. If you don't have time, skip it and go to the 2nd one. There is an error in the line below --- the date is old. It is left here for illustrative purposes. Don't forget those "L"s at the end; that indicates "long integer" in IDL and is essential. IDL> lsttoazRT,041226L,090800L,101955L,220855L,240. almanac for day 2004 12 26 JD(AST=0): 2453365.67 Computation via lsttoazRT (eggidl:mph) at: Mon Feb 7 11:57:08 2005 Source ra,dec: 101955 220855 10.331944 22.148611 Julian Date AST LMST az za RotAng Sep Slew Rotate dT 1 2453365.801739 03h14m30s 09h08m00s 253.9955 17.355 -60.2 1.80 0.0 0.0 0.0 2 2453365.805548 03h19m59s 09h13m30s 253.2318 16.096 -58.6 1.80 51.5 11.0 90.0 3 2453365.809357 03h25m28s 09h19m00s 252.2790 14.845 -56.7 1.81 51.3 11.1 90.0 4 2453365.813166 03h30m58s 09h24m30s 251.1123 13.601 -54.9 1.81 51.1 11.1 90.0 5 2453365.816975 03h36m27s 09h30m00s 249.6794 12.368 -52.9 1.82 50.8 11.2 90.0 6 2453365.820784 03h41m56s 09h35m30s 247.8823 11.146 -50.3 1.82 50.5 11.5 90.0 7 2453365.824593 03h47m25s 09h41m00s 245.5790 9.940 -47.2 1.84 50.2 11.9 90.0 8 2453365.828402 03h52m54s 09h46m30s 242.5790 8.758 -43.3 1.85 49.6 12.3 90.0 9 2453365.832211 03h58m23s 09h52m00s 238.6017 7.615 -38.5 1.87 48.6 12.8 90.0 10 2453365.836136 04h04m02s 09h57m40s 233.0097 6.498 -32.0 1.90 47.9 13.8 100.0 11 2453365.840060 04h09m41s 10h03m20s 225.1391 5.475 -23.5 1.95 45.6 15.0 100.0 12 2453365.844100 04h15m30s 10h09m10s 213.5452 4.588 -11.8 2.01 49.0 16.9 110.0 13 2453365.848486 04h21m49s 10h15m30s 195.7934 3.946 4.9 2.06 64.4 19.8 140.0 14 2453365.853449 04h28m58s 10h22m40s 171.1401 3.830 26.4 2.08 81.6 22.7 190.0 15 2453365.858297 04h35m57s 10h29m40s 149.6737 4.385 46.0 2.01 73.7 21.5 180.0 16 2453365.862453 04h41m56s 10h35m40s 136.5357 5.237 59.2 1.95 52.8 17.8 120.0 17 2453365.866262 04h47m25s 10h41m10s 128.1516 6.198 68.3 1.91 44.0 15.3 90.0 18 2453365.870071 04h52m54s 10h46m40s 122.1763 7.261 75.1 1.88 46.6 14.0 90.0 19 2453365.873880 04h58m23s 10h52m10s 117.8288 8.389 80.3 1.86 48.2 13.1 90.0 20 2453365.877689 05h03m52s 10h57m40s 114.5869 9.561 84.5 1.84 49.3 12.5 90.0 21 2453365.881498 05h09m21s 11h03m10s 112.1037 10.762 87.7 1.83 50.0 11.9 90.0 22 2453365.885307 05h14m50s 11h08m40s 110.1563 11.979 90.3 1.82 50.4 11.5 90.0 23 2453365.889116 05h20m20s 11h14m10s 108.6104 13.209 92.6 1.82 50.8 11.3 90.0 24 2453365.892925 05h25m49s 11h19m40s 107.3789 14.450 94.7 1.81 51.0 11.3 90.0 25 2453365.896734 05h31m18s 11h25m10s 106.3889 15.700 96.4 1.81 51.2 11.0 90.0 26 2453365.900543 05h36m47s 11h30m40s 105.5680 16.956 97.6 1.81 51.4 10.7 90.0 27 2453365.904352 05h42m16s 11h36m10s 104.8819 18.218 99.1 1.80 51.6 10.8 90.0 OPTANG beyond ALFA limit; consider alternatives Time runs out around 11h25 LST so you won't be able to do the last few drifts.