| 3 |
michaesp |
1 |
c ********************************************************************************
|
|
|
2 |
c * Transformation routine: LMSTOLM and PHSTOPH from library gm2em *
|
|
|
3 |
c ********************************************************************************
|
|
|
4 |
|
|
|
5 |
REAL FUNCTION LMSTOLM (PHIS, LAMS, POLPHI, POLLAM)
|
|
|
6 |
C
|
|
|
7 |
C**** LMSTOLM - FC:BERECHNUNG DER WAHREN GEOGRAPHISCHEN LAENGE FUER
|
|
|
8 |
C**** EINEN PUNKT MIT DEN KOORDINATEN (PHIS, LAMS)
|
|
|
9 |
C**** IM ROTIERTEN SYSTEM. DER NORDPOL DES SYSTEMS HAT
|
|
|
10 |
C**** DIE WAHREN KOORDINATEN (POLPHI, POLLAM)
|
|
|
11 |
C** AUFRUF : LAM = LMSTOLM (PHIS, LAMS, POLPHI, POLLAM)
|
|
|
12 |
C** ENTRIES : KEINE
|
|
|
13 |
C** ZWECK : BERECHNUNG DER WAHREN GEOGRAPHISCHEN LAENGE FUER
|
|
|
14 |
C** EINEN PUNKT MIT DEN KOORDINATEN (PHIS, LAMS)
|
|
|
15 |
C** IM ROTIERTEN SYSTEM. DER NORDPOL DIESES SYSTEMS HAT
|
|
|
16 |
C** DIE WAHREN KOORDINATEN (POLPHI, POLLAM)
|
|
|
17 |
C** VERSIONS-
|
|
|
18 |
C** DATUM : 03.05.90
|
|
|
19 |
C**
|
|
|
20 |
C** EXTERNALS: KEINE
|
|
|
21 |
C** EINGABE-
|
|
|
22 |
C** PARAMETER: PHIS REAL GEOGR. BREITE DES PUNKTES IM ROT.SYS.
|
|
|
23 |
C** LAMS REAL GEOGR. LAENGE DES PUNKTES IM ROT.SYS.
|
|
|
24 |
C** POLPHI REAL WAHRE GEOGR. BREITE DES NORDPOLS
|
|
|
25 |
C** POLLAM REAL WAHRE GEOGR. LAENGE DES NORDPOLS
|
|
|
26 |
C** AUSGABE-
|
|
|
27 |
C** PARAMETER: WAHRE GEOGRAPHISCHE LAENGE ALS WERT DER FUNKTION
|
|
|
28 |
C** ALLE WINKEL IN GRAD (NORDEN>0, OSTEN>0)
|
|
|
29 |
C**
|
|
|
30 |
C** COMMON-
|
|
|
31 |
C** BLOECKE : KEINE
|
|
|
32 |
C**
|
|
|
33 |
C** FEHLERBE-
|
|
|
34 |
C** HANDLUNG : KEINE
|
|
|
35 |
C** VERFASSER: D.MAJEWSKI
|
|
|
36 |
|
|
|
37 |
REAL LAMS,PHIS,POLPHI,POLLAM
|
|
|
38 |
|
|
|
39 |
DATA ZRPI18 , ZPIR18 / 57.2957795 , 0.0174532925 /
|
|
|
40 |
|
|
|
41 |
ZSINPOL = SIN(ZPIR18*POLPHI)
|
|
|
42 |
ZCOSPOL = COS(ZPIR18*POLPHI)
|
|
|
43 |
ZLAMPOL = ZPIR18*POLLAM
|
|
|
44 |
ZPHIS = ZPIR18*PHIS
|
|
|
45 |
ZLAMS = LAMS
|
|
|
46 |
IF(ZLAMS.GT.180.0) ZLAMS = ZLAMS - 360.0
|
|
|
47 |
ZLAMS = ZPIR18*ZLAMS
|
|
|
48 |
|
|
|
49 |
ZARG1 = SIN(ZLAMPOL)*(- ZSINPOL*COS(ZLAMS)*COS(ZPHIS) +
|
|
|
50 |
1 ZCOSPOL* SIN(ZPHIS)) -
|
|
|
51 |
2 COS(ZLAMPOL)* SIN(ZLAMS)*COS(ZPHIS)
|
|
|
52 |
ZARG2 = COS(ZLAMPOL)*(- ZSINPOL*COS(ZLAMS)*COS(ZPHIS) +
|
|
|
53 |
1 ZCOSPOL* SIN(ZPHIS)) +
|
|
|
54 |
2 SIN(ZLAMPOL)* SIN(ZLAMS)*COS(ZPHIS)
|
|
|
55 |
IF (ABS(ZARG2).LT.1.E-30) THEN
|
|
|
56 |
IF (ABS(ZARG1).LT.1.E-30) THEN
|
|
|
57 |
LMSTOLM = 0.0
|
|
|
58 |
ELSEIF (ZARG1.GT.0.) THEN
|
|
|
59 |
LMSTOLAM = 90.0
|
|
|
60 |
ELSE
|
|
|
61 |
LMSTOLAM = -90.0
|
|
|
62 |
ENDIF
|
|
|
63 |
ELSE
|
|
|
64 |
LMSTOLM = ZRPI18*ATAN2(ZARG1,ZARG2)
|
|
|
65 |
ENDIF
|
|
|
66 |
|
|
|
67 |
RETURN
|
|
|
68 |
END
|
|
|
69 |
|
|
|
70 |
|
|
|
71 |
REAL FUNCTION PHSTOPH (PHIS, LAMS, POLPHI, POLLAM)
|
|
|
72 |
C
|
|
|
73 |
C**** PHSTOPH - FC:BERECHNUNG DER WAHREN GEOGRAPHISCHEN BREITE FUER
|
|
|
74 |
C**** EINEN PUNKT MIT DEN KOORDINATEN (PHIS, LAMS) IM
|
|
|
75 |
C**** ROTIERTEN SYSTEM. DER NORDPOL DIESES SYSTEMS HAT
|
|
|
76 |
C**** DIE WAHREN KOORDINATEN (POLPHI, POLLAM)
|
|
|
77 |
C** AUFRUF : PHI = PHSTOPH (PHIS, LAMS, POLPHI, POLLAM)
|
|
|
78 |
C** ENTRIES : KEINE
|
|
|
79 |
C** ZWECK : BERECHNUNG DER WAHREN GEOGRAPHISCHEN BREITE FUER
|
|
|
80 |
C** EINEN PUNKT MIT DEN KOORDINATEN (PHIS, LAMS) IM
|
|
|
81 |
C** ROTIERTEN SYSTEM. DER NORDPOL DIESES SYSTEMS HAT
|
|
|
82 |
C** DIE WAHREN KOORDINATEN (POLPHI, POLLAM)
|
|
|
83 |
C** VERSIONS-
|
|
|
84 |
C** DATUM : 03.05.90
|
|
|
85 |
C**
|
|
|
86 |
C** EXTERNALS: KEINE
|
|
|
87 |
C** EINGABE-
|
|
|
88 |
C** PARAMETER: PHIS REAL GEOGR. BREITE DES PUNKTES IM ROT.SYS.
|
|
|
89 |
C** LAMS REAL GEOGR. LAENGE DES PUNKTES IM ROT.SYS.
|
|
|
90 |
C** POLPHI REAL WAHRE GEOGR. BREITE DES NORDPOLS
|
|
|
91 |
C** POLLAM REAL WAHRE GEOGR. LAENGE DES NORDPOLS
|
|
|
92 |
C** AUSGABE-
|
|
|
93 |
C** PARAMETER: WAHRE GEOGRAPHISCHE BREITE ALS WERT DER FUNKTION
|
|
|
94 |
C** ALLE WINKEL IN GRAD (NORDEN>0, OSTEN>0)
|
|
|
95 |
C**
|
|
|
96 |
C** COMMON-
|
|
|
97 |
C** BLOECKE : KEINE
|
|
|
98 |
C**
|
|
|
99 |
C** FEHLERBE-
|
|
|
100 |
C** HANDLUNG : KEINE
|
|
|
101 |
C** VERFASSER: D.MAJEWSKI
|
|
|
102 |
|
|
|
103 |
REAL LAMS,PHIS,POLPHI,POLLAM
|
|
|
104 |
|
|
|
105 |
DATA ZRPI18 , ZPIR18 / 57.2957795 , 0.0174532925 /
|
|
|
106 |
|
|
|
107 |
SINPOL = SIN(ZPIR18*POLPHI)
|
|
|
108 |
COSPOL = COS(ZPIR18*POLPHI)
|
|
|
109 |
ZPHIS = ZPIR18*PHIS
|
|
|
110 |
ZLAMS = LAMS
|
|
|
111 |
IF(ZLAMS.GT.180.0) ZLAMS = ZLAMS - 360.0
|
|
|
112 |
ZLAMS = ZPIR18*ZLAMS
|
|
|
113 |
ARG = COSPOL*COS(ZPHIS)*COS(ZLAMS) + SINPOL*SIN(ZPHIS)
|
|
|
114 |
|
|
|
115 |
PHSTOPH = ZRPI18*ASIN(ARG)
|
|
|
116 |
|
|
|
117 |
RETURN
|
|
|
118 |
END
|
|
|
119 |
|
|
|
120 |
|
|
|
121 |
REAL FUNCTION LMTOLMS (PHI, LAM, POLPHI, POLLAM)
|
|
|
122 |
C
|
|
|
123 |
C%Z% Modul %M%, V%I% vom %G%, extrahiert am %H%
|
|
|
124 |
C
|
|
|
125 |
C**** LMTOLMS - FC:UMRECHNUNG DER WAHREN GEOGRAPHISCHEN LAENGE LAM
|
|
|
126 |
C**** AUF EINEM PUNKT MIT DEN KOORDINATEN (PHIS, LAMS)
|
|
|
127 |
C**** IM ROTIERTEN SYSTEM. DER NORDPOL DES SYSTEMS HAT
|
|
|
128 |
C**** DIE WAHREN KOORDINATEN (POLPHI, POLLAM)
|
|
|
129 |
C** AUFRUF : LAM = LMTOLMS (PHI, LAM, POLPHI, POLLAM)
|
|
|
130 |
C** ENTRIES : KEINE
|
|
|
131 |
C** ZWECK : UMRECHNUNG DER WAHREN GEOGRAPHISCHEN LAENGE LAM AUF
|
|
|
132 |
C** EINEM PUNKT MIT DEN KOORDINATEN (PHIS, LAMS) IM
|
|
|
133 |
C** ROTIERTEN SYSTEM. DER NORDPOL DIESES SYSTEMS HAT
|
|
|
134 |
C** DIE WAHREN KOORDINATEN (POLPHI, POLLAM)
|
|
|
135 |
C** VERSIONS-
|
|
|
136 |
C** DATUM : 03.05.90
|
|
|
137 |
C**
|
|
|
138 |
C** EXTERNALS: KEINE
|
|
|
139 |
C** EINGABE-
|
|
|
140 |
C** PARAMETER: PHI REAL BREITE DES PUNKTES IM GEOGR. SYSTEM
|
|
|
141 |
C** LAM REAL LAENGE DES PUNKTES IM GEOGR. SYSTEM
|
|
|
142 |
C** POLPHI REAL GEOGR.BREITE DES N-POLS DES ROT. SYSTEMS
|
|
|
143 |
C** POLLAM REAL GEOGR.LAENGE DES N-POLS DES ROT. SYSTEMS
|
|
|
144 |
C** AUSGABE-
|
|
|
145 |
C** PARAMETER: WAHRE GEOGRAPHISCHE LAENGE ALS WERT DER FUNKTION
|
|
|
146 |
C** ALLE WINKEL IN GRAD (NORDEN>0, OSTEN>0)
|
|
|
147 |
C**
|
|
|
148 |
C** COMMON-
|
|
|
149 |
C** BLOECKE : KEINE
|
|
|
150 |
C**
|
|
|
151 |
C** FEHLERBE-
|
|
|
152 |
C** HANDLUNG : KEINE
|
|
|
153 |
C** VERFASSER: G. DE MORSIER
|
|
|
154 |
|
|
|
155 |
REAL LAM,PHI,POLPHI,POLLAM
|
|
|
156 |
|
|
|
157 |
DATA ZRPI18 , ZPIR18 / 57.2957795 , 0.0174532925 /
|
|
|
158 |
|
|
|
159 |
ZSINPOL = SIN(ZPIR18*POLPHI)
|
|
|
160 |
ZCOSPOL = COS(ZPIR18*POLPHI)
|
|
|
161 |
ZLAMPOL = ZPIR18*POLLAM
|
|
|
162 |
ZPHI = ZPIR18*PHI
|
|
|
163 |
ZLAM = LAM
|
|
|
164 |
IF(ZLAM.GT.180.0) ZLAM = ZLAM - 360.0
|
|
|
165 |
ZLAM = ZPIR18*ZLAM
|
|
|
166 |
|
|
|
167 |
ZARG1 = - SIN(ZLAM-ZLAMPOL)*COS(ZPHI)
|
|
|
168 |
ZARG2 = - ZSINPOL*COS(ZPHI)*COS(ZLAM-ZLAMPOL)+ZCOSPOL*SIN(ZPHI)
|
|
|
169 |
IF (ABS(ZARG2).LT.1.E-30) THEN
|
|
|
170 |
IF (ABS(ZARG1).LT.1.E-30) THEN
|
|
|
171 |
LMTOLMS = 0.0
|
|
|
172 |
ELSEIF (ZARG1.GT.0.) THEN
|
|
|
173 |
LMTOLMS = 90.0
|
|
|
174 |
ELSE
|
|
|
175 |
LMTOLMS = -90.0
|
|
|
176 |
ENDIF
|
|
|
177 |
ELSE
|
|
|
178 |
LMTOLMS = ZRPI18*ATAN2(ZARG1,ZARG2)
|
|
|
179 |
ENDIF
|
|
|
180 |
|
|
|
181 |
RETURN
|
|
|
182 |
END
|
|
|
183 |
|
|
|
184 |
|
|
|
185 |
REAL FUNCTION PHTOPHS (PHI, LAM, POLPHI, POLLAM)
|
|
|
186 |
C
|
|
|
187 |
C%Z% Modul %M%, V%I% vom %G%, extrahiert am %H%
|
|
|
188 |
C
|
|
|
189 |
C**** PHTOPHS - FC:UMRECHNUNG DER WAHREN GEOGRAPHISCHEN BREITE PHI
|
|
|
190 |
C**** AUF EINEM PUNKT MIT DEN KOORDINATEN (PHIS, LAMS)
|
|
|
191 |
C**** IM ROTIERTEN SYSTEM. DER NORDPOL DES SYSTEMS HAT
|
|
|
192 |
C**** DIE WAHREN KOORDINATEN (POLPHI, POLLAM)
|
|
|
193 |
C** AUFRUF : PHI = PHTOPHS (PHI, LAM, POLPHI, POLLAM)
|
|
|
194 |
C** ENTRIES : KEINE
|
|
|
195 |
C** ZWECK : UMRECHNUNG DER WAHREN GEOGRAPHISCHEN BREITE PHI AUF
|
|
|
196 |
C** EINEM PUNKT MIT DEN KOORDINATEN (PHIS, LAMS) IM
|
|
|
197 |
C** ROTIERTEN SYSTEM. DER NORDPOL DIESES SYSTEMS HAT
|
|
|
198 |
C** DIE WAHREN KOORDINATEN (POLPHI, POLLAM)
|
|
|
199 |
C** VERSIONS-
|
|
|
200 |
C** DATUM : 03.05.90
|
|
|
201 |
C**
|
|
|
202 |
C** EXTERNALS: KEINE
|
|
|
203 |
C** EINGABE-
|
|
|
204 |
C** PARAMETER: PHI REAL BREITE DES PUNKTES IM GEOGR. SYSTEM
|
|
|
205 |
C** LAM REAL LAENGE DES PUNKTES IM GEOGR. SYSTEM
|
|
|
206 |
C** POLPHI REAL GEOGR.BREITE DES N-POLS DES ROT. SYSTEMS
|
|
|
207 |
C** POLLAM REAL GEOGR.LAENGE DES N-POLS DES ROT. SYSTEMS
|
|
|
208 |
C** AUSGABE-
|
|
|
209 |
C** PARAMETER: ROTIERTE BREITE PHIS ALS WERT DER FUNKTION
|
|
|
210 |
C** ALLE WINKEL IN GRAD (NORDEN>0, OSTEN>0)
|
|
|
211 |
C**
|
|
|
212 |
C** COMMON-
|
|
|
213 |
C** BLOECKE : KEINE
|
|
|
214 |
C**
|
|
|
215 |
C** FEHLERBE-
|
|
|
216 |
C** HANDLUNG : KEINE
|
|
|
217 |
C** VERFASSER: G. DE MORSIER
|
|
|
218 |
|
|
|
219 |
REAL LAM,PHI,POLPHI,POLLAM
|
|
|
220 |
|
|
|
221 |
DATA ZRPI18 , ZPIR18 / 57.2957795 , 0.0174532925 /
|
|
|
222 |
|
|
|
223 |
ZSINPOL = SIN(ZPIR18*POLPHI)
|
|
|
224 |
ZCOSPOL = COS(ZPIR18*POLPHI)
|
|
|
225 |
ZLAMPOL = ZPIR18*POLLAM
|
|
|
226 |
ZPHI = ZPIR18*PHI
|
|
|
227 |
ZLAM = LAM
|
|
|
228 |
IF(ZLAM.GT.180.0) ZLAM = ZLAM - 360.0
|
|
|
229 |
ZLAM = ZPIR18*ZLAM
|
|
|
230 |
ZARG = ZCOSPOL*COS(ZPHI)*COS(ZLAM-ZLAMPOL) + ZSINPOL*SIN(ZPHI)
|
|
|
231 |
|
|
|
232 |
PHTOPHS = ZRPI18*ASIN(ZARG)
|
|
|
233 |
|
|
|
234 |
RETURN
|
|
|
235 |
END
|