I have a question concerning the follow of boundary condition write in the subroutine usclim.f90.
My simulation is a room with a heater.
The temperature of the heater depend on the temperature in a cell volume (which is calculate in the subroutine usproj.f90)
This is my code in the subroutine usclim.f90
Code: Select all
CALL GETFBR('x>=4.9 and x<=5.1 and y>=-2 and y<=0 and z>=0 and z<=1',NLELT,LSTELT)
!Température de consigne de la pièce
TempCon = 25.d0
If (Temppiece < TempCon) then
!Radiateur allume
TempRadia = 70.d0
Do ilelt = 1, nlelt
ifac = lstelt(ilelt)
ITYPFB(IFAC,IPHAS)=IPAROI
rcodcl(ifac,isca(1),1) = TempRadia
enddo
etatradia = 1.d0
else
! !Radiateur égal à la
température de la pièce
Do ilelt = 1, nlelt
TempRadia = Temppiece
ifac = lstelt(ilelt)
ITYPFB(IFAC,IPHAS)=IPAROI
rcodcl(ifac,isca(1),1) =
TempRadia
enddo
etatradia = 0.d0
endif
Best regards