Brest Water-Level for the HP-41

Overview
 

-This program computes the water-levels in Brest -France ( Latitude = 48° 23' 00"0' N, Longitude = 4° 30' 00"0' W )

-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

-The mean-sea-level given by "HCST" with these 2000 data is 4.135 m
-But it doesn't seem constant in Brest, so I've also used reference [2]

-Lines 19 to 22 and lines 1376-1377 gives Z₀ = 4098 + d / 184  millimeters where d is the number of days since 2000/01/01
-Change these lines if you want to employ other constants...
 

Program Listing
 
 

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

01 LBL "BREST"  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 35  15 +  16 COS  17 49  18 *  19 RCL 00  20 184  21 /  22 +  23 RCL 03  24 ST+ X  25 STO 05  26 104  27 +  28 COS  29 20  30 *  31 +  32 RCL 00  33 13.176396  34 *  35 STO 06  36 ST+ X  37 STO 12  38 116  39 -  40 STO 07  41 COS  42 RCL 00  43 .052954  44 *  45 235  46 +  47 STO 04  48 RCL 07  49 +  50 COS  51 .4  52 *  53 +  54 10  55 *  56 +  57 RCL 00  58 .111404  59 *  60 STO 08  61 RCL 12  62 -  63 ST+ X  64 RCL 03  65 +  66 RCL 00  67 360  68 *  69 STO 01  70 +  71 5  72 +  73 COS  74 LASTX  75 RCL 04  76 -  77 COS  78 5  79 /  80 +  81 4  82 *  83 +  84 RCL 01  85 RCL 06  86 ST- Y  87 RCL 03  88 -  89 STO 09  90 3  91 *  92 -  93 17  94 +  95 COS  96 LASTX  97 RCL 04  98 -  99 COS 100 5 101 / 102 + 103 4 104 * 105 + 106 RCL 01 107 RCL 06 108 - 109 RCL 12 110 - 111 RCL 03 112 + 113 RCL 08 114 + 115 76 116 + 117 COS 118 LASTX 119 RCL 04 120 - 121 COS 122 5.3 123 / 124 + 125 20 126 * 127 + 128 RCL 01 129 RCL 09 130 ST- Y 131 ST+ X 132 STO 10 133 - 134 RCL 08 135 - 136 104 137 + 138 COS 139 LASTX 140 RCL 04 141 - 142 COS 143 5 144 / 145 + 146 PI 147 * 148 + 149 RCL 01 150 RCL 03 151 + 152 RCL 12 153 - 154 STO 07 155 158 156 + 157 COS 158 LASTX 159 RCL 04 160 - 161 COS 162 5.3 163 / 164 + 165 66 166 * 167 + 168 RCL 07 169 RCL 05 170 + 171 36 172 + 173 COS 174 ST+ X 175 + 176 RCL 01 177 RCL 03 178 + 179 RCL 06 180 RCL 08 181 - 182 STO 11 183 - 184 146 185 - 186 COS 187 LASTX 188 RCL 04 189 + 190 COS 191 5 192 / 193 + 194 3 195 * 196 + 197 RCL 00

198 47 E-6 199 * 200 STO 00          201 RCL 05 202 - 203 RCL 01 204 + 205 46 206 - 207 COS 208 ST+ X 209 + 210 RCL 01 211 RCL 03 212 - 213 74 214 - 215 COS 216 22 217 * 218 + 219 RCL 01 220 163 221 + 222 COS 223 9 224 * 225 + 226 RCL 01 227 RCL 03 228 + 229 65 230 - 231 COS 232 LASTX 233 STO 07 234 RCL 04 235 ST+ 07 236 - 237 COS 238 50 239 / 240 - 241 RCL 07 242 COS 243 7.4 244 / 245 + 246 64 247 * 248 + 249 RCL 01 250 RCL 05 251 + 252 RCL 00 253 - 254 51 255 - 256 COS 257 ST+ X 258 + 259 RCL 01 260 RCL 03 261 + 262 RCL 05 263 + 264 92 265 + 266 COS 267 + 268 RCL 01 269 RCL 03 270 + 271 RCL 11 272 + 273 19 274 + 275 COS 276 3 277 * 278 + 279 RCL 01 280 RCL 12 281 + 282 RCL 03 283 + 284 144 285 - 286 COS 287 + 288 RCL 01 289 ST+ X 290 STO 02 291 RCL 10 292 3 293 * 294 - 295 STO 13 296 RCL 11 297 - 298 4 299 + 300 COS 301 ST+ X 302 + 303 RCL 13 304 RCL 10 305 + 306 STO 14 307 RCL 11 308 ST+ X 309 STO 15 310 - 311 71 312 - 313 COS 314 PI 315 * 316 + 317 RCL 13 318 125 319 + 320 COS 321 LASTX 322 RCL 04 323 - 324 COS 325 9 326 / 327 - 328 4 329 * 330 + 331 RCL 14 332 RCL 05 333 - 334 RCL 11 335 - 336 STO 07 337 145 338 + 339 COS 340 LASTX 341 RCL 04 342 - 343 COS 344 .4 345 * 346 + 347 5 348 * 349 + 350 RCL 07 351 RCL 08 352 ST+ X 353 + 354 60 355 + 356 COS 357 6 358 * 359 + 360 RCL 14 361 RCL 11 362 - 363 55 364 + 365 COS 366 LASTX 367 RCL 04 368 - 369 COS 370 13 371 / 372 - 373 19 374 * 375 + 376 RCL 13 377 RCL 11 378 + 379 49 380 + 381 COS 382 LASTX 383 RCL 04 384 - 385 COS 386 13 387 / 388 - 389 9 390 * 391 + 392 RCL 14 393 RCL 05 394 -

395 88 396 - 397 COS 398 LASTX 399 RCL 04          400 - 401 COS 402 5 403 / 404 + 405 11 406 * 407 + 408 RCL 02 409 RCL 10 410 - 411 RCL 15 412 - 413 167 414 + 415 COS 416 LASTX 417 RCL 04 418 - 419 COS 420 27 421 / 422 - 423 56 424 * 425 + 426 RCL 14 427 168 428 + 429 COS 430 LASTX 431 RCL 04 432 - 433 COS 434 26 435 / 436 - 437 85 438 * 439 + 440 RCL 14 441 RCL 03 442 + 443 RCL 00 444 - 445 174 446 + 447 COS 448 5 449 * 450 + 451 RCL 02 452 RCL 10 453 - 454 RCL 11 455 - 456 82.5 457 - 458 COS 459 LASTX 460 RCL 04 461 - 462 COS 463 27 464 / 465 - 466 416 467 * 468 + 469 RCL 14 470 RCL 11 471 + 472 45 473 - 474 COS 475 LASTX 476 RCL 04 477 - 478 COS 479 26 480 / 481 - 482 78 483 * 484 + 485 RCL 14 486 RCL 11 487 + 488 RCL 05 489 + 490 COS 491 + 492 RCL 02 493 RCL 12 494 - 495 117 496 + 497 COS 498 LASTX 499 RCL 04 500 - 501 COS 502 6 503 / 504 + 505 9 506 * 507 + 508 RCL 02 509 RCL 15 510 - 511 2 512 + 513 COS 514 7 515 * 516 + 517 RCL 02 518 RCL 12 519 - 520 RCL 03 521 + 522 RCL 00 523 + 524 2 525 - 526 COS 527 11 528 * 529 + 530 RCL 02 531 RCL 10 532 - 533 STO 07 534 27.64 535 + 536 COS 537 LASTX 538 RCL 04 539 - 540 COS 541 27 542 / 543 - 544 2049 545 * 546 + 547 RCL 07 548 RCL 03 549 + 550 RCL 00 551 - 552 26 553 - 554 COS 555 16 556 * 557 + 558 RCL 14 559 RCL 12 560 + 561 144 562 + 563 COS 564 7 565 * 566 + 567 RCL 02 568 RCL 11 569 - 570 20 571 - 572 COS 573 LASTX 574 RCL 04 575 - 576 COS 577 22 578 / 579 - 580 26 581 * 582 + 583 RCL 02 584 RCL 10 585 - 586 RCL 11 587 + 588 STO 07 589 12 590 - 591 COS

592 LASTX 593 RCL 04 594 - 595 COS 596 27 597 / 598 - 599 62 600 * 601 + 602 RCL 07          603 RCL 08 604 ST+ X 605 + 606 71 607 - 608 COS 609 LASTX 610 RCL 04 611 + 612 STO 07 613 COS 614 .44 615 * 616 + 617 RCL 07 618 RCL 04 619 + 620 COS 621 16 622 / 623 + 624 15 625 * 626 + 627 RCL 02 628 RCL 03 629 - 630 RCL 00 631 + 632 135 633 - 634 COS 635 42 636 * 637 + 638 RCL 02 639 148.1 640 - 641 COS 642 LASTX 643 RCL 04 644 - 645 COS 646 447 647 / 648 + 649 744 650 * 651 + 652 RCL 02 653 RCL 03 654 + 655 RCL 00 656 - 657 19 658 + 659 COS 660 5 661 * 662 + 663 RCL 02 664 RCL 05 665 + 666 54.7 667 + 668 COS 669 LASTX 670 STO 07 671 RCL 04 672 ST+ 07 673 - 674 COS 675 77 676 / 677 - 678 RCL 07 679 COS 680 .3 681 * 682 + 683 RCL 04 684 RCL 07 685 + 686 COS 687 31 688 / 689 + 690 211 691 * 692 + 693 RCL 02 694 RCL 11 695 + 696 162 697 - 698 COS 699 LASTX 700 STO 07 701 RCL 04 702 ST+ 07 703 - 704 COS 705 RCL 07 706 COS 707 + 708 27 709 / 710 - 711 14 712 * 713 + 714 RCL 02 715 RCL 11 716 + 717 RCL 05 718 + 719 147 720 + 721 COS 722 LASTX 723 RCL 04 724 + 725 COS 726 3 727 / 728 + 729 10 730 * 731 + 732 RCL 02 733 RCL 10 734 + 735 89 736 - 737 COS 738 LASTX 739 RCL 04 740 + 741 COS 742 27 743 / 744 - 745 16 746 * 747 + 748 RCL 02 749 RCL 12 750 + 751 121 752 + 753 COS 754 LASTX 755 RCL 04 756 + 757 COS 758 4 759 / 760 + 761 8 762 * 763 + 764 RCL 01 765 ST+ 02 766 ST+ 13 767 RCL 14 768 + 769 STO 14 770 RCL 03 771 - 772 121 773 - 774 COS 775 + 776 RCL 14 777 RCL 09 778 + 779 13 780 + 781 COS 782 LASTX 783 RCL 04 784 - 785 COS 786 18 787 / 788 -

789 20 790 * 791 + 792 RCL 14 793 RCL 10 794 + 795 RCL 03 796 + 797 STO 07          798 136 799 - 800 COS 801 3 802 * 803 + 804 RCL 07 805 RCL 06 806 + 807 172 808 + 809 COS 810 LASTX 811 RCL 04 812 + 813 COS 814 .4 815 * 816 + 817 4 818 * 819 + 820 RCL 02 821 RCL 03 822 - 823 97 824 + 825 COS 826 5 827 * 828 + 829 RCL 02 830 50 831 + 832 COS 833 PI 834 * 835 + 836 RCL 02 837 RCL 03 838 + 839 62 840 - 841 COS 842 LASTX 843 RCL 04 844 + 845 COS 846 7 847 / 848 + 849 6 850 * 851 + 852 RCL 01 853 ST+ 02 854 ST+ 14 855 RCL 13 856 + 857 STO 13 858 RCL 11 859 - 860 120 861 + 862 COS 863 ST+ X 864 + 865 RCL 14 866 RCL 15 867 - 868 5 869 + 870 COS 871 3 872 * 873 + 874 RCL 13 875 162 876 - 877 COS 878 3 879 * 880 + 881 RCL 14 882 RCL 11 883 - 884 85 885 + 886 COS 887 LASTX 888 RCL 04 889 - 890 COS 891 13 892 / 893 - 894 19 895 * 896 + 897 RCL 13 898 RCL 11 899 + 900 137 901 + 902 COS 903 5 904 * 905 + 906 RCL 14 907 RCL 05 908 - 909 146 910 + 911 COS 912 + 913 RCL 14 914 168 915 + 916 COS 917 LASTX 918 RCL 04 919 - 920 COS 921 13 922 / 923 - 924 54 925 * 926 + 927 RCL 02 928 RCL 10 929 - 930 RCL 11 931 - 932 126 933 - 934 COS 935 4 936 * 937 + 938 RCL 14 939 RCL 11 940 + 941 93 942 + 943 COS 944 4 945 * 946 + 947 RCL 14 948 RCL 06 949 + 950 RCL 08 951 + 952 62 953 + 954 COS 955 ST+ X 956 + 957 RCL 02 958 RCL 12 959 - 960 RCL 03 961 + 962 RCL 00 963 + 964 38 965 - 966 COS 967 ST+ X 968 + 969 RCL 02 970 RCL 10 971 - 972 43 973 - 974 COS 975 LASTX 976 RCL 04 977 - 978 COS 979 29 980 / 981 - 982 32 983 * 984 + 985 RCL 14

986 RCL 12          987 + 988 157 989 + 990 COS 991 LASTX 992 RCL 04 993 + 994 COS 995 .3 996 * 997 + 998 10 999 * 1000 + 1001 RCL 02 1002 RCL 11 1003 - 1004 163 1005 - 1006 COS 1007 + 1008 RCL 02 1009 RCL 10 1010 - 1011 RCL 11 1012 + 1013 132 1014 - 1015 COS 1016 4 1017 * 1018 + 1019 RCL 02 1020 120 1021 + 1022 COS 1023 + 1024 RCL 02 1025 RCL 05 1026 + 1027 85 1028 - 1029 COS 1030 + 1031 RCL 01 1032 ST+ X 1033 STO 01 1034 ST+ 02 1035 ST+ 14 1036 RCL 13 1037 + 1038 STO 13 1039 RCL 12 1040 - 1041 RCL 11 1042 - 1043 STO 07 1044 131 1045 + 1046 COS 1047 ST+ X 1048 + 1049 RCL 07 1050 RCL 05 1051 + 1052 STO 07 1053 24 1054 - 1055 COS 1056 LASTX 1057 RCL 04 1058 - 1059 COS 1060 7 1061 / 1062 - 1063 5 1064 * 1065 + 1066 RCL 07 1067 RCL 05 1068 + 1069 RCL 08 1070 ST+ X 1071 - 1072 20 1073 + 1074 COS 1075 ST+ X 1076 + 1077 RCL 13 1078 RCL 12 1079 - 1080 STO 07 1081 126 1082 - 1083 COS 1084 ST+ X 1085 + 1086 RCL 07 1087 RCL 08 1088 ST+ X 1089 + 1090 148 1091 - 1092 COS 1093 6 1094 * 1095 + 1096 RCL 13 1097 RCL 10 1098 - 1099 84 1100 + 1101 COS 1102 6 1103 * 1104 + 1105 RCL 13 1106 RCL 11 1107 - 1108 45 1109 - 1110 COS 1111 LASTX 1112 RCL 04 1113 - 1114 COS 1115 9 1116 / 1117 - 1118 19 1119 * 1120 + 1121 RCL 13 1122 RCL 10 1123 - 1124 RCL 11 1125 + 1126 6 1127 - 1128 COS 1129 LASTX 1130 RCL 04 1131 - 1132 COS 1133 9 1134 / 1135 - 1136 7 1137 * 1138 + 1139 RCL 14 1140 RCL 12 1141 - 1142 24 1143 + 1144 COS 1145 3 1146 * 1147 + 1148 RCL 13 1149 55 1150 + 1151 COS 1152 LASTX 1153 RCL 04 1154 - 1155 COS 1156 9 1157 / 1158 - 1159 30 1160 * 1161 + 1162 RCL 13 1163 RCL 05 1164 + 1165 80 1166 + 1167 COS 1168 ST+ X 1169 + 1170 RCL 14 1171 RCL 11 1172 - 1173 111 1174 + 1175 COS 1176 4 1177 * 1178 + 1179 RCL 13 1180 RCL 11 1181 + 1182 11

1183 - 1184 COS 1185 LASTX 1186 RCL 04 1187 - 1188 COS 1189 7 1190 / 1191 - 1192 5 1193 * 1194 + 1195 RCL 13 1196 RCL 06          1197 + 1198 RCL 08 1199 + 1200 42 1201 - 1202 COS 1203 LASTX 1204 RCL 04 1205 + 1206 COS 1207 .3 1208 * 1209 + 1210 ST+ X 1211 + 1212 RCL 14 1213 RCL 03 1214 - 1215 RCL 00 1216 + 1217 133 1218 - 1219 COS 1220 ST+ X 1221 + 1222 RCL 14 1223 140 1224 - 1225 COS 1226 LASTX 1227 RCL 04 1228 - 1229 COS 1230 14 1231 / 1232 - 1233 17 1234 * 1235 + 1236 RCL 14 1237 RCL 05 1238 + 1239 57 1240 + 1241 COS 1242 LASTX 1243 RCL 04 1244 + 1245 COS 1246 .3 1247 * 1248 + 1249 4 1250 * 1251 + 1252 RCL 02 1253 RCL 10 1254 - 1255 RCL 11 1256 - 1257 83 1258 + 1259 COS 1260 + 1261 RCL 14 1262 RCL 11 1263 + 1264 131 1265 + 1266 COS 1267 4 1268 * 1269 + 1270 RCL 02 1271 RCL 10 1272 - 1273 141 1274 + 1275 COS 1276 3 1277 * 1278 + 1279 RCL 14 1280 RCL 12 1281 + 1282 15 1283 + 1284 COS 1285 + 1286 RCL 01 1287 ST+ 02 1288 ST+ 14 1289 RCL 13 1290 + 1291 STO 13 1292 RCL 10 1293 - 1294 STO 15 1295 RCL 11 1296 - 1297 119 1298 - 1299 COS 1300 + 1301 RCL 15 1302 37 1303 - 1304 COS 1305 ST+ X 1306 + 1307 RCL 13 1308 RCL 11 1309 - 1310 COS 1311 + 1312 RCL 13 1313 93 1314 + 1315 COS 1316 4 1317 * 1318 + 1319 RCL 13 1320 RCL 11 1321 + 1322 COS 1323 + 1324 RCL 14 1325 129 1326 - 1327 COS 1328 ST+ X 1329 + 1330 RCL 01 1331 ST+ 13 1332 RCL 15 1333 + 1334 STO 15 1335 RCL 10 1336 - 1337 STO 14 1338 RCL 11 1339 - 1340 COS 1341 + 1342 RCL 14 1343 111 1344 + 1345 COS 1346 ST+ X 1347 + 1348 RCL 15 1349 RCL 11 1350 - 1351 176 1352 - 1353 COS 1354 ST+ X 1355 + 1356 RCL 15 1357 72 1358 - 1359 COS 1360 PI 1361 * 1362 + 1363 RCL 15 1364 RCL 11 1365 + 1366 25 1367 + 1368 COS 1369 + 1370 RCL 13 1371 106 1372 + 1373 COS 1374 ST+ X 1375 + 1376 4098 1377 + 1378 .1 1379 % 1380 END

 
      ( 1704 bytes / SIZE 016 )
 
 

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

 
Example:      2019 July 16th  16h41 UT

     16.072109  ENTER^       ( or 7.162019 in MDY format )
        16.41       XEQ "BREST"              >>>>                            WL = 6.317 m
 

Notes:

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

References:

[1] SHOM - http://www.shom.fr
[2] Guy Wöppelmann - THÈSE DE DOCTORAT DE L’UNIVERSITE DE LA ROCHELLE - "Trois cents ans de mesures marégraphiques en France :
     outils, méthodes et tendances des composantes du niveau de la mer au port de Brest"