CSURface in FEAP86

[Jacopo]

Dear FEAP community,
I am experiencing an issue with the CSURface command in FEAP86.
I can not apply a uniform normal traction distribution over a straight edge in a 2D problem.
Issuing the command results in a wavy distribution, instead of constant. 
Please find below a MWE that reproduces the problem. I tested the same toy example in FEAP83 and FEAPpv51,
where the results are fine. I hope I am not missing something trivial,
Thanks for your help and consideration,
Jacopo

[Prof. S. Govindjee]

There is a bug in the code; the distribution you see is an axis-symmetric distribution. 

To work around, remove your statements for global plane strain and you will get what you want (plus you problem will be plane strain  as that is the default in FEAP).

I'll post a true fix later.

[Jacopo]

Dear Professor,
Thank you for your reply, now it works properly.
Jacopo

[Prof. R.L. Taylor]

The bug is in ./program/pesurf.f at about line 238

Delete:   axifl = gtypfl