OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
section_init.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!|| section_init ../engine/source/tools/sect/section_init.F
25!||--- called by ------------------------------------------------------
26!|| resol_init ../engine/source/engine/resol_init.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| close_c ../common_source/tools/input_output/write_routines.c
30!|| cur_fil_c ../common_source/tools/input_output/write_routines.c
31!|| open_c ../common_source/tools/input_output/write_routines.c
32!|| read_r_c ../common_source/tools/input_output/write_routines.c
33!|| spmd_rbcast ../engine/source/mpi/generic/spmd_rbcast.F
34!||--- uses -----------------------------------------------------
35!|| inoutfile_mod ../common_source/modules/inoutfile_mod.f
36!|| message_mod ../engine/share/message_module/message_mod.F
37!|| names_and_titles_mod ../common_source/modules/names_and_titles_mod.F
38!||====================================================================
39 SUBROUTINE section_init (NSTRF,SECBUF,NOM_SECT,ISECTR,NSECT,IOLDSECT)
40C-----------------------------------------------
41C M o d u l e s
42C-----------------------------------------------
43 USE message_mod
44 USE inoutfile_mod
45 USE names_and_titles_mod, ONLY: nchartitle
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C C o m m o n B l o c k s
52C-----------------------------------------------
53#include "comlock.inc"
54#include "com08_c.inc"
55#include "chara_c.inc"
56#include "warn_c.inc"
57#include "task_c.inc"
58#include "com01_c.inc"
59C-----------------------------------------------
60C D u m m y A r g u m e n t s
61C-----------------------------------------------
62 INTEGER NSTRF(*),NOM_SECT(*),ISECTR,NSECT,IOLDSECT
63 my_real
64 . secbuf(*)
65C-----------------------------------------------
66C L o c a l V a r i a b l e s
67C-----------------------------------------------
68 INTEGER I,J, IR1, IR2, IFILNAM(2548),LROOTLEN,LEN,K0
69 my_real
70 . tt1, tt2, tt3, tmp(20)
71 CHARACTER(LEN=NCHARTITLE) :: FILNAM
72 CHARACTER LCHRUN*2,LCHRUN_P1*2,CH_IDSEC*10
73 LOGICAL FEXIST
74 real*4 r4
75 INTEGER :: LEN_TMP_NAME
76 CHARACTER(len=2048) :: TMP_NAME
77C-----------------------------------------------
78 tt1 = secbuf(2)
79 tt2 = secbuf(3)
80 tt3 = secbuf(4)
81 ir1=nstrf(4)
82 ir2=nstrf(5)
83C-----------------------------------------------
84C WRITE FILE
85C-----------------------------------------------
86 IF (ispmd==0) THEN
87 k0 = nstrf(25)
88 DO j=1,nsect
89 IF(nstrf(k0)>=1 .AND. nstrf(k0)<=10 )THEN
90 WRITE(lchrun,'(I2.2)')irun
91 lrootlen=0
92 DO i=1,500
93c
94 filnam(i:i)=char(nom_sect((j-1)*500+i))
95 IF(filnam(i:i)/=' ')LROOTLEN=LROOTLEN+1
96 ENDDO
97.AND. IF (LROOTLEN == 0 ABS(IOLDSECT) == 1) THEN
98 IOLDSECT = 1
99.AND. ELSEIF( LROOTLEN /= 0 (IOLDSECT >= 1)) THEN
100 IOLDSECT = 2
101 ENDIF
102 ENDIF
103 K0 = NSTRF(K0+24)
104 ENDDO
105.AND. IF(NSTRF(1)>=1 IOLDSECT == 1)THEN
106C FILNAM=ROOTNAM(1:ROOTLEN)//'SC'//CHRUN
107 WRITE(LCHRUN,'(i2.2)')IRUN
108 FILNAM=ROOTNAM(1:ROOTLEN)//'sc'//LCHRUN
109 LEN_TMP_NAME = OUTFILE_NAME_LEN + ROOTLEN + 4
110 TMP_NAME=OUTFILE_NAME(1:OUTFILE_NAME_LEN)//FILNAM(1:ROOTLEN+4)
111 DO I=1,LEN_TMP_NAME
112 IFILNAM(I)=ICHAR(TMP_NAME(I:I))
113 ENDDO
114 CALL CUR_FIL_C(42)
115 CALL OPEN_C(IFILNAM,LEN_TMP_NAME,0)
116 ELSEIF(NSTRF(1)>=1) THEN
117 K0 = NSTRF(25)
118 DO J=1,NSECT
119.AND. IF(NSTRF(K0)>=1 NSTRF(K0)<=10 )THEN
120 WRITE(LCHRUN,'(i2.2)')IRUN
121 LROOTLEN=0
122 DO I=1,500
123c
124 IF(CHAR(NOM_SECT((J-1)*500+I))/=' ')THEN
125 LROOTLEN=LROOTLEN+1
126 FILNAM(LROOTLEN:LROOTLEN)=CHAR(NOM_SECT((J-1)*500+I))
127 ENDIF
128 ENDDO
129 IF (LROOTLEN == 0) THEN
130 WRITE(CH_IDSEC,'(i10.10)')NSTRF(K0+23)
131 FILNAM=ROOTNAM(1:ROOTLEN)//CH_IDSEC//'sc'//LCHRUN
132 LEN_TMP_NAME = OUTFILE_NAME_LEN + ROOTLEN + 14
133 TMP_NAME=OUTFILE_NAME(1:OUTFILE_NAME_LEN)//FILNAM(1:ROOTLEN+14)
134 DO I=1,LEN_TMP_NAME
135 IFILNAM(I)=ICHAR(TMP_NAME(I:I))
136 ENDDO
137 CALL CUR_FIL_C(41+J)
138 CALL OPEN_C(IFILNAM,LEN_TMP_NAME,0)
139 ELSE
140 FILNAM=FILNAM(1:LROOTLEN)//'sc'//LCHRUN
141 LEN_TMP_NAME = OUTFILE_NAME_LEN + LROOTLEN + 4
142 TMP_NAME=OUTFILE_NAME(1:OUTFILE_NAME_LEN)//FILNAM(1:LROOTLEN+4)
143 DO I=1,LEN_TMP_NAME
144 IFILNAM(I)=ICHAR(TMP_NAME(I:I))
145 ENDDO
146 CALL CUR_FIL_C(41+J)
147 CALL OPEN_C(IFILNAM,LEN_TMP_NAME,0)
148 ENDIF
149 ENDIF
150 K0 = NSTRF(K0+24)
151 ENDDO
152 ENDIF
153C-----------------------------------------------
154C READ FILES
155C-----------------------------------------------
156 IF(NSTRF(2)>=1)THEN
157 LROOTLEN=0
158 DO I=1,500
159c
160 FILNAM(I:I)=CHAR(NOM_SECT((ISECTR-1)*500+I))
161 IF(FILNAM(I:I)/=' ')LROOTLEN=LROOTLEN+1
162 ENDDO
163 IF(TT==ZERO)THEN
164 TT1=ZERO
165 TT2=ZERO
166 FILNAM=FILNAM(1:LROOTLEN)//'sc01'
167 len_tmp_name = outfile_name_len + lrootlen + 4
168 tmp_name=outfile_name(1:outfile_name_len)//filnam(1:lrootlen+4)
169 INQUIRE(file=tmp_name(1:len_trim(tmp_name)),exist=fexist)
170
171 IF(.NOT.fexist) THEN
172 len_tmp_name = lrootlen + 4
173 tmp_name(1:len_tmp_name)=filnam(1:lrootlen+4)
174 INQUIRE(file=tmp_name(1:len_trim(tmp_name)),exist=fexist)
175 ENDIF
176
177 IF(fexist)THEN
178 CALL cur_fil_c(4)
179 DO i=1,len_tmp_name
180 ifilnam(i)=ichar(tmp_name(i:i))
181 ENDDO
182 CALL open_c(ifilnam,len_tmp_name,1)
183 CALL read_r_c(r4,1)
184 CALL close_c()
185 tt1=r4
186 ir1=1
187 ir2=1
188 ELSE
189 CALL ancmsg(msgid=188,anmode=aninfo,
190 . c1=filnam)
191 ierr=ierr+1
192 tstop = tt
193 ENDIF
194 filnam=filnam(1:lrootlen)//'SC02'
195 len_tmp_name = outfile_name_len + lrootlen + 4
196 tmp_name=outfile_name(1:outfile_name_len)//filnam(1:lrootlen+4)
197
198 INQUIRE(file=tmp_name(1:len_trim(tmp_name)),exist=fexist)
199
200 IF(.NOT.fexist) THEN
201 len_tmp_name = lrootlen + 4
202 tmp_name(1:len_tmp_name)=filnam(1:lrootlen+4)
203 INQUIRE(file=tmp_name(1:len_trim(tmp_name)),exist=fexist)
204 ENDIF
205
206 IF(fexist)THEN
207 ir2=2
208 CALL cur_fil_c(4)
209 DO i=1,len_tmp_name!LROOTLEN+4
210 ifilnam(i)=ichar(filnam(i:i))
211 ENDDO
212
213 CALL open_c(ifilnam,len_tmp_name,1)
214 CALL read_r_c(r4,1)
215 CALL close_c()
216 tt3=r4
217 ELSE
218 tt3=ep30
219 ENDIF
220 ELSE
221 tt1=zero
222 tt2=zero
223 WRITE(lchrun,'(I2.2)')irun
224 filnam=filnam(1:lrootlen)//'SC'//lchrun
225 len_tmp_name = outfile_name_len + len_trim(filnam)
226 tmp_name=outfile_name(1:outfile_name_len)//filnam(1:len_trim(filnam))
227
228 INQUIRE(file=tmp_name(1:len_trim(tmp_name)),exist=fexist)
229
230 IF(.NOT.fexist) THEN
231 len_tmp_name = len_trim(filnam)
232 tmp_name(1:len_tmp_name)=filnam(1:len_trim(filnam))
233 INQUIRE(file=tmp_name(1:len_trim(tmp_name)),exist=fexist)
234 ENDIF
235
236 IF(fexist)THEN
237 CALL cur_fil_c(4)
238 DO i=1,len_tmp_name
239 ifilnam(i)=ichar(tmp_name(i:i))
240 ENDDO
241
242 CALL open_c(ifilnam,len_tmp_name,1)
243 CALL read_r_c(r4,1)
244 CALL close_c()
245 tt1=r4
246 ir1=irun
247 ir2=irun
248 ELSE
249 CALL ancmsg(msgid=188,anmode=aninfo,
250 . c1=filnam)
251 ierr=ierr+1
252 tstop = tt
253 ENDIF
254 WRITE(lchrun_p1,'(I2.2)')irun+1
255 filnam=filnam(1:lrootlen)//'SC'//lchrun_p1
256 len_tmp_name = outfile_name_len + len_trim(filnam)
257 tmp_name=outfile_name(1:outfile_name_len)//filnam(1:len_trim(filnam))
258
259 INQUIRE(file=tmp_name,exist=fexist)
260
261 IF(.NOT.fexist) THEN
262 len_tmp_name = len_trim(filnam)
263 tmp_name(1:len_tmp_name)=filnam(1:len_tmp_name)
264 INQUIRE(file=tmp_name(1:len_trim(tmp_name)),exist=fexist)
265 ENDIF
266
267 IF(fexist)THEN
268 ir2=irun + 1
269 CALL cur_fil_c(4)
270 DO i=1,len_tmp_name
271 ifilnam(i)=ichar(tmp_name(i:i))
272 ENDDO
273
274 CALL open_c(ifilnam,len_tmp_name,1)
275 CALL read_r_c(r4,1)
276 CALL close_c()
277 tt3=r4
278 ELSE
279 tt3=ep30
280 ENDIF
281 ENDIF
282C
283 WRITE(lchrun,'(I2.2)')ir1
284 filnam=filnam(1:lrootlen)//'SC'//lchrun
285 len_tmp_name = outfile_name_len + len_trim(filnam)
286 tmp_name=outfile_name(1:outfile_name_len)//filnam(1:len_trim(filnam))
287
288 CALL cur_fil_c(4)
289 DO i=1,len_tmp_name
290 ifilnam(i)=ichar(tmp_name(i:i))
291 ENDDO
292
293 CALL open_c(ifilnam,len_tmp_name,1)
294 ENDIF
295 ENDIF
296C-----------------------------------------------
297 secbuf(2) = tt1
298 secbuf(3) = tt2
299 secbuf(4) = tt3
300C
301 nstrf(4) = ir1
302 nstrf(5) = ir2
303C
304C MAJ SPMD sur proc remote
305C
306 IF (nspmd > 1) THEN
307 IF(ispmd==0) THEN
308 tmp(1) = nstrf(4)
309 tmp(2) = nstrf(5)
310 tmp(3) = secbuf(2)
311 tmp(4) = secbuf(3)
312 tmp(5) = secbuf(4)
313 len = 5
314C 2*Len necessary for communication !!!
315 CALL spmd_rbcast(tmp,tmp,len,1,0,2)
316 ELSE
317 len = 5
318C 2*Len necessary for communication !!!
319 CALL spmd_rbcast(tmp,tmp,len,1,0,2)
320 nstrf(4) = tmp(1)
321 nstrf(5) = tmp(2)
322 secbuf(2)= tmp(3)
323 secbuf(3)= tmp(4)
324 secbuf(4)= tmp(5)
325 ENDIF
326 END IF
327C
328 RETURN
329 END
my_real
#define my_real
Definition cppsort.cpp:32
inoutfile_mod
Definition inoutfile_mod.F:96
inoutfile_mod::outfile_name
character(len=outfile_char_len) outfile_name
Definition inoutfile_mod.F:111
inoutfile_mod::outfile_name_len
integer outfile_name_len
Definition inoutfile_mod.F:109
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
section_init
subroutine section_init(nstrf, secbuf, nom_sect, isectr, nsect, ioldsect)
Definition section_init.F:40
spmd_rbcast
subroutine spmd_rbcast(tabi, tabr, n1, n2, from, add)
Definition spmd_rbcast.F:62
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
close_c
void close_c()
Definition write_routines.c:298
cur_fil_c
void cur_fil_c(int *nf)
Definition write_routines.c:477
read_r_c
void read_r_c(float *w, int *len)
Definition write_routines.c:756
open_c
void open_c(int *ifil, int *len, int *mod)
Definition write_routines.c:180