diff --git a/vlaplex.f90 b/vlaplex.f90 index a2304b2..6bc3423 100644 --- a/vlaplex.f90 +++ b/vlaplex.f90 @@ -14,10 +14,7 @@ real(dp), intent(in):: T real(dp):: Z - Z = 22.5 * (Temp_ref * T / eV_to_K /100.0)**0.6 - ! Z = max(Z, 1.0_dp) - ! Z = min(Z, 22.0_dp) - !Z = 12.0_dp + Z = 22.5 * (Temp_ref * T / eV_to_K / 100.0)**0.6 end function T_to_Z