DataMuseum.dk

top - metrics - download

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

└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ.  Detaljer om "HC8000" projekt.
    └─⟦0364f57e3⟧ 
        └─⟦7e928b248⟧ »algbib« 
            └─⟦this⟧ 

TextFile

;gosav
r=edit algcoriolis
l./set/,r/41/44/
l./The output/,i/
   7)  Two textstrings in < >, giving the names of two disc areas.
       In the first a input for ROTFIT is prepared, while a LIN-
       DA input is biuld up in the second (partial derivatives of
       d(ny) with respect to dW(vib), ga and Gb are given).
/
l.?close(Lu,true)?,l1
i/
readhead(in,head,1); i:=1;
open(Lu,4,string head(increase(i)),0);
readhead(in,head,1); i:=1;
open(Ll,4,string head(increase(i)),0);
/
l./stop:/
i?
for ku:= 1,2 do begin
for i:=1 step 1 until N do begin
   j:= qtn(i,1); k:= qtn(i,2);
   a:= trans(i,1); b:= trans(i,2); c:= trans(i,3); d:= trans(i,4);
   fu:= trans(i,13)+a-b; fl:= trans(i,14)+c-d;
   if trans(i,12+ku) > 1 then begin
   write(Lu,nl,1,<<ddd>,
   j shift(-16) extract 8,j shift(-8) extract 8,j extract 8,
   sp,2,k shift(-16) extract 8,k shift(-8) extract 8,k extract 8);
   k1:= abs((j shift(-8) extract 8)-(k shift(-8) extract 8));
   k2:= abs((j extract 8)-(k extract 8));
   k1:= if k1//2*2=k1 then 0 else 1;
   k2:= if k2//2*2=k2 then 0 else 1;
   Jm:= if k1=0 and k2=1 then 1 else
        if k1=1 and k2=0 then 3 else 2;
   write(Lu,<<ddd>,Jm,<<  dddddd.ddd>,case ku of (fu,fl));
   write(Ll,nl,1,<<   -d.dddddd>,
         trans(i,10+ku),trans(i,2+3*ku),trans(i,3+3*ku),
         sp,2,<:w:>,sp,2,<:1:>);
end end;
write(Lu,nl,1,<:-1:>);
end;
write(Lu,nl,1,<:<25>:>); write(Ll,nl,1,<:<25>:>);
?
l1,r/res)/res);/
l1,i/
close(Lu,true); close(Ll,true);
/,f
i r
▶EOF◀