This routine calculates the thrust, oblateness, and the thrust-, major- & minor axes. by calling LUTHRU.
Before calling this routine /LUJETS/ must have been filled, eg. by calling ZAUJSF. (NB ZAUCLU calls ZAUJSF).
none.
From /LUDAT/ : (filled by the call to LUSPHE) P(NLUND+1,...), major axis of sphericity tensor.
To /ZAEVT/ : THRUST, Thrust-,major- & minor- axes; THR, thrust; OBL, olateness, THTHRU, theta of thrust axis.