Index: trunk/MagicSoft/Mars/mastro/MAstro.cc
===================================================================
--- trunk/MagicSoft/Mars/mastro/MAstro.cc	(revision 8140)
+++ trunk/MagicSoft/Mars/mastro/MAstro.cc	(revision 8324)
@@ -446,4 +446,35 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Estimates the time at which a source culminates.
+//
+// ra: right ascension [rad]
+// elong: observers longitude [rad]
+// mjd: modified julian date (utc)
+//
+// return time in [-12;12]
+//
+Double_t MAstro::EstimateCulminationTime(Double_t mjd, Double_t elong, Double_t ra)
+{
+    // startime at 1.1.2000 for greenwich 0h
+    const Double_t gmt0 = 6.664520;
+
+    // difference of startime for greenwich for two calendar days [h]
+    const Double_t d0 = 0.06570982224;
+
+    // mjd of greenwich 1.1.2000 0h
+    const Double_t mjd0 = 51544;
+
+    // mjd today
+    const Double_t mjd1 = TMath::Floor(mjd);
+
+    // scale between star-time and sun-time
+    const Double_t scale = 1;//1.00273790926;
+
+    const Double_t UT = (ra-elong)*RadToHor() - (gmt0 + d0 * (mjd1-mjd0))/scale;
+
+    return fmod(2412 + UT, 24) - 12;
+}
 
 // --------------------------------------------------------------------------
Index: trunk/MagicSoft/Mars/mastro/MAstro.h
===================================================================
--- trunk/MagicSoft/Mars/mastro/MAstro.h	(revision 8140)
+++ trunk/MagicSoft/Mars/mastro/MAstro.h	(revision 8324)
@@ -31,6 +31,6 @@
     static Double_t HorToDeg()        { return 1./15; }
 
-    static Double_t RadToHor()        { return 15/TMath::TwoPi(); }
-    static Double_t HorToRad()        { return TMath::TwoPi()/15; }
+    static Double_t RadToHor()        { return 24/TMath::TwoPi(); }
+    static Double_t HorToRad()        { return TMath::TwoPi()/24; }
 
     // Angle/Time conversion functions
@@ -86,4 +86,6 @@
     static Int_t    GetMagicPeriod(Double_t mjd);
 
+    static Double_t EstimateCulminationTime(Double_t mjd, Double_t longit, Double_t ra);
+
     // Estimate some parameters around the sun
     static Double_t GetSunRaDec(Double_t mjd, Double_t &ra, Double_t &dec);
Index: trunk/MagicSoft/Mars/mastro/MObservatory.cc
===================================================================
--- trunk/MagicSoft/Mars/mastro/MObservatory.cc	(revision 8140)
+++ trunk/MagicSoft/Mars/mastro/MObservatory.cc	(revision 8324)
@@ -110,6 +110,6 @@
     *fLog << all;
     *fLog << underline << fObservatoryName << ":" << endl;
-    *fLog << "Latitude:  " << TMath::Abs(fLatitude*kRad2Deg)  << " deg " << (fLatitude  > 0 ? "W" : "E") << endl;
-    *fLog << "Longitude: " << TMath::Abs(fLongitude*kRad2Deg) << " deg " << (fLongitude < 0 ? "N" : "S") << endl;
+    *fLog << "Latitude:  " << TMath::Abs(fLatitude*kRad2Deg)  << " deg " << (fLatitude  > 0 ? "N" : "S") << endl;
+    *fLog << "Longitude: " << TMath::Abs(fLongitude*kRad2Deg) << " deg " << (fLongitude < 0 ? "E" : "W") << endl;
     *fLog << "Height:    " << fHeight << "m" << endl;
 }