source: PSPA/parmelaPSPA/trunk/change.f @ 405

      subroutine change
c---change card parameters are np, ne, value, where value is the
c---new value assigned to parameter number np for element number ne.
c-------------------------------------------------------------------
      save
c
      include 'param_sz.h'
      include 'constcom.h'
      include 'misccom.h'
      include 'syscom.h'
      include 'var_char.h'
      include 'ucom.h'
c
c--------------------------------------------------------------------------
c*
      np=vv(1)
      ne=vv(2)
      if(np.eq.1)then
      write(ndiag,*)' the length of and element can not be changed'
      call appendparm
      stop ' Abnormal stop change '
      endif
      el(np,ne)=vv(3)
      if(ntype(ne).eq.7.and.np.eq.4)then
c--this is a cell so change el parm 9 and 10
      phi=el(4,ne)*radian
      el(9,ne)=sin(phi)
      el(10,ne)=cos(phi)
      endif
      if(ntype(ne).eq.4)then
c--This is a bend so convert degress to radians
           if(np.ge.5.and.np.le.9)el(np,ne)=el(np,ne)*radian
           if(np.eq.4)then
               gr=el(4,ne)/erest
               bg=sqrt(gr*(2.+gr))
               el(10,ne)=bg
               el(4,ne)=el(1,ne)/el(5,ne)
           endif
           if(np.eq.5)el(4,ne)=el(1,ne)/el(5,ne)
      endif
      return
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*