DataMuseum.dk

Presents historical artifacts from the history of:

RC4000/8000/9000

This is an automatic "excavation" of a thematic subset of
artifacts from Datamuseum.dk's BitArchive.

See our Wiki for more about RC4000/8000/9000

Excavated with: AutoArchaeologist - Free & Open Source Software.


top - metrics - download

⟦9e52c10ec⟧ TextFile

    Length: 2304 (0x900)
    Types: TextFile
    Names: »setgaussctx«

Derivation

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

TextFile



;       proj proc             * page 1   27 11 77, 20.34;  

;  set_gauss_crv
;  *************

if listing.yes
char 10 12 10

set_gauss_crv = set 1

set_gauss_crv = algol

external real procedure set_gauss_crv
_____________________________________
_                (reg_lab, k_a);  
array             reg_lab, k_a;  

begin

  array           argm(-1:8);  

  integer         t;  

  real            mer_arc;  

  long            a0_m;  

  t := reg_lab(2) extract 12;  

if t = 2 <*geo*> then
begin

set_gauss_crv := 
k_a(-1)       := (extend 1) shift 46;
k_a(-2)       := 8;
 
chebychef_c(3, t, k_a(-2), argm(t),
_  0.5/(e_rad(long reg_lab(6), reg_lab(7), 3,
_  r_t_geo(pi*(0.25 + 0.25*argm(t)),
_  false add (1 shift 6))))**2, k_a(t));
 
end
 
else






  if t = 3 or t = 4 or t = 12 then
  <*  itm      utm      sb  *>
  begin

    comment scale of half major axis;  
    a0_m := (long reg_lab(6))*reg_lab(8);  
    a0_m := long reg_lab(6) - a0_m;  

    set_gauss_crv :=
    k_a(-1)       :=
    mer_arc       := m_arc(a0_m, reg_lab(7), 1, 
    _                      (extend 1) shift 46 <*90 dg*>);  
    k_a(-2)       := 6;

    chebychef_c(3, t, k_a(-2), argm(t), 
    _  0.5/(e_rad(a0_m, reg_lab(7), 3, 
    _  m_arc(a0_m, reg_lab(7), 2, 
    _  r_t_geo(0.5'-6*(argm(t)+1)*mer_arc, false add 512))))**2, 
    _  k_a(t));  

  end

  else

system(9)alarm:(t, <:<10>ill. setgc:>);

end set_gauss_crv;  

end

if warning.yes
(mode 0.yes
message set_gauss_crv not ok
lookup set_gauss_crv)
▶EOF◀