climel
C CLIMEL SOURCE PV 17/12/05 21:15:19 9646 C======================================================================= C COLlaborateur LIer MELeme C Restaure les pointeurs du maillage pMelem C======================================================================= integer nbsous, nbref integer iRef,iSous,nbnn,nbelem integer ipoPi,iPoint integer iNoeud,jNoeud -INC PPARAM -INC CCOPTIO -INC SMELEME -INC TMCOLAC pointeur pMelem.MELEME pointeur jcolac.ICOLAC pointeur pile.ITLACC C write(ioimp,*) 'Entre dans CLIMEL' if (pMelem.ne.0) then pile=jcolac.kcola(1) segact pMelem*mod nbsous=pMelem.lisous(/1) C write(ioimp,*) 'nbsous: ',nbsous nbref =pMelem.lisref(/1) C write(ioimp,*) 'nbref: ',nbref nbnn =pMelem.num(/1) nbelem=pMelem.num(/2) do iSous=1,nbsous C write(ioimp,*) 'Sous maillage :',iSous iPoPi=pMelem.lisous(iSous) C write(ioimp,*) 'Position dans la pile: ',iPoPi iPoint=pile.itlac(iPoPi) C write(ioimp,*) 'Pointeur: ',iPoint pMelem.lisous(iSous)=iPoint enddo do iRef=1,nbref C write(ioimp,*) 'Reference :',iRef iPoPi=pMelem.lisref(iRef) C write(ioimp,*) 'Position dans la pile: ',iPoPi iPoint=pile.itlac(iPoPi) C write(ioimp,*) 'Pointeur: ',iPoint C & pMelem.lisref(iRef) pMelem.lisref(iRef)=iPoint enddo c write(ioimp,*) 'Traduction' pile=jcolac.kcola(32) do jNoeud=1,nbelem do iNoeud=1,nbnn C write(ioimp,*) 'element',jNoeud,'noeud',iNoeud iPoPi=pMelem.num(iNoeud,jNoeud) pMelem.num(iNoeud,jNoeud)=pile.itlac(iPoPi) enddo enddo segdes pMelem else write(ioimp,*) 'Erreur: pointeur vers un objet MELEME nul' endif C write(ioimp,*) 'Sortie de CLIMEL' end
© Cast3M 2003 - Tous droits réservés.
Mentions légales