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