friction_factor_blasius_der Function

public function friction_factor_blasius_der(c_pt, Re) result(friction_der) bind(C)

Arguments

Type IntentOptional Attributes Name
type(c_ptr), value :: c_pt
real(kind=dp), value :: Re

Return Value real(kind=dp)


Source Code

  function friction_factor_blasius_der(c_pt,Re) result(friction_der) bind(C)
      type(c_ptr), value :: c_pt
      real(dp), value :: Re
      real(dp) :: friction_der
      real(dp) :: Remin
      type(blasius_cfg_t), pointer :: blasius_ptr
      call c_f_pointer(c_pt,blasius_ptr)      
  
      Remin=10.0_dp  
      if(Re>Remin) then
          friction_der=-blasius_ptr%alpha*blasius_ptr%beta*(Re**(-blasius_ptr%beta-1.0_dp))
      else
          friction_der=0.0_dp
      endif      
  end function friction_factor_blasius_der