36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
60
61
62
63#include "implicit_f.inc"
64
65
66
67#include "units_c.inc"
68
69
70
71
72 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
73 INTEGER,INTENT(IN)::IG
74 INTEGER,INTENT(IN)::IGTYP
75
76 TYPE(SUBMODEL_DATA),INTENT(IN)::LSUBMODEL(*)
77
78 INTEGER,INTENT(INOUT)::IGEO(*)
80 . INTENT(INOUT)::geo(*)
81
82
83
84 INTEGER I, J, IP, K, , IHBE, ISMSTR, ISORTH, ISHEAR
86 . pun
87 LOGICAL IS_AVAILABLE, IS_ENCRYPTED
88
89
90
91 DATA pun/0.1/
92
93 is_encrypted = .false.
94 is_available = .false.
95 ihbe=0
96 ismstr=0
97 isorth= 0
98
99 geo(3)=ismstr
100 IF(ismstr==3.AND.igtyp/=6.AND.igtyp/=14)THEN
101 geo(5)=ep06
102 ENDIF
103
104 igeo( 1)=ig
105 igeo(10)=ihbe
106 igeo(11)=igtyp
107 geo(12) =igtyp+pun
108 geo(171)=ihbe
109
110
111
112
114
115
116
117 CALL hm_get_floatv(
'THICK',geo(1),is_available,lsubmodel,unitab)
118
119 WRITE(iout,1010)ig
120 IF(.NOT. is_encrypted)THEN
121 WRITE(iout,1011)ig,geo(1)
122 ELSE
123 WRITE(iout,1099)ig
124 ENDIF
125
126 ishear = geo(37)
127 IF(ishear==0)THEN
128 geo(37)=0
129 ELSEIF(ishear==1)THEN
130 geo(37)=1
131 ELSEIF(ishear==2)THEN
132 geo(37)=0
133 ENDIF
134
135 igeo(1) =ig
136 igeo(11)=igtyp
137 igeo(17)=isorth
138 IF(geo(39)/=zero.AND.igeo( 9)== 0) igeo( 9)=nint(geo(39))
139 IF(geo(171)/=zero.AND.igeo(10)== 0) igeo(10)=nint(geo(171))
140
141 1010 FORMAT(
142 & 5x,'DUMMY PROPERTY SET'/,
143 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10//)
144 1011 FORMAT(
145 & 5x,'DUMMY PROPERTY SET'/,
146 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
147 & 5x,'THICKNESS . . . . . . . . . . . . . . .=',1pg20.13//)
148 1099 FORMAT(
149 & 5x,'DUMMY PROPERTY SET'/,
150 & 5x,'PROPERTY SET NUMBER . . . . . . . . . .=',i10/,
151 & 5x,'CONFIDENTIAL DATA'//)
152
153 RETURN
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
subroutine hm_option_is_encrypted(is_encrypted)