From 869e405d5a41c5c5d04c35c75bc69963176779f6 Mon Sep 17 00:00:00 2001 From: "E. Branlard" <1318316+ebranlard@users.noreply.github.com> Date: Wed, 23 Aug 2023 14:09:54 -0600 Subject: [PATCH] F/for 1733 AeroMap (#108) * Add example AeroMap case * Change extension of aero map input file was ".inp"; now ".dvr" --------- Co-authored-by: Bonnie Jonkman --- .../5MW_Land_AeroMap/5MW_Land_AeroMap.drv | 52 +++++++ .../5MW_Land_AeroMap/5MW_Land_AeroMap.outb | Bin 0 -> 9001 bytes .../5MW_Land_AeroMap/5MW_Land_DLL_WTurb.fst | 73 ++++++++++ .../NRELOffshrBsline5MW_Onshore_AeroDyn15.dat | 107 ++++++++++++++ .../NRELOffshrBsline5MW_Onshore_ElastoDyn.dat | 133 ++++++++++++++++++ ...ffshrBsline5MW_Onshore_ElastoDyn_Tower.dat | 54 +++++++ .../5MW_Land_AeroMap/plotFASTAeroMap.m | 47 +++++++ 7 files changed, 466 insertions(+) create mode 100644 glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_AeroMap.drv create mode 100644 glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_AeroMap.outb create mode 100644 glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_DLL_WTurb.fst create mode 100644 glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_AeroDyn15.dat create mode 100644 glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_ElastoDyn.dat create mode 100644 glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_ElastoDyn_Tower.dat create mode 100644 glue-codes/openfast/5MW_Land_AeroMap/plotFASTAeroMap.m diff --git a/glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_AeroMap.drv b/glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_AeroMap.drv new file mode 100644 index 000000000..43967453d --- /dev/null +++ b/glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_AeroMap.drv @@ -0,0 +1,52 @@ +------- OpenFAST AeroMap INPUT FILE ---------------------------------------------- +AeroMap generation for FAST Certification Test #18: NREL 5.0 MW Baseline Wind Turbine (Onshore) +---------------------- OpenFAST MODEL FILE --------------------------------------- +"5MW_Land_DLL_WTurb.fst" FstFile - Name of the primary OpenFAST input file (-) +---------------------- STEADY-STATE SIMULATION CONTROL -------------------------------------- +false Echo - Echo input data to .ech (flag) + 1e-4 Toler - Convergence tolerance for nonlinear solve residual equation [>0] (-) + 50 MaxIter - Maximum number of iteration steps for nonlinear solve [>0] (-) + 1 N_SSJac - Number of iteration steps to recalculate steady-state Jacobian (-) [1=every iteration step, 2=every other step] (Note: for large flexible blades, this almost always needs to be 1) + 1E+05 SSJacSclFact - Scaling factor used in steady-state Jacobians (-) [on order of blade mass in kg] +---------------------- STEADY-STATE CASES -------------------------------------- + 2 WindSpeedOrTSR - Choice of swept parameter (switch) { 1:wind speed; 2: TSR } + 36 NumCases - Number of cases to run +RotSpeed WndSpeedOrTSR Pitch +(rpm) (m/s or -) (deg) +8.000 3.000 0.000 +8.000 5.500 0.000 +8.000 8.000 0.000 +8.000 10.500 0.000 +8.000 13.000 0.000 +8.000 15.500 0.000 +8.000 3.000 5.000 +8.000 5.500 5.000 +8.000 8.000 5.000 +8.000 10.500 5.000 +8.000 13.000 5.000 +8.000 15.500 5.000 +8.000 3.000 10.000 +8.000 5.500 10.000 +8.000 8.000 10.000 +8.000 10.500 10.000 +8.000 13.000 10.000 +8.000 15.500 10.000 +8.000 3.000 15.000 +8.000 5.500 15.000 +8.000 8.000 15.000 +8.000 10.500 15.000 +8.000 13.000 15.000 +8.000 15.500 15.000 +8.000 3.000 20.000 +8.000 5.500 20.000 +8.000 8.000 20.000 +8.000 10.500 20.000 +8.000 13.000 20.000 +8.000 15.500 20.000 +8.000 3.000 25.000 +8.000 5.500 25.000 +8.000 8.000 25.000 +8.000 10.500 25.000 +8.000 13.000 25.000 +8.000 15.500 25.000 + diff --git a/glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_AeroMap.outb b/glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_AeroMap.outb new file mode 100644 index 0000000000000000000000000000000000000000..8113c2568d9dfad6fdeb221de3af20ffdf958b1d GIT binary patch literal 9001 zcmb`M33yG%`^QggjTDhY6)i?AjSv!9M3hWKYKw#<6oG99gEtCpc|Mz{*z1-%1c|5=8H_vnPnRn)$cfM!NymQXI zsjhe^iqZ-`M!6PlXO4@a46vBe%&As$R%Q+xZ?c%!C{w1%l591lv8+tyALC!3?&Ozdwpb^yG3LzCFza};HJ!1T#Dqu|pPORI%C(v^O)T1+ zVoA14=*ptTB2@h%$sb<9sj?9K$qh@_j8*8r(7;mxlzVuqUt=1oEK7MU&`h3z% zql~&kTSmtmTU*?K{{OG{;n<7^Q6APt$M_kBkF(mtvCe9V58}lOMmZ6=u6VG{F$9Cn~SKnY7>4DRV zkbLjqrARe@?0^Wg$2D$2+@-}M9ZExv#;o=&i8dW#x3@cv$HL*mMtHgrvtoWJ1Y26( zn~S#QIdU_#k3HLtxIXtJj$I;Uv-GWdDVc$jgfj;?><#@(Dlo1NliHP zLow`QS3Pk&77icb+OHZ??c{r)5U5@Eo`5!M@zHvuy5qcyM2@;xisSSkw~I(w|F`uN z=a0tzH2*)j^=bZgp@CI&wB|UzSHlKKD5OISol&30!%9=!`%Ug5dC`$wg+O`it5gmY zTW_pH^7XH_kw6h+IYr}MzIYkQyf3UECPvdb@Kv{J;{0+>r)E`j_;U15?n3t;t7t37 zKe=h!AEZOOv-2j!+w;rTLUgLLE*AY*o0dz7gP-@8U~37@kD(ubW;u-$)9)@oTi-sd z2F}ZOm)_?!?G{O23Zr#K~8`L>+p~kRT3}>(s?a zrF;88i99#=N77Lq9&pe3Kdze=vl^Qb#;q_~+-epAh z^FxsiI-KT<+R-b^QHSa|tQv=LvcR>Z5U3Q76=X-1evZ_yXVyWg8RBQ0SFN$21j*+Q z`y8qDb?t%kitWw4kh1R#${zv6y3;hapu*s(W5~%GCoWp_q+M+IK%M{=5a9;MxCL?cV ze~?c)A{Jfe&@N8wK8IAc1m8kB+!v?$qU37F9d(GkE$Cd(%a*Ps9TWc?C?Sq=^N%+< zK%{qjQ393XX4|n}{HGV?jci?WGTO@ICb|U3k|D($sGSP^6>ZSb5^X3>^F>@7{;}ge z6=z*2uk|w(*Dx+0a_0jMR1dzr70Lek#yHo@f2Z{j)zbHcV|&s&Ymj>2?7?^+lwQl5 z()|eOpJcl)b{$UhMYJoJ=%_=j=}+f@7X9O4)FG3`j^IFTc))U`tk$lL1oHh^`*B{S zx#v$vv8G}<@%6uhabBhE8+4xXFT=yC=x~}ZYGt4L6c2`S_~=7jD6X8~s=UVMEU+REjuXo-w#I~}PFPV0eWA7s*b ztN(oCSM*mFiH%79_1x!hUcSxih6KtrIjo8fr}?5SZ~K#@4)*1dGK?pRzS)I3nASW3 z<8kHMJ4od}jU$n4e(_!Olbbpk??v0UQ2y{SMa^m61J>KTo-Q?t<4#@GD{sSHh}e9`x~)AL1t>*j9Mt1k7V=ZkirA3a}m z-_IIFK)c?dfSxZ2x~u}^2X z*Fc-cy_DocH)_Ob zS7IW-mhY~n_jRrMxHz;ObU4iy9{5w4@jq_z`$p!RR-Y(gRsJ>>3ah|P(t$nq#4nGjqQhyvh#o5s zJL*vKZqj=2o?R<&owXj%SUFHR+Lrc*%u9Yjf~m#dd>SctSO1IFY}bYH~x`!YDpr*Ha!?u$6(-2~^?cF!QncKeL` zXim~mkJ4ok>Ck5UMfXMBGPkM@r}@HLZhX&ChZ?-`KE_j{Eu=$zwP!j9Tej?O=xErp zwgj7>Sl0$ArX)Uvln2NCM(Z-~DCv;rB1s49aEx?-3+)RSIg^9oScHS2HeZ;;(QZ6^ z5t1j(Hz8GdkoKkiW2AAedTyZaGver>E!fXWC$%L$`ce9H+(o{t;3fv!)^WiRK|6DcCL!du~KZcJEoFUfFLwQYpH;3-^`MXHF&U zub<=S9AWF-T~Mc-=W&Vlo6Vn%t(eB4U3}EE0bZZ>2!DRp7XM0iiK9OZ@kU#X&mtXm z8%Xc=JuMwC)#gfO4nEk{`%2%Q#nvJzd8`fn?^T~ z$f!%SZrtakBy3j~CoD(um2=uipkDapd=3}Q3uUt5FX~Q<#(w?O2nDZ0 ze}up2ZTPdL2`e1^VTd={4&Mvt{DLOBZuQy2wZg%`!ttXGEx~z2!BzVI1yt?Do7-y! zpJ?59bI}FIsh1K)OC0xRZYQL6tW9+eRKz{nZ+djCUJ^{+ea1|rUc7xso!%*F)3;gw##mxZ=!X3sp1yW?ynS{G-1C;ifZyWf7KRW-Qip4=nq5bJbTo- zY3|-07#9jXZiSG0VuIcST^Tl@-vvY~4YB->t|O6^5S_i??{%{5#w zp0xS1sM$H#uig)6@i>3+juD@))b~63!%+LWIn7gcBc>tF2Zh#+{kiOCoQKcdWc;2* zUG&HPJh8Vt*}836+mYjlOKnF=fVzHkBcz@eS&IY3#D8er)Pbu9U^|bUIti&yo=wja zHmqYF+St^I>$ec51i|tK!;rcn0|7zW&={D^zaVsi|LojtoN)%EL ze$RM6>}#x>c6z6;#P-NI<2~->uwb;cpWb?l1C>taeNh)v>eNPEEM?Fj4%C+9UF{B3 z+`7_vBW?b2a5+6Ulrg&+J#K37YrkCj{eQ7;YC(mcgm~Nv;Y7A literal 0 HcmV?d00001 diff --git a/glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_DLL_WTurb.fst b/glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_DLL_WTurb.fst new file mode 100644 index 000000000..95af4e6b5 --- /dev/null +++ b/glue-codes/openfast/5MW_Land_AeroMap/5MW_Land_DLL_WTurb.fst @@ -0,0 +1,73 @@ +------- OpenFAST EXAMPLE INPUT FILE ------------------------------------------- +FAST Certification Test #18: NREL 5.0 MW Baseline Wind Turbine (Onshore) +---------------------- SIMULATION CONTROL -------------------------------------- +False Echo - Echo input data to .ech (flag) +"FATAL" AbortLevel - Error level when simulation should abort (string) {"WARNING", "SEVERE", "FATAL"} + 60 TMax - Total run time (s) + 0.00625 DT - Recommended module time step (s) + 2 InterpOrder - Interpolation order for input/output time history (-) {1=linear, 2=quadratic} + 0 NumCrctn - Number of correction iterations (-) {0=explicit calculation, i.e., no corrections} + 99999 DT_UJac - Time between calls to get Jacobians (s) + 1000000 UJacSclFact - Scaling factor used in Jacobians (-) +---------------------- FEATURE SWITCHES AND FLAGS ------------------------------ + 1 CompElast - Compute structural dynamics (switch) {1=ElastoDyn; 2=ElastoDyn + BeamDyn for blades} + 1 CompInflow - Compute inflow wind velocities (switch) {0=still air; 1=InflowWind; 2=external from OpenFOAM} + 2 CompAero - Compute aerodynamic loads (switch) {0=None; 1=AeroDyn v14; 2=AeroDyn v15} + 1 CompServo - Compute control and electrical-drive dynamics (switch) {0=None; 1=ServoDyn} + 0 CompSeaSt - Compute sea state information (switch) {0=None; 1=SeaState} + 0 CompHydro - Compute hydrodynamic loads (switch) {0=None; 1=HydroDyn} + 0 CompSub - Compute sub-structural dynamics (switch) {0=None; 1=SubDyn; 2=External Platform MCKF} + 0 CompMooring - Compute mooring system (switch) {0=None; 1=MAP++; 2=FEAMooring; 3=MoorDyn; 4=OrcaFlex} + 0 CompIce - Compute ice loads (switch) {0=None; 1=IceFloe; 2=IceDyn} + 0 MHK - MHK turbine type (switch) {0=Not an MHK turbine; 1=Fixed MHK turbine; 2=Floating MHK turbine} +---------------------- ENVIRONMENTAL CONDITIONS -------------------------------- + 9.80665 Gravity - Gravitational acceleration (m/s^2) + 1.225 AirDens - Air density (kg/m^3) + 0 WtrDens - Water density (kg/m^3) + 1.464E-05 KinVisc - Kinematic viscosity of working fluid (m^2/s) + 335 SpdSound - Speed of sound in working fluid (m/s) + 103500 Patm - Atmospheric pressure (Pa) [used only for an MHK turbine cavitation check] + 1700 Pvap - Vapour pressure of working fluid (Pa) [used only for an MHK turbine cavitation check] + 0 WtrDpth - Water depth (m) + 0 MSL2SWL - Offset between still-water level and mean sea level (m) [positive upward] +---------------------- INPUT FILES --------------------------------------------- +"NRELOffshrBsline5MW_Onshore_ElastoDyn.dat" EDFile - Name of file containing ElastoDyn input parameters (quoted string) +"../5MW_Baseline/NRELOffshrBsline5MW_BeamDyn.dat" BDBldFile(1) - Name of file containing BeamDyn input parameters for blade 1 (quoted string) +"../5MW_Baseline/NRELOffshrBsline5MW_BeamDyn.dat" BDBldFile(2) - Name of file containing BeamDyn input parameters for blade 2 (quoted string) +"../5MW_Baseline/NRELOffshrBsline5MW_BeamDyn.dat" BDBldFile(3) - Name of file containing BeamDyn input parameters for blade 3 (quoted string) +"../5MW_Baseline/NRELOffshrBsline5MW_InflowWind_12mps.dat" InflowFile - Name of file containing inflow wind input parameters (quoted string) +"NRELOffshrBsline5MW_Onshore_AeroDyn15.dat" AeroFile - Name of file containing aerodynamic input parameters (quoted string) +"NRELOffshrBsline5MW_Onshore_ServoDyn.dat" ServoFile - Name of file containing control and electrical-drive input parameters (quoted string) +"unused" SeaStFile - Name of file containing sea state input parameters (quoted string) +"unused" HydroFile - Name of file containing hydrodynamic input parameters (quoted string) +"unused" SubFile - Name of file containing sub-structural input parameters (quoted string) +"unused" MooringFile - Name of file containing mooring system input parameters (quoted string) +"unused" IceFile - Name of file containing ice input parameters (quoted string) +---------------------- OUTPUT -------------------------------------------------- +True SumPrint - Print summary data to ".sum" (flag) + 5 SttsTime - Amount of time between screen status messages (s) + 99999 ChkptTime - Amount of time between creating checkpoint files for potential restart (s) +"default" DT_Out - Time step for tabular output (s) (or "default") + 0 TStart - Time to begin tabular output (s) + 4 OutFileFmt - Format for tabular (time-marching) output file (switch) {1: text file [.out], 2: binary file [.outb], 3: both 1 and 2, 4: uncompressed binary [.outb, 5: both 1 and 4} +True TabDelim - Use tab delimiters in text tabular output file? (flag) {uses spaces if false} +"ES10.3E2" OutFmt - Format used for text tabular output, excluding the time channel. Resulting field should be 10 characters. (quoted string) +---------------------- LINEARIZATION ------------------------------------------- +False Linearize - Linearization analysis (flag) +False CalcSteady - Calculate a steady-state periodic operating point before linearization? [unused if Linearize=False] (flag) + 3 TrimCase - Controller parameter to be trimmed {1:yaw; 2:torque; 3:pitch} [used only if CalcSteady=True] (-) + 0.001 TrimTol - Tolerance for the rotational speed convergence [used only if CalcSteady=True] (-) + 0.01 TrimGain - Proportional gain for the rotational speed error (>0) [used only if CalcSteady=True] (rad/(rad/s) for yaw or pitch; Nm/(rad/s) for torque) + 0 Twr_Kdmp - Damping factor for the tower [used only if CalcSteady=True] (N/(m/s)) + 0 Bld_Kdmp - Damping factor for the blades [used only if CalcSteady=True] (N/(m/s)) + 2 NLinTimes - Number of times to linearize (-) [>=1] [unused if Linearize=False] + 30, 60 LinTimes - List of times at which to linearize (s) [1 to NLinTimes] [used only when Linearize=True and CalcSteady=False] + 1 LinInputs - Inputs included in linearization (switch) {0=none; 1=standard; 2=all module inputs (debug)} [unused if Linearize=False] + 1 LinOutputs - Outputs included in linearization (switch) {0=none; 1=from OutList(s); 2=all module outputs (debug)} [unused if Linearize=False] +False LinOutJac - Include full Jacobians in linearization output (for debug) (flag) [unused if Linearize=False; used only if LinInputs=LinOutputs=2] +False LinOutMod - Write module-level linearization output files in addition to output for full system? (flag) [unused if Linearize=False] +---------------------- VISUALIZATION ------------------------------------------ + 0 WrVTK - VTK visualization data output: (switch) {0=none; 1=initialization data only; 2=animation; 3=mode shapes} + 1 VTK_type - Type of VTK visualization data: (switch) {1=surfaces; 2=basic meshes (lines/points); 3=all meshes (debug)} [unused if WrVTK=0] +true VTK_fields - Write mesh fields to VTK data files? (flag) {true/false} [unused if WrVTK=0] + 15 VTK_fps - Frame rate for VTK output (frames per second){will use closest integer multiple of DT} [used only if WrVTK=2 or WrVTK=3] diff --git a/glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_AeroDyn15.dat b/glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_AeroDyn15.dat new file mode 100644 index 000000000..582456ff2 --- /dev/null +++ b/glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_AeroDyn15.dat @@ -0,0 +1,107 @@ +------- AERODYN v15 for OpenFAST INPUT FILE ----------------------------------------------- +NREL 5.0 MW offshore baseline aerodynamic input properties. +====== General Options ============================================================================ +False Echo - Echo the input to ".AD.ech"? (flag) +"default" DTAero - Time interval for aerodynamic calculations {or "default"} (s) + 1 WakeMod - Type of wake/induction model (switch) {0=none, 1=BEMT, 2=DBEMT, 3=OLAF} [WakeMod cannot be 2 or 3 when linearizing] + 2 AFAeroMod - Type of blade airfoil aerodynamics model (switch) {1=steady model, 2=Beddoes-Leishman unsteady model} [AFAeroMod must be 1 when linearizing] + 1 TwrPotent - Type tower influence on wind based on potential flow around the tower (switch) {0=none, 1=baseline potential flow, 2=potential flow with Bak correction} + 0 TwrShadow - Calculate tower influence on wind based on downstream tower shadow? (switch) {0=none, 1=Powles model, 2=Eames model} +True TwrAero - Calculate tower aerodynamic loads? (flag) +False FrozenWake - Assume frozen wake during linearization? (flag) [used only when WakeMod=1 and when linearizing] +False CavitCheck - Perform cavitation check? (flag) [AFAeroMod must be 1 when CavitCheck=true] +False Buoyancy - Include buoyancy effects? (flag) +False CompAA - Flag to compute AeroAcoustics calculation [used only when WakeMod = 1 or 2] +"unused" AA_InputFile - AeroAcoustics input file [used only when CompAA=true] +====== Environmental Conditions =================================================================== +"default" AirDens - Air density (kg/m^3) +"default" KinVisc - Kinematic viscosity of working fluid (m^2/s) +"default" SpdSound - Speed of sound in working fluid (m/s) +"default" Patm - Atmospheric pressure (Pa) [used only when CavitCheck=True] +"default" Pvap - Vapour pressure of working fluid (Pa) [used only when CavitCheck=True] +====== Blade-Element/Momentum Theory Options ====================================================== [unused when WakeMod=0 or 3] + 2 SkewMod - Type of skewed-wake correction model (switch) {1=uncoupled, 2=Pitt/Peters, 3=coupled} [unused when WakeMod=0 or 3] +"default" SkewModFactor - Constant used in Pitt/Peters skewed wake model {or "default" is 15/32*pi} (-) [used only when SkewMod=2; unused when WakeMod=0 or 3] +True TipLoss - Use the Prandtl tip-loss model? (flag) [unused when WakeMod=0 or 3] +True HubLoss - Use the Prandtl hub-loss model? (flag) [unused when WakeMod=0 or 3] +True TanInd - Include tangential induction in BEMT calculations? (flag) [unused when WakeMod=0 or 3] +False AIDrag - Include the drag term in the axial-induction calculation? (flag) [unused when WakeMod=0 or 3] +False TIDrag - Include the drag term in the tangential-induction calculation? (flag) [unused when WakeMod=0,3 or TanInd=FALSE] +"Default" IndToler - Convergence tolerance for BEMT nonlinear solve residual equation {or "default"} (-) [unused when WakeMod=0 or 3] + 100 MaxIter - Maximum number of iteration steps (-) [unused when WakeMod=0] +====== Dynamic Blade-Element/Momentum Theory Options ============================================== [used only when WakeMod=2] + 2 DBEMT_Mod - Type of dynamic BEMT (DBEMT) model {1=constant tau1, 2=time-dependent tau1, 3=constant tau1 with continuous formulation} (-) [used only when WakeMod=2] + 4 tau1_const - Time constant for DBEMT (s) [used only when WakeMod=2 and DBEMT_Mod=1 or 3] +====== OLAF -- cOnvecting LAgrangian Filaments (Free Vortex Wake) Theory Options ================== [used only when WakeMod=3] +"unused" OLAFInputFileName - Input file for OLAF [used only when WakeMod=3] +====== Beddoes-Leishman Unsteady Airfoil Aerodynamics Options ===================================== [used only when AFAeroMod=2] + 3 UAMod - Unsteady Aero Model Switch (switch) {2=B-L Gonzalez, 3=B-L Minnema/Pierce, 4=B-L HGM 4-states, 5=B-L HGM+vortex 5 states, 6=Oye, 7=Boeing-Vertol} [used only when AFAeroMod=2] +True FLookup - Flag to indicate whether a lookup for f' will be calculated (TRUE) or whether best-fit exponential equations will be used (FALSE); if FALSE S1-S4 must be provided in airfoil input files (flag) [used only when AFAeroMod=2] + 0 UAStartRad - Starting radius for dynamic stall (fraction of rotor radius [0.0,1.0]) [used only when AFAeroMod=2; if line is missing UAStartRad=0] + 1 UAEndRad - Ending radius for dynamic stall (fraction of rotor radius [0.0,1.0]) [used only when AFAeroMod=2; if line is missing UAEndRad=1] +====== Airfoil Information ========================================================================= + 1 AFTabMod - Interpolation method for multiple airfoil tables {1=1D interpolation on AoA (first table only); 2=2D interpolation on AoA and Re; 3=2D interpolation on AoA and UserProp} (-) + 1 InCol_Alfa - The column in the airfoil tables that contains the angle of attack (-) + 2 InCol_Cl - The column in the airfoil tables that contains the lift coefficient (-) + 3 InCol_Cd - The column in the airfoil tables that contains the drag coefficient (-) + 4 InCol_Cm - The column in the airfoil tables that contains the pitching-moment coefficient; use zero if there is no Cm column (-) + 0 InCol_Cpmin - The column in the airfoil tables that contains the Cpmin coefficient; use zero if there is no Cpmin column (-) + 8 NumAFfiles - Number of airfoil files used (-) +"../5MW_Baseline/Airfoils/Cylinder1.dat" AFNames - Airfoil file names (NumAFfiles lines) (quoted strings) +"../5MW_Baseline/Airfoils/Cylinder2.dat" +"../5MW_Baseline/Airfoils/DU40_A17.dat" +"../5MW_Baseline/Airfoils/DU35_A17.dat" +"../5MW_Baseline/Airfoils/DU30_A17.dat" +"../5MW_Baseline/Airfoils/DU25_A17.dat" +"../5MW_Baseline/Airfoils/DU21_A17.dat" +"../5MW_Baseline/Airfoils/NACA64_A17.dat" +====== Rotor/Blade Properties ===================================================================== +True UseBlCm - Include aerodynamic pitching moment in calculations? (flag) +"../5MW_Baseline/NRELOffshrBsline5MW_AeroDyn_blade.dat" ADBlFile(1) - Name of file containing distributed aerodynamic properties for Blade #1 (-) +"../5MW_Baseline/NRELOffshrBsline5MW_AeroDyn_blade.dat" ADBlFile(2) - Name of file containing distributed aerodynamic properties for Blade #2 (-) [unused if NumBl < 2] +"../5MW_Baseline/NRELOffshrBsline5MW_AeroDyn_blade.dat" ADBlFile(3) - Name of file containing distributed aerodynamic properties for Blade #3 (-) [unused if NumBl < 3] +====== Hub Properties ============================================================================== [used only when Buoyancy=True] + 0 VolHub - Hub volume (m^3) + 0 HubCenBx - Hub center of buoyancy x direction offset (m) +====== Nacelle Properties ========================================================================== [used only when Buoyancy=True] + 0 VolNac - Nacelle volume (m^3) + 0, 0, 0 NacCenB - Position of nacelle center of buoyancy from yaw bearing in nacelle coordinates (m) +====== Tail Fin Aerodynamics ======================================================================= +False TFinAero - Calculate tail fin aerodynamics model (flag) +"unused" TFinFile - Input file for tail fin aerodynamics [used only when TFinAero=True] +====== Tower Influence and Aerodynamics ============================================================ [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or Buoyancy=True] + 12 NumTwrNds - Number of tower nodes used in the analysis (-) [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or Buoyancy=True] +TwrElev TwrDiam TwrCd TwrTI TwrCb ! TwrTI used only when TwrShadow=2; TwrCb used only when Buoyancy=True +(m) (m) (-) (-) (-) +0.0000000E+00 6.0000000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +8.5261000E+00 5.7870000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +1.7053000E+01 5.5740000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +2.5579000E+01 5.3610000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +3.4105000E+01 5.1480000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +4.2633000E+01 4.9350000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +5.1158000E+01 4.7220000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +5.9685000E+01 4.5090000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +6.8211000E+01 4.2960000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +7.6738000E+01 4.0830000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +8.5268000E+01 3.8700000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +8.7600000E+01 3.8700000E+00 1.0000000E+00 1.0000000E-01 0.0000000E+00 +====== Outputs ==================================================================================== +False SumPrint - Generate a summary file listing input options and interpolated properties to ".AD.sum"? (flag) + 0 NBlOuts - Number of blade node outputs [0 - 9] (-) + 1, 9, 19 BlOutNd - Blade nodes whose values will be output (-) + 0 NTwOuts - Number of tower node outputs [0 - 9] (-) + 1, 2, 6 TwOutNd - Tower nodes whose values will be output (-) + OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-) +RtAeroCp +RtAeroCt +RtAeroCq +RtAeroPwr +RtAeroFxh +RtAeroMxh +END of OutList section (the word "END" must appear in the first 3 columns of the last OutList line) +====== Outputs for all blade stations (same ending as above for B1N1.... =========================== [optional section] + 1 BldNd_BladesOut - Number of blades to output all node information at. Up to number of blades on turbine. (-) +"All" BldNd_BlOutNd - Future feature will allow selecting a portion of the nodes to output. Not implemented yet. (-) + OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-) +END (the word "END" must appear in the first 3 columns of this last OutList line in the optional nodal output section) +==================================================================================================== \ No newline at end of file diff --git a/glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_ElastoDyn.dat b/glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_ElastoDyn.dat new file mode 100644 index 000000000..4355b0e49 --- /dev/null +++ b/glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_ElastoDyn.dat @@ -0,0 +1,133 @@ +------- ELASTODYN for OpenFAST INPUT FILE ------------------------------------------- +NREL 5.0 MW Baseline Wind Turbine for Use in Offshore Analysis. Properties from Dutch Offshore Wind Energy Converter (DOWEC) 6MW Pre-Design (10046_009.pdf) and REpower 5M 5MW (5m_uk.pdf) +---------------------- SIMULATION CONTROL -------------------------------------- +False Echo - Echo input data to ".ech" (flag) + 3 Method - Integration method: {1: RK4, 2: AB4, or 3: ABM4} (-) +"DEFAULT" DT - Integration time step (s) +---------------------- DEGREES OF FREEDOM -------------------------------------- +True FlapDOF1 - First flapwise blade mode DOF (flag) +True FlapDOF2 - Second flapwise blade mode DOF (flag) +True EdgeDOF - First edgewise blade mode DOF (flag) +False TeetDOF - Rotor-teeter DOF (flag) [unused for 3 blades] +True DrTrDOF - Drivetrain rotational-flexibility DOF (flag) +True GenDOF - Generator DOF (flag) +True YawDOF - Yaw DOF (flag) +True TwFADOF1 - First fore-aft tower bending-mode DOF (flag) +True TwFADOF2 - Second fore-aft tower bending-mode DOF (flag) +True TwSSDOF1 - First side-to-side tower bending-mode DOF (flag) +True TwSSDOF2 - Second side-to-side tower bending-mode DOF (flag) +False PtfmSgDOF - Platform horizontal surge translation DOF (flag) +False PtfmSwDOF - Platform horizontal sway translation DOF (flag) +False PtfmHvDOF - Platform vertical heave translation DOF (flag) +False PtfmRDOF - Platform roll tilt rotation DOF (flag) +False PtfmPDOF - Platform pitch tilt rotation DOF (flag) +False PtfmYDOF - Platform yaw rotation DOF (flag) +---------------------- INITIAL CONDITIONS -------------------------------------- + 0 OoPDefl - Initial out-of-plane blade-tip displacement (meters) + 0 IPDefl - Initial in-plane blade-tip deflection (meters) + 0 BlPitch(1) - Blade 1 initial pitch (degrees) + 0 BlPitch(2) - Blade 2 initial pitch (degrees) + 0 BlPitch(3) - Blade 3 initial pitch (degrees) [unused for 2 blades] + 0 TeetDefl - Initial or fixed teeter angle (degrees) [unused for 3 blades] + 0 Azimuth - Initial azimuth angle for blade 1 (degrees) + 12.1 RotSpeed - Initial or fixed rotor speed (rpm) + 0 NacYaw - Initial or fixed nacelle-yaw angle (degrees) + 0 TTDspFA - Initial fore-aft tower-top displacement (meters) + 0 TTDspSS - Initial side-to-side tower-top displacement (meters) + 0 PtfmSurge - Initial or fixed horizontal surge translational displacement of platform (meters) + 0 PtfmSway - Initial or fixed horizontal sway translational displacement of platform (meters) + 0 PtfmHeave - Initial or fixed vertical heave translational displacement of platform (meters) + 0 PtfmRoll - Initial or fixed roll tilt rotational displacement of platform (degrees) + 0 PtfmPitch - Initial or fixed pitch tilt rotational displacement of platform (degrees) + 0 PtfmYaw - Initial or fixed yaw rotational displacement of platform (degrees) +---------------------- TURBINE CONFIGURATION ----------------------------------- + 3 NumBl - Number of blades (-) + 63 TipRad - The distance from the rotor apex to the blade tip (meters) + 1.5 HubRad - The distance from the rotor apex to the blade root (meters) + -2.5 PreCone(1) - Blade 1 cone angle (degrees) + -2.5 PreCone(2) - Blade 2 cone angle (degrees) + -2.5 PreCone(3) - Blade 3 cone angle (degrees) [unused for 2 blades] + 0 HubCM - Distance from rotor apex to hub mass [positive downwind] (meters) + 0 UndSling - Undersling length [distance from teeter pin to the rotor apex] (meters) [unused for 3 blades] + 0 Delta3 - Delta-3 angle for teetering rotors (degrees) [unused for 3 blades] + 0 AzimB1Up - Azimuth value to use for I/O when blade 1 points up (degrees) + -5.0191 OverHang - Distance from yaw axis to rotor apex [3 blades] or teeter pin [2 blades] (meters) + 1.912 ShftGagL - Distance from rotor apex [3 blades] or teeter pin [2 blades] to shaft strain gages [positive for upwind rotors] (meters) + -5 ShftTilt - Rotor shaft tilt angle (degrees) + 1.9 NacCMxn - Downwind distance from the tower-top to the nacelle CM (meters) + 0 NacCMyn - Lateral distance from the tower-top to the nacelle CM (meters) + 1.75 NacCMzn - Vertical distance from the tower-top to the nacelle CM (meters) + -3.09528 NcIMUxn - Downwind distance from the tower-top to the nacelle IMU (meters) + 0 NcIMUyn - Lateral distance from the tower-top to the nacelle IMU (meters) + 2.23336 NcIMUzn - Vertical distance from the tower-top to the nacelle IMU (meters) + 1.96256 Twr2Shft - Vertical distance from the tower-top to the rotor shaft (meters) + 87.6 TowerHt - Height of tower relative to ground level [onshore], MSL [offshore wind or floating MHK], or seabed [fixed MHK] (meters) + 0 TowerBsHt - Height of tower base relative to ground level [onshore], MSL [offshore wind or floating MHK], or seabed [fixed MHK] (meters) + 0 PtfmCMxt - Downwind distance from the ground level [onshore], MSL [offshore wind or floating MHK], or seabed [fixed MHK] to the platform CM (meters) + 0 PtfmCMyt - Lateral distance from the ground level [onshore], MSL [offshore wind or floating MHK], or seabed [fixed MHK] to the platform CM (meters) + 0 PtfmCMzt - Vertical distance from the ground level [onshore], MSL [offshore wind or floating MHK], or seabed [fixed MHK] to the platform CM (meters) + 0 PtfmRefzt - Vertical distance from the ground level [onshore], MSL [offshore wind or floating MHK], or seabed [fixed MHK] to the platform reference point (meters) +---------------------- MASS AND INERTIA ---------------------------------------- + 0 TipMass(1) - Tip-brake mass, blade 1 (kg) + 0 TipMass(2) - Tip-brake mass, blade 2 (kg) + 0 TipMass(3) - Tip-brake mass, blade 3 (kg) [unused for 2 blades] + 56780 HubMass - Hub mass (kg) + 115926 HubIner - Hub inertia about rotor axis [3 blades] or teeter axis [2 blades] (kg m^2) + 534.116 GenIner - Generator inertia about HSS (kg m^2) + 240000 NacMass - Nacelle mass (kg) + 2607890 NacYIner - Nacelle inertia about yaw axis (kg m^2) + 0 YawBrMass - Yaw bearing mass (kg) + 0 PtfmMass - Platform mass (kg) + 0 PtfmRIner - Platform inertia for roll tilt rotation about the platform CM (kg m^2) + 0 PtfmPIner - Platform inertia for pitch tilt rotation about the platform CM (kg m^2) + 0 PtfmYIner - Platform inertia for yaw rotation about the platform CM (kg m^2) +---------------------- BLADE --------------------------------------------------- + 17 BldNodes - Number of blade nodes (per blade) used for analysis (-) +"../5MW_Baseline/NRELOffshrBsline5MW_Blade.dat" BldFile(1) - Name of file containing properties for blade 1 (quoted string) +"../5MW_Baseline/NRELOffshrBsline5MW_Blade.dat" BldFile(2) - Name of file containing properties for blade 2 (quoted string) +"../5MW_Baseline/NRELOffshrBsline5MW_Blade.dat" BldFile(3) - Name of file containing properties for blade 3 (quoted string) [unused for 2 blades] +---------------------- ROTOR-TEETER -------------------------------------------- + 0 TeetMod - Rotor-teeter spring/damper model {0: none, 1: standard, 2: user-defined from routine UserTeet} (switch) [unused for 3 blades] + 0 TeetDmpP - Rotor-teeter damper position (degrees) [used only for 2 blades and when TeetMod=1] + 0 TeetDmp - Rotor-teeter damping constant (N-m/(rad/s)) [used only for 2 blades and when TeetMod=1] + 0 TeetCDmp - Rotor-teeter rate-independent Coulomb-damping moment (N-m) [used only for 2 blades and when TeetMod=1] + 0 TeetSStP - Rotor-teeter soft-stop position (degrees) [used only for 2 blades and when TeetMod=1] + 0 TeetHStP - Rotor-teeter hard-stop position (degrees) [used only for 2 blades and when TeetMod=1] + 0 TeetSSSp - Rotor-teeter soft-stop linear-spring constant (N-m/rad) [used only for 2 blades and when TeetMod=1] + 0 TeetHSSp - Rotor-teeter hard-stop linear-spring constant (N-m/rad) [used only for 2 blades and when TeetMod=1] +---------------------- DRIVETRAIN ---------------------------------------------- + 100 GBoxEff - Gearbox efficiency (%) + 97 GBRatio - Gearbox ratio (-) + 867637000 DTTorSpr - Drivetrain torsional spring (N-m/rad) + 6215000 DTTorDmp - Drivetrain torsional damper (N-m/(rad/s)) +---------------------- FURLING ------------------------------------------------- +False Furling - Read in additional model properties for furling turbine (flag) [must currently be FALSE) +"unused" FurlFile - Name of file containing furling properties (quoted string) [unused when Furling=False] +---------------------- TOWER --------------------------------------------------- + 20 TwrNodes - Number of tower nodes used for analysis (-) +"NRELOffshrBsline5MW_Onshore_ElastoDyn_Tower.dat" TwrFile - Name of file containing tower properties (quoted string) +---------------------- OUTPUT -------------------------------------------------- +False SumPrint - Print summary data to ".sum" (flag) + 1 OutFile - Switch to determine where output will be placed: {1: in module output file only; 2: in glue code output file only; 3: both} (currently unused) +True TabDelim - Use tab delimiters in text tabular output file? (flag) (currently unused) +"ES10.3E2" OutFmt - Format used for text tabular output (except time). Resulting field should be 10 characters. (quoted string) (currently unused) + 0 TStart - Time to begin tabular output (s) (currently unused) + 1 DecFact - Decimation factor for tabular output {1: output every time step} (-) (currently unused) + 0 NTwGages - Number of tower nodes that have strain gages for output [0 to 9] (-) + 10, 19, 28 TwrGagNd - List of tower nodes that have strain gages [1 to TwrNodes] (-) [unused if NTwGages=0] + 0 NBlGages - Number of blade nodes that have strain gages for output [0 to 9] (-) + 5, 9, 13 BldGagNd - List of blade nodes that have strain gages [1 to BldNodes] (-) [unused if NBlGages=0] + OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-) +"OoPDefl1" - Blade 1 out-of-plane and in-plane deflections and tip twist +"IPDefl1" - Blade 1 out-of-plane and in-plane deflections and tip twist +"RootMxb1" - In-plane bending, out-of-plane bending, and pitching moments at the root of blade 1 +"RootMyb1" - In-plane bending, out-of-plane bending, and pitching moments at the root of blade 1 +"RootMzb1" - In-plane bending, out-of-plane bending, and pitching moments at the root of blade 1 +END of OutList section (the word "END" must appear in the first 3 columns of the last OutList line) +====== Outputs for all blade stations (same ending as above for Spn1.... =========================== [optional section] + 1 BldNd_BladesOut - Number of blades to output all node information at. Up to number of blades on turbine. (-) +"All" BldNd_BlOutNd - Future feature will allow selecting a portion of the nodes to output. Not implemented yet. (-) + OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-) +END (the word "END" must appear in the first 3 columns of this last OutList line in the optional nodal output section) +---------------------------------------------------------------- + diff --git a/glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_ElastoDyn_Tower.dat b/glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_ElastoDyn_Tower.dat new file mode 100644 index 000000000..664185cbd --- /dev/null +++ b/glue-codes/openfast/5MW_Land_AeroMap/NRELOffshrBsline5MW_Onshore_ElastoDyn_Tower.dat @@ -0,0 +1,54 @@ +------- ELASTODYN V1.00.* TOWER INPUT FILE ------------------------------------- +NREL 5.0 MW offshore baseline tower input properties. +---------------------- TOWER PARAMETERS ---------------------------------------- + 11 NTwInpSt - Number of input stations to specify tower geometry + 1 TwrFADmp(1) - Tower 1st fore-aft mode structural damping ratio (%) + 1 TwrFADmp(2) - Tower 2nd fore-aft mode structural damping ratio (%) + 1 TwrSSDmp(1) - Tower 1st side-to-side mode structural damping ratio (%) + 1 TwrSSDmp(2) - Tower 2nd side-to-side mode structural damping ratio (%) +---------------------- TOWER ADJUSTMUNT FACTORS -------------------------------- + 1 FAStTunr(1) - Tower fore-aft modal stiffness tuner, 1st mode (-) + 1 FAStTunr(2) - Tower fore-aft modal stiffness tuner, 2nd mode (-) + 1 SSStTunr(1) - Tower side-to-side stiffness tuner, 1st mode (-) + 1 SSStTunr(2) - Tower side-to-side stiffness tuner, 2nd mode (-) + 1 AdjTwMa - Factor to adjust tower mass density (-) + 1 AdjFASt - Factor to adjust tower fore-aft stiffness (-) + 1 AdjSSSt - Factor to adjust tower side-to-side stiffness (-) +---------------------- DISTRIBUTED TOWER PROPERTIES ---------------------------- + HtFract TMassDen TwFAStif TwSSStif + (-) (kg/m) (Nm^2) (Nm^2) +0.0000000E+00 5.5908700E+03 6.1434300E+11 6.1434300E+11 +1.0000000E-01 5.2324300E+03 5.3482100E+11 5.3482100E+11 +2.0000000E-01 4.8857600E+03 4.6326700E+11 4.6326700E+11 +3.0000000E-01 4.5508700E+03 3.9913100E+11 3.9913100E+11 +4.0000000E-01 4.2277500E+03 3.4188300E+11 3.4188300E+11 +5.0000000E-01 3.9164100E+03 2.9101100E+11 2.9101100E+11 +6.0000000E-01 3.6168300E+03 2.4602700E+11 2.4602700E+11 +7.0000000E-01 3.3290300E+03 2.0645700E+11 2.0645700E+11 +8.0000000E-01 3.0530100E+03 1.7185100E+11 1.7185100E+11 +9.0000000E-01 2.7887500E+03 1.4177600E+11 1.4177600E+11 +1.0000000E+00 2.5362700E+03 1.1582000E+11 1.1582000E+11 +---------------------- TOWER FORE-AFT MODE SHAPES ------------------------------ + 0.7004 TwFAM1Sh(2) - Mode 1, coefficient of x^2 term + 2.1963 TwFAM1Sh(3) - , coefficient of x^3 term + -5.6202 TwFAM1Sh(4) - , coefficient of x^4 term + 6.2275 TwFAM1Sh(5) - , coefficient of x^5 term + -2.504 TwFAM1Sh(6) - , coefficient of x^6 term + -70.5319 TwFAM2Sh(2) - Mode 2, coefficient of x^2 term + -63.7623 TwFAM2Sh(3) - , coefficient of x^3 term + 289.737 TwFAM2Sh(4) - , coefficient of x^4 term + -176.513 TwFAM2Sh(5) - , coefficient of x^5 term + 22.0706 TwFAM2Sh(6) - , coefficient of x^6 term +---------------------- TOWER SIDE-TO-SIDE MODE SHAPES -------------------------- + 1.385 TwSSM1Sh(2) - Mode 1, coefficient of x^2 term + -1.7684 TwSSM1Sh(3) - , coefficient of x^3 term + 3.0871 TwSSM1Sh(4) - , coefficient of x^4 term + -2.2395 TwSSM1Sh(5) - , coefficient of x^5 term + 0.5357 TwSSM1Sh(6) - , coefficient of x^6 term + -121.21 TwSSM2Sh(2) - Mode 2, coefficient of x^2 term + 184.415 TwSSM2Sh(3) - , coefficient of x^3 term + -224.904 TwSSM2Sh(4) - , coefficient of x^4 term + 298.536 TwSSM2Sh(5) - , coefficient of x^5 term + -135.838 TwSSM2Sh(6) - , coefficient of x^6 term + + diff --git a/glue-codes/openfast/5MW_Land_AeroMap/plotFASTAeroMap.m b/glue-codes/openfast/5MW_Land_AeroMap/plotFASTAeroMap.m new file mode 100644 index 000000000..444bd6bfe --- /dev/null +++ b/glue-codes/openfast/5MW_Land_AeroMap/plotFASTAeroMap.m @@ -0,0 +1,47 @@ +function [Channels, ChannelNames, ChannelUnits] = plotFASTAeroMap(fileName) +% fileName = 'AeroMap.outb'; + + + [Channels, ChannelNames, ChannelUnits] = ReadFASTbinary(fileName); + %% + Y_indx = find(strcmp(ChannelNames,'Pitch'),1,"first"); + X_indx = find(strcmp(ChannelNames,'TSR'),1,"first"); + [X, Y, Z, indx] = getGriddedData(Channels, X_indx, Y_indx); +%% + figure; + tiledlayout(3,5) + for i=1:length(Z) + nexttile; + if strcmpi(ChannelNames(indx{i}),'RtAeroCp') + mn = max(-0.5,min(Channels(:,indx{i}))); + mx = max(Channels(:,indx{i})); + dx = min(mx,0.7)/20; + contourf(X,Y,Z{i},[mn 0:dx:0.7 mx]) + else + contourf(X,Y,Z{i},25); + end + xlabel([ChannelNames{X_indx} ' ' ChannelUnits{X_indx}]); + ylabel([ChannelNames{Y_indx} ' ' ChannelUnits{Y_indx}]) + title( [ChannelNames{indx{i}} ' ' ChannelUnits{indx{i}}]); + colorbar + end + +end + +function [x1, y1, Z, indx] = getGriddedData(Channels, X_indx, Y_indx) + + x = unique(Channels(:,X_indx), "sorted"); + y = unique(Channels(:,Y_indx), "sorted"); + %% + ic=0; + [x1,y1] = meshgrid(x, y); + for i=2:size(Channels,2) + if i~=X_indx && i~=Y_indx + ic = ic+1; + Z{ic}=griddata(Channels(:,X_indx), Channels(:,Y_indx), Channels(:,i), x1, y1,'nearest'); + indx{ic}=i; + end + end + + end +