function [ct_cm,Pmax_Pmin]=engine_list(num)
%returns the engine parameters for the selected engine number

switch num
 case 1
  %nothing defined
  ct=NaN*ones(1,5);
  cm=NaN*ones(1,5);
  Pmax=NaN;
  Pmin=NaN;
 case 2
  %NEXT - High Thrust EOL
  ct=[17.5451d-3,  11.7148d-3, 11.9373d-3, -1.9287d-3, 0.0871d-3];
  cm=[3.1007d-6, -2.4847d-6, 1.5187d-6, -0.2495d-6, 0.013d-6]; 
  Pmax=7.2517d0;
  Pmin=0.6165d0;
 case 3
  %NEXT - High Isp EOL
  ct=[-0.0895d-3,  52.323d-3, -14.021d-3, 2.9258d-3, -0.1889d-3];
  cm=[1.9438d-6, 0.1768d-6, -0.1898d-6, 0.0764d-6, -0.006d-6];
  Pmax=7.2517d0;
  Pmin=0.6165d0;
 case 4
  %BPT-4000
  ct=[49.466d-3,  37.751d-3, -14.282d-3, 6.2683d-3, -0.6574d-3];
  cm=[12.204d-6, -7.232d-6,2.712d-6, -0.2629d-6,0.0025d-6];
  Pmax=4.91d0;
  Pmin=1.15d0;
 case 5
  %HiVHAC Hall
  ct=[4.9665d-3,  42.908d-3, 8.6254d-3, -5.4533d-3, 0.7632d-3];
  cm=[1.2854d-6, 2.3871d-6, -0.5356d-6, 0.0385d-6, 0.0042d-6];
  Pmax=3.760d0;
  Pmin=0.335d0;
 case 6
  %NSTAR Engine-I
  ct=[8.6930d-3,  10.684d-3, 16.755d-3,-3.1460d-3, 0.d-3];
  cm=[1.29320d-6,-1.03390d-6, 1.12250d-6,-0.17760d-6, 0.d-6];
  Pmax=2.524d0;
  Pmin=0.650d0;
 case 7
  %NSTAR Engine-J
  ct=[-3.4318d-3,  37.365d-3, 0.d-3, 0.d-3, 0.d-3];
  cm=[0.74343d-6, 0.20951d-6, 0.25205d-6, 0.d-6, 0.d-6];
  Pmax=2.529d0;
  Pmin=0.649d0;
 case 8
  %NSTAR Engine-K
  ct=[17.780d-3, -20.834d-3, 53.105d-3,-19.841d-3, 2.68360d-3];
  cm=[2.23940d-6,-4.53150d-6, 5.41050d-6,-2.19890d-6, 0.323440d-6];
  Pmax=2.484d0;
  Pmin=0.557d0;
 case 9
  %NSTAR Engine-N
  ct=[23.085d-3, -33.902d-3, 64.478d-3,-23.529d-3, 3.03020d-3];
  cm=[2.48509d-6,-5.17626d-6, 5.93077d-6,-2.39744d-6, 0.356303d-6];
  Pmax=2.484d0;
  Pmin=0.557d0;
 case 10
  %NSTAR Engine-O
  ct=[26.652d-3, -50.827d-3, 88.404d-3,-35.727d-3, 5.09790d-3];
  cm=[2.55190d-6,-5.48960d-6, 6.39290d-6,-2.62270d-6, 0.391150d-6];
  Pmax=2.510d0;
  Pmin=0.524d0;
 case 11
  %Conventional TAL
  ct=[34.414d-3,  33.703d-3, 0.d-3, 0.d-3, 0.d-3];
  cm=[4.43720d-6, 0.62977d-6, 0.d-6, 0.d-6, 0.d-6];
  Pmax=6.390d0;
  Pmin=1.730d0;
 case 12
  %Direct Drive TAL
  ct=[-16.031d-3,  51.889d-3, 0.d-3, 0.d-3, 0.d-3];
  cm=[0.29695d-6, 1.88090d-6, 0.d-6, 0.d-6, 0.d-6];
  Pmax=4.290d0;
  Pmin=1.280d0;
 case 13
  %14 cm Xenon
  ct=[-1.5432d-3,  33.714d-3, 0.d-3, 0.d-3, 0.d-3];
  cm=[0.11003d-6, 0.85333d-6, 0.d-6, 0.d-6, 0.d-6];
  Pmax=1.206d0;
  Pmin=0.143d0;
 case 14
  %Advanced NSTAR
  ct=[0.2000d-3,  33.600d-3, 0.d-3, 0.d-3, 0.d-3];
  cm=[0.67650d-6, 0.72350d-6, 0.d-6, 0.d-6, 0.d-6];
  Pmax=4.400d0;
  Pmin=0.800d0;
 case 15
  %Hughes 13 cm XIPS
  ct=[-1.9905d-3,  34.869d-3, 0.d-3, 0.d-3, 0.d-3];
  cm=[0.22006d-6, 0.15156d-6, 0.61452d-6, 0.d-6, 0.d-6];
  Pmax=0.906d0;
  Pmin=0.270d0;
 case 16
  %NSTAR Engine-Q
  ct=[26.673d-3,-52.301d-3, 91.639d-3,-37.190d-3, 5.21200d-3];
  cm=[2.50570d-6,-5.35390d-6, 6.25090d-6,-2.53640d-6, 0.369830d-6];
  Pmax=2.567d0;
  Pmin=0.525d0;
 case 17
  %NSTAR [DS4]
  ct=[17.977d-3, -15.262d-3, 40.682d-3,-10.832d-3, 0.69466d-3];
  cm=[2.25370d-6,-4.67200d-6, 5.28080d-6,-2.04900d-6, 0.291610d-6];
  Pmax=2.490d0;
  Pmin=0.560d0;
 case 18
  %SPT-100 Hall
  ct=[-36.300d-3,  99.300d-3,-13.700d-3, 0.d-3, 0.d-3];
  cm=[-0.80000d-6, 6.00000d-6,-1.00000d-6, 0.d-6, 0.d-6];
  Pmax=2.300d0;
  Pmin=0.900d0;
 case 19
  %NSTAR ST4-IPS V4.0
  ct=[32.810d-3, -67.820d-3, 103.70d-3,-41.180d-3, 5.69200d-3];
  cm=[2.81200d-6,-6.13300d-6, 6.82880d-6,-2.70600d-6, 0.383500d-6];
  Pmax=2.567d0;
  Pmin=0.577d0;
 case 20
  %Hi Isp NSTAR V3.0
  ct=[19.973d-3, -12.577d-3, 33.159d-3,-9.5071d-3, 0.92663d-3];
  cm=[1.97290d-6,-3.09930d-6, 2.81730d-6,-0.84169d-6, 0.089495d-6];
  Pmax=3.445d0;
  Pmin=0.561d0;
 case 21
  %Hi Isp NSTAR V4.0
  ct=[19.594d-3, -17.998d-3, 44.562d-3,-14.171d-3, 1.50150d-3];
  cm=[1.75150d-6,-2.82040d-6, 2.90760d-6,-0.89542d-6, 0.095684d-6];
  Pmax=3.444d0;
  Pmin=0.561d0;
 case 22
  %Hi Isp NSTAR V1.1
  ct=[34.035d-3, -49.899d-3, 67.640d-3,-22.316d-3, 2.60460d-3];
  cm=[2.25790d-6,-3.75500d-6, 3.61600d-6,-1.20300d-6, 0.144680d-6];
  Pmax=3.168d0;
  Pmin=0.561d0;
 case 23
  %Hi Isp NSTAR V5.0
  ct=[17.175d-3,  6.1028d-3, 7.6391d-3,-1.0546d-3, 0.04500d-3];
  cm=[1.58190d-6,-1.22650d-6, 0.76473d-6,-0.14535d-6, 0.010120d-6];
  Pmax=5.390d0;
  Pmin=0.561d0;
 case 24
  %GRC 40 cm 4000 Isp
  ct=[-2.7480d-3, 33.7870d-3, 0.d-3, 0.d-3, 0.d-3];
  cm=[-0.01298d-6, 0.87544d-6, 0.d-6, 0.d-6, 0.d-6];
  Pmax=10.880d0;
  Pmin=2.350d0;
 case 25
  %GRC 40 cm 6500 Isp
  ct=[3.0487d-3, 20.6660d-3, 0.d-3, 0.d-3, 0.d-3];
  cm=[0.04250d-6, 0.32060d-6, 0.d-6, 0.d-6, 0.d-6];
  Pmax=29.130d0;
  Pmin=5.910d0;
 case 26
  %GRC Hi Isp 40 cm
  ct=[3.8237d-3, 43.5071d-3,-0.7603d-3, 0.d-3, 0.d-3];
  cm=[0.77840d-6, 1.25613d-6,-0.03955d-6, 0.d-6, 0.d-6];
  Pmax=10.013d0;
  Pmin=1.092d0;
 case 27
  %SPT-140 Hall
  ct=[-71.782d-3, 95.9890d-3,-4.9665d-3, 0.d-3, 0.d-3];
  cm=[-0.15895d-6, 3.74930d-6,-0.07550d-6, 0.d-6, 0.d-6];
  Pmax=5.000d0;
  Pmin=2.200d0;
 case 28
  %GRC 5 kW 3300 Isp
  ct=[6.3955d-3, 39.3810d-3, 0.d-3, 0.d-3, 0.d-3];
  cm=[0.91307d-6, 1.04115d-6, 0.d-6, 0.d-6, 0.d-6];
  Pmax=4.713d0;
  Pmin=1.092d0;
 case 29
  %Hi Isp NSTAR V5.1
  ct=[-2.1219d-3, 43.2765d-3,-5.8573d-3, 0.7410d-3,-0.03777d-3];
  cm=[0.61299d-6, 0.79780d-6,-0.10737d-6, 0.00392d-6, 0.000673d-6];
  Pmax=5.390d0;
  Pmin=0.561d0;
 case 30
  %Hi Isp NSTAR Max T/P
  ct=[14.318d-3,-21.8070d-3, 77.588d-3,-37.974d-3, 6.16280d-3];
  cm=[1.06940d-6,-1.40770d-6, 3.21620d-6,-1.62390d-6, 0.273350d-6];
  Pmax=2.523d0;
  Pmin=0.561d0;
 case 31
  %Hi Isp NSTAR Engine-Q
  ct=[34.967d-3,-75.4730d-3, 114.52d-3,-47.092d-3, 6.80200d-3];
  cm=[2.96250d-6,-6.79590d-6, 7.76930d-6,-3.22320d-6, 0.480590d-6];
  Pmax=2.523d0;
  Pmin=0.561d0;
 case 32
  %GRC - Unknown
  ct=[86.804d-3,-52.5270d-3, 28.639d-3,-4.1089d-3, 0.21726d-3];
  cm=[5.75140d-6,-4.34120d-6, 1.72110d-6,-0.25178d-6, 0.013626d-6];
  Pmax=7.252d0;
  Pmin=0.616d0;
 case 33
  %GRC - Unknown
  ct=[17.668d-3, 28.8940d-3,-.54734d-3, 0.d-3, 0.d-3];
  cm=[1.63790d-6, 0.43528d-6,-.012619d-6, 0.d-6, 0.d-6];
  Pmax=7.252d0;
  Pmin=0.616d0;
 case 34
 otherwise
  %nothing defined
  ct=NaN*ones(1,5);
  cm=NaN*ones(1,5);
  Pmax=NaN;
  Pmin=NaN;
end

ct_cm=[ct;cm];
Pmax_Pmin=[Pmax;Pmin];