OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_inigrav.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| hm_read_inigrav ../starter/source/initial_conditions/inigrav/hm_read_inigrav.F
25!||--- called by ------------------------------------------------------
26!|| lectur ../starter/source/starter/lectur.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
30!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
31!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
32!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
33!||--- uses -----------------------------------------------------
34!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
35!|| inigrav ../starter/share/modules1/inigrav_mod.F
36!|| message_mod ../starter/share/message_module/message_mod.F
37!|| submodel_mod ../starter/share/modules1/submodel_mod.f
38!||====================================================================
39 SUBROUTINE hm_read_inigrav(IGRV ,IBUF ,AGRV ,ITAB ,ITABM1 ,
40 . IGRPART ,NPC ,UNITAB ,ISKN ,
41 . ITAGND ,IGRSURF ,PLD ,BUFSF ,LSUBMODEL)
42C-----------------------------------------------
43C M o d u l e s
44C-----------------------------------------------
45 USE unitab_mod
46 USE message_mod
47 USE inigrav
48 USE groupdef_mod
49 USE submodel_mod
50 USE hm_option_read_mod
51 USE names_and_titles_mod , ONLY : nchartitle, ncharfield
52C-----------------------------------------------
53C I m p l i c i t T y p e s
54C-----------------------------------------------
55#include "implicit_f.inc"
56C-----------------------------------------------
57C C o m m o n B l o c k s
58C-----------------------------------------------
59#include "com04_c.inc"
60#include "units_c.inc"
61#include "param_c.inc"
62C-----------------------------------------------
63C D u m m y A r g u m e n t s
64C-----------------------------------------------
65 TYPE (unit_type_),INTENT(IN) ::unitab
66 INTEGER IGRV(NIGRV,*), IBUF(*), ITAB(*), ITABM1(*),NPC(*),
67 . iskn(liskn,*),itagnd(*)
68 my_real
69 . agrv(lfacgrv,*)
70 my_real pld(*), bufsf(*)
71 TYPE(submodel_data) LSUBMODEL(*)
72C-----------------------------------------------
73 TYPE (GROUP_) , DIMENSION(NGRPART) :: IGRPART
74 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
75C-----------------------------------------------
76C L o c a l V a r i a b l e s
77C-----------------------------------------------
78 my_real
79 . fcx,fcy,fac_m,fac_l,fac_t,ngx,ngy,ngz,dotprod
80 INTEGER I, NOD, NCUR, NOSKEW,NSKW,NN,IGS,UID,
81 . IAD,NS,IWA,J,K,ID,K1,K2,NCURS,N1,N2,NC,L,
82 . itag, iflagunit,flag_fmt,flag_fmt_tmp,ifix_tmp,iadpl
83 CHARACTER(LEN=NCHARFIELD) :: XYZ
84 CHARACTER X*1, Y*1, Z*1, XX*2, YY*2, ZZ*2, MESS*40
85 CHARACTER(LEN=NCHARTITLE)::TITR
86
87 INTEGER :: IGU,ISU,IGRAV,IBID, IG, IS, IDIR, PN1, ICURS, IIGRAV,IIG,IIS
88 my_real :: BX_,BY_,BZ_, GRAV0,NX,NY,NZ,NORM,PSURF,BID
89 LOGICAL :: lFOUND, lPLANAR_SURF, lUSER_SURF, lOUTP, lGRAV, lUNIQUE, IS_AVAILABLE
90 CHARACTER*2 :: CDIR
91
92 INTEGER :: M,ID_LIST(NINIGRAV)
93C-----------------------------------------------
94 DATA x/'X'/
95 DATA y/'Y'/
96 DATA z/'Z'/
97 DATA xx/'XX'/
98 DATA yy/'YY'/
99 DATA zz/'ZZ'/
100 DATA mess/'INITIAL GRAVITY LOADING DEFINITION '/
101C=======================================================================
102C
103 ! Initialization of variable
104 lplanar_surf = .false.
105 luser_surf = .false.
106 lgrav = .true.
107 loutp = .true.
108 is_available = .false.
109 ngx = zero
110 ngy = zero
111 ngz = zero
112C
113 ! Start reading /INIGRAV card
114 CALL hm_option_start('/INIGRAV')
115C
116 ! Loop over /INIGRAV
117 DO k=1,ninigrav
118C
119 ! Read title, ID and Unit ID
120 titr = ''
121 CALL hm_option_read_key(lsubmodel,
122 . option_id = id,
123 . unit_id = uid,
124 . option_titr = titr)
125C
126 ! Checking unit ID
127 iflagunit = 0
128 DO j=1,unitab%NUNITS
129 IF (unitab%UNIT_ID(j) == uid) THEN
130 iflagunit = 1
131 EXIT
132 ENDIF
133 ENDDO
134C
135 id_list(k)=id
136 lunique = .true.
137 DO m=1,k-1
138 IF(id==id_list(m))THEN
139 lunique=.false.
140 EXIT
141 ENDIF
142 ENDDO
143C
144 IF (uid /= 0 .AND. iflagunit == 0) THEN
145 CALL ancmsg(msgid=659,anmode=aninfo,msgtype=msgerror,
146 . i2=uid,i1=id,
147 . c1='INITIAL GRAVITY LOADING',
148 . c2='INITIAL GRAVITY LOADING',
149 . c3= titr)
150 loutp = .false.
151 ENDIF
152C
153 ! Reading 1st card : ids
154 CALL hm_get_intv('grpart_ID',igu,is_available,lsubmodel)
155 CALL hm_get_intv('surf_ID' ,isu,is_available,lsubmodel)
156 CALL hm_get_intv('grav_id' ,IGRAV,IS_AVAILABLE,LSUBMODEL)
157C
158 ! Reading 2nd card : pressure, etc
159 CALL HM_GET_FLOATV('pref' ,PSURF,IS_AVAILABLE, LSUBMODEL, UNITAB)
160 ! Only if ISU == 0
161 IF (ISU == 0) THEN
162 CALL HM_GET_FLOATV('bx',BX_,IS_AVAILABLE, LSUBMODEL, UNITAB)
163 CALL HM_GET_FLOATV('by',BY_,IS_AVAILABLE, LSUBMODEL, UNITAB)
164 CALL HM_GET_FLOATV('bz',BZ_,IS_AVAILABLE, LSUBMODEL, UNITAB)
165 ENDIF
166C
167 ! Checking Gravity ID
168 lFOUND = .FALSE.
169 IIGRAV = 0
170 GRAV0 = ZERO
171 DO IG=1,NGRAV
172 IF (IGRAV == IGRV(5,IG)) THEN
173 lFOUND = .TRUE.
174 IIGRAV = IG
175 ICURS = IGRV(3,IG)
176 IF (ICURS > 0) THEN
177 PN1 = NPC(ICURS)
178 GRAV0 = AGRV(1,IG)*PLD(PN1+1)
179 ELSE
180 GRAV0 = AGRV(1,IG)
181 ENDIF
182 IDIR = MOD(IGRV(2,IG),10)
183 NGX = ZERO
184 NGY = ZERO
185 NGZ = ZERO
186 SELECT CASE (IDIR)
187 CASE(1)
188 CDIR(1:2) =' x'
189 NGX = ONE
190 CASE(2)
191 CDIR(1:2) =' y'
192 NGY = ONE
193 CASE(3)
194 CDIR(1:2) =' z'
195 NGZ = ONE
196 END SELECT
197 CDIR(1:1)="+"
198 IF (GRAV0 < ZERO) THEN
199 CDIR(1:1)="-"
200 NGX = -NGX
201 NGY = -NGY
202 NGZ = -NGZ
203 ENDIF
204 EXIT
205 ENDIF
206 ENDDO
207 ! Wrong gravity ID
208.NOT. IF (lFOUND) THEN
209 CALL ANCMSG(MSGID=73,ANMODE=ANINFO,MSGTYPE=MSGERROR,
210 . I1=ID,
211 . C1=TITR,
212 . I2= IGRAV,
213 . C2='does not refer to a valid /grav id')
214 lOUTP = .FALSE.
215 lGRAV = .FALSE.
216 ENDIF
217 ! Inigrav ID duplicated
218.NOT. IF (lUNIQUE) THEN
219 CALL ANCMSG(MSGID=73,ANMODE=ANINFO,MSGTYPE=MSGERROR,
220 . I1=ID,
221 . C1=TITR,
222 . I2= ID,
223 . C2='identifier is duplicated')
224 ENDIF
225C
226 ! Checking surface ID
227 IIS = 0
228 IF (ISU > 0) THEN
229 lfound = .FALSE.
230 DO IS=1,NSURF
231 IF (ISU == IGRSURF(IS)%ID)THEN
232 SELECT CASE(IGRSURF(IS)%TYPE)
233 CASE(0)
234 IIS = IS
235 lUSER_SURF=.TRUE.
236 IADPL = IGRSURF(IS)%IAD_BUFR
237 lfound = .TRUE.
238 BX_ = ZERO
239 BY_ = ZERO
240 BZ_ = ZERO
241 NX = ZERO
242 NY = ZERO
243 NZ = ZERO
244 NX = ZERO
245 NY = ZERO
246 NZ = ZERO
247 EXIT
248 CASE(200)
249 IIS = IS
250 lPLANAR_SURF=.TRUE.
251 IADPL = IGRSURF(IS)%IAD_BUFR
252 lfound = .TRUE.
253 BX_ = BUFSF(IADPL+1)
254 BY_ = BUFSF(IADPL+2)
255 BZ_ = BUFSF(IADPL+3)
256 NX = BUFSF(IADPL+4)- BUFSF(IADPL+1)
257 NY = BUFSF(IADPL+5)- BUFSF(IADPL+2)
258 NZ = BUFSF(IADPL+6)- BUFSF(IADPL+3)
259 NORM = SQRT(NX*NX+NY*NY+NZ*NZ)
260 NX = NX / NORM
261 NY = NY / NORM
262 NZ = NZ / NORM
263 EXIT
264 END SELECT
265 ENDIF
266 ENDDO
267 ! Wrong surface ID
268.NOT. IF (lFOUND) THEN
269 CALL ANCMSG(MSGID=73,ANMODE=ANINFO,MSGTYPE=MSGERROR,
270 . I1=ID,
271 . C1=TITR,
272 . I2= ISU,
273 . C2='does not refer to a valid /surf id')
274 lOUTP = .FALSE.
275 ENDIF
276 ELSE
277 NX = NGX
278 NY = NGY
279 NZ = NGZ
280 ENDIF
281C
282 ! Checking GRPART ID
283 lfound = .FALSE.
284 IIG = 0
285 IAD = NGRNOD+NGRBRIC+NGRQUAD+NGRSHEL+NGRSH3N+NGRTRUS+NGRBEAM+NGRSPRI
286 IF (IGU > 0) THEN
287 DO IG=1,NGRPART
288 IF (IGU == IGRPART(IG)%ID) THEN
289 IIG = IG
290 lfound = .TRUE.
291 EXIT
292 ENDIF
293 ENDDO
294 ! Wrong GRPART ID
295.NOT. IF (lFOUND) THEN
296 CALL ANCMSG(MSGID=73,ANMODE=ANINFO,MSGTYPE=MSGERROR,
297 . I1=ID,
298 . C1=TITR,
299 . I2= IGU,
300 . C2='does not refer to a valid grpart ')
301 lOUTP = .FALSE.
302 ENDIF
303 ENDIF
304C
305 ! Checking the normal
306.AND. IF (lPLANAR_SURF lGRAV) THEN
307 DOTPROD = NX*NGX+NY*NGY+NZ*NGZ
308 IF(ABS(DOTPROD)<=EM20)THEN
309 CALL ANCMSG(MSGID=73,ANMODE=ANINFO,MSGTYPE=MSGERROR,
310 . I1=ID,
311 . C1=TITR,
312 . I2= IGRAV,
313 . C2='refer to a gravity direction colinear to the input surface')
314 lOUTP = .FALSE.
315 ENDIF
316 ENDIF
317C
318 ! Storing IDs in INIGRV table
319 INIGRV(1,K) = IIG
320 INIGRV(2,K) = IIS
321 INIGRV(3,K) = IIGRAV
322 INIGRV(4,K) = ID
323C
324 ! Storing real data in LINIGRAV table
325 LINIGRAV(01,K) = BX_
326 LINIGRAV(02,K) = BY_
327 LINIGRAV(03,K) = BZ_
328 LINIGRAV(04,K) = NX
329 LINIGRAV(05,K) = NY
330 LINIGRAV(06,K) = NZ
331 LINIGRAV(07,K) = GRAV0
332 LINIGRAV(08,K) = NGX
333 LINIGRAV(09,K) = NGY
334 LINIGRAV(10,K) = NGZ
335 LINIGRAV(11,K) = PSURF
336C
337 ! Printout data
338 IF (lPLANAR_SURF) THEN
339 WRITE (IOUT,2000)
340 WRITE (IOUT,FMT='(a)') ''
341 WRITE (IOUT,3000) IGU,ISU,IGRAV,BX_,BY_,BZ_, PSURF
342 WRITE (IOUT,3001) CDIR(2:2)
343 WRITE (IOUT,3002) GRAV0
344 IF(lPLANAR_SURF) WRITE (IOUT,3003) NX,NY,NZ
345 ELSEIF(lUSER_SURF)THEN
346 WRITE (IOUT,2001)
347 WRITE (IOUT,FMT='(a)') ''
348 WRITE (IOUT,3005) IGU,ISU,IGRAV, PSURF
349 WRITE (IOUT,3001) CDIR(2:2)
350 WRITE (IOUT,3002) GRAV0
351 IF(lUSER_SURF)WRITE (IOUT,3004)
352 ENDIF
353 ENDDO !next K
354C-----------
355 RETURN
356C-----------
357 2000 FORMAT(//
358 .' initial gravity loading '/
359 .' ----------------------- '/
360 .' grpart_id surf_id grav_id bx by bz psurf ')
361
362 2001 FORMAT(//
363 .' initial gravity loading '/
364 .' ----------------------- '/
365 .' grpart_id surf_id grav_id psurf ')
366
367 3000 FORMAT(2X,I10,2X,I10,2X,I10,2X,E12.4,2X,E12.4,2X,E12.4,2X,E12.4)
368 3005 FORMAT(2X,I10,2X,I10,2X,I10,3X,E12.4)
369
370 3001 FORMAT(' gravity orientation : ',1X,A2)
371 3002 FORMAT(' gravity VALUE : ',2X,E12.4)
372 3003 FORMAT(' surface orientation : ',2X,E12.4,2X,E12.4,2X,E12.4)
373 3004 FORMAT(' user defined surface')
374C-----------
375 RETURN
376 END
my_real
#define my_real
Definition cppsort.cpp:32
the
end diagonal values have been computed in the(sparse) matrix id.SOL
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:251
hm_read_inigrav
subroutine hm_read_inigrav(igrv, ibuf, agrv, itab, itabm1, igrpart, npc, unitab, iskn, itagnd, igrsurf, pld, bufsf, lsubmodel)
Definition hm_read_inigrav.F:42
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:410
groupdef_mod
Definition groupdef_mod.F:662
hm_option_read_mod
Definition hm_option_read_mod.F:408
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
names_and_titles_mod::ncharfield
integer, parameter ncharfield
Definition names_and_titles_mod.F:1005
submodel_mod
Definition submodel_mod.F:712
unitab_mod
Definition unitab_mod.F:513
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
starter
program starter
Definition starter.F:39
submodel_mod::submodel_data
Definition submodel_mod.F:723
unitab_mod::unit_type_
Definition unitab_mod.F:520