34
35
36
39
40
41
42#include "implicit_f.inc"
43
44
45
46#include "impl1_c.inc"
47
48
49
50 INTEGER, INTENT(IN) :: MTN
51 INTEGER JFT, JLT ,IR,IS,IT
52
54 . ht(*)
55
56
57
58 INTEGER I,MX,J,K
59
60
61 TYPE(L_KTBUFEP_) , POINTER :: LBUF
62 TYPE(MLAW_TAG_) , POINTER :: MTAG
63
65 lbuf =>
ktbuf_str(ng_imp)%KTBUFEP(ir,is,it)
66 IF (mtag%L_A_KT>0) THEN
67 DO i=jft,jlt
68 ht(i)= lbuf%A_KT(i)
69 IF (ht(i)==one) ht(i)= zero
70 ENDDO
71 ELSE
72 DO i=jft,jlt
73 ht(i)= zero
74 ENDDO
75 END IF
76
77 RETURN
type(ktbuf_struct_), dimension(:), allocatable, target ktbuf_str