Commit 6f4aaa00 authored by Paul Stryck's avatar Paul Stryck
Browse files

set temperature on inlet

parent ecda90a5
......@@ -83,19 +83,20 @@ CONTAINS
DO iface=nIntFaces, nIntFaces+nBndFaces
IF (family_f(iface) == topology_get_ifamily(trim('YMIN1'))) THEN
IF (topo_family(family_f(iface))%velocity_userCode) THEN
!IF (topo_family(family_f(iface))%velocity_userCode) THEN
! Define Inlet Velocity
! ...
END IF
!END IF
IF (topo_family(family_f(iface))%temperature_userCode) THEN
! Define Inlet Temperature
! ...
x = cent_c_3(cell2_f(iface), 1)
t_c(1 + TANH(slopeFact * (1 + 2*x)))
END IF
END IF
END DO
IF ( controls_initial%Velocity_userCode )THEN
DO icell=1,nIntCells
DO icell=1,nIntCells+nBndCells
! Define (initial) field velocities
x = cent_c_3(icell, 1)
y = cent_c_3(icell, 2)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment