function [datos,nev] = BDPropFis(prop,comp)
%Banco de datos de propiedades físicas
%Datos de banco de datos de ASPEN PLUS PURE 10

% Por filas se devuelven las propiedades especificadas en el vector prop. 
% Las columnas son para cada componente del vector comp.
% Si la propiedad es un vector de parámetros en lugar de un escalar, todos
% aparecen seguidos. En el vector nev se recogen las posiciones que ocupa cada propiedad

% Ejemplo de uso:
%[dat,nev]=BDPropfis([2,45,6,47], [4,5])

%dat =
%  1.0e+007 *
%         0         0
%    4.9507    7.6350
%    0.0000    0.0000
%         0         0
%         0         0
%         0         0
%    0.0000    0.0000
%    0.0001    0.0001
%         0         0
%    0.0000    0.0000
%   -0.0007   -0.0010
%         0         0
%         0         0
%    0.0000    0.0000
%    0.0000    0.0000
%    0.0000    0.0000
%    0.0000    0.0000
%    0.0001    0.0001

%nev =
%     1     7     1     9

%---------------------------------------------------------------------------------

% Componentes puros (escalares). Propiedades 1 a 43

%  H=H2         M=CH4       B=C6H6       T=C7H8      F=C12H10
ComPuE=[
   340         340          28.5         30.8        10.9505    %  API
   0           0            0            0           0          %  CHARGE
   0           0            0            0           0          %  CHI
   1363020     459454000    1337.33      45026.8     37626.4    %  DCPLS       J/KMOL-K
   0           -50490000    129600000    122200000   276300000  %  DGFORM      J/KMOL
   0           0            0            0           0          %  DGFVK       J/KMOL
   0           0            0            0           0          %  DGSFRM      J/KMOL
   -4200000    -89040000    0            0           0          %  DHAQFM      J/KMOL
   0           -74520000    82880000      50170000   178490000  %  DHFORM      J/KMOL
   0           0            0            0           0          %  DHFVK       J/KMOL
   0           0            0            0           0          %  DHSFRM      J/KMOL
   896545      8171280      30803600     33363100    48661000   %  DHVLB       J/KMOL
   1           1            1            1           1          %  DLWC
   1           1            1            1           1          %  DVBLNC
   0.572       1.152        2.4          2.968       4.264      %  GMUQQ
   0.572       1.152        2.4          2.968       4.264      %  GMUQQ1
   0.416612    1.12393      3.19051      3.92287     6.04351    %  GMUQR
   -241820000  -802620000   -3136000000  -3734000000 -6248000   %  HCOM        J/KMOL
   117100      941400       9866000      6636000     18580000   %  HFUS        J/KMOL
   0           0            0            1.13763E-25 0          %  MUP         (J*CUM)**.5
   2.01588     16.04276     78.11364     92.14052    154.2114   %  MW
   -0.215993   0.0115478    0.2103       0.264012    0.402873   %  OMEGA
   -0.215993   0.0115478    0.2103       0.264012    0.402873   %  OMGPR
   0           0.01         0.212        0.26        0.402873   %  OMRKSS
   1313000     4599000      4895000      4108000     3380000    %  PC          N/SQM
   1313000     4599000      4895000      4108000     3380000    %  PCPR        N/SQM
   2047000     4617000      4924000      4053000     3380000    %  PCRKSS      N/SQM
   0.3187      0.2882       0.2694       0.2647      0.2747     %  RKTZRA
   130680      267100.4     426482.4     562902.8    0          %  S025E       J/KMOL-K
   0.3         0.3          0.8844       0.8718      0.993328   %  SG
   20.39       111.66       353.24       383.78      528.15     %  TB          K
   33.19       190.564      562.05       591.75      773        %  TC          K
   33.19       190.564      562.05       591.75      773        %  TCPR        K
   43.6        190.63       562.6        593.95      773        %  TCRKSS      K
   13.95       90.694       278.68       178.18      342.2      %  TFP         K
   13.95       90.694       278.68       178.18      342.2      %  TPT         K
   298.15      298.15       298.15       298.15      298.15     %  TREFHS      K
   2.01588     16.04276     78.11364     92.14052    154.2114   %  TRUEMW
   0.0285681   0.0379694    0.0958294    0.118095    0.183313   %  VB          CUM/KMOL
   0.064147    0.0986       0.256        0.316       0.497      %  VC          CUM/KMOL
   0.064147    0.0986       0.256        0.316       0.497      %  VCRKT       CUM/KMOL
   0.0535578   0.0535578    0.0885091    0.105918    0.155639   %  VLSTD       CUM/KMOL
   0.305       0.286        0.268        0.264       0.261      %  ZC
];


% Calor específico del gas ideal. Propiedad 44

%CPIGDP
%K, J/KMOL-K

%H       M       B       T       F
CPIGDP=[
27617   33298   44767   58140   107590
9560    79933   230850  286300  421050
2466    2086.9  1479.2  1440.6  1904.1
3760    41602   168360  189800  417850
567.6   991.96  677.66  -650.43 828.81
250     50      200     200     200
1500    1500    1500    1500    1500
];


% Calor de vaporización. Propiedad 45  
%DHVLDP
%K, J/KMOL

%H       M       B       T       F

DHVLDP=[
1012700 10194000   45346000  49507000  76350000
0.698   0.26087    0.39053   0.37742   0.39182
-1.817  -0.14694   0         0         0
1.447   0.22154    0         0         0
0       0          0         0         0
13.95   90.69      278.68    178.18    342.2
33.19   190.56     562.05    591.75    773
];

% Densidad del líquido. Propiedad 46
%DNLDIP
%K, KMOL/M3

%H       M       B       T       F
DNLDIP=[
5.414   2.9214  1.0259  0.8792  0.52257
0.34893 0.28976 0.26666 0.27136 0.25833
33.19   190.56  562.05  591.75  773
0.2706  0.28881 0.28394 0.29241 0.27026
0       0       0       0       0
13.95   90.69   278.68  178.18  342.2
33.19   190.56  562.05  591.75  773
];

%Presión de vapor. Propiedad 47
%PLXANT
%K,Pa

%H       M       B       T       F
PLXANT=[
12.69   39.205  83.107  76.945  77.314
-94.896 -1324.4 -6486.2 -6729.8 -9910.4
0       0       0       0       0
0       0       0       0       0
1.1125  -3.4366 -9.2194 -8.179  -7.5079
0.00032915      0.000031019     0.0000069844    0.0000053017    2.2385E-18
2       2       2       2       6
13.95   90.69   278.68  178.18  342.2
33.19   190.56  562.05  591.75  773
];

% Tensión superficial. Propiedad 48
% SIGDIP
% K, N/m

%H               M              B               T               F
SIGDIP=[
0.005345        0.036557        0.071815        0.066779        0.072185
1.0646          1.1466          1.2362          1.2442          1.2267
0               0               0               0               0
0               0               0               0               0
0               0               0               0               0
13.95           90.69           278.68          178.18          342.2
33.19           190.56          562.05          591.75          773
];

n=0;
for k=1:length(prop)
   
      switch prop(k)
      % Si la propiedad depende de T   
      case 44 
         nd=7;
         
         for i=1:length(comp)
            datos(n+1:n+nd,i)=CPIGDP(:,comp(i));
         end

      case 45
         nd=7;
         for i=1:length(comp)
            datos(n+1:n+nd,i)=DHVLDP(:,comp(i));
         end

      case 46
         nd=7;
         for i=1:length(comp)
            datos(n+1:n+nd,i)=DNLDIP(:,comp(i));
         end

      case 47
         nd=9;
         for i=1:length(comp)
            datos(n+1:n+nd,i)=PLXANT(:,comp(i));
         end
      case 48
         nd=7;
         for i=1:length(comp)
            datos(n+1:n+nd,i)=SIGDIP(:,comp(i));
         end

      otherwise % Si la propiedad es un escalar
         nd=1;
         for i=1:length(comp)
            datos(n+1,i)=ComPuE(prop(k),comp(i));
         end

      end
   nev(k)=nd;   
   n=n+nd;
end


