OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_funct.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_funct ../starter/source/tools/curve/hm_read_funct.F
25!||--- called by ------------------------------------------------------
26!|| lectur ../starter/source/starter/lectur.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| fretitl ../starter/source/starter/freform.F
30!|| hm_get_float_array_index ../starter/source/devtools/hm_reader/hm_get_float_array_index.F
31!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
32!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
33!|| hm_option_count ../starter/source/devtools/hm_reader/hm_option_count.F
34!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
35!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
36!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
37!||--- uses -----------------------------------------------------
38!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
39!|| message_mod ../starter/share/message_module/message_mod.F
40!|| submodel_mod ../starter/share/modules1/submodel_mod.F
41!|| table_mod ../starter/share/modules1/table_mod.F
42!||====================================================================
43 SUBROUTINE hm_read_funct(NPC ,PLD ,NFUNCT ,TABLE ,NPTS_ALLOC,
44 . NOM_OPT ,FUNCRYPT ,UNITAB ,LSUBMODEL)
45C-----------------------------------------------
46C M o d u l e s
47C-----------------------------------------------
48 USE table_mod
49 USE message_mod
50 USE submodel_mod
51 USE hm_option_read_mod
52 USE unitab_mod
53 USE names_and_titles_mod , ONLY : nchartitle
54C-----------------------------------------------
55C I m p l i c i t T y p e s
56C-----------------------------------------------
57#include "implicit_f.inc"
58C-----------------------------------------------
59C C o m m o n B l o c k s
60C-----------------------------------------------
61#include "scr17_c.inc"
62#include "units_c.inc"
63C-----------------------------------------------
64C D u m m y A r g u m e n t s
65C-----------------------------------------------
66 INTEGER NFUNCT, NPTS_ALLOC
67 INTEGER NPC(*),FUNCRYPT(*)
68 my_real pld(npts_alloc)
69 TYPE(ttable) TABLE(*)
70 INTEGER NOM_OPT(LNOPT1,*)
71 TYPE(submodel_data), DIMENSION(NSUBMOD), INTENT(IN) :: LSUBMODEL
72 TYPE(unit_type_), INTENT(IN) :: UNITAB
73C-----------------------------------------------
74C L o c a l V a r i a b l e s
75C-----------------------------------------------
76 INTEGER I,L,FUNC_ID,NPTS,STAT,N,ISMOOTH
77 INTEGER :: IPYTHON ! is it a FUNCT_PYTHON
78 my_real time, funct
79 my_real xscale,yscale,xshift,yshift
80 CHARACTER(LEN=NCHARTITLE) :: TITR
81 CHARACTER MESS*40,KEY*20
82 DATA mess/' FUNCTION & TABLES DEFINITION '/
83 LOGICAL :: IS_ENCRYPTED, IS_AVAILABLE
84 INTEGER :: NB_FUNCT, NB_FUNCT_SMOOTH, IPT, NPT
85C--------------------------------------------------
86C B e g i n n i n g o f S u b r o u t i n e
87C--------------------------------------------------
88 IF (nfunct == 0) RETURN
89
90 stat = 0
91
92 WRITE (iout,2000) nfunct
93
94! Initialization
95 is_encrypted = .false.
96 is_available = .false.
97
98 npc(1)=1
99 l =0
100C--------------------------------------------------
101C READING /FUNCT ( ISMOOTH = 0, IPYTHON = 0)
102C--------------------------------------------------
103 CALL hm_option_count('/FUNCT', nb_funct)
104 CALL hm_option_count('/FUNCT_SMOOTH', nb_funct_smooth)
105 IF (nb_funct > 0) THEN
106 CALL hm_option_start('/FUNCT')
107 DO i = 1, nb_funct
108c
109 CALL hm_option_read_key(lsubmodel,
110 . option_titr = titr,
111 . option_id = func_id,
112 . keyword1 = key)
113c
114 ismooth = 0
115 ipython = 0
116 IF(key(6:12) == '_SMOOTH') ismooth = 1
117 IF(key(6:12) == '_PYTHON') ipython = 1
118
119 IF(ismooth == 0 .AND. ipython == 0 ) THEN
120 CALL hm_option_is_encrypted(is_encrypted)
121 l = l + 1
122 CALL fretitl(titr,nom_opt(lnopt1-ltitr+1,l),ltitr)
123 nom_opt(1, l) = func_id
124 npc(nfunct + l + 1) = func_id
125 npc(2 * nfunct + l + 1) = ismooth
126 npc(l + 1) = npc(l)
127 npts = 0
128 WRITE(iout, 2100) func_id
129! Number of points
130 CALL hm_get_intv('numberofpoints', npt, is_available, lsubmodel)
131c
132 DO ipt = 1, npt
133 CALL hm_get_float_array_index('points', time, 2 * ipt - 1, is_available, lsubmodel, unitab)
134 CALL hm_get_float_array_index('points', funct, 2 * ipt, is_available, lsubmodel, unitab)
135 IF (.NOT. is_encrypted) THEN
136 WRITE(iout,'(3X,1PG20.13,2X,1G20.13)') time,funct
137 ENDIF
138 npts = npts + 1
139 pld(npc(l + 1)) = time
140 IF (npts > 1) THEN
141 IF (pld(npc(l+1)) <= pld(npc(l+1)-2)) THEN
142! Decreasing time line
143 CALL ancmsg(msgid = 156, msgtype = msgerror, anmode = aninfo_blind_1,
144 . i1 = func_id, c1 = titr, i2 = npts, i3 = npts-1)
145 ENDIF
146 ENDIF
147 npc(l + 1) = npc(l + 1) + 1
148 pld(npc(l + 1)) = funct
149 npc(l + 1) = npc(l + 1) + 1
150 ENDDO
151c
152 IF (npt < 2) THEN
153 CALL ancmsg(msgid=1874, msgtype=msgwarning, anmode=aninfo_blind_1,
154 . i1=func_id,
155 . c1=titr)
156 END IF
157c
158! build table structure
159 table(l)%NOTABLE=func_id
160 table(l)%NDIM =1
161!
162 ALLOCATE(table(l)%X(1),stat=stat)
163 IF (stat /= 0) CALL ancmsg(msgid=268,anmode=aninfo,
164 . msgtype=msgerror,
165 . c1='TABLE')
166 ALLOCATE(table(l)%X(1)%VALUES(npts),stat=stat)
167 IF (stat /= 0) CALL ancmsg(msgid=268,anmode=aninfo,
168 . msgtype=msgerror,
169 . c1='TABLE')
170 ALLOCATE(table(l)%Y,stat=stat) ! Y or Y(1) ?
171 IF (stat /= 0) CALL ancmsg(msgid=268,anmode=aninfo,
172 . msgtype=msgerror,
173 . c1='TABLE')
174 ALLOCATE(table(l)%Y%VALUES(npts),stat=stat)
175 IF (stat /= 0) CALL ancmsg(msgid=268,anmode=aninfo,
176 . msgtype=msgerror,
177 . c1='table')
178!
179 DO N=1,NPTS
180 TABLE(L)%X(1)%VALUES(N)=PLD(NPC(L)+2*N-2)
181 TABLE(L)%Y%VALUES(N) =PLD(NPC(L)+2*N-1)
182 ENDDO
183 IF (IS_ENCRYPTED) THEN
184 WRITE(IOUT,'(a)')'confidential data'
185 FUNCRYPT(L) = 1
186 ENDIF
187 ENDIF
188 ENDDO
189 ENDIF
190C--------------------------------------------------
191C READING /FUNCT_SMOOTH ( ISMOOTH = 1)
192C--------------------------------------------------
193 IF (NB_FUNCT_SMOOTH > 0) THEN
194 CALL HM_OPTION_START('/funct_smooth')
195 DO I = 1, NB_FUNCT_SMOOTH
196c
197 CALL HM_OPTION_READ_KEY(LSUBMODEL,
198 . OPTION_TITR = TITR,
199 . OPTION_ID = FUNC_ID,
200 . KEYWORD1 = KEY)
201c
202 CALL HM_OPTION_IS_ENCRYPTED(IS_ENCRYPTED)
203 ISMOOTH = 1
204 L = L + 1
205 CALL FRETITL(TITR,NOM_OPT(LNOPT1-LTITR+1,L),LTITR)
206 NOM_OPT(1, L) = FUNC_ID
207 NPC(NFUNCT + L + 1) = FUNC_ID
208 NPC(2 * NFUNCT + L + 1) = ISMOOTH
209 NPC(L + 1) = NPC(L)
210 NPTS = 0
211 WRITE(IOUT, 2200) FUNC_ID
212c
213 CALL HM_GET_FLOATV('a_scale_x' ,XSCALE ,IS_AVAILABLE ,LSUBMODEL ,UNITAB)
214 CALL HM_GET_FLOATV('f_scale_y' ,YSCALE ,IS_AVAILABLE ,LSUBMODEL ,UNITAB)
215 CALL HM_GET_FLOATV('a_shift_x' ,XSHIFT ,IS_AVAILABLE ,LSUBMODEL ,UNITAB)
216 CALL HM_GET_FLOATV('f_shift_y' ,YSHIFT ,IS_AVAILABLE ,LSUBMODEL ,UNITAB)
217 IF (XSCALE == ZERO) XSCALE = ONE
218 IF (YSCALE == ZERO) YSCALE = ONE
219
220 WRITE (IOUT,2300)
221.NOT. IF ( IS_ENCRYPTED)
222 . WRITE(IOUT,'(3x,1pg20.13,3(2x,1g20.13))') XSCALE,YSCALE,XSHIFT,YSHIFT
223 WRITE (IOUT,2400)
224
225! Number of points
226 CALL HM_GET_INTV('numberofpoints', NPT, IS_AVAILABLE, LSUBMODEL)
227c
228 DO IPT = 1, NPT
229 CALL HM_GET_FLOAT_ARRAY_INDEX('points', TIME, 2 * IPT - 1, IS_AVAILABLE, LSUBMODEL, UNITAB)
230 CALL HM_GET_FLOAT_ARRAY_INDEX('points', FUNCT, 2 * IPT, IS_AVAILABLE, LSUBMODEL, UNITAB)
231c
232 TIME = TIME * XSCALE + XSHIFT
233 FUNCT = FUNCT * YSCALE + YSHIFT
234c
235.NOT. IF ( IS_ENCRYPTED) THEN
236 WRITE(IOUT,'(3x,1pg20.13,2x,1g20.13)') TIME,FUNCT
237 ENDIF
238 NPTS = NPTS + 1
239 PLD(NPC(L + 1)) = TIME
240 IF (NPTS > 1) THEN
241 IF (PLD(NPC(L+1)) <= PLD(NPC(L+1)-2)) THEN
242! Decreasing time line
243 CALL ANCMSG(MSGID = 156, MSGTYPE = MSGERROR, ANMODE = ANINFO_BLIND_1,
244 . I1 = FUNC_ID, C1 = TITR, I2 = NPTS, I3 = NPTS-1)
245 ENDIF
246 ENDIF
247 NPC(L + 1) = NPC(L + 1) + 1
248 PLD(NPC(L + 1)) = FUNCT
249 NPC(L + 1) = NPC(L + 1) + 1
250 ENDDO
251c
252 IF (NPT < 2) THEN
253 CALL ANCMSG(MSGID=1874, MSGTYPE=MSGWARNING, ANMODE=ANINFO_BLIND_1,
254 . I1=FUNC_ID,
255 . C1=TITR)
256 END IF
257c
258! build table structure
259 TABLE(L)%NOTABLE=FUNC_ID
260 TABLE(L)%NDIM =1
261!
262 ALLOCATE(TABLE(L)%X(1),STAT=stat)
263 IF (STAT /= 0) CALL ANCMSG(MSGID=268,ANMODE=ANINFO,
264 . MSGTYPE=MSGERROR,
265 . C1='table')
266 ALLOCATE(TABLE(L)%X(1)%VALUES(NPTS),STAT=stat)
267 IF (STAT /= 0) CALL ANCMSG(MSGID=268,ANMODE=ANINFO,
268 . MSGTYPE=MSGERROR,
269 . C1='table')
270 ALLOCATE(TABLE(L)%Y,STAT=stat)
271 IF (STAT /= 0) CALL ANCMSG(MSGID=268,ANMODE=ANINFO,
272 . MSGTYPE=MSGERROR,
273 . C1='table')
274 ALLOCATE(TABLE(L)%Y%VALUES(NPTS),STAT=stat)
275 IF (STAT /= 0) CALL ANCMSG(MSGID=268,ANMODE=ANINFO,
276 . MSGTYPE=MSGERROR,
277 . C1='table')
278!
279 DO N=1,NPTS
280 TABLE(L)%X(1)%VALUES(N)=PLD(NPC(L)+2*N-2)
281 TABLE(L)%Y%VALUES(N) =PLD(NPC(L)+2*N-1)
282 ENDDO
283 IF (IS_ENCRYPTED) THEN
284 WRITE(IOUT,'(a)')'confidential data'
285 FUNCRYPT(L) = 1
286 ENDIF
287 ENDDO
288 ENDIF
289
290C
291 RETURN
292C-----------------------------------------------------------------
2932000 FORMAT(//
294 . ' load curves'/
295 . ' -----------'/
296 . ' number of load curves. . . . . . . . =',I10/)
2972100 FORMAT(/' load curve id . . . . . . . . . . . =',I10//
298 . ' x y ')
2992200 FORMAT(/' load smooth curve id . . . . . . . =',I10)
3002300 FORMAT(/' xscale yscale xshift
301 . yshift ')
3022400 FORMAT(/' x y ')
303 END
my_real
#define my_real
Definition cppsort.cpp:32
hm_get_float_array_index
subroutine hm_get_float_array_index(name, rval, index, is_available, lsubmodel, unitab)
Definition hm_get_float_array_index.F:105
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:528
hm_option_count
subroutine hm_option_count(entity_type, hm_option_number)
Definition hm_option_count.F:105
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:288
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:253
hm_read_funct
subroutine hm_read_funct(npc, pld, nfunct, table, npts_alloc, nom_opt, funcrypt, unitab, lsubmodel)
Definition hm_read_funct.F:45
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:414
id
initmumps id
Definition multiplerhs_example.m:4
hm_option_read_mod
Definition hm_option_read_mod.F:412
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008
submodel_mod
Definition submodel_mod.F:717
table_mod
Definition table_mod.F:113
unitab_mod
Definition unitab_mod.F:517
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:895
fretitl
subroutine fretitl(titr, iasc, l)
Definition freform.F:615
submodel_mod::submodel_data
Definition submodel_mod.F:728
table_mod::ttable
Definition table_mod.F:126
unitab_mod::unit_type_
Definition unitab_mod.F:524