Long-arc LARES analysis
- import and screen observations (screening orbit according to V_MINUS, V_PLUS)
- generate apriori 1d-orbit
- generate satellite(-group) specific 1d-NEQ
- long-arc stacking of seven 1d-NEQs and satellite(-group) specific 7d-solution
- reconstruct 1d-orbits and compute residuals
- perform full validation of 7d-solution
- generate orbit prediction

PID
SCRIPT
PARAM1
PARAM2
PARAM3
PARAM4
PARAM5
PARAM6
Input
Output
Copy required files and adapt MEANPOLE in ERP
001
L_COPY_L1

$(NPDIR)/$(NPDAT)_$yyyymmdd.npt,np2

$(STAINF).STA

$(STAINF).ABB

$(BLQ).BLQ

$(CRDINF).CRD,VEL,PSD

$(CRDINF).FIX

$(CRDINF).SIG

$(RGBINF)

$(COM).SLR

$(ATL)

$(PCVINF)

$(SINEX).SKL

$(EOP).ERP

AOD1B$YD+-.AOD

Extract ERP-file covering screening arc
003
U_POLUPD__
PRED
$(EOP).ERP
$(L0)$(SAT)$YD+0.ERP
Generate FIX-file for RGB estimation
006
U_RGBUPD_V

$(RGBINF).FIX_intervals

$(RGBINF).FIX

RGB$(SAT)$YD+0.FIX
Choice of apriori orbit (C7Dd, L7Dd, EXDd, CPF)
009
L_CHKORB_2
015
010
Case CPF: create PRE
010
U_GETCPF_V
L60$YD+-.CPF
011
U_CPFSP3_V
L60$YD+-.CPF
L60_$YD+-.PRE
012
U_CCPREO_V
L60_$YD+-.PRE
AP$(SAT)$YD+-.PRE
Generate screening orbit (depending on CHKORBDyyddd.FLG)
015
U_ORBGEN_S

$(L0)$(SAT)$YD+0.PRE

or according to

CHKORB$(SAT)$YD+0.FLG

$(L0)$(SAT)$YD+0.STD

$(L0)$(SAT)$YD+0.RPR

Import observation and meteo files
018
U_CHKORX_V
sel$YD+0_$(SAT).orx
020
U_NPTOBS_V
sel$YD+0_$(SAT).orx

????$YD+0$(SAT).RZH,RZO

????$YD+00.MET

Station selection
022
U_STASEL_V
A
F
D
stalst.$YD+0$(SAT)
Create a priori coordinate file
032
U_COOVEL__
$(CRDINF).CRD,VEL,PSD
$(ID)$YD+0.CRD
Observation screening
040
U_SLRED0_V

GPSEST

????$YD+0$(SAT).RZH,RZO

????$YD+00.MET

$(ID)$YD+0.CRD

$(L0)$(SAT)$YD+0.STD

$(L0)$(SAT)$YD+0.RPR

QLRSUM

$(SAT)$YD+00_$(IT).RES

SATMRK

$(SAT)$YD+00R$(IT).EDT

SLREDT_0.$YD+0$(SAT)

GPSEST

$(SAT)$YD+00_$(IT).OUT,RES






QLRSUM

$(SAT)$YD+00R$(IT).SUM,EDT

SATMRK

????$YD+0$(SAT).RZH,RZO

045
U_SLRED2_V
20

GPSEST

????$YD+0$(SAT).RZH,RZO

????$YD+00.MET

$(ID)$YD+0.CRD

$(L0)$(SAT)$YD+0.STD

$(L0)$(SAT)$YD+0.RPR

QLRSUM

$(SAT)$YD+0_$(IT).RES

SATMRK

$(SAT)$YD+0R$(IT).EDT

SLREDT.$YD+0$(SAT)

LAST_EDT$(SAT)$YD+0.FLG

GPSEST

$(SAT)$YD+0_$(IT).OUT,RES




QLRSUM

$(SAT)$YD+0R$(IT).SUM,EDT

SATMRK

????$YD+0$(SAT).RZH,RZO

Cut out 1d observation files
050
L_OBSSPL_V
1
h
????$YD+0$(SAT).RZH,RZO
????$YD+0$(SAT).RZH,RZO
055
L_OBSSPL_V
????$YD+0$(SAT).RZH,RZO
????$YD+0$(SAT).RZH,RZO
Generate 1d apriori orbit
060
U_ORBGEN_S

$(L0)$(SAT)$YD+0.PRE

or according to

CHKORB$(SAT)$YD+0.FLG

$(L0)$(SAT)$(L1)$YD+0.STD

$(L0)$(SAT)$(L1)$YD+0.RPR

$(L0)$(SAT)$(L1)$YD+0.GRV

Generate 1d normal equation (solution skipped)
070
U_GPSEST_S

????$YD+0$(SAT).RZH,RZO

????$YD+00.MET

$(ID)$YD+0.CRD

$(L0)$(SAT)$(L1)$YD+0.STD

$(L0)$(SAT)$(L1)$YD+0.RPR

$(L1)$(SAT)$YD+0.NQ0

$(L1)$(SAT)$YD+0.RP2

Stack seven daily LARES NEQs (long-arc) and generate weekly solution (if requested)
100
U_SKIPSL_V
900
110
SKIPSOL_$(SAT).FLG
LARES 7d-NEQ for combination, dynamic orbit parameters eliminated after stacking
110
L_ADDNEQ_S
-6
0

neqsel_1d_$YD+0.txt

neqsel_1d_$YD+0.txt: $(L1)$(SAT)YD+-.NQ0

$(AP)$(SAT)$(L1)$YD+-.STD,RPR,RP2

$(ID)$YD+0.CRD

$(EOP).ERP

$(LX)$(SAT)$(SID)$YD+0.NQ0

$(LX)$(SAT)$(SID)$YD+0.CRD

$(LX)$(SAT)$(SID)$YD+0.ERP

$(LX)$(SAT)$(SID)$YD+0.SLR

$(LX)$(SAT)$(SID)$YD+0.OUT

empirical 1/rev-cross-track parameters fixed to 0
111
L_ADDNEQ_S
-6
0
w

neqsel_1d_$YD+0.txt: $(L1)$(SAT)YD+-.NQ0

$(AP)$(SAT)$(L1)$YD+-.STD,RPR,RP2

$(ID)$YD+0.CRD

$(EOP).ERP

$(LX)$(SAT)$(SID)$(SID2)$YD+0.NQ0

$(LX)$(SAT)$(SID)$(SID2)$YD+0.CRD

$(LX)$(SAT)$(SID)$(SID2)$YD+0.ERP

$(LX)$(SAT)$(SID)$(SID2)$YD+0.SLR

$(LX)$(SAT)$(SID)$(SID2)$YD+0.OUT

empirical 1/rev-cross-track and -along-track parameters fixed to 0
112
L_ADDNEQ_S
-6
0
s

neqsel_1d_$YD+0.txt: $(L1)$(SAT)YD+-.NQ0

$(AP)$(SAT)$(L1)$YD+-.STD,RPR,RP2

$(ID)$YD+0.CRD

$(EOP).ERP

$(LX)$(SAT)$(SID)$(SID2)$YD+0.NQ0

$(LX)$(SAT)$(SID)$(SID2)$YD+0.CRD

$(LX)$(SAT)$(SID)$(SID2)$YD+0.ERP

$(LX)$(SAT)$(SID)$(SID2)$YD+0.SLR

$(LX)$(SAT)$(SID)$(SID2)$YD+0.OUT

Orbit reconstruction (if requested)
115
L_SKIPOR_V
900
120
SKIPORB_$(SAT).FLG
Generate ELE-files, long-arc solution resubstituted
120
L_ADDNEQ_S
-6
0

$(L1)$(SAT)YD+-.NQ0

$(LX)$(SAT)$(SID)$YD+0.CRD

$(LX)$(SAT)$(SID)$YD+0.ERP

$(L1)$(SAT)$(SID)o$YD+-.ELE

$(LX)$(SAT)$(SID)o$YD+0.SLR

$(LX)$(SAT)$(SID)o$YD+0.OUT

Update SLR-file (re-estimated in second ADDNEQ2 run)
125
L_SNXSLR_S

$(LX)$(SAT)$(SID)o$YD+0.SLR

$(COM).SLR

$(LX)$(SAT)$(SID)o$YD+0M.SLR
Compute 1d-orbits and residuals: days MJD - 6 ... MJD
130
L_ORBGEN_T
_
6543210
131
L_ORBGEN_S
PARALLEL=130

$(L1)$(SAT)$(SID)o$(DAT).ELE

$(LX)$(SAT)$(SID)$YD+0.ERP

$(L1)$(SAT)$(SID)$(DAT).STD
140
L_STDPRE_T
_
6543210
141
L_STDPRE_S
PARALLEL=140

$(L1)$(SAT)$(SID)$(DAT).STD

$(LX)$(SAT)$(SID)$YD+0.ERP

$(LX)$(SAT)$(SID)$YD+0.CRD

$(L1)$(SAT)$(SID)$(DAT).PRE
150
L_GPSEST_T
_
6543210
151
L_GPSEST_S
PARALLEL=150

????$(DAT)$(SAT).RZH,RZO

????$(DAT)$(SAT).MET

$(LX)$(SAT)$(SID)$YD+0.CRD

$(L1)$(SAT)$(SID)$(DAT).STD

$(LX)$(SAT)$(SID)$YD+0.ERP

$(LX)$(SAT)$(SID)o$YD+0M.SLR

$(L1)$(SAT)$YD+0$(SID)$(DAT).RES,OUT
160
L_QLRSUM_T
_
6543210
161
L_QLRSUM_S
PARALLEL=160

$(L1)$(SAT)$YD+0$(SID)$(DAT).RES

$(L1)$(SAT)$(SID)$(DAT).STD

$(LX)$(SAT)$(SID)$YD+0.CRD

$(LX)$(SAT)$(SID)$YD+0.ERP

RES$YD+0$(L1)$(SAT)$(SID)$(DAT).OUT,SUM,PLT
Perform full validation (if requested)
200
L_SKIPVL_V
300
202
SKIPVAL_$(SAT).FLG
Generate list of daily PRE-files
202
L_CHKORB_S
-6
0
L1
sel$(SID)$YD+0_$(SAT).orb
Fit 7-day orbit (7 arcs) to seven continuous 1-day PRE-Files
205
L_ORBGEN_S
-6
0
0

sel${sid}$YD+0_$(SAT).orb: $(L1)$(SAT)${sid}$YD+-.PRE

$(LX)$(SAT)$(SID)$YD+0.ERP

$(LX)$(SAT)$(SID)$YD+0.STD,RPR
Generate list of 1d observation files for 7d solution
207
L_CHKOBS_S
-6
0
sel$YD+0_$(SAT).obs
Merge seven 1d meteo files
208
L_MERGEM_V
-6
0
sel$YD+0_$(SAT).obs: ????$YD+-$(SAT).RZH,RZO
????$YD+0X.MET
Generate residuals of 7-day orbit
210
U_GPSEST_S

sel$YD+0_$(SAT).obs: ????$YD+-$(SAT).RZH,RZO

????$YD+0X.MET

$(LX)$(SAT)$(SID)$YD+0.CRD

$(LX)$(SAT)$(SID)$YD+0.STD

$(LX)$(SAT)$(SID)$YD+0.ERP

$(LX)$(SAT)$(SID)o$YD+0.SLR

$(LX)$(SAT)$(SID)$YD+0.RES

$(LX)$(SAT)$(SID)RES$YD+0.OUT

215
U_QLRSUM_S

$(LX)$(SAT)$(SID)$YD+0.RES

$(LX)$(SAT)$(SID)$YD+0.STD

$(LX)$(SAT)$(SID)$YD+0.ERP

$(LX)$(SAT)$(SID)$YD+0.CRD

RES$YD+0$(LX)$(SAT)$(SID).OUT,SUM,PLT
Generate list of used stations
217
L_STALST_V
L1
$(LX)$(SAT)$(SID)$YD+0.OUT
stalst.$YD+0$(SAT)
Update FIX file
220
L_FIXUPD_C
LXS
ONE

$(CRDINF).FIX

RES$YD+0$(LX)$(SAT)$(SID).OUT

$(ID)$(SAT)$YD+0.FIX
Compare orbit with previous prediction
230
L_ORBCMP_S
1
L1
EXT

$(L1)$(SAT)$(SID)$YD+0.PRE

$(EXT)$(SAT)$(SID)$YD-1.PRE

$(LX)$(SAT)$(SID)$YD+0.ERP

$(L1)$(SAT)$(SID)$YD+0.STD

CMP$YD+0M1_$(L1)$(SAT)$(SID)_$(EXT)$(SAT)$(SID).LST,OUT
Compare station coordinates
240
U_HELMR1_S

$(ID)$YD+0.CRD

$(LX)$(SAT)$(SID)$YD+0.CRD

HLM$YD+0$(LX)$(SAT)$(SID)_ALL.OUT
250
U_HELMR1_S

$(ID)$YD+0.CRD

$(LX)$(SAT)$(SID)$YD+0.CRD

$(ID)$(SAT)$YD+0.FIX

HLM$YD+0$(LX)$(SAT)$(SID)_FIX.OUT
260
L_COMPAR_2
LX
10
$(LX)$(SAT)$(SID)?????.CRD
COMP$YD+0_$(LX)$(SAT)$(SID).PLT,OUT
Generate orbit prediction (if requested)
300
L_SKIPPR_V
900
302
SKIPPRED_$(SAT).FLG
Extend ERP-file beyond solution period with apriori values
302
L_PMERGE_S
L0
LX

$(L0)$(SAT)$YD+0.ERP

$(LX)$(SAT)$(SID)$YD+0.ERP

$(EXT)$(SAT)$(SID)$YD+0.ERP
Fit seven continuous PRE-files by one arc (including propagation: V_PRED)
305
L_ORBGEN_S
-6
0
PRED

sel$(SID)$YD+0_$(SAT).orb: $(L1)$(SAT)$(SID)$YD+-.PRE

$(EXT)$(SAT)$(SID)$YD+0.ERP

$(EXT)$(SAT)$(SID)$YD+0.STD
Convert predicted orbits: STD -> PRE
310
L_STDPRE_S
-6
0
PRED

$(EXT)$(SAT)$(SID)$YD+0.STD

$(EXT)$(SAT)$(SID)$YD+0.ERP

$(LX)$(SAT)$(SID)$YD+0.CRD

$(EXT)$(SAT)$(SID)$YD+0.PRE
End of BPE
900
L_ILRS_PRC
ALL
LX
L1

CHKORB$(SAT)$YD+0.FLG

SLREDT_0.$YD+0$(SAT)

SLREDT.$YD+0$(SAT)

$(SAT)$YD+00R0.OUT

LAST_EDT$(SAT)$YD+0.FLG

$(SAT)$YD+0R${nIter}.OUT

$(LX)$(SAT)$(SID)$YD+0.OUT

RES$YD+0$(LX)$(SAT)$(SID).OUT

CMP$YD+0M1_$(L1)$(SAT)$(SID)_$(EXT)$(SAT)$(SID).OUT

HLM$YD+0$(LX)$(SAT)$(SID)_ALL.OUT

HLM$YD+0$(LX)$(SAT)$(SID)_FIX.OUT

COMP$YD+0_$(LX)$(SAT)$(SID).OUT

ILRS7D_$(LX)$(SAT)$(SID).$Y_$+0
990
BPE_CLN2
35
999
BPE_BYE
-25
SLR