Astrometric Library
Benutzung
Die altbekannte ARI Programmbibliothek kann mit Hilfe des Compilerparameters -lmari eingebunden werden, wobei auch noch der Suchpfad -L/usr/local/lib mit anzugeben ist. Z.B.:
g77 -L/usr/local/lib -lmari myprog.f
Liste der Unterprogramme mit Beschreibung
Eine ausführliche Beschreibungen kann man den Programmquellen /usr/local/src/libmari/*.f entnehmen.
ANGLE - BERECHNUNG VON ALPHA UND DELTA AUS DEM VEKTOR DER
RICHTUNGSKOSIN.
AUFRUF: CALL ANGLE (V,ALP,DEL)
-----------------------------------------------------------------------
BEJD - BESSELSCHE EPOCHEN IN JULIANISCHE DATEN UMWANDELN.
AUFRUF: CALL BEJD(BE,TJD)
-----------------------------------------------------------------------
BEJE - UMWANDLUNG VON BESSELSCHER IN JULIANISCHER EPOCHE.
AUFRUF CALL BEJE(XBEP,XJEP)
-----------------------------------------------------------------------
BEPOCH - BERECHNUNG VON BESSELSCHER EPOCHE AUS JULIANISCHEM DATUM.
AUFRUF: CALL BEPOCH(XJ,BE)
-----------------------------------------------------------------------
CHARAD - UMWANDLUNG EINES WINKELS, ANGEGEBEN IN GRAD, BOGEN- BZW.
ZEITSTUNDEN, MINUTEN ODER SEKUNDEN, IN BOGENMASS.
AUFRUF: CALL CHARAD(C1,C2,RDG,IER,LUNERB)
-----------------------------------------------------------------------
DIPOWI - ANTRAGEN VON ABSTAND UND POSITIONSWINKEL AN REKTASZENSION
UND DEKLINATION.
AUFRUF: CALL DIPOWI(AALT,DALT,DIST,POWI,ANEU,DNEU)
-----------------------------------------------------------------------
DIRCOS - BERECHNUNG VON KARTESISCHEN KOORDINATEN (RICHTUNGSKOSINUSSEN
AUS LAENGEN- UND BREITENKOORDINATEN) Z.B. ALPHA UND DELTA.
AUFRUF: CALL DIRCOS(ALP,DEL,V)
-----------------------------------------------------------------------
DLDA - ABLEITUNG DER RICHTUNGSKOSINUSSE (KARTES. VEKTORKOMPONENTEN)
EINES ORTES NACH EINER LAENGENKOORDINATE (Z.B. ALPHA).
AUFRUF: CALL DLDA(ALP,DEL,V)
-----------------------------------------------------------------------
DLDD - ABLEITUNG DER RICHTUNGSKOSINUSSE (KARTES. VEKTORKOMPONENTEN)
EINES ORTES NACH EINER BREITENKOORDINATE (Z.B. DELTA).
AUFRUF: CALL DLDD(ALP,DEL,V)
-----------------------------------------------------------------------
DREMTX - BERECHNUNG DER DREHMATRIX UM DIE ACHSE MIT DREHWINKEL
AUFRUF: CALL DREMTX(IACHS,WINKEL,T)
-----------------------------------------------------------------------
DSCONV - CONVERTIERT DEN BISHERIGEN DOPPELSTERNCODE IN DAS NEUE
DS-SYSTEM. (FUER DIE IDENTIFIZIERUNG)
AUFRUF: CALL DSCONV(DATSANA,SP1,SP2,IDSNEU,IERR)
-----------------------------------------------------------------------
DSCONVER- CONVERTIERT DEN BISHERIGEN DOPPELSTERNCODE IN DAS NEUE
DS-SYSTEM. (FUER DIE IDENTIFIZIERUNG)
AUFRUF: CALL DSCONV(DATSANA,SP1,SP2,IDSNEU,IERR)
-----------------------------------------------------------------------
ELAB - BERECHNET KORREKTUR DER A,D, WEGEN ELIPTISCHER ABERRATION
IM KOORDINATENSYSTEM 1950.0.
AUFRUF: CALL ELAB(A,D,AN,DN,EPA,EPD)
-----------------------------------------------------------------------
ELABER - BEHANDLUNG DER ELLIPTISCHEN ABERRATION.
AUFRUF: CALL ELABER(IVERS,XEQJ,ALPHA1,DELTA1,ALPHA2,DELTA2)
-----------------------------------------------------------------------
ELBA - EINRECHNEN DER ELLIPTISCHEN ABERRATION IN STERNOERTERN.
(UMKEHRUNG DER SUBROUTINE ELAB)
AUFRUF: CALL ELBA(A,D,AN,DN,EPA,EPD)
-----------------------------------------------------------------------
EQUECL - EKLIPTIKALE KOORDINATEN IN SEQUATORIALE KOORDINATEN
UMRECHNEN UND UMGEKEHRT.
AUFRUF: CALL EQUECL(X1,X2,Y1,Y2,EPS,L)
-----------------------------------------------------------------------
EQUGAL - TRANSFORMATION VON AEQUATORIALEN KOORDIN. UND EIGEN-
BEWEGUNGEN IN GALAKTISCHE UND UMGEKEHRT.
AUFRUF: CALL EQUGAL(L,ALP,DEL,PMA,PMD,ALPNEU,DELNEU,
PMANEU,PMDNEU)
-----------------------------------------------------------------------
EQUIN - PRAEZEDIEREN VON AEQUATORIALEN KOORDINATEN UND EIGEN-
BEWEGUNGEN. (ACHTUNG ALTE NEWCOMBSCHE PRAEZESSION)
AUFRUF: CALL EQUIN(EQU,EQUNEU,ALP,DEL,PMA,PMD,ALPNEU,DELNEU,
PMANEU,PMDNEU)
-----------------------------------------------------------------------
EQUIN5 - PRAEZEDIEREN VON AEQUATORIALEN KOORDINATEN UND EIGEN-
BEWEGUNGEN. (ACHTUNG NEUE (IAU 1976) PRAEZESSION)
AUFRUF: CALL EQUIN5(TJD,TJDNEU,ALP,DEL,PMA,PMD,ALPNEU,
DELNEU,PMANEU,PMDNEU)
-----------------------------------------------------------------------
ETERM - BERECHNUNG DER KORREKTIONEN, DIE VON DER ELLIPTIZITAET
DER ERDBAHN HERRUEHREN (LINEARER TEIL).
AUFRUF: CALL ETERM(ALF,DEL C0,D0,EPS,C1,C2)
-----------------------------------------------------------------------
EULMTX - BERECHNUNG DER DREHMATRX AUS DEN DREI EULERSCHEN WINKELN.
AUFRUF: CALLEULMTX(PHT, THETA,PSI,T)
-----------------------------------------------------------------------
GMADD - ADDITION ZWEIER MATRIZEN.
AUFRUF: CALL GMADD(A,B,R,M,N)
-----------------------------------------------------------------------
GMPRD - MULTIPLIKATION ZWEIER REAL*8 MATRIZEN.
AUFRUF: CALL GMPRD(A,B,R,M,N,L)
-----------------------------------------------------------------------
GNOINV - UMKEHRUNG DER GNOMONISCHEN ABBILDUNG.
AUFRUF: CALL GNOINV(ALPHA0,DELTA0,ALPHA,DELTA,X,Y)
-----------------------------------------------------------------------
GNOMAP - GMOMISCHE ABBILDUNG DER EINHEITSKUGEL AUF DIE
TANGENDIALEBENE,
AUFRUF: CALL GNOMAP(ALPHA0,DELTA0,ALPHA,DELTA,X,Y)
-----------------------------------------------------------------------
HADE - BERECHNUNG VON STUNDENWINKEL UND DEKLINATION AUS
GEOGR. BREITE, ZENITDISTANZ UND AZIMUT.
AUFRUF: CALL HADE(PHI,ZD,AZ,HA,DELTA)
-----------------------------------------------------------------------
HISTO - ERSTELLEN EINES HISTOGRAMMS. ZEICHNEN EINES HISTOGRAMMS.
BERECHNUNG VON MITTELWERT UNS STREUUNG.
AUFRUF: CALL HISTO(IHTYP,LUN,XWERT,XMIN,XMAX,NINT,X,Y
XQUER,SIGMA,IER,TEXTX,TEXTY,TITEL)
-----------------------------------------------------------------------
HXPRNT - ABDRUCKEN DES HEX-CODES EINER DATEI FESTER SATZLAENGE.
AUFRUF: CALL HXPRNT(LENI,LENO,ZEILI,ZEILO,DSN)
-----------------------------------------------------------------------
IJ - HILFSPROGRAMM FUER DAS PROGRAMMPAKET SLSS.
AUFRUF: CALL IJ(I,J,MODE,NS)
-----------------------------------------------------------------------
INTEP - SPLINE-INTERPOLATION EINES 1-DIM. ARRAYS (GLATT)
AUFRUF: CALL INTEP(X,Y,XARRAY,YARRAY,IMAX,IERR)
-----------------------------------------------------------------------
JDBE - UMWANDLUNG DES JULIANISCHEN DATUMS IN BESSELSCHE EPOCHE
AUFRUF: CALL JDBE(JD,BE)
-----------------------------------------------------------------------
JEBE - UMWANDLUNG VON BESSELSCHER IN JULIANISCHER EPOCHE.
AUFRUF: CALL JEBE(XJEP,XBEP)
-----------------------------------------------------------------------
JEJD - UMWANDLUNG DER JULIANISCHEN EPOCHE IN JULIANISCHES DATUM.
AUFRUF: CALL JEJD(JE,JD)
-----------------------------------------------------------------------
JULDAT - BERECHNUNG DES JULIANISCHEN DATUMS UND DER
BESSELSCHEN EPOCHE AUF BUERGEL. DATUM, MITTL. ORTSZEIT
UND GEOGR. LAENGE DES BEOBACHTUNGSORTS.
AUFRUF: CALL JULDAT(IJ,IM,IT,IH,JM,SEK,SEKL,XL,BE)
-----------------------------------------------------------------------
JULKAL - JULIANISCHES DATUM IN BUERGLICHES DATUM + UT UMWANDELN.
AUFRUF: CALL JULKAL(TJD,TJAHR,TMONAT,TAG,UTH,UTM,UTS)
-----------------------------------------------------------------------
KALBE - BERECHNUNG DER BESSELSCHEN EPOCHE AUS EPOCHE, JAHR,
MONAT,TAG, STD., MIN., SEC. IN BUERGERLICHER ZEIT ( UT ).
AUFRUF: CALL KALBE(JAHR,MONAT,TAG,ST,MIN,SEC,BE )
-----------------------------------------------------------------------
KALJUL - BERECHNUNG DES JULIANISCHEN DATUMS FUER DIE EPOCH, JAHR,
MONAT, TAG, STD., MIN., SEC. IN BUERGERLICHER ZEIT ( UT ).
AUFRUF: CALL KALBE(JAHR,MONAT,TAG,ST,MIN,SEC,JD)
-----------------------------------------------------------------------
KPRZ - BESTIMMUNG DER PRAEZESSIONSWINKEL NACH VERSCHIEDENEN
KONSTANTEN JE NACH WAHL DES PARAMETERS IP.
AUFRUF: CALL KPRZ(IP,T0,T1,ZET,TETA,ZETA)
-----------------------------------------------------------------------
KSTETR - BESTIMMUNG DES 'ELLIPTISCHEN TEILS DER ABERATION
DAY NUMBERS'.
AUFRUF: CALL KSTETR(BY,C0,D0,EPS)
-----------------------------------------------------------------------
KSTNUT - BERECHNUNG DEL LANG- UND KURZPERIODISCHEN NUTATION
IN LAENGE UND SCHIEFE (BEZOGEN AUF ECLIPTIC OF DATE).
AUFRUF CALL KSTNUT(T,DPSTL,DPSIS,DEPSL,DEPSS,EPS)
-----------------------------------------------------------------------
KSTPR5 - STANDARTDREHWINKEL FUER PRAEZESSION ERZEUGEN.
(ACHTUNG: NEUE PRAEZESSION TAU 1976)
AUFRUF: CALL KSTPR5(TJD1,TJD2,ZET,TETA,ZETA)
-----------------------------------------------------------------------
KSTPRZ - STANDARTDREHWINKEL FUER PRAEZESSION ERZEUGEN.
(ACHTUNG: NEWCOMB]S PRAEZESSION 'ALTE KONSTANTE')
AUFRUF: CALL KSTPRZ(TJD1,TJD2,ZET,TETA,ZETA)
-----------------------------------------------------------------------
LINFIT - ES WIRD EINE AUSGLEICHSGERADE DURCH DIE PUNKTE
T(I) (ARGUMENTE) X(I) (WERTE) MIT I=1, ... , N
BESTIMMT.
AUFRUF CALL LINFIT(N,T,X,W,MT,NX,M,X0,SMX,SM.SX0,RHO,S0)
-----------------------------------------------------------------------
MTXVEK - MULTIPLIKATION EINER REAL*8-MATRIX MIT EINEM VEKTOR.
AUFRUF: CALL MTXVEK(A,VIN,VOUT,N,M)
-----------------------------------------------------------------------
NUTA - BERECHNUNG DER ELEMENTE UND DER DREHMATRIX DER NUTATION.
AUFRUF: CALL NUTA(T,RN,EP, PS,O,F,D,RL,RS)
-----------------------------------------------------------------------
NUTMTX - BERECHNUNG DER DREHMATRIX DER NUTATION (AUS PRAMETERN,
DIE MIT KSTNUT BESTIMMT WURDEN).
AUFRUF: CALL NUTMTX(DPSI,DEPS,EPS,XN)
-----------------------------------------------------------------------
PLATE - BESTIMMEN DER AC-PLATTEN, AUF DENEN SICH EIN STERN
MIT GEGEBENEN ALPHA, DELTA BEFINDET.
AUFRUF: CALL PLATE(A,D,IPLAT)
-----------------------------------------------------------------------
POLDEC - UMRECHNUNG VON POLDISTANZEN IN DEKLINATION.
AUFRUF: CALL POLDEC(IVERS,IPG,IPM,PS,DEC)
-----------------------------------------------------------------------
POWI - POSITIONSWINKEL AUS ZWEI PAAREN VON SPHAERISCHEN
KOORDINATEN (LAENGE,BREITE) BERECHNEN.
AUFRUF: CALL POWI(ALP1,DEL1,ALP2,DEL2,PW,IER)
-----------------------------------------------------------------------
PRCPOS - PRAEZESSIONSUEBERGABE VON OERTERN
AUFRUF: CALL PRCPOS(ALPOLD,DELOLD,ALPNEW,DELNEW,PX)
-----------------------------------------------------------------------
PRECES - ANBRINGEN EINER DREHMATRIX AN KOORDINATEN UND/ODER
EIGENBEWEGUNGEN (Z.B. PRAEZESSIONMATRIX)
AUFRUF: CALL PRECES(ALP,DEL,PMA,PMD.PX,ALPNEU,DELNEU,
PMANEU,PMDNEU)
-----------------------------------------------------------------------
PREFK5 - KORREKTUR VON POSITIONEN UND EIGENBEWEGUNGEN DES
FK4 NACH FRICKE.
AUFRUF: CALL PREFK5(AOLD,DOLD,AMUOLD,DMUOLD,EPOCH,
ANEW,DNEW,AMUNEW,DMUNES)
-----------------------------------------------------------------------
PRZMTX - DREHMATRIX AUS PRAEZESSIONSWINKEL BERECHNEN.
AUFRUF: CALL PRZMTX(ZET,TETA,ZETA,PX)
-----------------------------------------------------------------------
RADCHA - BEQUEME UMRECHNUNG VON WINKELN, DIE IM BOGENMASS ANGEGEBEN
SIND, IN EINHEITEN VON STUNDEN ODER GRAD, MINUTEN UND
SEKUNDEN ODER ODER IN DEREN SINNVOLLEN KOMBINATION.
AUFRUF: CALL RADCHA(RAD,FORMAT,WERT,IER,LUNERR)
-----------------------------------------------------------------------
REFRAT - BERECHNUNG DER ASTRONOMISCHEN REFRAKTION
SCHEINBARE ZENITDISTANZ IN WAHRE ZENITDISTANZ.
AUFRUF: CALL REFRAT(ZAPP,ZTRUE)
-----------------------------------------------------------------------
REFRTA - BERECHNUNG DER ASTRONOMISCHEN REFRAKTION
WAHRE ZENITDISTANZ IN SCHEINBARE ZENITDISTANZ.
AUFRUF: CALL REFRTA(ZTRUE,ZAPP)
-----------------------------------------------------------------------
SCLPRD - SKALARPRODUKT ZWEIER REAL*8-VEKTOREN.
AUFRUF: CALL SCLPRD(A,B,R,M)
-----------------------------------------------------------------------
SMPY - MULTIPLIKATION EINER MATRIZ MIT EINEM SKALAR.
AUFRUF: SMPY(A,C,R,N,M)
-----------------------------------------------------------------------
SPACEV - BERECHNUNG VON RAUMGESCHWINDIGKEIT UND DREIDIM. ORT
AUS POSITION, PARALLAXE, EIGENBEWEGUNG UND RADIAL-
GESCHWINDIGKEIT.
AUFRUF: CALL SPACEV(ALP,DEL,PMA,PMD,PAR,RVEL,X,V)
-----------------------------------------------------------------------
SPRING - UEBERSPRINGEN VON SAETZEN (RECONDS) IN EINER UNFORMATIERTEN
SEQUENTIELLEN DATEI.
AUFRUF: CALL SPRING(LUN,N)
-----------------------------------------------------------------------
STHMS - BERECHNUNG DER STUNDEN, MINUTEN UND SEKUNDEN (AUF 3
DEZIMALEN) AUS DEN GEGEBENEN STUNDEN (MIT DEZIMALEN).
AUFRUF: CALL STHMS(X,XH,XM,XS)
-----------------------------------------------------------------------
STUT - BERECHNUNG DER MITTLEREN WELTZEIT IND STUNDEN, MINUTEN UND
SEKUNDEN SOWIE DES JULIANISCHEN DATUMS FUER DEN SIDERISCHEN
TAG GSD, MITTLERE SID. ZEIT STH,STM,STS NACH AENA.
AUFRUF: CALL STUT(GSD,STH,STM,STS,UTH,UTM,UTS,TJD)
-----------------------------------------------------------------------
TPMPM - EPOCHENUEBERTRAGUNG VON OERTERN UND EIGENGEWEGUNGEN.
AUFRUF: CALL TPMPM(ALP,DEL,PMA,PMD,FORESH,TJD,TJDNEU,
ALPNEU,DELNEU,PMANEU,PMDNEU)
-----------------------------------------------------------------------
TPMPOS - EPOCHENUEBERTRAGUNG VON MITTLEREN OERTERN.
('EIGENBEWEGUNG ANBRINGEN')
AUFRUF: CALL TPMPOS(ALP,DEL,PMA,PMD,FORESH,TJD,TJDNEU,
ALPNEU,DELNEU,PMANEU,PMDNEU)
-----------------------------------------------------------------------
TPMSPC - EPOCHENUEBERTRAGUNG VON OERTERN UND EIGENBEWEGUNGEN
UNTER DER ANNAHME KONSTANTER RAUMGESCHWINDIGKEIT.
AUFRUF: CALL TPMSPC(T0,T1,A0,D0,XMY0,XMYS0,PAR0,VEL0,
A1,D1,XMY1,XMYS1,PAR1,VEL1,IER)
-----------------------------------------------------------------------
TSPV2 - EPOCHEN- UND AEQUINOX-UEBERTRAGUNG VON OERTERN UND
EIGENBEWEGUNGEN
EPOCHENUEBERTRAGUNG MIT BENUTZUNG DER RAUMGESCHW.
AUFRUF: CALL TSPV2(IP,XEQ,XEP,YEQ,YEP,XA,XD,XM,
XMS,XP0,XVEL0,YA,YD,YM,YMS,YP0,
YVEL0,IER)
-----------------------------------------------------------------------
TSPVEL - EPOCHEN- UND AEQUINOX-UEBERTRAGUNG VON OERTERN UND
EIGENBEWEGUNGEN
EPOCHENUEBERTRAGUNG MIT BENUTZUNG DER RAUMGESCHW.
AUFRUF: CALL TSPVEL(XEQ,XEP,YEQ,YEP,XA,XD,XM,
XMS,XP0,XVEL0,YA,YD,YM,YMS,YP0,
YVEL0,IER)
-----------------------------------------------------------------------
UELWEG - AUFSPALTEN 'UEBERLOCHTER' ZIFFERN IN ZIFFER UND
UEBERLOCH.
AUFRUF: CALL UELWEG(SPALTE,ZIFFER,UELOCH,IER)
-----------------------------------------------------------------------
UTST - BERECHNUNG DER MITTLEREN SIDERISCHEN ZEIT IN STUNDE,
MINUTEN UND SEKUNDEN SOWIE DER GREENWICH SIDERIAL DAY
NUMBER FUER EIN JULIANISCHES DATUM (JD) UND EINE MITTLERE
WELTZEIT (UT)
AUFRUF: CALL UTST(TJD,UTH,UTM,UTS,STH,STM,STS,GSD)
-----------------------------------------------------------------------
UTST2 - BERECHNUNG DER GREENWICH SIDERIAL DAY NUMBER
(UND DAMIT DER MITTELREN SIDERISCHEN ZEIT)
FUER EIN JULIANISCHES DATUM (JD).
AUFRUF: CALL UTST2(TJD,GSD)
-----------------------------------------------------------------------
VEKPRD - DAS GEWOEHNLICHER VEKTORPRODUKT BERECHNEN.
AUFRUF: CALL VEKPRD(V1,V2,V3)
-----------------------------------------------------------------------
ZDAZ - BERECHNUNG VON ZENITDISTANZ UND ASIMUT AUS
GEOGR. BREITE, STUNDENWINKEL UND DEKLINATION.
AUFRUF: CALL ZDAZ(PHI,HA,DELTA,ZD,AZ)
-----------------------------------------------------------------------
