87/12/17. 
- 21 - 
HARANB 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
PROGRAM HARANB(INPUT,OUTPUT,TAPE60,TAPE70,TAPE80,TAPE62) 
PANSCH, HARMONISCHE ANALYSE VON GEZEITEN UNO GEZE I TENSTR 0EME N 
MIT GANG-DARSTELLUNG DURCH LINEARES GANGGLIED UND NG 
TERME EINER FOURIERREIHE, 
GERUFEN VON CCL-PROCEDUR HARANB. 
UNTERPROGRAMMSATZ FUER HARMONISCHE ANALYSEN *HAMANB' WIRD GERUFEN, 
SOURCE-PROGRAMME AUF FILE *PANPROG/UN = HYPAN ‘ , 
UEBERSETZTER UNTER-PROGRAMMSATZ AUF REL-FILE *HYLIB/UN=HYLIB* 
UEBERSETZTES HAUPTPROGRAMM AUF FILE *PAN RELB/UN = HYPAN *. 
FUER DIE BEOBACHTUNGSDATEN SIND 2 STRUKTUREN MOEGLICH, 
(SIEHE SUBR. INPUT): 
1. GENORMTE DATENSAETZE FUER GLE I CHABSTAEND I GE BEOBACHTUNGEN, 
KLEINSTMOEGLICHES DELTA-T (DT) IST 1 MINUTE, 
EINLESEPROGRAMM IST SUBR. RDJR, 
BESCHREIBUNGEN DER DATENSAETZE SIEHE DORT. 
SOLCHE DATENSAETZE ENTSTEHEN Z.B. BEI DER DIGITALISIERUNG VON 
PEGEL80EGEN AM * MOP-D I G I PL AN * GERAET UEBER DAS PROGRAMMSYSTEM 
‘ J AHRPG * , * SUBR. STDOUT*. 
2. SONDERFORMAT FUER BEOBACHTUNGEN AUS SPANIEN, 
(SIEHE SUBR. INPUT) : 
UMSCHALTUNG AUF DIE RICHTIGE DATENSTRUKTUR GESCHIEHT 
MITTELS PARAMETER ’ F M‘ VON DER CCL-PROC. ‘HARANB* 
HARANB IST SO KONSTRUIERT, DASS BEOBACHTUNGEN FUER HOECHSTENS 1 
KALENDERJAHR EINGELESEN WERDEN. WENN BEOBACHTUNGSREIHEN VORLIEGEN, 
IN DENEN EIN JAHRESWECHSEL VORKOMMT, IST ES ZWECKMAESS I G , DEN EIN 
GABETEIL VON HARANB (Z.B. UEBER SUBR. INPUT) ENTSPRECHEND ZU 
AENDERN ODER SUBR. HAMANB DIREKT ZU BENUTZEN. 
PARAMETER WIE: 
PEGELNR., NAME DER EINGABEDATEI, DEREN USER-NR., ANFANGSDATUM, 
FORMATPARAMETER, SEITENZAHLEN DER J/V-TAFEL FUER DIE PARTIALTIDEN 
WERDEN UEBER DEN FILE ‘INPUT*, AM BESTEN MITTELS DER 
CCL-PROC. ‘HARANB* VOM PROGRAMM UEBERNOMMEN. 
HIERFUER ENTHAELT DAS PROGR. FOLGENDE READ-STATEMENTS, DIE VON 
DER CCL-PROC. ERFUELLT WERDEN (NFMT IST DER FORMAT-SCHAL TER): 
FILES, DESSEN USER-NR., UND TOLERANZC 
T I 0 N 
READ*, NG.ISTN, XMT , JF , UN , T0L , NF, 
READ*, ISTN, X M T,J F,U N 
READ *, NA, NE 
READ ANFANGSDATUM DER MESSUNG 
READ * ,JA,MA, ITA, ISA, IMINA 
N = N + 1 
READ 110, ENDE, IPG(N), NFMT 
FORMAT(A3,6X,I3,2 3 X,12) 
IF(ENDE.NE.TENDE) GOTO 105 
105 
110 
8EDEUTUNG EINIGER PARAMETER: 
NMAX = ANZAHL DER EINGELESENEN BEOBACHTUNGEN 
NMAX MUSS UNGERADE SEIN 
NA, NE = ANFANG UND ENDE DES AUSDRUCKS 
NUR WENN NA.GT.NE IST, GIBT ES EINE FEHLERVERTEILUNG, 
WOBEI DANN JEDOCH DER ARRAY DER BEOBACHTUNGEN ZER- 
STOERT WIRD. 
DT = ABSTAND DER MESSUNGEN IN MINUTEN 
IPG(N) = TIDEN UEBER SEITENZAHL DER J/V TAFEL IDENTIFIZIERT 
NG = ANZAHL DER GEFORDERTEN GANGTERME 
TOL, NF(RE AL) , F SIND GROESSEN, DIE FUER STEUERUNG DER ITERATION 
DER AUSGLEICHUNG BENOETIGT WERDEN, WERTE SIEHE CCL.-PROC. HARANB, 
DIE DORT ANGEGEBENEN HABEN SICH BEWAEHRT. 
WEITERE DATEN ZUR STEUERUNG DES INPUTS WERDEN GESTEUERT 
DURCH PARAMETER NFMT ( = FM IN PROCEDUR HARANB), VGL. 
SUBROUTINE INPUT. 
DIE DIMENSIONIERUNG REICHT FUER: 
61 PARTIALTIDEN, PARTIALTIDEN ♦ GANGTERME LT.70, 
18500 BEOBACHTUNGEN 
!!!!!! KEINE SICHERUNG GEGEN UEBERLAUF !!!!!!! 
HARANB 
1 
HARANB 
2 
HARAN8 
3 
HARANB 
4 
HARANB 
5 
HARANB 
6 
HARAN8 
7 
HARANB 
8 
HARAN8 
9 
HARANB 
10 
HARANB 
11 
HARANB 
12 
HARANB 
13 
HARANB 
14 
HARANB 
15 
HARANB 
16 
HARANB 
17 
HARANB 
18 
HARANB 
19 
HARANB 
20 
HARAN8 
21 
HARANB 
22 
HARANB 
23 
HARANB 
24 
HARANB 
25 
HARANB 
26 
HARANB 
27 
HARANB 
28 
HARANB 
29 
HARANB 
30 
HARAN8 
31 
HARANB 
32 
HARANB 
33 
HARAN8 
34 
HARANB 
35 
HARANB 
36 
HARANB 
37 
HARANB 
38 
EHARANB 
39 
^HARANB 
40 
HARANB 
41 
HARAN8 
42 
HARAN8 
43 
HARANB 
44 
HARANB 
45 
HARAN8 
4 6 
HARANB 
47 
HARANB 
48 
HARANB 
49 
HARAN8 
50 
HARAN8 
51 
HARANB 
52 
HARANB 
53 
HARANB 
54 
H AR AN 8 
55 
HARANB 
56 
HARANB 
57 
HARANB 
58 
HARANB 
59 
HARANB 
60 
HARANB 
61 
HARANB 
62 
HARANB 
63 
HARANB 
64 
HARANB 
65 
HARAN8 
66 
HARANB 
67 
HARANB 
68 
HARANB 
69 
HARANB 
70 
HARANB 
71 
HARANB 
72 
HARANB 
73