IAPWS_IF97

该程序为计算水和水蒸汽性质的程序，采用的是IAPWS_IF97公式

function out = IAPWS_IF97（funin1in2）
% IAPWS_IF97（funin1in2）
%   102 water functions of water properties based on the International Association on Properties of Water and Steam
%   Industrial Formulation 1997 （IAPWS-IF97） IAPWS-IF97-S01 IAPWS-IF97-S03rev IAPWS-IF97-S04 IAPWS-IF97-S05 Revised
%   Advisory Note No. 3 Thermodynamic Derivatives from IAPWS Formulations 2008 Release on the IAPWS Formulation 2008
%   for the Viscosity of Ordinary Water Substance 2008 Revised Release on the IAPWS Formulation 1985 for the Thermal
%   Conductivity of Ordinary Water Substance.
%
%   fun is the desired function that may take 1 input in1 or 2 inputs in1 and in2. in1 and in2 should be scalar
%   column vector or matrix and in1 and in2 should be the same size. If a row vector is entered it is transposed. If a
%   scalar is entered for one input and the other input is a vector or matrix then the scalar is repeated to form a
%   vector or matrix of the same size as the other input.
if nargin<2 out = NaN; return end
if ~any（strcmpi（fun{‘k_pT‘‘k_ph‘‘mu_pT‘‘mu_ph‘‘dmudh_ph‘‘dmudp_ph‘‘dhLdp_p‘‘dhVdp_p‘‘dvdp_ph‘‘dvdh_ph‘‘dTdp_ph‘‘cp_ph‘‘h_pT‘‘v_pT‘ ...
        ‘vL_p‘‘vV_p‘‘hL_p‘‘hV_p‘‘T_ph‘‘v_ph‘‘psat_T‘‘Tsat_p‘‘h1_pT‘‘h2_pT‘‘h3_rhoT‘‘v1_pT‘‘v2_pT‘‘cp1_pT‘ ...
        ‘cp2_pT‘‘cp3_rhoT‘‘cv3_rhoT‘‘alphav1_pT‘‘alphav2_pT‘‘alphap3_rhoT‘‘betap3_rhoT‘‘kappaT1_pT‘‘kappaT2_pT‘...
        ‘dgammadtau1_pT‘‘dgammadpi1_pT‘‘dgammadtautau1_pT‘‘dgammadpipi1_pT‘‘dgammadpitau1_pT‘‘dgammadtau2_pT‘ ...
        ‘dgammadpi2_pT‘‘dgammadtautau2_pT‘‘dgammadpipi2_pT‘‘dgammadpitau2_pT‘‘dphidtau3_rhoT‘‘dphiddelta3_rhoT‘ ...
        ‘dphidtautau3_rhoT‘‘dphiddeltatau3_rhoT‘‘dphiddeltadelta3_rhoT‘‘dTsatdpsat_p‘‘T1_ph‘‘T2a_ph‘‘T2b_ph‘‘T2c_ph‘ ...
        ‘T3a_ph‘‘T3b_ph‘‘v3a_ph‘‘v3b_ph‘‘h2bc_p‘‘h3ab_p‘‘TB23_p‘‘pB23_T‘‘p3sat_h‘‘v3a_pT‘‘v3b_pT‘‘v3c_pT‘ ...
        ‘v3d_pT‘‘v3e_pT‘‘v3f_pT‘‘v3g_pT‘‘v3h_pT‘‘v3i_pT‘‘v3j_pT‘‘v3k_pT‘‘v3l_pT‘‘v3m_pT‘‘v3n_pT‘‘v3o_pT‘ ...
        ‘v3p_pT‘‘v3q_pT‘‘v3r_pT‘‘v3s_pT‘‘v3t_pT‘‘v3u_pT‘‘v3v_pT‘‘v3w_pT‘‘v3x_pT‘‘v3y_pT‘‘v3z_pT‘‘T3ab_p‘ ...
        ‘T3cd_p‘‘T3ef_p‘‘T3gh_p‘‘T3ij_p‘‘T3jk_p‘‘T3mn_p‘‘T3op_p‘‘T3qu_p‘‘T3rx_p‘‘T3uv_p‘‘T3wx_p‘}））
    out = NaN;return
end
if nargin==2
    dim = size（in1）;
    if length（dim）>2out = NaN;returnend
    if dim（1）==1 && dim（2）>1in1 = in1‘;end
    out = feval（funin1）;
end
if nargin==3
    dim1 = size（in1）;dim2 = size（in2）;
    if length（dim1）>2 || length（dim2）>2out = NaN;returnend
    if any（dim1~=dim2）;
        if dim1==ones（12）in1 = in1*ones（dim2）;dim1=dim2;
        elseif dim2==ones（12）in2 = in2*ones（dim1）;
        elseif dim1==fliplr（dim2）in1 = in1‘;dim1=dim2;
        else out = NaN;return
        end
    end
    if dim1（1）==1 && dim1（2）>1in1 = in1‘;in2 = in2‘;end
    out = feval（funin1in2）;
end
end
%% stand alone functions
function k = k_pT（pT）
% 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      188538  2012-03-19 01:38  IAPWS_IF97.m
     文件        1334  2012-03-26 15:36  license.txt