Home    Search    Reference Manuals    Return    Track Group Documentation





Users Manual for Program TRC_IPSD


Table of Contents

Introduction
Input data commands
Examples


Introduction

Program TRC_IPSD creates track irregularities which fulfills a given PSD-spectra. In command RNDFIG the user can choose a random number which will be used when selecting the phase angles between different waves in the track irregularity spectra. As an alternative the user can in command INFIL feed an old existent track into program TRC_IPSD and the old track will work as a base when generating the new track which fulfills required PSD-spectra.

See also:

TRC_PSDPL Plot the PSD-spectra of the track
PSD Create PSD-spectra of a column in the track file from the autocorrelation
PSD_FOUR Create PSD-spectra of a column in the track file from the Fourier spectra

Input data commands

Input data is read in free format, valid separators between the input values are <space>, <comma>, <tab>, <equal sign> or <carriage return>. The commands can be written both in lower and upper case letters. The operation of the program is controlled by the commands described below; some of the commands also need arguments.


DRAW_ERRI_LEVELS
Draw high and low ERRI-levels in the PSD-diagrams.
The lines are drawn in a turquoise color. DRAW_ERRI_LEVELS can be set equal to "y" for yes or "n" for no.
Declared= Character*1    Default= "y"

GAUGE
Average gauge of the generated track.
Declared= Real*8    Default= 1.435 [m]

INFIL
Existing track data file written in trac, trax or trax_wdesign format. Phase angles between different waves from an existing track will be used when creating the new track UTFIL. In case you do not have measured track irregularities, phase angles can be chosen randomly if you set INFIL equal to "no".
Declared= Character*132    Default= "no"

PSD_LAT
Definition of PSD-spectra in lateral direction. The PSD-spectra is given as node coordinates in the following way: Y1,X1, Y2,X2, Y3,,, Yn,Xn. The Y values are amplitudes given in the unit [m*m/1/m]. The X values are spatial frequencies and the values are given in [1/m]. Whole range from infinite long waves down to the shortest possible waves must be defined in the spectra, i.e. X1 must start at 0. and Xn must be greater than 1/(2*DX)
Declared= Real*8(100)    Default= 100*0.

PSD_VERT
Definition of PSD-spectra in vertical direction. The PSD-spectra is given as node coordinates in the following way: Y1,X1, Y2,X2, Y3,,, Yn,Xn. The Y values are amplitudes given in the unit [m*m/1/m]. The X values are spatial frequencies and the values are given in [1/m]. Whole range from infinite long waves down to the shortest possible waves must be defined in the spectra, i.e. X1 must start at 0. and Xn must be greater than 1/(2*DX)
Declared= Real*8(100)    Default= 100*0.

PSD_GAUGE
Definition of PSD-spectra in gauge direction. The PSD-spectra is given as node coordinates in the following way: Y1,X1, Y2,X2, Y3,,, Yn,Xn. The Y values are amplitudes given in the unit [m*m/1/m]. The X values are spatial frequencies and the values are given in [1/m]. Whole range from infinite long waves down to the shortest possible waves must be defined in the spectra, i.e. X1 must start at 0. and Xn must be greater than 1/(2*DX)
Declared= Real*8(100)    Default= 100*0.

PSD_FI
Definition of PSD-spectra in cant direction. The PSD-spectra is given as node coordinates in the following way: Y1,X1, Y2,X2, Y3,,, Yn,Xn. The Y values are amplitudes given in the unit [m*m/1/m]. The X values are spatial frequencies and the values are given in [1/m]. Whole range from infinite long waves down to the shortest possible waves must be defined in the spectra, i.e. X1 must start at 0. and Xn must be greater than 1/(2*DX)
Declared= Real*8(100)    Default= 100*0.

RNDFIG
Initial value when calculating the random series numbers.
The random numbers are used as phase angles between different waves in the track irregularity spectra. RNDFIG must be greater than 0.001 and less than 0.499.
Declared= Real*8    Default= 0.1

UTFIL
File containing the output results.
The file contains a new track fulfilling the spectra given in PSD_LAT, PSD_VERT, PSD_GAUGE and PSD_FI. If the file extension of UTFIL equals "trax_wdesign", the output file will be written in trax_wdesign format. Otherwise the output will be written in trax-format.
Declared= Character*132    Default= 'trc_ipsd.trax'

XSTART
Start coordinate when reading track irregularity data from INFIL.
The track read from INFIL should start at least 250[m] before XSTART to avoid disturbances from the end. If INFILE starts too close to XSTART a warning message will be written.
Declared= Real*8   
Default= -1.E99 (Automatically set XSTART 250[m] after beginning of INFIL)

XSTOP
End coordinate when reading track irregularity data from INFIL.
The track read from INFIL should end at least 250[m] after XSTOP to avoid disturbances from the end. If INFILE starts too close to XSTOP a warning message will be written.
Declared= Real*8   
Default= 1.E99 (Automatically set XSTOP 250[m] before end of INFIL)

DX
Equidistant step between the points in the X-axis.
Declared= Real*8    Default= 0.5 [m]

Examples:


Example #1

Create a PSD spectra with constant amplitude for all frequencies

##
##      Input data for program TRC_IPSD
##

 RNDFIG= 0.1
 INFIL = $gentrc/K0_4a21-n-k.trax  # Read phase angles from an existed trac-file.
 UTFIL = trc_ipsd.trax

 DX    = 0.5
 XSTART=-1.e99  XSTOP= 1.e99  GAUGE= 1.435

#
#        White PSD spectra
#
  PSD_LAT  = 1.E-6,     0,    1.E-6,     1.      # 1[mm] in lateral direction
  PSD_VERT = 1.E-6,     0,    1.E-6,     1.      # 1[mm] in vertical direction
  PSD_GAUGE= 1.E-6,     0,    1.E-6,     1.      # 1[mm] in gauge
  PSD_FI   = .444444e-6 0,    .444444e-6 1.      # 1[mm] in cant



Example #2

Create a PSD spectra according to ERRI QUESTION B176

In the ERRI report B176 the following PSD-spectras was presented:


  Sa=               Aa*0.8246^2              
      (omega^2+0.0206^2) * (omega^2+0.8246^2)

  Sv=               Av*0.8246^2              
      (omega^2+0.0206^2) * (omega^2+0.8246^2)

  Sc=               Av/bo^2*0.8246^2*omega^2                      
      (omega^2+0.0206^2) * (omega^2+0.8246^2) * (omega^2+0.4380^2)

Where:
Sa = PSD-spectra in lateral direction (alignment)
Sv = PSD-spectra in vertical direction
Sc = PSD-spectra of cant (expressed in radians)

The coefficients Aa and Av are as follows:

Low level of irregularity:

Aa = 2.119e-7 [ m rad ]
Av = 4.032e-7 [ m rad ]

High level of irregularity:

Aa = 6.125e-7 [ m rad ]
Av = 1.080e-6 [ m rad ]

##
##      Input data for program TRC_IPSD
##

 RNDFIG= 0.1
 INFIL = $gentrc/K0_4a21-n-k.trax  # Read phase angles from an existed trac-file.
 UTFIL = trc_ipsd.trax

 DX    = 0.5
 XSTART=-1.e99  XSTOP= 1.e99  GAUGE= 1.435

#
#        ERRI B176 level high
#
# Sa= 2*pi*6.125e-7*0.8246**2 / ((2*pi*freq)**2+0.0206**2) / ((2*pi*freq)**2+0.8246**2)   Lateral
# Sv= 2*pi*1.080e-6*0.8246**2 / ((2*pi*freq)**2+0.0206**2) / ((2*pi*freq)**2+0.8246**2)   Vertical
# Sc= 2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*freq)**2 / ((2*pi*freq)**2+0.0206**2) / ((2*pi*freq)**2+0.8246**2) / ((2*pi*freq)**2+0.4380**2)
#
  PSD_LAT  =
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.      )**2+0.0206**2) / ((2*pi*  0.      )**2+0.8246**2)`   0.             # Infinite wave length;  Unit [m^2/(1/m)]
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.003183)**2+0.0206**2) / ((2*pi*  0.003183)**2+0.8246**2)`   0.003183       # 314.16m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.010000)**2+0.0206**2) / ((2*pi*  0.010000)**2+0.8246**2)`   0.010000       # 100.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.020000       #  50.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.050000)**2+0.0206**2) / ((2*pi*  0.050000)**2+0.8246**2)`   0.050000       #  20.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.100000)**2+0.0206**2) / ((2*pi*  0.100000)**2+0.8246**2)`   0.100000       #  10.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.200000)**2+0.0206**2) / ((2*pi*  0.200000)**2+0.8246**2)`   0.200000       #   5.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.500000)**2+0.0206**2) / ((2*pi*  0.500000)**2+0.8246**2)`   0.500000       #   2.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  1.000000)**2+0.0206**2) / ((2*pi*  1.000000)**2+0.8246**2)`   1.000000       #   1.00m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*  2.000000)**2+0.0206**2) / ((2*pi*  2.000000)**2+0.8246**2)`   2.000000       #   0.50m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi* 10.000000)**2+0.0206**2) / ((2*pi* 10.000000)**2+0.8246**2)`  10.000000       #   0.10m waves
     `2*pi*6.125e-7*0.8246**2 / ((2*pi*100.000000)**2+0.0206**2) / ((2*pi*100.000000)**2+0.8246**2)` 100.000000       #   0.01m waves
#
  PSD_VERT =
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.      )**2+0.0206**2) / ((2*pi*  0.      )**2+0.8246**2)`   0.             # Infinite wave length;  Unit [m^2/(1/m)]
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.003183)**2+0.0206**2) / ((2*pi*  0.003183)**2+0.8246**2)`   0.003183       # 314.16m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.010000)**2+0.0206**2) / ((2*pi*  0.010000)**2+0.8246**2)`   0.010000       # 100.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.020000       #  50.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.050000)**2+0.0206**2) / ((2*pi*  0.050000)**2+0.8246**2)`   0.050000       #  20.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.100000)**2+0.0206**2) / ((2*pi*  0.100000)**2+0.8246**2)`   0.100000       #  10.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.200000)**2+0.0206**2) / ((2*pi*  0.200000)**2+0.8246**2)`   0.200000       #   5.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  0.500000)**2+0.0206**2) / ((2*pi*  0.500000)**2+0.8246**2)`   0.500000       #   2.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  1.000000)**2+0.0206**2) / ((2*pi*  1.000000)**2+0.8246**2)`   1.000000       #   1.00m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*  2.000000)**2+0.0206**2) / ((2*pi*  2.000000)**2+0.8246**2)`   2.000000       #   0.50m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi* 10.000000)**2+0.0206**2) / ((2*pi* 10.000000)**2+0.8246**2)`  10.000000       #   0.10m waves
     `2*pi*1.080e-6*0.8246**2 / ((2*pi*100.000000)**2+0.0206**2) / ((2*pi*100.000000)**2+0.8246**2)` 100.000000       #   0.01m waves
#
  PSD_GAUGE=
     `0.05*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.             # Infinite wave length;  Unit [m^2/(1/m)]
     `0.05*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.003183       # 314.16m waves
     `0.05*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.010000       # 100.00m waves
     `0.05*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2)`   0.020000       #  50.00m waves
     `0.25*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.050000)**2+0.0206**2) / ((2*pi*  0.050000)**2+0.8246**2)`   0.050000       #  20.00m waves
     `0.50*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.100000)**2+0.0206**2) / ((2*pi*  0.100000)**2+0.8246**2)`   0.100000       #  10.00m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.200000)**2+0.0206**2) / ((2*pi*  0.200000)**2+0.8246**2)`   0.200000       #   5.00m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*  0.500000)**2+0.0206**2) / ((2*pi*  0.500000)**2+0.8246**2)`   0.500000       #   2.00m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*  1.000000)**2+0.0206**2) / ((2*pi*  1.000000)**2+0.8246**2)`   1.000000       #   1.00m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*  2.000000)**2+0.0206**2) / ((2*pi*  2.000000)**2+0.8246**2)`   2.000000       #   0.50m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi* 10.000000)**2+0.0206**2) / ((2*pi* 10.000000)**2+0.8246**2)`  10.000000       #   0.10m waves
     `1.00*2*pi*6.125e-7*0.8246**2 / ((2*pi*100.000000)**2+0.0206**2) / ((2*pi*100.000000)**2+0.8246**2)` 100.000000       #   0.01m waves
#
  PSD_FI   =
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.      )**2 / ((2*pi*  0.      )**2+0.0206**2) / ((2*pi*  0.      )**2+0.8246**2) / ((2*pi*  0.      )**2+0.4380**2)`   0.             # Infinite wave length;  Unit [rad^2/(1/m)] 
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.010000)**2 / ((2*pi*  0.010000)**2+0.0206**2) / ((2*pi*  0.010000)**2+0.8246**2) / ((2*pi*  0.010000)**2+0.4380**2)`   0.010000       # 100.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.020000)**2 / ((2*pi*  0.020000)**2+0.0206**2) / ((2*pi*  0.020000)**2+0.8246**2) / ((2*pi*  0.020000)**2+0.4380**2)`   0.020000       #  50.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.050000)**2 / ((2*pi*  0.050000)**2+0.0206**2) / ((2*pi*  0.050000)**2+0.8246**2) / ((2*pi*  0.050000)**2+0.4380**2)`   0.050000       #  20.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.100000)**2 / ((2*pi*  0.100000)**2+0.0206**2) / ((2*pi*  0.100000)**2+0.8246**2) / ((2*pi*  0.100000)**2+0.4380**2)`   0.100000       #  10.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.200000)**2 / ((2*pi*  0.200000)**2+0.0206**2) / ((2*pi*  0.200000)**2+0.8246**2) / ((2*pi*  0.200000)**2+0.4380**2)`   0.200000       #   5.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  0.500000)**2 / ((2*pi*  0.500000)**2+0.0206**2) / ((2*pi*  0.500000)**2+0.8246**2) / ((2*pi*  0.500000)**2+0.4380**2)`   0.500000       #   2.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  1.000000)**2 / ((2*pi*  1.000000)**2+0.0206**2) / ((2*pi*  1.000000)**2+0.8246**2) / ((2*pi*  1.000000)**2+0.4380**2)`   1.000000       #   1.00m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*  2.000000)**2 / ((2*pi*  2.000000)**2+0.0206**2) / ((2*pi*  2.000000)**2+0.8246**2) / ((2*pi*  2.000000)**2+0.4380**2)`   2.000000       #   0.50m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi* 10.000000)**2 / ((2*pi* 10.000000)**2+0.0206**2) / ((2*pi* 10.000000)**2+0.8246**2) / ((2*pi* 10.000000)**2+0.4380**2)`  10.000000       #   0.10m waves
     `2*pi*1.080e-6/bo**2 * 0.8246**2 * (2*pi*100.000000)**2 / ((2*pi*100.000000)**2+0.0206**2) / ((2*pi*100.000000)**2+0.8246**2) / ((2*pi*100.000000)**2+0.4380**2)` 100.000000       #   0.01m waves