C 6/22/07 Initial release C VMOLWPT CALL HUNT(TEMP,N,T,JL,JH) FRACJ = MIN(1.0,MAX(0.0,(T - TEMP(JL))/(TEMP(JH) - TEMP(JL)))) PX = P CALL HUNTP(PRES,M,PX,KL,KH) FRACK = MIN(1.0,MAX(0.0,SNGL((PX - PRES(KL))/(PRES(KH) - PRES(KL))))) WMOL = (1. - FRACJ)*(1. - FRACK)*VMW(JL,KL) + . (1. - FRACK)*FRACJ*VMW(JH,KL) + . FRACK*FRACJ*VMW(JH,KH) + . (1. - FRACJ)*FRACK*VMW(JL,KH)