32 . PM_STACK ,GEO_STACK ,IGEO_STACK)
41#include "implicit_f.inc"
51 INTEGER,
INTENT(IN) :: IGEO(NPROPGI,NUMGEO),IGEO_STACK(4* NPT_STACK+2,NS_STACK)
52 my_real,
INTENT(IN) :: geo(npropg,numgeo),geo_stack(6*npt_stack+1,ns_stack),
53 . pm_stack(20,ns_stack),bufgeo(*)
57 INTEGER I, MY_ID, MY_PID,
58 . iadbuf, nuparam, iadmat, njmat, iadpid, njpid, iadfun, njfun, iadtab, njtab
59 CHARACTER(LEN=NCHARTITLE) :: TITR
60 CHARACTER (LEN=255) :: VARNAME
61 DOUBLE PRECISION TEMP_DOUBLE
67 CALL fretitl2(titr,igeo(npropgi-ltitr+1,my_pid),ltitr)
69 IF(len_trim(titr)/=0)
THEN
70 CALL qaprint(titr(1:len_trim(titr)),igeo(1,my_pid),0.0_8)
72 CALL qaprint(
'A_PID_FAKE_NAME',igeo(1,my_pid),0.0_8)
75 IF(igeo(i,my_pid)/=0)
THEN
78 WRITE(varname,
'(A,I0)')
'IGEO_',i
79 CALL qaprint(varname(1:len_trim(varname)),igeo(i,my_pid),0.0_8)
83 IF(geo(i,my_pid)/=zero)
THEN
86 WRITE(varname,
'(A,I0)')
'GEO_',i
87 temp_double = geo(i,my_pid)
88 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
96 iadbuf =igeo(57,my_pid)
97 nuparam=igeo(52,my_pid)
99 IF(bufgeo(iadbuf+i-1)/=zero)
THEN
102 WRITE(varname,
'(A,I0)')
'BUFGEO_IADBUF_',i
103 temp_double = bufgeo(iadbuf+i-1)
104 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
107 iadfun =igeo(58,my_pid)
108 njfun =igeo(53,my_pid)
110 IF(bufgeo(iadfun+i-1)/=zero)
THEN
113 WRITE(varname,
'(A,I0)')
'BUFGEO_IADFUN_',i
114 temp_double = bufgeo(iadfun+i-1)
115 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
118 iadmat =igeo(59,my_pid)
119 njmat =igeo(54,my_pid)
121 IF(bufgeo(iadmat+i-1)/=zero)
THEN
124 WRITE(varname,
'(A,I0)')
'BUFGEO_IADMAT_',i
125 temp_double = bufgeo(iadmat+i-1)
126 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
129 iadpid =igeo(60,my_pid)
130 njpid =igeo(55,my_pid)
132 IF(bufgeo(iadpid+i-1)/=zero)
THEN
135 WRITE(varname,
'(A,I0)')
'BUFGEO_IADPID_',i
136 temp_double = bufgeo(iadpid+i-1)
137 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
140 iadtab =igeo(61,my_pid)
141 njtab =igeo(56,my_pid)
143 IF(bufgeo(iadtab+i-1)/=zero)
THEN
146 WRITE(varname,
'(A,I0)')
'BUFGEO_IADTAB_',i
147 temp_double = bufgeo(iadtab+i-1)
148 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
155 CALL qaprint(
'STACK_NAME_NO', my_id,0.0_8)
158 IF (pm_stack(i,my_pid) /= zero)
THEN
159 WRITE(varname,
'(A,I0)')
'STACK_PM_',i
160 temp_double = pm_stack(i,my_pid)
161 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
166 IF (geo_stack(i,my_pid) /= zero)
THEN
167 WRITE(varname,
'(A,I0)')
'STACK_GEO_',i
168 temp_double = geo_stack(i,my_pid)
169 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
174 IF (igeo_stack(i,my_pid) /= 0)
THEN
175 WRITE(varname,
'(A,I0)')
'STACK_IGEO_',i
176 CALL qaprint(varname(1:len_trim(varname)),igeo_stack(i,my_pid),0.0_8)