Error found in element therm2d
In Feap-pv 5.1, element therm2d.f, the write statements below causes a program crash since gradt contains three values and the format 2003 only accepts two.
write(iow,2003) n_el,ma,xx,flux(1:2,l),gradt
write(*,2003) n_el,ma,xx,flux(1:2,l),gradt
This is easily fixed:
write(iow,2003) n_el,ma,xx,flux(1:2,l),gradt(1:2)
write(*,2003) n_el,ma,xx,flux(1:2,l),gradt(1:2)
However, after the fix, output of element fluxes for gauss point 1 seems OK but for point 2-4 is zero
Sample output from Feap 5.1
elmt matl 1-coord 2-coord 1-flux 2-flux 1-grad 2-grad
350 1 33.634 78.634 -0.786E-10 -0.819E-08 0.786E-02 0.819E+00
350 1 35.366 78.634 0.000E+00 0.000E+00 0.786E-02 0.813E+00
350 1 35.366 80.366 0.000E+00 0.000E+00 0.211E-02 0.813E+00
350 1 33.634 80.366 0.000E+00 0.000E+00 0.211E-02 0.819E+00
Run with version 4.1 I gives the following result
350 1 33.634 78.634 -0.786E-10 -0.819E-08 0.786E-02 0.819E+00
350 1 35.366 78.634 -0.786E-10 -0.813E-08 0.786E-02 0.813E+00
350 1 35.366 80.366 -0.211E-10 -0.813E-08 0.211E-02 0.813E+00
350 1 33.634 80.366 -0.211E-10 -0.819E-08 0.211E-02 0.819E+00
Kind regards
Björn Gustafsson
