#include "implicit_f.inc"
#include "vect01_c.inc"
#include "param_c.inc"
#include "scr17_c.inc"
Go to the source code of this file.
|
| subroutine | m24insph (pm, ang, sf, vk0, rob, ipart, ipartsp, nel) |
◆ m24insph()
| subroutine m24insph |
( |
| pm, |
|
|
| ang, |
|
|
| sf, |
|
|
| vk0, |
|
|
| rob, |
|
|
integer, dimension(lipart1,*) | ipart, |
|
|
integer, dimension(*) | ipartsp, |
|
|
integer | nel ) |
Definition at line 29 of file m24insph.F.
31
32
33
34#include "implicit_f.inc"
35
36
37
38#include "vect01_c.inc"
39#include "param_c.inc"
40#include "scr17_c.inc"
41
42
43
44 INTEGER IPART(LIPART1,*),IPARTSP(*),NEL
45 my_real pm(npropm,*), ang(nel,6), sf(nel,3), vk0(*), rob(*)
46
47
48
49 INTEGER I, IPRT, MX, IG, ISK
50
51 DO i=lft,llt
52 ang(i,1)=one
53 ang(i,2)=zero
54 ang(i,3)=zero
55 ang(i,4)=zero
56 ang(i,5)=one
57 ang(i,6)=zero
58 ENDDO
59
60 DO i=lft,llt
61 sf(i,1)=-one
62 sf(i,2)=-one
63 sf(i,3)=-one
64 ENDDO
65
66 DO i=lft,llt
67 iprt=ipartsp(i+nft)
68 mx=ipart(1,iprt)
69 rob(i)=pm(30,mx)
70 vk0(i)=pm(46,mx)
71 ENDDO
72
73 RETURN