Changeset 9864 in ntrip


Ignore:
Timestamp:
Oct 31, 2022, 9:49:03 PM (4 weeks ago)
Author:
stuerze
Message:

minor changes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/scripts/compareSP3Plot.pl

    r9836 r9864  
    11#!/usr/bin/env perl
    22# ========================================================================
    3 # cmbPlot.pl
     3# compareSP3Plot.pl
    44# ========================================================================
    55#
     
    2525use File::Basename;
    2626use Date::Manip;
     27use Time::Piece 1.30;
    2728use List::MoreUtils qw( minmax );
    2829use experimental 'smartmatch';
    29 
    30 use Gps_Date; # Bernese
    3130
    3231# -----------------------------------------------------------------------------
     
    6160# generate data sets for gnuplot
    6261# -----------------------------------------------------------------------------
    63 my $old_epochSec = 0;
    64 my $epochSec = 0;
    65 my $epochDiff = 0;
     62my $epoSec = 0;
    6663# for pdf generation
    6764my ($png, $page, $headline, $headline_text);
     
    7370    my (%COMPARISON);
    7471    my (@array);
    75     my ($sys, $sat, $epo, $mjd, $radial, $along, $out, $drange, $clk, $clkRed);
     72    my ($sys, $sat, $epo, $radial, $along, $out, $drange, $clk, $clkRed);
    7673    # -----------------------------------------------------------------------------
    7774    # create pdf for plot results
     
    8986        if ( $_ !~ /^!/) { #print "$_\n";
    9087            @array = split(/\s+/, $_);
    91             $mjd = $array[0];
    92             $epo  = gps_date("-mjd","$mjd","-o %Y-%m-%d_%H:%M:%S");
    93             $sys = substr $array[1], 0, 1;
     88            $epo   = $array[0];   #print "$epo\n";
     89            $epo   =~s/(\.\d*)//; #print "$epo\n";
     90            $sys    = substr $array[1], 0, 1;
    9491            $sat    = $array[1];#print "$sat\n";
    9592            $radial = $array[2];#print "$radial\n";
     
    213210    if (("ALL" ~~ @plotTypes ) || ("CLK" ~~ @plotTypes)) {
    214211        print "Plot CLOCK and CLOCK REDUCED BY RADIAL COMPONENT ...\n";
    215         my $yrange = 10.0;
     212        my $yrange = 0.6;
    216213        $page = $pdf->page(); $page->mediabox('A4');
    217214        $headline = sprintf("Clock residuals and reduced clock residuals(%s)", $inputFilename);
     
    229226            my $chart_c = Chart::Gnuplot->new(
    230227                          output => $pngName_c, terminal => 'png',title => "Clock only",
    231                           ylabel => "Clock Residuals [m]",# yrange => ["-$yrange", "$yrange"],
     228                          ylabel => "Clock Residuals [m]", yrange => ["-$yrange", "$yrange"],
    232229                          xlabel => "Time [h]", timeaxis => 'x', xtics => {labelfmt => '%H:%M', rotate => '-270',},
    233230                          legend => { position => "outside right",},
     
    235232            my $chart_r = Chart::Gnuplot->new(
    236233                          output => $pngName_r, terminal => 'png',title => "Clock minus radial component",
    237                           ylabel => "Clock Residuals [m]", #yrange => ["-$yrange", "$yrange"],
     234                          ylabel => "Clock Residuals [m]", yrange => ["-$yrange", "$yrange"],
    238235                          xlabel => "Time [h]", timeaxis => 'x', xtics => {labelfmt => '%H:%M', rotate => '-270',},
    239236                          legend => { position => "outside right",},
     
    286283    $pdf->end();
    287284    system("rm *png");
    288     #system("okular $inputDir/$pdf_name&");
     285    #system("evince $inputDir/$pdf_name&");
    289286}
    290287
    291288
    292 #foreach my $t(@array) {print"$t \n ";}
    293 #print Dumper \%AMB;
    294289#########################################
    295290sub help {
     
    306301DESCRIPTION:
    307302
    308       Script to pplot BNC's PPP logfiles
     303      Script to plot BNC's SP3 Comparison Results
    309304
    310305EOI_HILFE
Note: See TracChangeset for help on using the changeset viewer.