Film cooling is an important technique to ensure safe operation and performance fulfillment of turbines. Its ultimate goal is to protect the axial turbine blades from high gas temperatures. An appropriate study is necessary in order to obtain a reliable representation of the flow characteristics involved in such phenomena. Because of the high computational cost of high-fidelity simulations, the low-fidelity simulation method Reynolds Averaged Navier Stokes (RANS) is commonly used in practical configurations. However, the majority of the current turbulent heat flux models fail to accurately predict heat transfer in film cooling flows. Recent work suggests the use of machine learning models to improve turbulent closure in these flows. In the present work, a machine learning model for spatially varying turbulent Prandtl number previously described in the literature is applied to a transverse film cooling flow consisting of a jet square channel. The results obtained in the present work were compared to adiabatic effectiveness experimental data available in the literature to assess the performance of the machine learning model. The results shown that for low blowing ratios (BR = 0.2 and BR = 0.4) the proposed machine learning model has poor performance. However, for the case with the highest blowing ratio (BR = 0.8), the proposed model presented better results. These results are then explained in terms of the resulting turbulent Prandtl number field and suggest that the training set is not appropriate for capturing the turbulent heat flux in fully attached jets in crossflow.