function delta(Te) REAL(rp), INTENT(IN) :: Te REAL(rp) :: delta delta = VTe(Te)*cparams_ss%deltao end function delta