Changeset 10383 in ntrip for trunk/BNC/scripts/Bnc.pm
- Timestamp:
- Mar 11, 2024, 12:52:02 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/scripts/Bnc.pm
r10259 r10383 493 493 my $epochSec = 0; 494 494 my $epochDiff = 0; 495 my ( @hlp, @N, @E, @U, %SATNUM, @TRP s, @CLKs, @OFF_GLOs, @OFF_GALs, @OFF_BDSs);496 my ( @EPOCH s, @EPOCHs_CLK, @EPOCHs_OFF_GLO, @EPOCHs_OFF_GAL, @EPOCHs_OFF_BDS );495 my ( @hlp, @N, @E, @U, %SATNUM, @TRP, @CLK, @OFF_GLO, @OFF_GAL, @OFF_BDS); 496 my ( @EPOCH, @EPOCH_CLK, @EPOCH_OFF_GLO, @EPOCH_OFF_GAL, @EPOCH_OFF_BDS ); 497 497 my ( %AMB, %RES, %ELE, %ION, %BIA ); 498 498 my ( $station, $lki, $sys, $sat, $amb ); … … 530 530 531 531 if ( $ln =~ /\bdN\b/ ) { 532 push ( @EPOCH s, $epochSec ); # besser $epo ?532 push ( @EPOCH, $epochSec ); # besser $epo ? 533 533 $old_epochSec = $epochSec; 534 534 535 #2015-08-27_13:59:50.000 DIEP1 X = 3842152.9054 +- 0.0242 Y = 563402.0331 +- 0.0176 Z = 5042888.5182 +- 0.0319 dN = 0.0130 +- 0.0193 dE = -0.0032 +- 0.0178 dU = -0.0248 +- 0.0349535 #2015-08-27_13:59:50.000 DIEP1 X = 3842152.9054 +- 0.0242 Y = 563402.0331 +- 0.0176 Z = 5042888.5182 +- 0.0319 dN = 0.0130 +- 0.0193 dE = -0.0032 +- 0.0178 dU = -0.0248 +- 0.0349 536 536 $station = $hlp[1]; 537 537 … … 585 585 } 586 586 elsif ( $ln =~ /\bREC_CLK\b/ ) { 587 push ( @EPOCH s_CLK, $epochSec );588 push ( @CLK s, $hlp[2] + $hlp[3] );587 push ( @EPOCH_CLK, $epochSec ); 588 push ( @CLK, $hlp[2] + $hlp[3] ); 589 589 } 590 590 elsif ( $ln =~ /\bOFF_GLO\b/ ) { # 2015-08... OFF_GLO 52.6806 -3.8042 +- 9.0077 591 push ( @EPOCH s_OFF_GLO, $epochSec );592 push ( @OFF_GLO s, $hlp[2] + $hlp[3] );591 push ( @EPOCH_OFF_GLO, $epochSec ); 592 push ( @OFF_GLO, $hlp[2] + $hlp[3] ); 593 593 } 594 594 elsif ( $ln =~ /\bOFF_GAL\b/ ) { # 2015-08... OFF_GAL 52.6806 -3.8042 +- 9.0077 595 push ( @EPOCH s_OFF_GAL, $epochSec );596 push ( @OFF_GAL s, $hlp[2] + $hlp[3] );595 push ( @EPOCH_OFF_GAL, $epochSec ); 596 push ( @OFF_GAL, $hlp[2] + $hlp[3] ); 597 597 } 598 598 elsif ( $ln =~ /\bOFF_BDS\b/ ) { # 2015-08... OFF_BDS 52.6806 -3.8042 +- 9.0077 599 push ( @EPOCH s_OFF_BDS, $epochSec );600 push ( @OFF_BDS s, $hlp[2] + $hlp[3] );599 push ( @EPOCH_OFF_BDS, $epochSec ); 600 push ( @OFF_BDS, $hlp[2] + $hlp[3] ); 601 601 } 602 602 elsif ( $ln =~ /\bSATNUM\b/ ) { # 2015-09... SATNUM G 8 … … 604 604 } 605 605 elsif ( $ln =~ /\bTRP\b/ ) { # 2015-08... TRP 2.3803 +0.1009 +- 0.0324 606 push ( @TRP s, $hlp[2] + $hlp[3] );606 push ( @TRP, $hlp[2] + $hlp[3] ); 607 607 } 608 608 } # ----- next line ----- … … 611 611 close $fh; 612 612 613 my $nof_epochs = scalar @EPOCH s;613 my $nof_epochs = scalar @EPOCH; 614 614 DEBUG( "$station: epochs:$nof_epochs, North displac.: " 615 615 . scalar @N … … 618 618 . ", Up displac.: " 619 619 . scalar @U 620 . ", TRP s:"621 . scalar @TRP s620 . ", TRP:" 621 . scalar @TRP 622 622 ); 623 if ( $nof_epochs != scalar @N ) 624 if ( $nof_epochs != scalar @TRP s ) { LOGDIE "number of epochs and TRPsnot equal\n" }625 if ( @CLK s && scalar @EPOCHs_CLK != scalar @CLKs ) { LOGDIE "number of epochs and CLKsnot equal\n" }626 if ( @OFF_GLO s && scalar @EPOCHs_OFF_GLO != scalar @OFF_GLOs ) { LOGDIE "number of epochs and OFF_GLOsnot equal\n" }627 if ( @OFF_GAL s && scalar @EPOCHs_OFF_GAL != scalar @OFF_GALs ) { LOGDIE "number of epochs and OFF_GALsnot equal\n" }628 if ( @OFF_BDS s && scalar @EPOCHs_OFF_BDS != scalar @OFF_BDSs ) { LOGDIE "number of epochs and OFF_BDSsnot equal\n" }623 if ( $nof_epochs != scalar @N ) { LOGDIE "number of epochs and residuals not equal\n" } 624 if ( $nof_epochs != scalar @TRP ) { LOGDIE "number of epochs and TRP not equal\n" } 625 if ( @CLK && scalar @EPOCH_CLK != scalar @CLK ) { LOGDIE "number of epochs and CLK not equal\n" } 626 if ( @OFF_GLO && scalar @EPOCH_OFF_GLO != scalar @OFF_GLO ) { LOGDIE "number of epochs and OFF_GLO not equal\n" } 627 if ( @OFF_GAL && scalar @EPOCH_OFF_GAL != scalar @OFF_GAL ) { LOGDIE "number of epochs and OFF_GAL not equal\n" } 628 if ( @OFF_BDS && scalar @EPOCH_OFF_BDS != scalar @OFF_BDS ) { LOGDIE "number of epochs and OFF_BDS not equal\n" } 629 629 630 630 if ( !$station ) { WARN "could not grep stationname from file: $file\n" } 631 631 632 632 my %data = ( 633 EPOCH S => \@EPOCHs,634 N 635 E 636 U 637 SATNUM 638 TRP s => \@TRPs,639 CLK s => \@CLKs,640 OFF_GLO s => \@OFF_GLOs,641 OFF_GAL s => \@OFF_GALs,642 OFF_BDS s => \@OFF_BDSs,643 RES 644 AMB 645 ELE 646 ION 647 BIA 633 EPOCH => \@EPOCH, 634 N => \@N, 635 E => \@E, 636 U => \@U, 637 SATNUM => \%SATNUM, 638 TRP => \@TRP, 639 CLK => \@CLK, 640 OFF_GLO => \@OFF_GLO, 641 OFF_GAL => \@OFF_GAL, 642 OFF_BDS => \@OFF_BDS, 643 RES => \%RES, 644 AMB => \%AMB, 645 ELE => \%ELE, 646 ION => \%ION, 647 BIA => \%BIA, 648 648 ); 649 649
Note:
See TracChangeset
for help on using the changeset viewer.