FEAP User Forum

FEAP => nurbFEAP => Topic started by: resammc on July 28, 2017, 02:00:30 AM

Title: OUTMesh command (bug)
Post by: resammc on July 28, 2017, 02:00:30 AM

Dear FEAP community,

The OUTMesh command in nurbfeap does not produce the mesh output properly (it lacks correct side information) most of the times. For instance, the side data for the OUTMesh outcome of the attached input file is:

Quote

NSIDes        6
  side   1   4   1     0     0     0     0
  side   2   4   1     0     0     0     0
  side   3   4   1     0     0     0     0
  side   4   3   2     0     0     0
  side   5   3   2     0     0     0
  side   6   3   2     0     0     0

which lacks control points data.

The problem appears to be in poutm.f file. I think

Quote

      include   'cdata.h'

      logical    setvar,palloc

      ...

      if(np(79).eq.0) then
        setvar = palloc( 79,'IPOS ',numnp,  1)
        call pseqn(mr(np(79)),numnp)
      endif
 
      ...

must be added before calling poutsid subroutine in outm_nurbs subroutine.

Quote

NSIDes        6
  side   1   4   1     1     2     3     4
  side   2   4   1     5     6     7     8
  side   3   4   1     9    10    11    12
  side   4   3   2     1     5     9
  side   5   3   2     2     6    10
  side   6   3   2     3     7    11

or, the user must use a command which activates the np(79) array (e.g. TIE) before calling OUTMesh.

Kind regards,

RM