Subversion Repositories lagranto.um

Rev

Rev 3 | Only display areas with differences | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

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