CoolProp: R1233zd does not work in the saturation region close to the bubble point

Try the following minimal working example

#include "CoolProp.h"
#include <stdio.h>

using namespace CoolProp;

int main()
{
    printf("%f\n",PropsSI("T", "P", 136730, "H", 260000.00, "R1233zd(E)"));
    printf("%f\n",PropsSI("T", "P", 136730, "H", 263277.27, "R1233zd(E)"));
    printf("%f\n",PropsSI("T", "P", 136730, "H", 263277.35, "R1233zd(E)"));
    printf("%f\n",PropsSI("T", "P", 136730, "H", 268000.00, "R1233zd(E)"));
    printf("%f\n",PropsSI("T", "P", 136730, "H", 275000.00, "R1233zd(E)"));
    printf("%f\n",PropsSI("T", "P", 136730, "H", 293000.00, "R1233zd(E)"));
    printf("%f\n",PropsSI("T", "P", 136730, "H", 300000.00, "R1233zd(E)"));
    return 0;
}

On my machine ti produces the following results

$ g++ -o t -DNO_FMTLIB -Iinclude test.cpp -L. -lCoolProp -ldl
$ ./t
296.914485
299.609517
inf
inf
inf
299.609875
299.609875

The expected behavior is to see 299.609875 in place of inf

@jowr could you please have a look, I’m keeping @casella and @albertoleva in the loop.

About this issue

  • Original URL
  • State: closed
  • Created 2 years ago
  • Comments: 22 (11 by maintainers)

Most upvoted comments

Deleting the enthalpy ancillary fixes the issue and the expected values are returned.