Skip to content
Snippets Groups Projects
Commit fd4c2408 authored by Behrens, Prof. Dr. Jörn's avatar Behrens, Prof. Dr. Jörn
Browse files

still small bug fixes

parent faeaa570
No related branches found
No related tags found
No related merge requests found
......@@ -127,7 +127,7 @@ SUBROUTINE plot_vtu(p_mesh, c_filename, &
INTEGER( KIND = GRID_SI) :: i_numberofvtupoints
INTEGER( KIND = GRID_SI) :: i_numberofvtufaces
INTEGER( KIND = GRID_SI) :: i_numberofvtucells
REAL ( KIND = GRID_SR ) :: r_heightfactor
REAL ( KIND = GRID_SR ) :: r_heightfactor, r_aux
TYPE(t_vtu_data), DIMENSION(6) :: p_stcnodedata
REAL (KIND = GRID_SR), DIMENSION(:,:), POINTER :: r_stcval
INTEGER( KIND = GRID_SI) :: i_stcbeg
......@@ -394,7 +394,7 @@ SUBROUTINE plot_vtu(p_mesh, c_filename, &
END IF
DO i_cnt = 1, i_nodedata
IF (p_nodedata(i_cnt)%c_name(1:5) == 'layer') THEN
READ(i_hcnt,*) p_nodedata(i_cnt)%c_name(6:32)
READ(i_hcnt,'(I8)') p_nodedata(i_cnt)%c_name(6:32)
i_stcbeg = p_mesh%i_nnumber*(i_hcnt-1)+1
i_stcend = p_mesh%i_nnumber*i_hcnt
r_stcval(1:1,i_stcbeg:i_stcend) = p_nodedata(i_cnt)%p_vdata(1:1,:)
......@@ -443,7 +443,8 @@ SUBROUTINE plot_vtu(p_mesh, c_filename, &
DO i_hcnt=1,i_numberofcells
i_stcbeg = (i_hcnt-1)*(i_lay-1)+1
i_stcend = i_hcnt*(i_lay-1)
r_stcval(1:1,i_stcbeg:i_stcend) = p_celldata(i_cnt)%p_vdata(1:1,i_hcnt)
r_aux = p_celldata(i_cnt)%p_vdata(1,i_hcnt)
r_stcval(1,i_stcbeg:i_stcend) = r_aux
END DO
p_stcnodedata(1)%c_name = p_celldata(i_cnt)%c_name
p_stcnodedata(1)%i_size = p_celldata(i_cnt)%i_size
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment