anchor.diff:      if( ndof==2 ) mdof = 3
       if( ndof==3 ) mdof = 6
       if( ndof==6 ) mdof = 12
  -    if( ndof==6 ) shell33_coeff = 2
  -
  -    do it = 1, hecMESH%material%n_mat
  -      tmp = fstrSOLID%materials(it)%totallyr
  -      if (tmp >= ntot_lyr)then
  -        ntot_lyr = tmp
  -      endif
  -    enddo
  -    do it=1,hecMESH%n_elem_type
  -      j =  hecMESH%elem_type_item(it)
  -      if (j == 641)then
  -        if (mixedflag == 0) then
  -          mixedflag = 2
  -        endif
  -      elseif (j == 781 .or. j == 761)then
  -        mixedflag = 1
  -        shell33_coeff=2
  -      endif
  -      if (mixedflag == 1)then
  -        mdof = 12
  -      endif
  -    enddo
  -    n_lyr = 1
  +
  +    if (is_33shell == 1)then
  +      mdof = 12
  +    endif
   
       call hecmw_nullify_result_data( fstrRESULT )
       ncomp = 0
  
positive.code:     type (hecmwST_result_data):: fstrRESULT
      real(kind=kreal), pointer :: tnstrain(:), testrain(:)
      character(len=4) :: char_n_lyr
      character(len=4) :: char_num
      real(kind=kreal), allocatable   ::unode781(:),rnode781(:)
  
      integer(kind=kint) :: n_lyr, ntot_lyr, mixedflag, tmp, it, shell33_coeff
      integer(kind=kint) :: i, j, k, ndof, mdof, ncomp, nitem, iitem, ecomp, eitem, jitem, nn, mm
  
      tnstrain => fstrSOLID%TNSTRAIN
      testrain => fstrSOLID%TESTRAIN
  
      tmp = 0
      ntot_lyr = 1
      mixedflag = 0
      j = 0
      shell33_coeff=1;
  
      mm = hecMESH%n_node
      if( hecMESH%n_elem>hecMESH%n_node ) mm = hecMESH%n_elem
  
negative.code:     ecomp = 0
      eitem = 0
  
  ! --- COUNT SUM OF ALL NITEM
  ! --- DISPLACEMENT
          if( fstrSOLID%output_ctrl(4)%outinfo%on(1) ) then
            ncomp = ncomp + 1
            nitem = nitem + n_comp_valtype( fstrSOLID%output_ctrl(4)%outinfo%vtype(1), ndof )
          endif
  ! --- ROTATION (Only for 781 shell)
          if( fstrSOLID%output_ctrl(4)%outinfo%on(18) .and. is_33shell == 1 ) then
            ncomp = ncomp + 1
            nitem = nitem + n_comp_valtype( fstrSOLID%output_ctrl(4)%outinfo%vtype(18), ndof )
          endif
  