1!     CalculiX - A 3-dimensional finite element program
2!              Copyright (C) 1998-2021 Guido Dhondt
3!
4!     This program is free software; you can redistribute it and/or
5!     modify it under the terms of the GNU General Public License as
6!     published by the Free Software Foundation(version 2);
7!
8!
9!     This program is distributed in the hope that it will be useful,
10!     but WITHOUT ANY WARRANTY; without even the implied warranty of
11!     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12!     GNU General Public License for more details.
13!
14!     You should have received a copy of the GNU General Public License
15!     along with this program; if not, write to the Free Software
16!     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
17!
18      subroutine materialdata_dvi(shcon,nshcon,imat,dvi,t1l,ntmat_,
19     &  ithermal)
20!
21      implicit none
22!
23!     determines the dynamic viscosity
24!
25      integer imat,ntmat_,id,nshcon(*),four,ithermal
26!
27      real*8 t1l,shcon(0:3,ntmat_,*),dvi
28!
29      four=4
30!
31      if(ithermal.eq.0) then
32         dvi=shcon(2,1,imat)
33      else
34         call ident2(shcon(0,1,imat),t1l,nshcon(imat),four,id)
35         if(nshcon(imat).eq.0) then
36            continue
37         elseif(nshcon(imat).eq.1) then
38            dvi=shcon(2,1,imat)
39         elseif(id.eq.0) then
40            dvi=shcon(2,1,imat)
41         elseif(id.eq.nshcon(imat)) then
42            dvi=shcon(2,id,imat)
43         else
44            dvi=shcon(2,id,imat)+
45     &           (shcon(2,id+1,imat)-shcon(2,id,imat))*
46     &           (t1l-shcon(0,id,imat))/
47     &           (shcon(0,id+1,imat)-shcon(0,id,imat))
48         endif
49      endif
50!
51      return
52      end
53
54
55
56
57
58
59
60