Hello,
I have an object that uniformly extends up through the flow layers. It's represented by applying head losses at the nodes on the vertical Z(I3D).
It has a max height (1.5m). Head losses should not be applied above this height. However, when I have flows that exceed this height head losses aren't applied on any of the nodes above or below the max height.
My code is below. Does anyone know why this is happening?
The code is simple, and should say that head losses are applied for nodes only below the max height:
DO IPLAN = NPLAN-29,NPLAN !for layers at that node
I3D = I1+NPOIN2*(IPLAN-1) !get the 3D node IDs for nodes I1
!height application
IF (Z(I3D).LE.SEMX)THEN !if node height < max height (SEMX) apply h.loss
NORM =SQRT(UN3%R(I3D)**2+VN3%R(I3D)**2+WN3%R(I3D)**2)
S1U%R(I3D)=0.5D0*A*CDV*RF*NORM
S1V%R(I3D)=0.5D0*A*CDV*RF*NORM
S1W%R(I3D)=0.5D0*A*CDV*RF*NORM
END IF
END DO
END IF
END DO