import to data.dat
#
# import <fileformat> <filename>
#
# file-format string should be composed from the following letters:
#    t ... time
#    o ... observed
#    a ... adjusted
#    c ... calculated
#    w ... point weight 
#    u ... point error
#    r ... residuals(observed)
#    R ... residuals(adjusted)
#    1 ... attribute #1
#    2 ... attribute #2
#    3 ... attribute #3
#    4 ... attribute #4
#    5 ... weight(attribute #1)
#    6 ... weight(attribute #2)
#    7 ... weight(attribute #3)
#    8 ... weight(attribute #4)
#    i ... ignore
#
# note: 
#    if you import point errors/weights ('u', 'w') the weights 
#    will be automatically activated for use in calculations!
#

loadproject myproject.p04
#
# loadproject <filename>
#
#      load project from file <filename>
#

fourier 0 nyq o y
#
# fourier <from> <to> <data> <compact-mode>
#         
#     <from> ... lower frequency limit
#     <to>   ... higher frequency limit or 'nyq' (=computes spectrum up to Nyquist freq)
#     <data> ... choose one of these letters:
#                o ... observed
#                a ... adjusted
#                r ... residuals(observed)
#                R ... residuals(adjusted)
#                1 ... spectral window 1
#                2 ... spectral window 2
#     <compact-mode> ... choose y or n
#

fit o 1
fit o 2
#
# fit <data> <mode>
#
#     <data> ... choose one of these letters:
#                o ... observed
#                a ... adjusted
#     <mode> ... calculation mode
#     	     	 1 ... improve amplitudes & phases
#     	     	 2 ... improve frequencies, amplitudes & phases
#

addharmonics 15
#
# addharmonics <maxnumber>
#
# Adds <maxnumber> harmonics of the frequency that was found 
# in the previous step. This command only defines the additional frequencies, 
# no least squares fit is performed at this stage.
#

addsubharmonic 2
#
# addsubharmonic <number>
#
# Add a specific subharmonic (1/<number>) of the frequency that was found 
# in the previous step. This command only registers the subharmonic, 
# no least squares fit is performed at this stage.
#

savefreqs data.per
#
# savefreqs <file>
#
# save the current frequencies, amplitudes & phases in <file>
#

importfreqs data.per
#
# importfreqs <file>
#
# import frequencies, amplitudes & phases from <file>
#

savefourier 0 nyq r n data.fou
#
# savefourier <from> <to> <data> <compact-mode> <file>
#
#     compute fourier spectrum and save it in <file>
#     <from> ... any number
#     <to>   ... any number or 'nyq' (=computes spectrum up to Nyquist freq)
#     <data> ... choose one of these letters:
#                o ... observed
#                a ... adjusted
#                r ... residuals(observed)
#                R ... residuals(adjusted)
#                1 ... spectral window 1
#                2 ... spectral window 2
#     <compact-mode> ... choose y or n
#
#

saveproject myproject.p04
#
# saveproject <filename>
#
#      save project in <filename>
#

createartificial d 1000.0 1500.0 1.0 0.0 artificial.dat a
# 
# createartificial d <starttime> <endtime> <step> <leading> <outputfile> <appendmode>
#
#      generate artificial data based on the current frequency solution
#      <starttime>  ... desired start time 
#      <endtime>    ... desired end time 
#      <step>       ... time step
#      <leading>    ... additional leading/trailing time
#      <outputfile> ... output file containing the generated data
#      <appendmode> ... choose one of these letters:
#                       a ... append to outputfile if it already exists
#                       n ... create new file
# 
# createartificial f <inputfile> <outputfile> <appendmode>
#
#      generate artificial data based on the current frequency solution
#      <inputfile>  ... input file with one column containing the times used for data generation
#      <outputfile> ... output file containing the generated data
#      <appendmode> ... choose one of these letters:
#                       a ... append to outputfile if it already exists
#                       n ... create new file
#
#

calcsnr 2.0 0 r snr.txt
#
# calcsnr <box size> <step rate quality> <data> <outputfile> 
#
#     compute signal to noise ratio for active frequencies
#
#     <box size> ... size of the frequency range to be used for the noise calculation: 
#                    [frequency-box/2, frequency+box/2]
#     <step rate quality> ... defines the step rate
#                0 ... high
#                1 ... medium
#                2 ... low
#     <data> ... choose one of these letters:
#                o ... observed
#                a ... adjusted
#                r ... residuals(observed)
#                R ... residuals(adjusted)
#                1 ... spectral window 1
#                2 ... spectral window 2# 
#      <outputfile> ... output file containing the results
#


calcerrors a uncertainties.txt
calcerrors ls uncertainties.txt o 2 t0
calcerrors mc uncertainties.txt 1000 t0 r
#
# calcerrors a <outputfile> 
#
#     calculate analytically computed uncertainties
#     <outputfile> ... output file containing the results
#
# calcerrors ls <outputfile> <data> <mode>
# calcerrors ls <outputfile> <data> <mode> t0
#
#     calculate uncertaintes derived from least-squares fit
#     <outputfile> ... output file containing the results
#     <data> ... choose one of these letters:
#                o ... observed
#                a ... adjusted
#     <mode> ... calculation mode
#     	     	 1 ... improve amplitudes & phases
#     	     	 2 ... improve frequencies, amplitudes & phases
#     t0 ... use original zeropoint in time
#
# calcerrors mc <outputfile> <num_experiments> 
# calcerrors mc <outputfile> <num_experiments> t0
# calcerrors mc <outputfile> <num_experiments> r
#
#     calculate uncertaintes derived by monte carlo simulation
#     <outputfile> ... output file containing the results
#     <num_experiments> ... number of monte carlo experiments
#     t0 ... use original zeropoint in time
#     r ... use system time to initialize random generator
#


exit
#
# exit Period04
#