Hello :) as the title already states I have issues with using the VFILM subroutine. 

I want to use that subroutine in an explicit simulation. For simplicity I tested it with just a square workpiece.

That workpiece has an initial temperature of 300°C and the top surface is cooled via Surface film condition with fictional values h=1000mW/mm²K and T_sink=20°C.

The problem is, when using the funktion in the GUI, I have a result of HFL~5.3e+4 with a temperature around 80°C.

When I am using the subroutine with the exact same values, my result is HFL~1.6e+3 and a temperature around 295°C with the surface being defined as Surface (Mesh) and HFL~6.9e+2 and a temperature at nealy 300°C with the surface defined as Surface (Geometry).

I appended the .inp, where I changed the "F" to "FNU" in
*Sfilm
Surf_Contact, FNU, 20., 0.

and the .for file (using Abaqus 2025)
      SUBROUTINE VFILM(nblock,ndim,kStep,kIncr,stepTime,totalTime,
    * jElemUid,jNodeUid,numExfv,temp,field,curCoords,areaNode,jltyp,
    * sname,h,sink)
     INCLUDE 'vaba_param.inc'
C
     INTEGER nblock, ndim, kStep, kIncr, numExfv
     INTEGER jElemUid(nblock), jNodeUid(nblock), jltyp(nblock)
     INTEGER km
     DOUBLE PRECISION stepTime(2), totalTime
     DOUBLE PRECISION temp(nblock), field(nblock,numExfv)
     DOUBLE PRECISION curCoords(nblock,ndim), areaNode(nblock)
     DOUBLE PRECISION h(nblock), sink(nblock)
     CHARACTER*80 sname
C
     DO 100 km = 1, nblock
        h(km)    = 1000.0D0
        sink(km) = 20.0D0
 100 CONTINUE
     RETURN
     END
 

If someone has an idea, I appreciate the help :)