| 1 | subroutine vprops
|
|---|
| 2 | C-----------------------------------------------------------------------
|
|---|
| 3 | C
|
|---|
| 4 | C Set material properties
|
|---|
| 5 | C
|
|---|
| 6 | C Material type: 0 for default (PARAM and PCOND/PRHOCP)
|
|---|
| 7 | C 1 for constant props;
|
|---|
| 8 | C 2 for fortran function;
|
|---|
| 9 | C
|
|---|
| 10 | C-----------------------------------------------------------------------
|
|---|
| 11 | include 'SIZE'
|
|---|
| 12 | include 'INPUT'
|
|---|
| 13 | include 'SOLN'
|
|---|
| 14 | include 'TSTEP'
|
|---|
| 15 | LOGICAL IFKFLD,IFEFLD
|
|---|
| 16 | C
|
|---|
| 17 | if (nio.eq.0.and.loglevel.gt.2)
|
|---|
| 18 | $ write(6,*) 'vprops', ifield
|
|---|
| 19 |
|
|---|
| 20 | NXYZ1 = lx1*ly1*lz1
|
|---|
| 21 | NEL = NELFLD(IFIELD)
|
|---|
| 22 | NTOT1 = NXYZ1*NEL
|
|---|
| 23 | C
|
|---|
| 24 | IF (ISTEP.EQ.0) THEN
|
|---|
| 25 | C
|
|---|
| 26 | C First time around, set defaults
|
|---|
| 27 | C
|
|---|
| 28 | ifvarp(ifield) = .false.
|
|---|
| 29 | if (iflomach) ifvarp(ifield) = .true.
|
|---|
| 30 |
|
|---|
| 31 | if (.not.ifvarp(ifield)) then ! check all groups
|
|---|
| 32 | do iel=1,nel
|
|---|
| 33 | igrp = igroup(iel)
|
|---|
| 34 | itype = matype(igrp,ifield)
|
|---|
| 35 | if(itype.ne.0) ifvarp(ifield) = .true.
|
|---|
| 36 | enddo
|
|---|
| 37 | endif
|
|---|
| 38 |
|
|---|
| 39 | itest = 0 ! test against all processors
|
|---|
| 40 | if (ifvarp(ifield)) itest = 1
|
|---|
| 41 | itest = iglmax(itest,1)
|
|---|
| 42 | if (itest.gt.0) ifvarp(ifield) = .true.
|
|---|
| 43 |
|
|---|
| 44 | endif
|
|---|
| 45 | C
|
|---|
| 46 | C Fill up property arrays every time step
|
|---|
| 47 | DO 1000 IEL=1,NEL
|
|---|
| 48 | C
|
|---|
| 49 | IGRP=IGROUP(IEL)
|
|---|
| 50 |
|
|---|
| 51 | if (ifuservp) then
|
|---|
| 52 | C
|
|---|
| 53 | C User specified fortran function (pff 2/13/01)
|
|---|
| 54 | CALL NEKUVP (IEL)
|
|---|
| 55 | DIFMIN = VLMIN(VDIFF(1,1,1,IEL,IFIELD),NXYZ1)
|
|---|
| 56 | IF (DIFMIN .LE. 0.0) THEN
|
|---|
| 57 | WRITE (6,100) DIFMIN,IFIELD,IGRP
|
|---|
| 58 | CALL EXITT
|
|---|
| 59 | endif
|
|---|
| 60 | C
|
|---|
| 61 | ELSE IF(MATYPE(IGRP,IFIELD).EQ.1)THEN
|
|---|
| 62 | C
|
|---|
| 63 | C Constant property within groups of elements
|
|---|
| 64 | C
|
|---|
| 65 | CDIFF = CPGRP(IGRP,IFIELD,1)
|
|---|
| 66 | CTRANS = CPGRP(IGRP,IFIELD,2)
|
|---|
| 67 | CALL CFILL(VDIFF (1,1,1,IEL,IFIELD),CDIFF,NXYZ1)
|
|---|
| 68 | CALL CFILL(VTRANS(1,1,1,IEL,IFIELD),CTRANS,NXYZ1)
|
|---|
| 69 | IF (CDIFF.LE.0.0) THEN
|
|---|
| 70 | WRITE(6,100) CDIFF,IFIELD,IGRP
|
|---|
| 71 | 100 FORMAT(2X,'ERROR: Non-positive diffusivity ('
|
|---|
| 72 | $ ,G12.3,') specified for field',I2,', group',I2
|
|---|
| 73 | $ ,' element',I4,'.'
|
|---|
| 74 | $ ,/,'ABORTING in VPROPS',//)
|
|---|
| 75 | CALL EXITT
|
|---|
| 76 | endif
|
|---|
| 77 | C
|
|---|
| 78 | ELSE IF(MATYPE(IGRP,IFIELD).EQ.2)THEN
|
|---|
| 79 | C
|
|---|
| 80 | C User specified fortran function
|
|---|
| 81 | C
|
|---|
| 82 | CALL NEKUVP (IEL)
|
|---|
| 83 | C
|
|---|
| 84 | if(optlevel.le.2) then
|
|---|
| 85 | DIFMIN = VLMIN(VDIFF(1,1,1,IEL,IFIELD),NXYZ1)
|
|---|
| 86 | IF (DIFMIN .LE. 0.0) THEN
|
|---|
| 87 | WRITE (6,100) DIFMIN,IFIELD,IGRP
|
|---|
| 88 | CALL EXITT
|
|---|
| 89 | endif
|
|---|
| 90 | endif
|
|---|
| 91 | C
|
|---|
| 92 | ELSE IF(MATYPE(IGRP,IFIELD).EQ.0)THEN
|
|---|
| 93 | C
|
|---|
| 94 | C Default constant property
|
|---|
| 95 | C
|
|---|
| 96 | CDIFF = CPFLD(IFIELD,1)
|
|---|
| 97 | CTRANS = CPFLD(IFIELD,2)
|
|---|
| 98 | c write(6,*) 'vdiff:',ifield,cdiff,ctrans
|
|---|
| 99 | CALL CFILL(VDIFF (1,1,1,IEL,IFIELD),CDIFF,NXYZ1)
|
|---|
| 100 | CALL CFILL(VTRANS(1,1,1,IEL,IFIELD),CTRANS,NXYZ1)
|
|---|
| 101 | IF (CDIFF.LE.0.0) THEN
|
|---|
| 102 | WRITE(6,200) CDIFF,IFIELD
|
|---|
| 103 | 200 FORMAT(2X,'ERROR: Non-positive diffusivity ('
|
|---|
| 104 | $ ,G12.3,') specified for field',I2,'.',/
|
|---|
| 105 | $ ,'ABORTING in VPROPS',//)
|
|---|
| 106 | CALL EXITT
|
|---|
| 107 | endif
|
|---|
| 108 | endif
|
|---|
| 109 | C
|
|---|
| 110 | 1000 CONTINUE
|
|---|
| 111 | C
|
|---|
| 112 | return
|
|---|
| 113 | end
|
|---|
| 114 | C
|
|---|
