Hong-Kong Water-Level for the HP-41

Overview
 

-This program computes the water-levels in Hong-Kong - China  ( Latitude 22° 17' 00"0 N, Longitude 114° 11' 00"0 E )

-I've calculated the harmonic constituents by the routine listed in "Harmonic Constituents for the free42" ( "HCST" )
 and 2000 water-levels given by the SHOM - http://www.shom.fr between March 2018 and June 2019
 

Program Listing
 
 

Data Registers:   R00 thru R15: temp
Flags: /
Subroutines: /
 
 
 

01 LBL "HKONG"  02 DEG  03 HR  04 24  05 /  06 X<>Y  07 1.012  08 DDAYS  09 -  10 STO 00  11 .985647  12 *  13 STO 03  14 64  15 +  16 COS  17 100  18 *  19 RCL 03  20 ST+ X  21 STO 05  22 166  23 +  24 COS  25 57  26 *  27 +  28 RCL 00  29 13.176396  30 *  31 STO 06  32 RCL 00  33 .111404  34 *  35 STO 08  36 -  37 STO 09  38 131  39 +  40 STO 07  41 COS  42 RCL 00  43 .052954  44 *  45 235  46 +  47 STO 04  48 RCL 07  49 X<>Y  50 ST+ 07  51 -  52 COS  53 RCL 07  54 COS  55 +  56 15  57 /  58 -  59 5  60 *  61 +  62 RCL 06  63 ST+ X  64 STO 10  65 RCL 05  66 -  67 STO 11  68 127  69 +  70 COS  71 LASTX  72 STO 07  73 RCL 04  74 ST+ 07  75 -  76 COS  77 14  78 /  79 +  80 RCL 07  81 COS  82 15  83 /  84 -  85 31  86 *  87 +  88 RCL 10  89 10  90 +  91 COS  92 LASTX  93 RCL 04  94 +  95 STO 07  96 COS  97 2.4  98 /  99 + 100 RCL 04 101 RCL 07 102 + 103 COS 104 26 105 / 106 + 107 18 108 * 109 + 110 RCL 00 111 360 112 * 113 STO 01 114 RCL 09

115 - 116 RCL 10            117 - 118 RCL 03 119 + 120 120 121 - 122 COS 123 LASTX 124 RCL 04 125 - 126 COS 127 5.3 128 / 129 + 130 54 131 * 132 + 133 RCL 01 134 RCL 10 135 - 136 RCL 03 137 + 138 11.5 139 - 140 COS 141 LASTX 142 RCL 04 143 - 144 COS 145 5.3 146 / 147 + 148 288 149 * 150 + 151 RCL 01 152 RCL 05 153 - 154 RCL 00 155 47 E-6 156 * 157 STO 00 158 + 159 97 160 - 161 COS 162 + 163 RCL 01 164 RCL 03 165 - 166 175.5 167 + 168 COS 169 108 170 * 171 + 172 RCL 01 173 171 174 + 175 COS 176 ST+ X 177 + 178 RCL 01 179 RCL 03 180 + 181 168.7 182 - 183 COS 184 LASTX 185 STO 07 186 RCL 04 187 ST+ 07 188 - 189 COS 190 50 191 / 192 - 193 RCL 07 194 COS 195 7.4 196 / 197 + 198 RCL 04 199 RCL 07 200 + 201 COS 202 343 203 / 204 - 205 358 206 * 207 + 208 RCL 01 209 RCL 03 210 + 211 RCL 09 212 + 213 47 214 - 215 COS 216 LASTX 217 STO 07 218 RCL 04 219 ST+ 07 220 - 221 COS 222 34 223 / 224 - 225 RCL 07 226 COS 227 5 228 /

229 + 230 15 231 * 232 + 233 RCL 01            234 ST+ X 235 STO 02 236 RCL 09 237 RCL 11 238 + 239 ST+ X 240 - 241 174 242 - 243 COS 244 ST+ X 245 + 246 RCL 02 247 RCL 11 248 3 249 * 250 - 251 125 252 + 253 COS 254 LASTX 255 RCL 04 256 - 257 COS 258 9 259 / 260 - 261 16 262 * 263 + 264 RCL 02 265 RCL 09 266 - 267 RCL 11 268 - 269 RCL 10 270 - 271 STO 07 272 56 273 + 274 COS 275 ST+ X 276 + 277 RCL 07 278 RCL 08 279 ST+ X 280 STO 13 281 + 282 173 283 + 284 COS 285 ST+ X 286 + 287 RCL 02 288 RCL 11 289 ST+ X 290 - 291 RCL 09 292 - 293 160 294 - 295 COS 296 4 297 * 298 + 299 RCL 02 300 RCL 10 301 - 302 RCL 11 303 - 304 STO 07 305 103 306 - 307 COS 308 + 309 RCL 07 310 RCL 13 311 + 312 129 313 - 314 COS 315 10 316 * 317 + 318 RCL 02 319 RCL 11 320 ST+ X 321 - 322 STO 14 323 80 324 - 325 COS 326 LASTX 327 RCL 04 328 - 329 COS 330 26 331 / 332 - 333 19 334 * 335 + 336 RCL 14 337 RCL 03 338 + 339 RCL 00 340 - 341 122 342 -

343 COS 344 ST+ X 345 + 346 RCL 02 347 RCL 11           348 - 349 RCL 09 350 - 351 15 352 - 353 COS 354 LASTX 355 RCL 04 356 - 357 COS 358 27 359 / 360 - 361 84 362 * 363 + 364 RCL 14 365 RCL 09 366 + 367 41 368 + 369 COS 370 LASTX 371 RCL 04 372 - 373 COS 374 26 375 / 376 - 377 11 378 * 379 + 380 RCL 02 381 RCL 10 382 - 383 STO 07 384 107 385 - 386 COS 387 LASTX 388 RCL 04 389 - 390 COS 391 6 392 / 393 + 394 10 395 * 396 + 397 RCL 07 398 RCL 13 399 + 400 COS 401 9 402 * 403 + 404 RCL 02 405 RCL 10 406 - 407 RCL 03 408 + 409 RCL 00 410 + 411 9 412 + 413 COS 414 ST+ X 415 + 416 RCL 02 417 RCL 11 418 - 419 STO 15 420 100.5 421 + 422 COS 423 LASTX 424 RCL 04 425 - 426 COS 427 27 428 / 429 - 430 403 431 * 432 + 433 RCL 15 434 RCL 03 435 + 436 RCL 00 437 - 438 166 439 - 440 COS 441 + 442 RCL 02 443 RCL 09 444 - 445 169 446 - 447 COS 448 + 449 RCL 15 450 RCL 09 451 + 452 STO 07 453 10 454 + 455 COS 456 8

457 * 458 + 459 RCL 07 460 RCL 13           461 + 462 51 463 - 464 COS 465 LASTX 466 RCL 04 467 + 468 COS 469 2 470 / 471 + 472 ST+ X 473 + 474 RCL 02 475 RCL 03 476 - 477 RCL 00 478 + 479 136 480 - 481 COS 482 12 483 * 484 + 485 RCL 02 486 58.8 487 - 488 COS 489 LASTX 490 RCL 04 491 - 492 COS 493 447 494 / 495 + 496 160 497 * 498 + 499 RCL 02 500 RCL 03 501 + 502 RCL 00 503 - 504 26 505 + 506 COS 507 ST+ X 508 + 509 RCL 02 510 RCL 05 511 + 512 144 513 + 514 COS 515 LASTX 516 STO 07 517 RCL 04 518 ST+ 07 519 - 520 COS 521 77 522 / 523 - 524 RCL 07 525 COS 526 .3 527 * 528 + 529 RCL 04 530 RCL 07 531 + 532 COS 533 31 534 / 535 + 536 43 537 * 538 + 539 RCL 02 540 RCL 11 541 + 542 65 543 - 544 COS 545 ST+ X 546 + 547 RCL 01 548 ST+ 02 549 ST+ 15 550 RCL 14 551 + 552 STO 14 553 RCL 03 554 - 555 56 556 - 557 COS 558 5 559 * 560 + 561 RCL 02 562 RCL 11 563 - 564 RCL 06 565 RCL 03 566 - 567 STO 12 568 - 569 41 570 +

571 COS 572 LASTX 573 RCL 04 574 - 575 COS 576 18 577 / 578 - 579 16 580 * 581 + 582 RCL 15           583 RCL 03 584 + 585 115 586 + 587 COS 588 LASTX 589 STO 07 590 RCL 04 591 ST+ 07 592 - 593 COS 594 17 595 / 596 - 597 RCL 07 598 COS 599 7 600 / 601 + 602 15 603 * 604 + 605 RCL 02 606 RCL 03 607 - 608 173 609 - 610 COS 611 4 612 * 613 + 614 RCL 02 615 RCL 03 616 + 617 46 618 - 619 COS 620 LASTX 621 RCL 04 622 + 623 COS 624 7 625 / 626 + 627 8 628 * 629 + 630 RCL 01 631 ST+ 02 632 ST+ 15 633 RCL 14 634 + 635 STO 14 636 RCL 11 637 ST+ X 638 - 639 72 640 - 641 COS 642 ST+ X 643 + 644 RCL 14 645 RCL 11 646 - 647 58 648 - 649 COS 650 ST+ X 651 + 652 RCL 14 653 RCL 09 654 - 655 24 656 - 657 COS 658 LASTX 659 RCL 04 660 - 661 COS 662 13 663 / 664 - 665 8 666 * 667 + 668 RCL 14 669 44 670 + 671 COS 672 LASTX 673 RCL 04 674 - 675 COS 676 13 677 / 678 - 679 25 680 * 681 + 682 RCL 02 683 RCL 11 684 - 685 RCL 09

686 - 687 130 688 + 689 COS 690 4 691 * 692 + 693 RCL 15  694 RCL 03 695 - 696 RCL 00            697 + 698 10 699 - 700 COS 701 ST+ X 702 + 703 RCL 15 704 130 705 - 706 COS 707 LASTX 708 RCL 04 709 - 710 COS 711 29 712 / 713 - 714 17 715 * 716 + 717 RCL 15 718 RCL 05 719 + 720 50 721 + 722 COS 723 LASTX 724 RCL 04 725 + 726 COS 727 .3 728 * 729 + 730 5 731 * 732 + 733 RCL 02 734 114 735 + 736 COS 737 ST+ X 738 + 739 RCL 01 740 RCL 02 741 + 742 RCL 03 743 + 744 12 745 + 746 COS 747 + 748 RCL 01 749 ST+ 01 750 ST+ X 751 ST+ 02 752 ST+ 15 753 RCL 14 754 + 755 STO 14 756 RCL 11 757 - 758 STO 15 759 RCL 09 760 - 761 162 762 + 763 COS 764 4 765 * 766 + 767 RCL 15 768 122 769 - 770 COS 771 4 772 * 773 + 774 RCL 14 775 48 776 + 777 COS 778 4 779 * 780 + 781 RCL 01 782 RCL 15 783 + 784 STO 15 785 56 786 - 787 COS 788 + 789 RCL 15 790 RCL 05 791 + 792 136 793 + 794 COS 795 + 796 .1 797 % 798 1.38 799 + 800 END

 
      ( 1041 bytes / SIZE 016 )
 
 

STACK	INPUTS	OUTPUTS
Y	Date	/
X	Time ( UT )	Water-Level (m)

 
Example:      2019 July 07th  16h41 UT

     7.072109  ENTER^
        16.41     XEQ "HKONG"              >>>>                            WL = 1.312 m
 

Notes:

-Compared to the 2000 water-levels used to calculate the harmonic constituents, the root-mean-square error is about 0.6 cm and the maximum error is 3.4 cm
-But nothing proves the errors are never larger on other intervals of time !
 
 

Reference:

[1] SHOM - http://www.shom.fr