echo on;
% --- BIGD3HILBERT loads files and plots their 3-D Hilbert spectrum.
% Data loaded is 'cp71301'-'cp713100'.
%
% Non MATLAB Library routines used are: NSPAB, FSPECIAL.
%
% See marked lines in the code to change cases.
%
% N. E. Huang (NASA GSFC)	08.Feb. 2001 Initial
% --- bigd3hilbert.m --- Version 08.Feb.2001 ----- %

clear;					% Start Fresh

load cp71301; cp71301=cp71301';		% Load HHT Results and Rotate
load cp71302; cp71302=cp71302';
load cp71303; cp71303=cp71303';
load cp71304; cp71304=cp71304';
load cp71305; cp71305=cp71305';
load cp71306; cp71306=cp71306';
load cp71307; cp71307=cp71307';
load cp71308; cp71308=cp71308';
load cp71309; cp71309=cp71309';
load cp713010; cp713010=cp713010'; 
load cp713011; cp713011=cp713011';
load cp713012; cp713012=cp713012';
load cp713013; cp713013=cp713013';
load cp713014; cp713014=cp713014';
load cp713015; cp713015=cp713015';
load cp713016; cp713016=cp713016';
load cp713017; cp713017=cp713017';
load cp713018; cp713018=cp713018';
load cp713019; cp713019=cp713019';
load cp713020; cp713020=cp713020';
load cp713021; cp713021=cp713021';					% Load HHT Results and Rotate
load cp713022; cp713022=cp713022';
load cp713023; cp713023=cp713023';
load cp713024; cp713024=cp713024';
load cp713025; cp713025=cp713025';
load cp713026; cp713026=cp713026';
load cp713027; cp713027=cp713027';
load cp713028; cp713028=cp713028';
load cp713029; cp713029=cp713029';
load cp713030; cp713030=cp713030'; 
load cp713031; cp713031=cp713031';
load cp713032; cp713032=cp713032';
load cp713033; cp713033=cp713033';
load cp713034; cp713034=cp713034';
load cp713035; cp713035=cp713035';
load cp713036; cp713036=cp713036';
load cp713037; cp713037=cp713037';
load cp713038; cp713038=cp713038';
load cp713039; cp713039=cp713039';
load cp713040; cp713040=cp713040';
load cp713041; cp713041=cp713041';					% Load HHT Results and Rotate
load cp713042; cp713042=cp713042';
load cp713043; cp713043=cp713043';
load cp713044; cp713044=cp713044';
load cp713045; cp713045=cp713045';
load cp713046; cp713046=cp713046';
load cp713047; cp713047=cp713047';
load cp713048; cp713048=cp713048';
load cp713049; cp713049=cp713049';
load cp713050; cp713050=cp713050'; 
load cp713051; cp713051=cp713051';
load cp713052; cp713052=cp713052';
load cp713053; cp713053=cp713053';
load cp713054; cp713054=cp713054';
load cp713055; cp713055=cp713055';
load cp713056; cp713056=cp713056';
load cp713057; cp713057=cp713057';
load cp713058; cp713058=cp713058';
load cp713059; cp713059=cp713059';
load cp713060; cp713060=cp713060';
load cp713061; cp713061=cp713061';					% Load HHT Results and Rotate
load cp713062; cp713062=cp713062';
load cp713063; cp713063=cp713063';
load cp713064; cp713064=cp713064';
load cp713065; cp713065=cp713065';
load cp713066; cp713066=cp713066';
load cp713067; cp713067=cp713067';
load cp713068; cp713068=cp713068';
load cp713069; cp713069=cp713069';
load cp713070; cp713070=cp713070'; 
load cp713071; cp713071=cp713071';
load cp713072; cp713072=cp713072';
load cp713073; cp713073=cp713073';
load cp713074; cp713074=cp713074';
load cp713075; cp713075=cp713075';
load cp713076; cp713076=cp713076';
load cp713077; cp713077=cp713077';
load cp713078; cp713078=cp713078';
load cp713079; cp713079=cp713079';
load cp713080; cp713080=cp713080';
load cp713081; cp713081=cp713081';					% Load HHT Results and Rotate
load cp713082; cp713082=cp713082';
load cp713083; cp713083=cp713083';
load cp713084; cp713084=cp713084';
load cp713085; cp713085=cp713085';
load cp713086; cp713086=cp713086';
load cp713087; cp713087=cp713087';
load cp713088; cp713088=cp713088';
load cp713089; cp713089=cp713089';
load cp713090; cp713090=cp713090'; 
load cp713091; cp713091=cp713091';
load cp713092; cp713092=cp713092';
load cp713093; cp713093=cp713093';
load cp713094; cp713094=cp713094';
load cp713095; cp713095=cp713095';
load cp713096; cp713096=cp713096';
load cp713097; cp713097=cp713097';
load cp713098; cp713098=cp713098';
load cp713099; cp713099=cp713099';
load cp7130100; cp7130100=cp7130100';

[n1, x, k]=nspab(cp71301(:, 1:4),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n2, x, k]=nspab(cp71302(:, 1:5),128, 0,5,0, 26.54); 
[n3, x, k]=nspab(cp71303(:, 1:6),128, 0,5,0, 26.54); 
[n4, x, k]=nspab(cp71304(:, 1:5),128, 0,5,0, 26.54); 
[n5, x, k]=nspab(cp71305(:, 1:6),128, 0,5,0, 26.54); 
[n6, x, k]=nspab(cp71306(:, 1:5),128, 0,5,0, 26.54); 
[n7, x, k]=nspab(cp71307(:, 1:5),128, 0,5,0, 26.54); 
[n8, x, k]=nspab(cp71308(:, 1:6),128, 0,5,0, 26.54); 
[n9, x, k]=nspab(cp71309(:, 1:5),128, 0,5,0, 26.54); 
[n10, x, k]=nspab(cp713010(:, 1:5),128, 0,5,0, 26.54); 
[n11, x, k]=nspab(cp713011(:, 1:5),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n12, x, k]=nspab(cp713012(:, 1:5),128, 0,5,0, 26.54); 
[n13, x, k]=nspab(cp713013(:, 1:6),128, 0,5,0, 26.54); 
[n14, x, k]=nspab(cp713014(:, 1:6),128, 0,5,0, 26.54); 
[n15, x, k]=nspab(cp713015(:, 1:6),128, 0,5,0, 26.54); 
[n16, x, k]=nspab(cp713016(:, 1:6),128, 0,5,0, 26.54); 
[n17, x, k]=nspab(cp713017(:, 1:5),128, 0,5,0, 26.54); 
[n18, x, k]=nspab(cp713018(:, 1:5),128, 0,5,0, 26.54); 
[n19, x, k]=nspab(cp713019(:, 1:5),128, 0,5,0, 26.54); 
[n20, x, k]=nspab(cp713020(:, 1:5),128, 0,5,0, 26.54); 
[n21, x, k]=nspab(cp713021(:, 1:5),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n22, x, k]=nspab(cp713022(:, 1:5),128, 0,5,0, 26.54); 
[n23, x, k]=nspab(cp713023(:, 1:5),128, 0,5,0, 26.54); 
[n24, x, k]=nspab(cp713024(:, 1:5),128, 0,5,0, 26.54); 
[n25, x, k]=nspab(cp713025(:, 1:6),128, 0,5,0, 26.54); 
[n26, x, k]=nspab(cp713026(:, 1:6),128, 0,5,0, 26.54); 
[n27, x, k]=nspab(cp713027(:, 1:6),128, 0,5,0, 26.54); 
[n28, x, k]=nspab(cp713028(:, 1:6),128, 0,5,0, 26.54); 
[n29, x, k]=nspab(cp713029(:, 1:6),128, 0,5,0, 26.54); 
[n30, x, k]=nspab(cp713030(:, 1:5),128, 0,5,0, 26.54); 
[n31, x, k]=nspab(cp713031(:, 1:6),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n32, x, k]=nspab(cp713032(:, 1:6),128, 0,5,0, 26.54); 
[n33, x, k]=nspab(cp713033(:, 1:6),128, 0,5,0, 26.54); 
[n34, x, k]=nspab(cp713034(:, 1:5),128, 0,5,0, 26.54); 
[n35, x, k]=nspab(cp713035(:, 1:5),128, 0,5,0, 26.54); 
[n36, x, k]=nspab(cp713036(:, 1:6),128, 0,5,0, 26.54); 
[n37, x, k]=nspab(cp713037(:, 1:5),128, 0,5,0, 26.54); 
[n38, x, k]=nspab(cp713038(:, 1:5),128, 0,5,0, 26.54); 
[n39, x, k]=nspab(cp713039(:, 1:5),128, 0,5,0, 26.54); 
[n40, x, k]=nspab(cp713040(:, 1:6),128, 0,5,0, 26.54); 
[n41, x, k]=nspab(cp713041(:, 1:5),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n42, x, k]=nspab(cp713042(:, 1:6),128, 0,5,0, 26.54); 
[n43, x, k]=nspab(cp713043(:, 1:6),128, 0,5,0, 26.54); 
[n44, x, k]=nspab(cp713044(:, 1:6),128, 0,5,0, 26.54); 
[n45, x, k]=nspab(cp713045(:, 1:5),128, 0,5,0, 26.54); 
[n46, x, k]=nspab(cp713046(:, 1:6),128, 0,5,0, 26.54); 
[n47, x, k]=nspab(cp713047(:, 1:6),128, 0,5,0, 26.54); 
[n48, x, k]=nspab(cp713048(:, 1:5),128, 0,5,0, 26.54); 
[n49, x, k]=nspab(cp713049(:, 1:5),128, 0,5,0, 26.54); 
[n50, x, k]=nspab(cp713050(:, 1:5),128, 0,5,0, 26.54); 
[n51, x, k]=nspab(cp713051(:, 1:5),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n52, x, k]=nspab(cp713052(:, 1:5),128, 0,5,0, 26.54); 
[n53, x, k]=nspab(cp713053(:, 1:5),128, 0,5,0, 26.54); 
[n54, x, k]=nspab(cp713054(:, 1:6),128, 0,5,0, 26.54); 
[n55, x, k]=nspab(cp713055(:, 1:6),128, 0,5,0, 26.54); 
[n56, x, k]=nspab(cp713056(:, 1:6),128, 0,5,0, 26.54); 
[n57, x, k]=nspab(cp713057(:, 1:6),128, 0,5,0, 26.54); 
[n58, x, k]=nspab(cp713058(:, 1:6),128, 0,5,0, 26.54); 
[n59, x, k]=nspab(cp713059(:, 1:5),128, 0,5,0, 26.54); 
[n60, x, k]=nspab(cp713060(:, 1:6),128, 0,5,0, 26.54); 
[n61, x, k]=nspab(cp713061(:, 1:6),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n62, x, k]=nspab(cp713062(:, 1:6),128, 0,5,0, 26.54); 
[n63, x, k]=nspab(cp713063(:, 1:6),128, 0,5,0, 26.54); 
[n64, x, k]=nspab(cp713064(:, 1:5),128, 0,5,0, 26.54); 
[n65, x, k]=nspab(cp713065(:, 1:6),128, 0,5,0, 26.54); 
[n66, x, k]=nspab(cp713066(:, 1:6),128, 0,5,0, 26.54); 
[n67, x, k]=nspab(cp713067(:, 1:6),128, 0,5,0, 26.54); 
[n68, x, k]=nspab(cp713068(:, 1:6),128, 0,5,0, 26.54); 
[n69, x, k]=nspab(cp713069(:, 1:6),128, 0,5,0, 26.54); 
[n70, x, k]=nspab(cp713070(:, 1:6),128, 0,5,0, 26.54); 
[n71, x, k]=nspab(cp713071(:, 1:5),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n72, x, k]=nspab(cp713072(:, 1:6),128, 0,5,0, 26.54); 
[n73, x, k]=nspab(cp713073(:, 1:5),128, 0,5,0, 26.54); 
[n74, x, k]=nspab(cp713074(:, 1:5),128, 0,5,0, 26.54); 
[n75, x, k]=nspab(cp713075(:, 1:5),128, 0,5,0, 26.54); 
[n76, x, k]=nspab(cp713076(:, 1:5),128, 0,5,0, 26.54); 
[n77, x, k]=nspab(cp713077(:, 1:5),128, 0,5,0, 26.54); 
[n78, x, k]=nspab(cp713078(:, 1:5),128, 0,5,0, 26.54); 
[n79, x, k]=nspab(cp713079(:, 1:6),128, 0,5,0, 26.54); 
[n80, x, k]=nspab(cp713080(:, 1:6),128, 0,5,0, 26.54); 
[n81, x, k]=nspab(cp713081(:, 1:6),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n82, x, k]=nspab(cp713082(:, 1:6),128, 0,5,0, 26.54); 
[n83, x, k]=nspab(cp713083(:, 1:6),128, 0,5,0, 26.54); 
[n84, x, k]=nspab(cp713084(:, 1:5),128, 0,5,0, 26.54); 
[n85, x, k]=nspab(cp713085(:, 1:6),128, 0,5,0, 26.54); 
[n86, x, k]=nspab(cp713086(:, 1:6),128, 0,5,0, 26.54); 
[n87, x, k]=nspab(cp713087(:, 1:6),128, 0,5,0, 26.54); 
[n88, x, k]=nspab(cp713088(:, 1:5),128, 0,5,0, 26.54); 
[n89, x, k]=nspab(cp713089(:, 1:5),128, 0,5,0, 26.54); 
[n90, x, k]=nspab(cp713090(:, 1:6),128, 0,5,0, 26.54); 
[n91, x, k]=nspab(cp713091(:, 1:6),128, 0,5,0, 26.54);		% Do Hilbert Spectra
[n92, x, k]=nspab(cp713092(:, 1:6),128, 0,5,0, 26.54); 
[n93, x, k]=nspab(cp713093(:, 1:6),128, 0,5,0, 26.54); 
[n94, x, k]=nspab(cp713094(:, 1:5),128, 0,5,0, 26.54); 
[n95, x, k]=nspab(cp713095(:, 1:5),128, 0,5,0, 26.54); 
[n96, x, k]=nspab(cp713096(:, 1:6),128, 0,5,0, 26.54); 
[n97, x, k]=nspab(cp713097(:, 1:5),128, 0,5,0, 26.54); 
[n98, x, k]=nspab(cp713098(:, 1:6),128, 0,5,0, 26.54); 
[n99, x, k]=nspab(cp713099(:, 1:6),128, 0,5,0, 26.54); 
[n100, x, k]=nspab(cp7130100(:, 1:6),128, 0,5,0, 26.54); 

q=fspecial('ga', 7, 0.7);		  							% Make Filter

n1=filter2(q, n1);n1=filter2(q, n1);n1=filter2(q, n1);	% Apply Filter
n2=filter2(q, n2);n2=filter2(q, n2);n2=filter2(q, n2);
n3=filter2(q, n3);n3=filter2(q, n3);n3=filter2(q, n3);
n4=filter2(q, n4);n4=filter2(q, n4);n4=filter2(q, n4);
n5=filter2(q, n5);n5=filter2(q, n5);n5=filter2(q, n5);
n6=filter2(q, n6);n6=filter2(q, n6);n6=filter2(q, n6);
n7=filter2(q, n7);n7=filter2(q, n7);n7=filter2(q, n7);
n8=filter2(q, n8);n8=filter2(q, n8);n8=filter2(q, n8);
n9=filter2(q, n9);n9=filter2(q, n9);n9=filter2(q, n9);
n10=filter2(q, n10);n10=filter2(q, n10);n10=filter2(q, n10);
n11=filter2(q, n11);n11=filter2(q, n11);n11=filter2(q, n11);	% Apply Filter
n12=filter2(q, n12);n12=filter2(q, n12);n12=filter2(q, n12);
n13=filter2(q, n13);n13=filter2(q, n13);n13=filter2(q, n13);
n14=filter2(q, n14);n14=filter2(q, n14);n14=filter2(q, n14);
n15=filter2(q, n15);n15=filter2(q, n15);n15=filter2(q, n15);
n16=filter2(q, n16);n16=filter2(q, n16);n16=filter2(q, n16);
n17=filter2(q, n17);n17=filter2(q, n17);n17=filter2(q, n17);