#include "implicit_f.inc"
#include "mvsiz_p.inc"
Go to the source code of this file.
|
| subroutine | i9imp3 (iloc, cst, nor, xface, n1, n2, n3, ssc, ttc, thk, ans, lft, llt, nft) |
◆ i9imp3()
| subroutine i9imp3 |
( |
integer, dimension(*) | iloc, |
|
|
| cst, |
|
|
| nor, |
|
|
intent(in) | xface, |
|
|
intent(in) | n1, |
|
|
intent(in) | n2, |
|
|
intent(in) | n3, |
|
|
intent(in) | ssc, |
|
|
intent(in) | ttc, |
|
|
intent(in) | thk, |
|
|
intent(inout) | ans, |
|
|
integer, intent(inout) | lft, |
|
|
integer, intent(inout) | llt, |
|
|
integer, intent(inout) | nft ) |
Definition at line 28 of file i9imp3.F.
33
34
35
36#include "implicit_f.inc"
37
38
39
40#include "mvsiz_p.inc"
41
42
43
44 INTEGER, INTENT(INOUT) :: LFT
45 INTEGER, INTENT(INOUT) :: LLT
46 INTEGER, INTENT(INOUT) :: NFT
47 INTEGER ILOC(*)
49 my_real,
DIMENSION(MVSIZ),
INTENT(IN) :: xface,ssc,ttc,thk
50 my_real,
DIMENSION(MVSIZ),
INTENT(IN) :: n1,n2,n3
51 my_real,
DIMENSION(MVSIZ),
INTENT(INOUT) ::ans
52
53
54
55
56
57
58 INTEGER I, IL
59
60
61
62 DO i=lft,llt
63 il=i+nft
64 ans(i)=
min(zero,(ans(i)*xface(i)-thk(i)))
65 IF(ans(i)*xface(i)==zero)iloc(il)=-iloc(il)
66 ENDDO
67
68 DO i=lft,llt
69 il=i+nft
70 nor(1,il)=n1(i)
71 nor(2,il)=n2(i)
72 nor(3,il)=n3(i)
73 cst(1,il)=ssc(i)
74 cst(2,il)=ttc(i)
75 ENDDO
76
77 RETURN