DataMuseum.dk

top - metrics - download

    Length: 1536 (0x600)
    Types: TextFile
    Names: »gravi67tx«

└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ.  Detaljer om "HC8000" projekt.
    └─⟦0364f57e3⟧ 
        └─⟦80900d603⟧ »giprocfile« 
└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
    └─⟦b2ec5d50f⟧ 
        └─⟦80900d603⟧ »giprocfile« 
            └─⟦this⟧ 

TextFile

; gravi67_tx     * page  ;



;  
;  g r a v i _ 6 7
;  _______________
;  

gravi_67=set 1
gravi_67=algol message.no

external real procedure gravi_67(fi);  
____________________________________

value fi;  
real  fi;  

begin

  comment

  The procedure computes the normal gravity at a 
  point  with the geodetic latitude fi using the

  _        GRAVITY FORMULA 1967

  from     System Geodesique de Rererence 1967, 
  _        Bureau Central de
  _        L'Association Internationale de Geodesie.
  _        Paris, 1971.

  The angle is in units of radians and the computed
  gravity in units of gal (cm/s**2).

  G.I. prog index no. 76001.
  Willy Lehmann Weng, march 1976.
 
  ;  
  real sin_2_fi;  
  sin_2_fi := sin(fi)**2;  
  gravi_67 := 978.031 845 58 *
  _          (1 + 0.001 931 663 383 21 * sin_2_fi) /
  _           sqrt(1 - 0.006 694 605 328 56 * sin_2_fi);  
end gravi_67;  
end




if ok.no
mode warning.yes


if warning.yes
(mode 0.yes
message gravi_67 not ok
lookup gravi_67)


end


finis
▶EOF◀