usermat/init_dyn21.f at master · declmal/usermat · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
      subroutine usercomm
c
c******************************************************************
c|  Livermore Software Technology Corporation  (LSTC)             |
c|  ------------------------------------------------------------  |
c|  Copyright 1987-2008 Livermore Software Tech. Corp             |
c|  All rights reserved                                           |
c******************************************************************
c
c     to define threadprivate common (aux14loc) for the history variable
c     and any new threadprivate common blocks by the user
c
c     I/O unit 241~245 are reversed for user material usage
c
c     the maximum number of history variables (NHISVAR) is defined in
c     the include file (nhisparm.inc)
c     the user can increase the number of history variables by
c     editing NHISVAR in the include (nhisparm.inc)
c
c Note:
c  For SMP or MPP version, all threadprivate common blocks must be loaded
c  and initialized in the beginning of job
c
c  If the user adds any new threadprivate common blocks, the following
c  three lines must be added:
c   1. common /xxxloc/xxxzz(nlq) # define common block
c   2. c$omp threadprivate (/xxxloc/)  # declare it as threadprivate common
c   3. xxxzz(1)=0.0              # initailze the first variable
c
c For example, to add a new threadprivate common /userloc/
c
c 123456789012345678901234567890123456789012345678901234567890
c      common/user1loc/user1zz1(nlq),user1zz2(nlq)
c c$omp threadprivate (/user1loc/)
c      user1zz1(1)=0.0
c
c---------------------------------------------------------------
cdir$ OPTIMIZE:0
c     this routine must be compiler with -O0 for Intel compiler
c     do not remove the above compiler option line
c
      include 'nlqparm'
      include 'nhisparm.inc'
c
      parameter (NBEL8A=max(100,2*(NHISVUE+45)))
      parameter (NBEL7 =max(2601,
     . 1180+576*NXDOFUE+64*NXDOFUE*NXDOFUE+12*5))
      common/aux14loc/ax14zz1(nlq,7+NHISVAR)
      common/bel8loc/bel8lz1(nlq,2*(NHISVUE+45))
      common/bel8aloc/bel8lz2(nlq,NBEL8A)
      common/bel7loc/bel7lz(nlq,NBEL7)
c
c$omp threadprivate (/aux14loc/)
c$omp threadprivate (/bel8loc/)
c$omp threadprivate (/bel8aloc/)
c$omp threadprivate (/bel7loc/)
      ax14zz1(1,1)=0.0
      bel8lz1(1,1)=0.0
      bel8lz2(1,1)=0.0
      bel7lz (1,1)=0.0
c
      return
      end
      subroutine usercomm1
c
c******************************************************************
c|  Livermore Software Technology Corporation  (LSTC)             |
c|  ------------------------------------------------------------  |
c|  Copyright 1987-2008 Livermore Software Tech. Corp             |
c|  All rights reserved                                           |
c******************************************************************
c
      include 'nlqparm'
      include 'nhisparm.inc'
      parameter (NBEL8A=max(100,2*(NHISVUE+45)))
      parameter (NBEL7 =max(2601,
     . 1180+576*NXDOFUE+64*NXDOFUE*NXDOFUE+12*5))
      common/usercom/iusercomm(4)
      parameter (nconst_max=max(1474,(NHISVAR+26)*8+41+1*24))
      common/strhx/ideh,mxh,ixhx(8),xe(3,8),ve(3,8),strnh(6),
     . ressh,hisvar(nconst_max)
      ideh=0
c
c     define the lenght of 4 common blocks
      iusercomm(1)=nlq*(7+NHISVAR)
      iusercomm(2)=nlq*2*(NHISVUE+45)
      iusercomm(3)=nlq*NBEL8A
      iusercomm(4)=nlq*NBEL7
c
      return
      end