xref: /phasta/phSolver/common/gensav.f (revision 513954ef803c300cddba2bb96b4a5dac0b93489a)
159599516SKenneth E. Jansen        subroutine gensav (ientmp, mattmp, ien,    mater)
259599516SKenneth E. Jansenc
359599516SKenneth E. Jansenc----------------------------------------------------------------------
459599516SKenneth E. Jansenc
559599516SKenneth E. Jansenc  This routine saves the element block data.
659599516SKenneth E. Jansenc
759599516SKenneth E. Jansenc input:
859599516SKenneth E. Jansenc  ientmp (npro,nshl)   : nodal connectivity
959599516SKenneth E. Jansenc  mattmp (npro)        : material type flag
1059599516SKenneth E. Jansenc
1159599516SKenneth E. Jansenc output:
1259599516SKenneth E. Jansenc  ien    (npro,nshl)   : nodal connectivity
1359599516SKenneth E. Jansenc  mater  (npro)        : material type flag
1459599516SKenneth E. Jansenc
1559599516SKenneth E. Jansenc
1659599516SKenneth E. Jansenc Zdenek Johan, Winter 1992.
1759599516SKenneth E. Jansenc----------------------------------------------------------------------
1859599516SKenneth E. Jansenc
19*513954efSKenneth E. Jansen        use readarrays
20*513954efSKenneth E. Jansen        use fncorpmod
2159599516SKenneth E. Jansen        include "common.h"
2259599516SKenneth E. Jansenc
2359599516SKenneth E. Jansen        dimension   ientmp(npro,nshl),
2459599516SKenneth E. Jansen     &              mattmp(npro),           ien(npro,nshl),
2559599516SKenneth E. Jansen     &              mater(npro)
26*513954efSKenneth E. Jansen       integer*8    ienG(npro,nshl)
2759599516SKenneth E. Jansenc
2859599516SKenneth E. Jansenc.... save the element data
2959599516SKenneth E. Jansenc
3059599516SKenneth E. Jansen        do i = 1, nshl
3159599516SKenneth E. Jansen          ien(:,i) = ientmp(:,i)
3259599516SKenneth E. Jansen        enddo
33*513954efSKenneth E. Jansen        if(usingpetsc.eq.1) then
34*513954efSKenneth E. Jansen          do i = 1, nshl
35*513954efSKenneth E. Jansen            if(numpe .ne. 1) then
36*513954efSKenneth E. Jansen              ienG(:,i) = fncorp(abs(ientmp(:,i)))
37*513954efSKenneth E. Jansen            else
38*513954efSKenneth E. Jansen              ienG(:,i) = abs(ientmp(:,i))
39*513954efSKenneth E. Jansen            endif
40*513954efSKenneth E. Jansen          enddo
41*513954efSKenneth E. Jansen        endif
4259599516SKenneth E. Jansenc
4359599516SKenneth E. Jansen        mater = mattmp
4459599516SKenneth E. Jansenc
4559599516SKenneth E. Jansenc.... end
4659599516SKenneth E. Jansenc
4759599516SKenneth E. Jansen        return
4859599516SKenneth E. Jansen        end
49