[inner7549]

Hello everyone，I'm a beginner of flow-3D，and I'm trying to put a heat source which is able to catch the free surface and heating the fluid at the catched surface from the z-axis direction.
I added these codes in the heatfl.F90,
----------------------------------------------------------------------------------------
if (xi(i).lt.(1.0e-4+5e-3*t).and.xi(i).gt.(0.8e-4+5e-3*t).and. &
& yj(j).lt.2.6e-4.and.yj(j).gt.2.4e-4.and. &
& zk(k).lt.1.8e-4.and.zk(k).gt.1.0e-4) then
if (abs(fn(ijk)-0.5).lt.1.and.fn(ijkp).lt.0.6.and.fn(ijkm).gt.0) heatfl1=5e10
endif
----------------------------------------------------------------------------------------
and it did add a heat source to the surface, but it just continue finding another mesh cell to add heat source which did not exist since the laser just heat the fluid surface at the top level.
It seems that only ijk, instead of any other cell index can be passed to fn subroutine. So I just want to get a skeleton of flow-3D and make change to other subroutine that can influence the whole process to work as a counter for heating the surface. Or there are any other way to do things like whar I said? Thank you.

[SUTATO]

i note that free surface elevation(hyel) is the Z-coordinate of free surface. but i can't add the "hyel" into the subroutine.