hp41programs

BREST

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 Z0 = 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"