31
32
33
36
37
38
39#include "implicit_f.inc"
40
41
42
43#include "com04_c.inc"
44
45
46
48 INTEGER, INTENT(IN) :: ILAS(*)
49
50
51
52 LOGICAL :: OK_QA
53 CHARACTER (LEN=255) :: VARNAME
54 INTEGER K,J,DET_IID,I1,I2,I3,I4,I5,I6
55 DOUBLE PRECISION TEMP_DOUBLE
56 INTEGER :: NL,K1,K2,IFUNC,ITYP, IFK, MFK, IAFUNC, NCELL
57 DOUBLE PRECISION :: CHALEUR, FSI, XK0, HNUK, DAR_HNUK_2,ONE_TIMESCAL, FI, ALPHA
58
59
60
61
63
64 k1 = 1
65 k2 = 1
66
67 DO k=1,nlaser
68
70 ifunc = ilas(k1+3)
71 iafunc = ilas(k1+4)
72 ityp = ilas(k1+5)
73
74 chaleur = xlas(k2)
75 fi = xlas(k2+1)
77
78 xk0 = xlas(k2+3)
79 hnuk = xlas(k2+4)
80 dar_hnuk_2 = xlas(k2+5)
81 one_timescal = xlas(k2+6)
82
83
84 WRITE(varname,'(A,I0,A,I0,A)') 'LASER(',k ,')%ILAS(',k1,')=NL :'
85 CALL qaprint(varname(1:len_trim(varname)),ilas(k1),0.0_8)
86 WRITE(varname,'(A,I0,A,I0,A)') 'LASER(',k ,')%ILAS(',k1+3,')=IFUNC :'
87 CALL qaprint(varname(1:len_trim(varname)),ilas(k1+3),0.0_8)
88 WRITE(varname,'(A,I0,A,I0,A)') 'LASER(',k ,')%ILAS(',k1+4,')=IAFUNC :'
89 CALL qaprint(varname(1:len_trim(varname)),ilas(k1+4),0.0_8)
90 WRITE(varname,'(A,I0,A,I0,A)') 'LASER(',k ,')%ILAS(',k1+5,')=ITYP :'
91 CALL qaprint(varname(1:len_trim(varname)),ilas(k1+5),0.0_8)
92
93 ncell=0
95 ncell=ncell+1
96 WRITE(varname,'(A,I0,A,I0,A,I0,A)') 'LASER(',k ,')%ILAS(',k1+5+j,')=ELEM(',ncell,') :'
97 CALL qaprint(varname(1:len_trim(varname)),ilas(k1+5+j),0.0_8)
98 WRITE(varname,'(A,I0,A,I0,A)') 'LASER(',k ,')%ILAS(',k1+5+j+1,')=NG :'
99 CALL qaprint(varname(1:len_trim(varname)),ilas(k1+5+j+1),0.0_8)
100 ENDDO
101
102 WRITE(varname,'(A,I0,A,I0,A,I0,A)') 'LASER(',k ,')%XLAS(',k2,')=CHALEUR :'
103 CALL qaprint(varname(1:len_trim(varname)),0,chaleur
104 WRITE(varname,'(A,I0,A,I0,A,I0,A)') 'laser(',K ,')%XLAS(',K2+1,')=fi :'
105 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),0,FI)
106 WRITE(VARNAME,'(a,i0,a,i0,a,i0,a)
') 'laser(
',K ,')%XLAS(
',K2+2,')=
alpha :
'
107 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),0,ALPHA)
108
109 WRITE(VARNAME,'(a,i0,a,i0,a,i0,a)') 'laser(',K ,')%XLAS(',K2+3,')=xk0 :'
110 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),0,XK0)
111 WRITE(VARNAME,'(a,i0,a,i0,a,i0,a)') 'laser(',K ,')%XLAS(',K2+4,')=hnuk :'
112 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),0,HNUK)
113 WRITE(VARNAME,'(a,i0,a,i0,a,i0,a)') 'laser(',K ,')%XLAS(',K2+5,')=dar_hnuk_2 :'
114 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),0,DAR_HNUK_2)
115 WRITE(VARNAME,'(a,i0,a,i0,a,i0,a)') 'laser(',K ,')%XLAS(',K2+6,')=un_timescal :'
116 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),0,ONE_TIMESCAL)
117
118 IFK = 6 + 2*(NL+1)
119 MFK = 7
120 K1 = K1 + IFK
121 K2 = K2 + MFK
122
123 ENDDO !next K
124 ENDIF
125
126 RETURN
logical function myqakey(value)
@purpose Check if a given value is part of the values set by env variable Useful to make a condition ...
subroutine qaprint(name, idin, value)
@purpose print one entry to QA extract file example of call for real print CALL QAPRINT('MY_LABEL',...
character *2 function nl()