anchor.diff:            enddo
             iitem = iitem + nn
           endif
  -      endif
  -
   
   ! --- STRESS @node
  -      if(fstrSOLID%output_ctrl(4)%outinfo%on(4)) then
  -
  -        if(ndof==6 .or. mixedflag==1 ) then
  -          nn = 6
  -          ncomp = ncomp + 1
  -          fstrRESULT%nn_dof(ncomp) = nn
  -          fstrRESULT%node_label(ncomp) = 'NodalSTRESS_Layer'//trim(char_n_lyr)//'+'
  -          do i = 1, hecMESH%n_node
  -            do j = 1, 6
  -              fstrRESULT%node_val_item(nitem*(i-1)+j+iitem) = fstrSOLID%STRESS(12*ntot_lyr*(i-1)+12*(n_lyr-1)+j)
  -            enddo
  -          enddo
  -          iitem = iitem + nn
  -
  -          ncomp = ncomp + 1
  -          fstrRESULT%nn_dof(ncomp) = nn
  -          fstrRESULT%node_label(ncomp) = 'NodalSTRESS_Layer'//trim(char_n_lyr)//'-'
  -          do i = 1, hecMESH%n_node
  -            do j = 1, 6
  -              fstrRESULT%node_val_item(nitem*(i-1)+j+iitem) = fstrSOLID%STRESS(12*ntot_lyr*(i-1)+12*(n_lyr-1)+j+6)
  -            enddo
  -          enddo
  -          iitem = iitem + nn
  -
  -        else
  +        if(fstrSOLID%output_ctrl(4)%outinfo%on(4)) then
             ncomp = ncomp + 1
             nn = n_comp_valtype( fstrSOLID%output_ctrl(4)%outinfo%vtype(4), ndof )
             fstrRESULT%nn_dof(ncomp) = nn
  
positive.code:     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
  
      call hecmw_nullify_result_data( fstrRESULT )
      ncomp = 0
  
negative.code:           enddo
            iitem = iitem + nn
          endif
  
  ! --- Princ STRESS Vector @node
          if(fstrSOLID%output_ctrl(4)%outinfo%on(23)) then
            call make_principal_stress_value(fstrSOLID, hecMESH, 1, ndof)
            do k=1,3
              write(cnum, '(i0)') k
              ncomp = ncomp + 1
              nn = n_comp_valtype( fstrSOLID%output_ctrl(4)%outinfo%vtype(23), ndof )
              fstrRESULT%nn_dof(ncomp) = nn
              fstrRESULT%node_label(ncomp) = 'NodalPrincipalSTRESSVector'//'_'//trim(cnum)
              do i = 1, hecMESH%n_node
                do j = 1, nn
                  fstrRESULT%node_val_item(nitem*(i-1)+j+iitem) = fstrSOLID%PSTRESS_VECT((nn)*(i-1)+j,k)
                enddo
              enddo
              iitem = iitem + nn
            end do
          endif
  
  ! --- Princ STRAIN @node
          if( fstrSOLID%output_ctrl(4)%outinfo%on(21)) then
            call make_principal_strain_value(fstrSOLID, hecMESH, 1, ndof)
            nn = n_comp_valtype( fstrSOLID%output_ctrl(4)%outinfo%vtype(21), ndof )
            ncomp = ncomp + 1
            fstrRESULT%nn_dof(ncomp) = nn
  