44
45
46
47
48
49
50
57 USE format_mod , ONLY : lfield
58
59
60
61#include "implicit_f.inc"
62
63
64
65#include "com04_c.inc"
66#include "units_c.inc"
67
68
69
70 INTEGER,INTENT(IN) :: ICODE(NUMNOD), ISKEW(*), ITAB(NUMNOD), ITABM1(*), IKINE(*),IBCSLAG(5,*)
71 INTEGER,INTENT(IN) :: LAG_NCF,LAG_NKF,LAG_NHF,IKINE1LAG(*)
72 INTEGER,INTENT(INOUT) :: LINALE(*)
73 TYPE(), DIMENSION(NSUBMOD), INTENT(IN) :: LSUBMODEL
74 TYPE(UNIT_TYPE_), INTENT(IN) :: UNITAB
75
76
77
78 INTEGER GRID_DOF, BID, IFORM, NODE_ID1, NODE_ID2,POS,IAD0,NNOD
79 INTEGER I,JJ(12), IC, NC, N, NUSR, IS, IC1, IC2, IC3, IC4,
80 . NOSYS, J,J10(10),IGR,IGRS,ISU,IBCALE,J3(3),K,
81 . IC0, IC01, IC02, IC03, IC04, ID ,ILAGM, NBCSLAG,
82 . FLAG_FMT,FLAG_FMT_TMP,IFIX_TMP,IKINE1(3*NUMNOD),IGRNS
83 CHARACTER(LEN=NCHARKEY) :: KEY,KEY2
84 CHARACTER(LEN=NCHARFIELD) :: STRING, CHAR_XYZ
85 CHARACTER(LEN=NCHARTITLE) :: TITR
86 CHARACTER :: OPT*8,MESS*40
88 LOGICAL :: IS_AVAILABLE
89 INTEGER :: WX, WY, WZ
90
91 TYPE (GROUP_) ,TARGET, DIMENSION(NGRNOD) :: IGRNOD
92
93
94
95 INTEGER USR2SYS,MY_OR,CHECK_NEW,NGR2USR
96
97 INTEGER, DIMENSION(:), POINTER :: INGR2USR
98
99
100
101 pos=0
102 igrns=0
103 tstart=zero
104 tstop=zero
105
106 WRITE(iout,1000)
108
109 DO i = 1, nalelk
111 CALL hm_get_intv(
'node_ID1', node_id1, is_available, lsubmodel)
112 CALL hm_get_intv(
'node_ID2', node_id2, is_available, lsubmodel)
113 CALL hm_get_intv(
'grnod_ID', igr, is_available, lsubmodel)
114 CALL hm_get_intv(
'Wx', wx, is_available, lsubmodel)
115 CALL hm_get_intv(
'Wy', wy, is_available, lsubmodel)
116 CALL hm_get_intv(
'Wz', wz, is_available, lsubmodel)
117 CALL hm_get_intv(
'Iform', iform, is_available, lsubmodel)
118 CALL hm_get_floatv(
'Tstart', tstart, is_available, lsubmodel, unitab)
119 CALL hm_get_floatv(
'Tstop', tstop, is_available, lsubmodel, unitab)
120 IF (tstop == zero) tstop = ep30
121 ingr2usr => igrnod(1:ngrnod)%ID
122 igrs=
ngr2usr(igr,ingr2usr,ngrnod)
123 IF(igrs==0)THEN
124 WRITE(iout,*) ' -- NODE GROUP',igr,' NOT FOUND'
125 WRITE(istdo,*) ' -- NODE GROUP',igr,' NOT FOUND'
126 ENDIF
127
128 j3(1) = wx
129 j3(2) = wy
130 j3(3) = wz
131 ic = j3(1) * 4 + j3(2) * 2 + j3(3)
132 char_xyz = ' '
133 k=lfield
134 IF(j3(3)==1 )THEN
135 char_xyz(k:k)='Z'
136 k=k-1
137 ENDIF
138 IF(j3(2)==1 )THEN
139 char_xyz(k:k)='Y'
140 k=k-1
141 ENDIF
142 IF(j3(1)==1 )THEN
143 char_xyz(k:k)='X'
144 ENDIF
145
146 igrns = igrnod(igrs)%SORTED
147
148 IF(iform==0.AND.igrns /= 1)THEN
149 CALL ancmsg(msgid=271,anmode=aninfo,msgtype=msgerror,i1=n,i2=iform)
150 ELSE
151 WRITE(iout,1100) n,trim(titr),node_id1,node_id2,igr,char_xyz,iform,tstart,tstop
152 ENDIF
153
154 linale(pos+1)=n
155 linale(pos+2)=node_id1
156 linale(pos+3)=node_id2
157 linale(pos+4)=-1
158 linale(pos+5)=ic
159 linale(pos+6)=iform
160 linale(pos+7)=igrs
161
162 pos = pos+1+6
163 mess=''
164 mess(1:10)='ALE LINKS'
165 bid=
usr2sys(node_id1,itabm1,mess,n)
166 bid=
usr2sys(node_id2,itabm1,mess,n)
167 ENDDO
168
1691000 FORMAT(/
170 . ' ALE LINKS DEFINITIONS '/
171 . ' ---------------------- '/)
172
1731100 FORMAT( /5x,'ALE LINK ID ',i10,': ',a,
174 . /10x,'MAIN NODE 1 . . . . . . . . . . . . . ',i10
175 . /10x,'
main node 2 . . . . . . . . . . . . .
',I10
176 . /10X,'group identifier
for secondary nodes. . . .
',I10
177 . /10X,'grid velocity directions to link. . . . ',A
178 . /10X,'formulation . . . . . . . . . . . . . . ',I10
179 . /10X,'start time. . . . . . . . . . . . . . . ',1PG20.13
180 . /10X,'stop time . . . . . . . . . . . . . . . ',1PG20.13/)
181
182 RETURN
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
subroutine hm_option_start(entity_type)
for(i8=*sizetab-1;i8 >=0;i8--)
integer, parameter nchartitle
integer, parameter ncharkey
integer, parameter ncharfield
integer function ngr2usr(iu, igr, ngr)
int main(int argc, char *argv[])
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)
integer function usr2sys(iu, itabm1, mess, id)