sortri
C SORTRI SOURCE PV 17/12/05 21:17:22 9646 IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) C --------------------------------------------------------------------- C C CAS DES OBJETS RIGIDITES C ET DES SUPERELEMNETS DONT ON DEMANDE LE SAUVETAGE C LE POINTEUR EST MIS NEGATIF (PILE 3) C C PROGRAMME PAR FARVACQUE - REPRIS PAR LENA C APPELE PAR: SAUV C APPELLE: C======================================================================= C TABLEAU KCOLA : C 1 MELEME 2 CHPOIN 3 MRIGID 4 MCHAFF 5 MCHELM 6 C 7 8 MSOLUT 9 MSTRUC 10 11 MAFFEC 12 MSOSTU C 13 IMATRI 14 MJONCT 15 MATTAC 16 MMATRI 17 MDEFOR 18 MLREEL C 19 MLENTI 20 MCHARG 21 MODELE 22 MEVOLL 23 MSUPER C======================================================================= C -INC SMRIGID -INC SMSUPER -INC PPARAM -INC CCOPTIO -INC TMCOLAC C C C **** CAS DES OBJETS RIGIDITES: ON NE SAUVE QUE LES MMATRI DES OBJETS C **** SPECIFIES PAR L'UTILISATEUR. POUR LES RECONNAITRE ON MET LEUR C **** POINTEUR NEGATIF SEGACT ICOLAC ITLACC=KCOLA(3) IF (ITLACC.LE.0) GO TO 2 N=ITLAC(/1) IF(N.EQ.0) GO TO 1 ideb = kcolac(3)+1 DO 6 IEL=ideb,N MRIGID=ITLAC(IEL) if(mrigid.eq.0) go to 6 SEGACT MRIGID*MOD ICHOLE=-ABS(ICHOLE) SEGDES MRIGID 6 CONTINUE 1 CONTINUE 2 CONTINUE C ------MEME TRAVAIL POUR LES SUPER ELEMENTS-------------- ITLACC=KCOLA(23) IF (ITLACC.LE.0) GO TO 20 N=ITLAC(/1) IF(N.EQ.0) GO TO 10 DO 11 IEL=1,N MSUPER=ITLAC(IEL) if(msuper.eq.0) go to 11 SEGACT MSUPER MRIGID=MRIGTO SEGACT MRIGID*MOD ICHOLE=-ABS(ICHOLE) SEGDES MRIGID MRIGID=MSURAI SEGACT MRIGID*MOD ICHOLE=-ABS(ICHOLE) SEGDES MRIGID MRIGID=MSUMAS IF(MRIGID.NE.0) THEN SEGACT MRIGID*MOD ICHOLE=-ABS(ICHOLE) SEGDES MRIGID ENDIF SEGDES MSUPER 11 CONTINUE 10 CONTINUE 20 CONTINUE SEGDES ICOLAC RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales