Geodæsi øvelse 4

 

 

(login på ikos )

setenv DISPLAY ???:0

 

 

Programmet "geocol12" køres interaktivt (se vedlagt eksempel!)

/disk1/cct/geod/E4/geocol12

ex4h.out (outputfil som indeholder geoide højder )

 

/disk1/cct/geod/E4/geocol12

ex4g.out (outputfil som indeholder tyngde anomalier )

 

Plot geoide højder + tyngde anomalier ved hjælp af GMT

/disk1/cct/geod/E4/ovelse4g.gmt konturplot af tyngde anomalier

/disk1/cct/geod/E4/ovelse4h.gmt konturplot af geoide højder

/disk1/cct/geod/E4/ovelse4g_c.gmt "image" af tyngde anomalier

/disk1/cct/geod/E4/ovelse4h_c.gmt "image" af geoide højder

/disk1/cct/geod/E4/run_3.job (tyngde anomalier ved 0 km)

/disk1/cct/geod/E4/run_4.job (tyngde anomalier ved 30 km)

/disk1/cct/geod/E4/ceylon.gmt "image" af tyngde anomalier

ved 0 km og 30 km

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eksempel på anvendelse af "geocol12"

/disk1/cct/geod/E4$ geocol12

GEODETIC COLLOCATION, VERSION 17 JUNE 1997, RELEASE 12a (UNIX)

Wed Oct 10 11:58:01 2001

NOTE THAT THE FUNCTIONALS ARE IN SPHERICAL APPROXIMATION

MEAN RADIUS = RE = 6371 KM AND MEAN GRAVITY 981 KGAL USED.

MAX NUMBER OF OBS= 5600, MAX NUMBER OF PARAMETERS=239

MAX NUMBER OF OBS IN GIVEN REF. FRAME =5600

SIZE OF NORMAL EQ. BLOCKS=79600, SIZE OF POT.COFF. BLOCK= 130322

INTERACTIVE INPUT (T/F)

t

INPUT: LTRAN, TRUE IF NON-STANDARD REF. SYSTEM IS USED

LPOT, TRUE IF SPHERICAL HARMONIC EXPANSION IS USED

LTEST, TRUE IF TEST-OUTPUT IS NEEDED

LLEG, TRUE IF LEGEND IS TO BE OUTPUT

LPARAM,TRUE IF PARAMETERS ARE TO BE DETERMINED

LNCOL, TRUE IF COLLOCATION IS NOT USED

LIOSOL,TRUE IF SOLUTION IS STORED OR RECOVERED

f t f f f t f

ARE ALL PARAMETERS OK ?

t

INPUT CODE FOR BASIC REFERENCE SYSTEM:

0: USER DEFINED, 1: ED50 NORTH SEA, 2: ED50/EDOC,

3: NAD1927 /NEW MEXICO, 4: GRS67, 5: GRS80, 6: NWL9D,

7: BEST CURRENT, 8: BEST CUR. FAROE ISL, 9: ED50 FOR SF,

10: IAG-75, 11: KRASSOWSKY, DDR, 12: GERMAN DHDN, BESS.

5

REFERENCE SYSTEM:

+ GRS1980.

A = 6378137.00 M

1/F = 298.2572221

GM= 0.3986005000E+15

REF.GRAVITY AT EQUATOR = 978032.6772 MGAL

POTENTIAL AT REF.ELL. = 62636860.8500 M**2/SEC**2

INPUT NAME OF POT.COEFF. SET

EGM96 TO DEGREE 360

SOURCE OF THE POTENTIAL COEFFICIENTS USED:

EGM TO DEGREE 360

INPUT: GM, SEMI-MAJOR AXIS (M), C(2,0), MAX. DEGREE

LFM, TRUE IF COEFF. IN INPUT STREEM AND *1.0D6

LBIN, TRUE IF ON BINARY FORM

LFORM, TRUE IF FORMAT IS INPUT

LINT, TRUE IF STORED AS INTEGERS

3.986005D14 6378136.2 0.0 360 f f t f

GM A COFF(5) MAX.DEGREE

0.39860050E+15 6378136.2 0.0000 360

INPUT FORMAT (2I4,2D18.0) F.EX.

(2i4,2d20.12)

INPUT NAME OF FILE HOLDING COEFF.

/disk1/cct/geod/E4/EGM96

NAME OF FILE HOLDING COEFFICIENTS: EGM96

COEFFICIENTS UP TO N=5

2 0 -0.484165372E-03 0.000000000E+00

2 1 -0.186987636E-09 0.119528012E-08

2 2 0.243914352E-05 -0.140016684E-05

3 0 0.957254174E-06 0.000000000E+00

3 1 0.202998882E-05 0.248513159E-06

3 2 0.904627769E-06 -0.619025944E-06

3 3 0.721072657E-06 0.141435627E-05

4 0 0.539873864E-06 0.000000000E+00

4 1 -0.536321617E-06 -0.473440266E-06

4 2 0.350694106E-06 0.662671573E-06

4 3 0.990771804E-06 -0.200928369E-06

4 4 -0.188560803E-06 0.308853169E-06

INPUT: LGRID - TRUE IF COMPUTATIONS IN A GRID

LERR - TRUE IF ERROR ESTIMATES ARE TO BE COMPUTED

OR REPRODUCED IN OUTPUT

LCOMP- TRUE IF COMPUTED VALUES ARE SUBTRACTED FROM OBSERVED

t f f

INPUT GRID SPECIFICATION

MIN, MAX LATITUDE, MIN, MAX LONGITUDE, STEP IN LAT AND LONG

FUNCTIONAL TYPE (CODE), NEG. VALUE THEN SPH.EXP. SUBTRACTED

COORD.SYSTEM CODE (-1 THEN GLOBAL SYSTEM)

HEIGHT OF GRID POINTS (M)

LMAP - PRIMITIVE MAP OUTPUT

LPUNCH- OUTPUT TO UNIT 17

LMEAN - MEAN VALUES OUTPUT

LMAP&LPUNCH SIMULTANEOUS TRUE, ONLY MAP OUTPUT TO UNIT 17

54.5 57.5 7.0 13.0 0.5 0.5 11 -1 0.0 f t f (11 ændres til 13 for tyngde anomalier!)

GRID CONSIST OF 7 * 13 POINTS

INPUT NAME OF FILE TO HOLD RESULT

ex4h.out (ændres til ex4g.out for tyngde anomalier!)

SIMULTANEOUS OUTPUT TO FILE: ex4h.out

ALL SPECIFICATIONS OK ?

t

SELECTED GEOCENTRIC SYSTEM USED.

0 NO LATITUDE LONGITUDE H ZETA (M)

DEGREES DEGREES M

POT

1 57.500000 7.000000 0.0 41.69

2 57.500000 7.500000 0.0 41.34

3 57.500000 8.000000 0.0 40.87

4 57.500000 8.500000 0.0 40.26

5 57.500000 9.000000 0.0 39.64

.

.

89 54.500000 12.000000 0.0 38.39

90 54.500000 12.500000 0.0 37.84

91 54.500000 13.000000 0.0 37.10

STOP ?

t

GEOCOL TERMINATED AT:

Wed Oct 10 12:02:01 2001

gaia:/disk1/abbas/E04$