Changeset 9864 in ntrip for trunk/BNC/scripts

Timestamp:

Oct 31, 2022, 9:49:03 PM (2 years ago)

Author:

stuerze

Message:

minor changes

File:

• trunk/BNC/scripts/compareSP3Plot.pl (modified) (10 diffs)

trunk/BNC/scripts/compareSP3Plot.pl

@@ -1 +1 @@
 #!/usr/bin/env perl
 # ========================================================================
-# cmbPlot.pl
+# compareSP3Plot.pl
 # ========================================================================
 #
@@ -25 +25 @@
 use File::Basename;
 use Date::Manip;
+use Time::Piece 1.30;
 use List::MoreUtils qw( minmax );
 use experimental 'smartmatch';
-
-use Gps_Date; # Bernese
 
 # -----------------------------------------------------------------------------
@@ -61 +60 @@
 # generate data sets for gnuplot
 # -----------------------------------------------------------------------------
-my $old_epochSec = 0;
-my $epochSec = 0;
-my $epochDiff = 0;
+my $epoSec = 0;
 # for pdf generation
 my ($png, $page, $headline, $headline_text);
@@ -73 +70 @@
     my (%COMPARISON);
     my (@array);
-    my ($sys, $sat, $epo, $mjd, $radial, $along, $out, $drange, $clk, $clkRed);
+    my ($sys, $sat, $epo, $radial, $along, $out, $drange, $clk, $clkRed);
     # -----------------------------------------------------------------------------
     # create pdf for plot results
@@ -89 +86 @@
         if ( $_ !~ /^!/) { #print "$_\n";
             @array = split(/\s+/, $_);
-            $mjd = $array[0];
-            $epo  = gps_date("-mjd","$mjd","-o %Y-%m-%d_%H:%M:%S");
-            $sys = substr $array[1], 0, 1;
+            $epo   = $array[0];   #print "$epo\n";
+            $epo   =~s/(\.\d*)//; #print "$epo\n";
+            $sys    = substr $array[1], 0, 1;
             $sat    = $array[1];#print "$sat\n";
             $radial = $array[2];#print "$radial\n";
@@ -213 +210 @@
     if (("ALL" ~~ @plotTypes ) || ("CLK" ~~ @plotTypes)) {
         print "Plot CLOCK and CLOCK REDUCED BY RADIAL COMPONENT ...\n";
-        my $yrange = 10.0;
+        my $yrange = 0.6;
         $page = $pdf->page(); $page->mediabox('A4');
         $headline = sprintf("Clock residuals and reduced clock residuals(%s)", $inputFilename);
@@ -229 +226 @@
             my $chart_c = Chart::Gnuplot->new(
                           output => $pngName_c, terminal => 'png',title => "Clock only",
-                          ylabel => "Clock Residuals [m]",# yrange => ["-$yrange", "$yrange"],
+                          ylabel => "Clock Residuals [m]", yrange => ["-$yrange", "$yrange"],
                           xlabel => "Time [h]", timeaxis => 'x', xtics => {labelfmt => '%H:%M', rotate => '-270',},
                           legend => { position => "outside right",},
@@ -235 +232 @@
             my $chart_r = Chart::Gnuplot->new(
                           output => $pngName_r, terminal => 'png',title => "Clock minus radial component",
-                          ylabel => "Clock Residuals [m]", #yrange => ["-$yrange", "$yrange"],
+                          ylabel => "Clock Residuals [m]", yrange => ["-$yrange", "$yrange"],
                           xlabel => "Time [h]", timeaxis => 'x', xtics => {labelfmt => '%H:%M', rotate => '-270',},
                           legend => { position => "outside right",},
@@ -286 +283 @@
     $pdf->end();
     system("rm *png");
-    #system("okular $inputDir/$pdf_name&");
+    #system("evince $inputDir/$pdf_name&");
 }
 
 
-#foreach my $t(@array) {print"$t \n ";}
-#print Dumper \%AMB;
 #########################################
 sub help {
@@ -306 +301 @@
 DESCRIPTION:
 
-      Script to pplot BNC's PPP logfiles
+      Script to plot BNC's SP3 Comparison Results
 
 EOI_HILFE