Nutation for the HP-41
Overview
-"NUT" calculates the nutation in longitude and the nutation in obliquity
( IAU 2000_R06 expressions )
-The terms smaller than 0"00017 are neglected.
Program Listing
Data Registers:
• R00 = T = time from J2000.0 ( 2000/01/01 12h TT )
( Register R00 is to be initialized before executing "NUT" )
unit = 1000 years = 365250 days
R01 = Dpsi = Nutation in longitude ( in decimal degree )
R03-R04-R05-R11-R12: temp
R02 = Deps = Nutation in obliquity ( in decimal degree )
R06 to R10 = Delaunay arguments ( l l' F D Omega )
Flags: /
Subroutines: /
| 01 LBL "NUT"
02 DEG 03 RCL 00 04 .886 05 RCL 00 06 70 07 / 08 + 09 * 10 4771988.676 11 + 12 * 13 134.963 14 + 15 STO 06 16 CLX 17 65 18 / 19 CHS 20 359990.503 21 + 22 * 23 2.471 24 - 25 STO 07 26 CLX 27 .354 28 * 29 CHS 30 4832020.175 31 + 32 * 33 93.272 34 + 35 STO 08 36 CLX 37 546 38 / 39 .177 40 - 41 * 42 4452671.114 43 + 44 * 45 62.15 46 - 47 STO 09 48 CLX 49 467 50 / 51 .208 52 + 53 * 54 19341.363 55 - 56 * 57 125.045 58 + 59 STO 10 60 1 61 STO 12 62 P-R 63 STO 11 64 883 65 RCL 00 66 * 67 920523 68 + 69 * 70 X<>Y 71 STO 03 72 154 73 * 74 + 75 STO 02 76 RCL 11 77 334 78 * 79 RCL 00 80 17419 81 * 82 1720642 83 + 84 RCL 03 85 * 86 - 87 STO 01 88 RCL 08 89 RCL 09 90 - 91 RCL 10 92 + 93 STO 13 94 ST+ X 95 STO 05 96 RCL 12 97 P-R 98 STO 11 99 303 100 RCL 00 101 * 102 57303 103 - 104 * 105 X<>Y 106 STO 03 107 46 108 * 109 + 110 ST- 02 111 RCL 03 112 164 113 RCL 00 114 * 115 131709 116 + 117 * 118 RCL 11 119 137 120 * 121 + 122 ST- 01 123 RCL 07 124 RCL 12 125 P-R 126 STO 11 127 RCL 00 128 18 129 * 130 739 131 - 132 * 133 X<>Y |
134 STO 03
135 19 136 * 137 + 138 ST- 02 139 RCL 03 140 364 141 RCL 00 142 * 143 14759 144 - 145 * 146 RCL 11 147 118 148 * 149 - 150 ST- 01 151 RCL 08 152 RCL 10 153 + 154 ST+ X 155 STO 04 156 RCL 12 157 P-R 158 STO 11 159 RCL 00 160 49 161 * 162 9785 163 - 164 * 165 ST- 02 166 CLX 167 RCL 00 168 23 169 * 170 22764 171 + 172 * 173 RCL 11 174 28 175 * 176 - 177 ST- 01 178 RCL 10 179 ST+ X 180 RCL 12 181 P-R 182 RCL 00 183 47 184 * 185 8975 186 - 187 * 188 ST+ 02 189 CLX 190 RCL 00 191 20 192 * 193 20746 194 + 195 * 196 ST+ 01 197 RCL 07 198 RCL 05 199 - 200 RCL 12 201 P-R 202 RCL 00 203 30 204 * 205 959 206 - 207 * 208 ST+ 02 209 CLX 210 RCL 00 211 49 212 * 213 2158 214 - 215 * 216 ST+ 01 217 RCL 05 218 RCL 07 219 + 220 RCL 12 221 P-R 222 68 223 RCL 00 224 * 225 2244 226 - 227 * 228 ST- 02 229 CLX 230 123 231 RCL 00 232 * 233 5168 234 - 235 * 236 RCL 04 237 RCL 10 238 - 239 RCL 12 240 P-R 241 2007 242 * 243 ST+ 02 244 CLX 245 RCL 00 246 37 247 * 248 3873 249 + 250 * 251 - 252 RCL 04 253 RCL 06 254 + 255 RCL 12 256 P-R 257 1290 258 * 259 ST+ 02 260 CLX 261 3015 262 * 263 - 264 RCL 06 265 RCL 12 266 P-R |
267 68
268 * 269 ST- 02 270 CLX 271 7112 272 * 273 + 274 RCL 06 275 RCL 04 276 - 277 RCL 12 278 P-R 279 533 280 * 281 ST- 02 282 CLX 283 1235 284 * 285 - 286 RCL 06 287 RCL 10 288 + 289 RCL 12 290 P-R 291 332 292 * 293 ST- 02 294 CLX 295 631 296 * 297 + 298 RCL 05 299 RCL 10 300 - 301 RCL 12 302 P-R 303 690 304 * 305 ST- 02 306 CLX 307 1282 308 * 309 + 310 RCL 06 311 RCL 09 312 ST+ X 313 STO 03 314 - 315 SIN 316 1570 317 * 318 - 319 RCL 03 320 SIN 321 634 322 * 323 + 324 RCL 06 325 RCL 04 326 - 327 RCL 03 328 - 329 RCL 12 330 P-R 331 255 332 * 333 ST+ 02 334 CLX 335 596 336 * 337 + 338 RCL 06 339 RCL 10 340 - 341 RCL 12 342 P-R 343 314 344 * 345 ST+ 02 346 CLX 347 580 348 * 349 + 350 RCL 04 351 RCL 06 352 + 353 RCL 10 354 - 355 RCL 12 356 P-R 357 264 358 * 359 ST+ 02 360 CLX 361 516 362 * 363 - 364 RCL 06 365 RCL 09 366 - 367 ST+ X 368 SIN 369 477 370 * 371 + 372 RCL 06 373 RCL 08 374 - 375 ST+ X 376 STO 14 377 RCL 10 378 - 379 RCL 12 380 P-R 381 242 382 * 383 ST- 02 384 CLX 385 459 386 * 387 - 388 RCL 03 389 RCL 04 390 + 391 RCL 12 392 P-R 393 164 394 * 395 ST+ 02 396 CLX 397 386 398 * 399 - |
400 RCL 07
401 ST+ X 402 RCL 05 403 - 404 RCL 12 405 P-R 406 139 407 * 408 ST- 02 409 CLX 410 325 411 * 412 - 413 RCL 04 414 RCL 06 415 ST+ X 416 + 417 RCL 12 418 P-R 419 132 420 * 421 ST+ 02 422 CLX 423 310 424 * 425 - 426 RCL 05 427 RCL 06 428 + 429 RCL 12 430 P-R 431 123 432 * 433 ST- 02 434 CLX 435 286 436 * 437 + 438 RCL 06 439 ST+ X 440 SIN 441 292 442 * 443 + 444 RCL 08 445 ST+ X 446 SIN 447 259 448 * 449 + 450 RCL 08 451 RCL 09 452 - 453 ST+ X 454 SIN 455 218 456 * 457 - 458 RCL 06 459 RCL 04 460 - 461 RCL 10 462 + 463 RCL 12 464 P-R 465 108 466 * 467 ST- 02 468 CLX 469 204 470 * 471 - 472 RCL 07 473 ST+ X 474 SIN 475 167 476 * 477 + 478 RCL 05 479 RCL 07 480 ST+ X 481 + 482 RCL 12 483 P-R 484 69 485 * 486 ST+ 02 487 CLX 488 158 489 * 490 - 491 RCL 06 492 RCL 10 493 + 494 RCL 03 495 - 496 RCL 12 497 P-R 498 70 499 * 500 ST+ 02 501 CLX 502 129 503 * 504 - 505 RCL 07 506 RCL 10 507 - 508 RCL 12 509 P-R 510 64 511 * 512 ST+ 02 513 CLX 514 127 515 * 516 + 517 RCL 06 518 RCL 03 519 - 520 RCL 04 521 - 522 RCL 10 523 + 524 RCL 12 525 P-R 526 52 527 * 528 ST+ 02 529 CLX 530 102 531 * 532 + |
533 RCL 07
534 RCL 10 535 + 536 RCL 12 537 P-R 538 86 539 * 540 ST+ 02 541 CLX 542 141 543 * 544 - 545 RCL 06 546 RCL 10 547 - 548 RCL 03 549 - 550 RCL 12 551 P-R 552 80 553 * 554 ST- 02 555 CLX 556 152 557 * 558 - 559 RCL 14 560 SIN 561 110 562 * 563 + 564 RCL 03 565 RCL 04 566 + 567 RCL 06 568 + 569 RCL 12 570 P-R 571 33 572 * 573 ST+ 02 574 CLX 575 77 576 * 577 - 578 RCL 04 579 RCL 07 580 + 581 RCL 12 582 P-R 583 33 584 * 585 ST- 02 586 CLX 587 76 588 * 589 + 590 RCL 06 591 RCL 07 592 + 593 RCL 03 594 - 595 SIN 596 74 597 * 598 - 599 RCL 07 600 RCL 04 601 - 602 RCL 12 603 P-R 604 31 605 * 606 ST+ 02 607 CLX 608 71 609 * 610 + 611 RCL 03 612 RCL 04 613 + 614 RCL 10 615 - 616 RCL 12 617 P-R 618 34 619 * 620 ST+ 02 621 CLX 622 66 623 * 624 - 625 RCL 03 626 RCL 06 627 + 628 SIN 629 66 630 * 631 + 632 RCL 05 633 RCL 06 634 ST+ X 635 + 636 SIN 637 64 638 * 639 + 640 RCL 03 641 RCL 10 642 + 643 RCL 12 644 P-R 645 33 646 * 647 ST+ 02 648 CLX 649 63 650 * 651 - 652 RCL 04 653 RCL 06 654 + 655 RCL 03 656 - 657 RCL 10 658 - 659 RCL 12 660 P-R 661 30 662 * 663 ST- 02 664 CLX 665 58 |
666 *
667 + 668 RCL 06 669 RCL 09 670 - 671 ST+ X 672 RCL 10 673 - 674 RCL 12 675 P-R 676 30 677 * 678 ST+ 02 679 CLX 680 58 681 * 682 + 683 RCL 06 684 RCL 08 685 + 686 ST+ X 687 RCL 10 688 + 689 RCL 12 690 P-R 691 27 692 * 693 ST+ 02 694 CLX 695 53 696 * 697 - 698 81 699 RCL 00 700 359586 701 * 702 - 703 RCL 13 704 + 705 STO 11 706 31 707 + 708 SIN 709 32 710 * 711 ST+ 02 712 X<> 11 713 59 714 - 715 SIN 716 60 717 * 718 - 719 RCL 03 720 RCL 10 721 - 722 RCL 12 723 P-R 724 27 725 * 726 ST+ 02 727 CLX 728 49 729 * 730 + 731 RCL 07 732 RCL 05 733 - 734 RCL 10 735 + 736 RCL 12 737 P-R 738 27 739 * 740 ST+ 02 741 CLX 742 48 743 * 744 + 745 RCL 06 746 RCL 07 747 - 748 SIN 749 47 750 * 751 + 752 RCL 07 753 RCL 03 754 - 755 SIN 756 43 757 * 758 - 759 RCL 09 760 SIN 761 42 762 * 763 - 764 RCL 06 765 RCL 09 766 - 767 ST+ X 768 RCL 10 769 + 770 RCL 12 771 P-R 772 22 773 * 774 ST- 02 775 CLX 776 41 777 * 778 + 779 RCL 06 780 RCL 08 781 ST+ X 782 - 783 SIN 784 41 785 * 786 + 787 RCL 06 788 RCL 09 789 - 790 SIN 791 40 792 * 793 - 794 RCL 05 795 RCL 07 796 + 797 RCL 10 798 - |
799 RCL 12
800 P-R 801 19 802 * 803 ST- 02 804 CLX 805 36 806 * 807 + 808 RCL 06 809 RCL 07 810 + 811 SIN 812 34 813 * 814 - 815 RCL 06 816 RCL 08 817 ST+ X 818 + 819 SIN 820 RCL 06 821 RCL 07 822 - 823 RCL 09 824 - 825 SIN 826 - 827 33 828 * 829 + 830 RCL 06 831 ST+ X 832 RCL 04 833 - 834 SIN 835 31 836 * 837 + 838 RCL 06 839 3 840 * 841 RCL 04 842 + 843 SIN 844 RCL 06 845 RCL 07 846 - 847 RCL 04 848 + 849 SIN 850 + 851 29 852 * 853 - 854 RCL 07 855 RCL 04 856 - 857 RCL 03 858 - 859 STO 11 860 RCL 06 861 + 862 SIN 863 28 864 * 865 + 866 RCL 11 867 SIN 868 26 869 * 870 + 871 RCL 06 872 RCL 07 873 + 874 RCL 04 875 + 876 SIN 877 25 878 * 879 + 880 RCL 06 881 ST+ X 882 RCL 10 883 - 884 SIN 885 23 886 * 887 + 888 RCL 06 889 ST+ X 890 RCL 10 891 + 892 SIN 893 22 894 * 895 + 896 44 897 RCL 00 898 44462 899 * 900 - 901 SIN 902 21 903 * 904 + 905 RCL 06 906 RCL 05 907 - 908 RCL 10 909 + 910 SIN 911 RCL 07 912 RCL 13 913 - 914 COS 915 - 916 RCL 06 917 RCL 10 918 ST+ X 919 + 920 SIN 921 - 922 20 923 * 924 + 925 ST+ 01 926 36 E7 927 ST/ 01 928 ST/ 02 929 RCL 02 930 RCL 01 931 END |
( 1348 bytes / SIZE 015 )
| STACK | INPUTS | OUTPUTS |
| Y | / | Deps |
| X | / | Dpsi |
where Dpsi = nutation in longitude & Deps = nutation in obliquity, expressed in decimal degrees.
Example: T = 1.2 i-e 3200/01/10 12h TT
1.2 STO 00 XEQ "NUT"
>>>> Dpsi = 0°003626735 = R01
---Execution time = 85s---
X<>Y Deps = 0°001710644 = R02
or, in arcseconds:
Dpsi = 13"056
Deps = 6"158
Notes:
-The precision is of the order of one milliarcsecond over the interval
[1000,3000]
-If you prefer reckoning the time from 2000/01/01 0h TT, replace lines
13-23-33-45-57-698 by
128.431 2.964 86.657 68.245 125.071 82 respectively.
-In this case, store 1.2 + 0.5 / 365250 into R00 for
the example above.
References:
[1] Table 5.3a: Nutation in longitude (IAU 2000_R06 expression)
derived from the IAU 2000A
lunisolar and planetary components
with slight IAU 2006 adjustments
[2] Table 5.3b: Nutation in obliquity (IAU 2000_R06 expression)
derived from the IAU 2000A
lunisolar and planetary components
with slight IAU 2006 adjustments