First commit

Benjamin Renard
Showing 1172 changed files with 134559 additions and 0 deletions Show diff stats
.gitignore
.project
amda_plus/README
amda_plus/amdalib/a2000.pro
amda_plus/amdalib/a2000omni.pro
amda_plus/amdalib/average.pro
amda_plus/amdalib/averagelap.pro
amda_plus/amdalib/averagemip.pro
amda_plus/amdalib/ba2000_0_init.pro
amda_plus/amdalib/ba2000_1_init.pro
amda_plus/amdalib/ba2000_2_init.pro
amda_plus/amdalib/ba2000_init.pro
amda_plus/amdalib/bcain_mso_0_init.pro
amda_plus/amdalib/bcain_mso_1_init.pro
amda_plus/amdalib/bcain_mso_2_init.pro
amda_plus/amdalib/bcain_mso_init.pro
amda_plus/amdalib/bcain_tot_init.pro
amda_plus/amdalib/bmorsch_mso_0_init.pro
amda_plus/amdalib/bmorsch_mso_1_init.pro
amda_plus/amdalib/bmorsch_mso_2_init.pro
+++ a/.gitignore
@@ -0,0 +1,9 @@
+AMDA_IHM/data/
+AMDA_IHM/generic_data/EPNResources/*.xml
+AMDA_IHM/generic_data/LocalData/
+AMDA_IHM/generic_data/Param/*.xml
+AMDA_IHM/generic_data/RemoteData/
+AMDA_IHM/generic_data/Statistics/
+AMDA_IHM/generic_data/newKernelDDBase/
+AMDA_IHM/js/lib/ext
+AMDA_IHM/help/simu/simu*
+++ a/.project
@@ -0,0 +1,11 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<projectDescription>
+	<name>AMDA-NG4</name>
+	<comment></comment>
+	<projects>
+	</projects>
+	<buildSpec>
+	</buildSpec>
+	<natures>
+	</natures>
+</projectDescription>
+++ a/amda_plus/README
@@ -0,0 +1,12 @@
+These idl codes should be added to old AMDA installation in order AMDA-NG to work.
+
+Copy 
+
+     amdalib/* to /AMDALIB/lib
+     ddlib/*   to /DDLIB/lib
+     predefined/*  to  pro specific for amda-ng don't copy; keep as it is 
+     pro/* to   pro specific for amda-ng don't copy; keep as it is 
+     templates/* - pro specific for amda-ng don't copy; keep as it is
+
+Normally there should be no interference between old AMDA and AMDA-NG (TO CHECK)
+So AMDA installation could be used with old AMDA as well as with AMDA-NG
 \ No newline at end of file
+++ a/amda_plus/amdalib/a2000.pro
@@ -0,0 +1,64 @@
+pro a2000, Time, Orbit
+
+  StepDefaut = 3600.0
+  FillValue = -1.e31
+  Number = N_elements(Time)
+  Bgsm = fltarr(3)
+  xyz = fltarr(3)
+  
+
+  num31 = where(Orbit[0,*] ne FillValue, Index31)
+  if (Index31 le 0) then return
+
+  Time = Time[num31]
+  Orbit = Orbit[*,num31]
+
+ ; StepAve = (Time[Number-1] - Time[0])/20.0 > StepDefaut
+  StepAve = StepDefaut;
+  TimeIntDst = StepAve > 3600.;
+  TimeIntDstS = call_external(!DDLIB + 'DD_idl.so','GetTime', double(TimeIntDst), /S_Value)
+  StepAveS = call_external(!DDLIB + 'DD_idl.so','GetTime', double(StepAve), /S_Value)
+      
+      
+      bm = fltarr(3);
+      bdd = fltarr(3,7);
+      UT = lonarr(7);
+ 
+   p = [2., 1., 1., 1., 1., 1., 1.];
+   result = call_external(!PROJLIB + 'a2012fa_idl.so', 'pstatus', p);
+   for TimeCurr = Time[0], Time[Number-1], double(StepAve) do begin
+
+        numOrb = where((Time - TimeCurr LT StepAve) AND (Time - TimeCurr GE 0.0), IndexOrb)
+        if (IndexOrb LE 0) then continue
+        
+        params = sw_monitor_a99_(TimeCurr,StepAve,Orbit[*,numOrb[IndexOrb/2]])
+        STime = call_external(!DDLIB + 'DD_idl.so','GetTime', double(TimeCurr+StepAve/2.0), /S_Value)
+       
+         AE = ae_(STime, StepAveS, TimeDst, RetSize, Sampling, /CLOSE);
+         AL = RetSize GT 0 ? mean(AE[1,*]) : -100.0;
+  
+        Dst = dst_(STime, TimeIntDstS, TimeDst, RetSize, Sampling, /CLOSE);
+        Dst_av = RetSize GT 0 ? mean(Dst) : -10.;
+        ro = params[0];
+        v = params[1];
+        bimf =params[2:4];
+
+        for i = 0, IndexOrb - 1 do begin 
+          res = call_external(!DDLIB + 'DD_idl.so','GetGavTime', UT, double(TimeCurr+StepAve/2.0));
+         
+	  x = Orbit[*,numOrb[i]];      
+  
+          result = call_external(!PROJLIB + 'a2012fa_idl.so', 'model',  $
+                                                UT[0], UT[1], UT[2], float(UT[3]),$
+						ro, v, bimf, Dst_av, AL, $
+					        x, bm, bdd);
+          if bm[0] eq 0.0 then bm(*) = FillValue ;  
+          Orbit[*,numOrb[i]] = bm;        
+        endfor         
+  endfor
+
+ return
+end
+
+  
+  
 \ No newline at end of file
+++ a/amda_plus/amdalib/a2000omni.pro
@@ -0,0 +1,64 @@
+pro a2000omni, Time, Orbit
+
+  StepDefaut = 3600.0
+  FillValue = -1.e31
+  Number = N_elements(Time)
+  Bgsm = fltarr(3)
+  xyz = fltarr(3)
+  
+
+  num31 = where(Orbit[0,*] ne FillValue, Index31)
+  if (Index31 le 0) then return
+
+  Time = Time[num31]
+  Orbit = Orbit[*,num31]
+
+ ; StepAve = (Time[Number-1] - Time[0])/20.0 > StepDefaut
+  StepAve = StepDefaut;
+  TimeIntDst = StepAve > 3600.;
+  TimeIntDstS = call_external(!DDLIB + 'DD_idl.so','GetTime', double(TimeIntDst), /S_Value)
+  StepAveS = call_external(!DDLIB + 'DD_idl.so','GetTime', double(StepAve), /S_Value)
+      
+      
+      bm = fltarr(3);
+      bdd = fltarr(3,7);
+      UT = lonarr(7);
+ 
+   p = [2., 1., 1., 1., 1., 1., 1.];
+   result = call_external(!PROJLIB + 'a2012fa_idl.so', 'pstatus', p);
+   for TimeCurr = Time[0], Time[Number-1], double(StepAve) do begin
+
+        numOrb = where((Time - TimeCurr LT StepAve) AND (Time - TimeCurr GE 0.0), IndexOrb)
+        if (IndexOrb LE 0) then continue
+        
+        params = sw_monitor_a99_omni_(TimeCurr,StepAve,Orbit[*,numOrb[IndexOrb/2]])
+        STime = call_external(!DDLIB + 'DD_idl.so','GetTime', double(TimeCurr+StepAve/2.0), /S_Value)
+       
+         AE = ae_(STime, StepAveS, TimeDst, RetSize, Sampling, /CLOSE);
+         AL = RetSize GT 0 ? mean(AE[1,*]) : -100.0;
+  
+        Dst = dst_(STime, TimeIntDstS, TimeDst, RetSize, Sampling, /CLOSE);
+        Dst_av = RetSize GT 0 ? mean(Dst) : -10.;
+        ro = params[0];
+        v = params[1];
+        bimf =params[2:4];
+
+        for i = 0, IndexOrb - 1 do begin 
+          res = call_external(!DDLIB + 'DD_idl.so','GetGavTime', UT, double(TimeCurr+StepAve/2.0));
+         
+	  x = Orbit[*,numOrb[i]];      
+  
+          result = call_external(!PROJLIB + 'a2012fa_idl.so', 'model',  $
+                                                UT[0], UT[1], UT[2], float(UT[3]),$
+						ro, v, bimf, Dst_av, AL, $
+					        x, bm, bdd);
+          if bm[0] eq 0.0 then bm(*) = FillValue ;  
+          Orbit[*,numOrb[i]] = bm;        
+        endfor         
+  endfor
+
+ return
+end
+
+  
+  
 \ No newline at end of file
+++ a/amda_plus/amdalib/average.pro
@@ -0,0 +1,149 @@
+; $Id: average.pro,v 1.6 2012/07/06 13:24:01 budnik Exp $
+;====================================================================
+;--------------- Averaging for 1D & 2D arrays (second dim -> Time) 
+;====================================================================
+
+pro AVERAGE, AverStructure, T, Val 
+common GraphC, Graph, GraphN, Item
+
+ FillValue = Graph[GraphN].FillValue[Item]
+ Yinfo = SIZE(Val) 
+ if (Yinfo[0] GT 2) then RETURN
+ RetSize = Yinfo[0] eq 2 ? Yinfo[2] : Yinfo[1]
+
+ if (Yinfo[0] eq 2) then index =  Yinfo[1] eq 1 ? 0 : Yinfo[1]/2 ;  
+ 
+ specialCase = 0; 
+ 
+ if (Graph[GraphN].Step GT AverStructure.Sampling*2.0) then begin
+           if (finite(FillValue)) then begin
+
+              if (Yinfo[0] eq 1) then num31 = where(Val ne FillValue, Index31) $
+              else begin
+                 num31 = where(Val[index,0:RetSize-1] ne FillValue, Index31);
+ ; artificial truc for different spectra definitions....
+                if (index gt 1) then begin                 
+                  num31_1 = where(Val[1,0:RetSize-1] ne FillValue, Index31_1);
+                  if (Index31_1 ne Index31) then begin                       
+                        specialCase = 1;                                                                     
+                  endif else begin
+                        if (Index31 le 0) then begin
+                            numTest =  where(Val[*,0:RetSize-1] ne FillValue, IndexTest);
+                            if (IndexTest gt Yinfo[1]) then specialCase = 1;  
+                        endif
+                  endelse
+              endif
+              endelse 
+           endif else begin
+             if (Yinfo[0] eq 1) then num31 = where(finite(Val), Index31) $
+             else begin
+                num31 = where(finite(Val[index,0:RetSize-1]), Index31); 
+; artificial truc for different spectra definitions....
+                if (index gt 1) then begin                 
+                  num31_1 = where(finite(Val[1,0:RetSize-1]), Index31_1);
+                  if (Index31_1 ne Index31) then begin                       
+                        specialCase = 1;                                                                     
+                  endif else begin
+                        if (Index31 le 0) then begin
+                            numTest =  where(finite(Val[*,0:RetSize-1]), IndexTest);
+                            if (IndexTest gt Yinfo[1]) then specialCase = 1;  
+                        endif
+                  endelse
+              endif
+            endelse    
+        endelse
+   
+         if (specialCase) then Index31 = n_elements(T);
+ 
+         if (Index31 le 0) then begin
+               AverStructure.Nres = -1L
+             return
+         endif
+  
+       if (specialCase) then begin          
+            Val = REFORM(Val[*,*],Yinfo[1]*Index31)
+       endif else begin 
+            T = T[num31]
+            Val = Yinfo[0] eq 2 ? REFORM(Val[*,num31],Yinfo[1]*Index31) : Val[num31];
+       endelse
+   
+      MF =  AverStructure.Nres eq -1L ? Val : [*(AverStructure.LastVal),VAL]
+      Time = AverStructure.Nres eq -1L ? T : [*(AverStructure.LastTime),T]
+      N_Time = N_elements(Time)
+      NAve = long(((Time[N_Time-1] < Graph[GraphN].TotalTime) - Time[0]) / Graph[GraphN].Step + 1)> 1L
+     
+      if Yinfo[0] eq 2 then MF = REFORM(MF, Yinfo[1], N_Time) 
+
+       numStart = where(abs((*AverStructure.Time)-Time[0]) LE Graph[GraphN].Step/2.0, IndexStart)
+       if (IndexStart le 0) then begin
+            AverStructure.Nres = -1L
+            return
+       endif
+       NStart = numStart[0]
+             
+      for i = NStart, NStart+NAve-1L  do begin   
+         num = where(abs((*AverStructure.Time)[i]-Time) LE Graph[GraphN].Step/2.0, Index)
+         if (Index LE 0) AND (Graph[GraphN].Sampling[Item] LT Graph[GraphN].MaxSampling[Item]) then $
+              num = where(abs((*AverStructure.Time)[i]-Time) LE Graph[GraphN].MaxSampling[Item]/2.0, Index)
+         if (Index GT 0) then begin
+          if (Yinfo[0] eq 2) then begin
+             if (specialCase) then begin
+                 for k = 0, Yinfo[1]-1 do begin
+                    if (Index GT 1) then begin
+                        temp = 0.0;
+                        numK = 0;
+                        for kk = 0, Index-1 do begin
+                           if  (finite(FillValue)) then begin
+                                if (finite(MF[k,num[kk]])) then begin
+                                    temp += MF[k,num[kk]];
+                                    numK++;
+                                endif
+                          endif else begin
+                                if (MF[k,num[kk]] ne FillValue ) then begin
+                                    temp += MF[k,num[kk]];
+                                    numK++;
+                                endif
+                          endelse
+                         endfor                          
+                        (*AverStructure.Val)[i*Yinfo[1]+k]  = temp/ numK; 
+                    endif else (*AverStructure.Val)[i*Yinfo[1]+k] = MF[k,num[0]]
+                 endfor
+             endif else (*AverStructure.Val)[i*Yinfo[1]:i*Yinfo[1]+Yinfo[1]-1] = Index GT 1 ? total(MF[*,num],2)/Index : MF[*,num[0]] 
+          endif else (*AverStructure.Val)[i] = Index GT 1 ? total(MF[num])/Index : MF[num[0]]
+         endif  
+      endfor  
+      AverStructure.Nres = Index GT 0 ? (num[Index-1]+1) < (N_Time-1) : -1L  
+      
+      if (Graph[GraphN].CurrentSec LT (Graph[GraphN].Sections-1)) AND (AverStructure.Nres ne -1) then begin
+            *(AverStructure.LastTime) = Time[AverStructure.NRes:*]
+            *(AverStructure.LastVal) = Yinfo[0] eq 2 ? REFORM(MF[*,AverStructure.NRes:N_Time-1], Yinfo[1]*(N_Time-AverStructure.Nres)) :   MF[AverStructure.NRes:N_Time-1] 
+      endif
+      
+      if (Graph[GraphN].CurrentSec eq (Graph[GraphN].Sections-1)) then begin
+        if ptr_valid(AverStructure.LastTime) then ptr_free, AverStructure.LastTime
+        if ptr_valid(AverStructure.LastVal) then ptr_free, AverStructure.LastVal
+      endif
+      
+   endif else begin ;if (Step GT Sampling*2.0)
+
+     records = N_elements(*AverStructure.Time);
+ 
+     if Yinfo[0] eq 2 then Val = reform(Val, Yinfo[1]*RetSize);
+     if (records eq 0) then begin
+         (*AverStructure.Time) = T; : [*AverStructure.Time,T]
+         (*AverStructure.Val) =  Val;
+     endif else begin 
+         numBefore = where((*AverStructure.Time)[records-1] GT T, IndexBefore);
+         if (IndexBefore eq 0) then begin
+                (*AverStructure.Time) =  [*AverStructure.Time,T];
+                (*AverStructure.Val) =   [*AverStructure.Val,Val];
+         endif else begin
+                (*AverStructure.Time) =  [*AverStructure.Time,T[IndexBefore-1:RetSize-1]];
+                if (Yinfo[0] eq 1) then  (*AverStructure.Val) =  [*AverStructure.Val,Val[IndexBefore-1:RetSize-1]]
+                if (Yinfo[0] eq 2) then  (*AverStructure.Val) =  [*AverStructure.Val,Val[Yinfo[1]*(IndexBefore-1):Yinfo[1]*RetSize-1]]  
+         endelse
+     endelse
+    endelse  
+
+return
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/averagelap.pro
@@ -0,0 +1,88 @@
+; $Id:   $
+;====================================================================
+;---------- Special LAP averaging  2D array
+;====================================================================
+
+pro AVERAGELAP, AverStructure, T, Val 
+common GraphC, Graph, GraphN, Item
+ 
+ Yinfo = SIZE(Val)  
+ RetSize =  Yinfo[2];
+ Val = reform(Val, Yinfo[1]*RetSize);
+
+ if (Graph[GraphN].Step GT AverStructure.Sampling*2.0) then begin
+  
+      MF =  AverStructure.Nres eq -1L ? Val : [*(AverStructure.LastVal),VAL]
+      Time = AverStructure.Nres eq -1L ? T : [*(AverStructure.LastTime),T]
+      N_Time = N_elements(Time)
+      NAve = long(((Time[N_Time-1] < Graph[GraphN].TotalTime) - Time[0]) / Graph[GraphN].Step + 1)> 1L
+ 
+      MF = REFORM(MF, Yinfo[1], N_Time) 
+
+       numStart = where(abs((*AverStructure.Time)-Time[0]) LE Graph[GraphN].Step/2.0, IndexStart)
+       if (IndexStart le 0) then begin
+            AverStructure.Nres = -1L
+            return
+       endif
+
+      NStart = numStart[0]
+             
+      for i = NStart, NStart+NAve-1L  do begin 
+  
+         num = where(abs((*AverStructure.Time)[i]-Time) LE Graph[GraphN].Step/2.0, Index)
+         if (Index LE 0) AND (Graph[GraphN].Sampling[Item] LT Graph[GraphN].MaxSampling[Item]) then $
+              num = where(abs((*AverStructure.Time)[i]-Time) LE Graph[GraphN].MaxSampling[Item]/2.0, Index)
+
+         if (Index GT 0) then begin         
+                 for k = 0, Yinfo[1]-2 do begin
+                    if (Index GT 1) then begin
+                        temp = 0.0;
+                        numK = 0;
+                        for kk = 0, Index-1 do begin                         
+                                if (finite(MF[k,num[kk]])) then begin
+                                    temp += MF[k,num[kk]];
+                                    numK++;
+                                endif  ; if (finite                      
+                         endfor   ; for kk                        
+                        (*AverStructure.Val)[i*Yinfo[1]+k]  = temp/ numK; 
+                    endif else (*AverStructure.Val)[i*Yinfo[1]+k] = MF[k,num[0]] ; if (Index GT 1)
+               endfor ;   for k = 0,  
+                (*AverStructure.Val)[i*Yinfo[1]+Yinfo[1]-1] = MF[Yinfo[1]-1,num[0]]               
+         endif ; if (Index GT 0)
+  
+      endfor  ; for i = NStart
+
+      AverStructure.Nres = Index GT 0 ? (num[Index-1]+1) < (N_Time-1) : -1L  
+      
+      if (Graph[GraphN].CurrentSec LT (Graph[GraphN].Sections-1)) AND (AverStructure.Nres ne -1) then begin
+ 
+            *(AverStructure.LastTime) = Time[AverStructure.NRes:*]
+            *(AverStructure.LastVal) = REFORM(MF[*,AverStructure.NRes:N_Time-1], Yinfo[1]*(N_Time-AverStructure.Nres)) 
+      endif
+      
+      if (Graph[GraphN].CurrentSec eq (Graph[GraphN].Sections-1)) then begin
+        if ptr_valid(AverStructure.LastTime) then ptr_free, AverStructure.LastTime
+        if ptr_valid(AverStructure.LastVal) then ptr_free, AverStructure.LastVal
+      endif
+      
+   endif else begin ;if (Step GT Sampling*2.0)
+
+     records = N_elements(*AverStructure.Time);
+     
+     if (records eq 0) then begin
+         (*AverStructure.Time) = T; : [*AverStructure.Time,T]
+         (*AverStructure.Val) =  Val;
+     endif else begin 
+         numBefore = where((*AverStructure.Time)[records-1] GT T, IndexBefore);
+         if (IndexBefore eq 0) then begin
+                (*AverStructure.Time) =  [*AverStructure.Time,T];
+                (*AverStructure.Val) =   [*AverStructure.Val,Val];
+         endif else begin
+                (*AverStructure.Time) =  [*AverStructure.Time,T[IndexBefore-1:RetSize-1]];              
+                (*AverStructure.Val) =  [*AverStructure.Val,Val[Yinfo[1]*(IndexBefore-1):Yinfo[1]*RetSize-1]]  
+         endelse
+     endelse
+    endelse  
+
+return
+end
+++ a/amda_plus/amdalib/averagemip.pro
@@ -0,0 +1,90 @@
+; $Id:   $
+;====================================================================
+;---------- Special MIP averaging  2D array
+;====================================================================
+
+pro AVERAGEMIP, AverStructure, T, Val 
+common GraphC, Graph, GraphN, Item
+
+ Yinfo = SIZE(Val)  
+
+ if (Yinfo[0] eq 2) then begin
+      Val = reform(Val, Yinfo[1]*Yinfo[2]);
+      RetSize = Yinfo[2];
+endif else RetSize = Yinfo[1];
+ 
+ if (Graph[GraphN].Step GT AverStructure.Sampling*2.0) then begin 
+      MF =  AverStructure.Nres eq -1L ? Val : [*(AverStructure.LastVal),VAL]
+      Time = AverStructure.Nres eq -1L ? T : [*(AverStructure.LastTime),T]
+      N_Time = N_elements(Time)
+      NAve = long(((Time[N_Time-1] < Graph[GraphN].TotalTime) - Time[0]) / Graph[GraphN].Step + 1)> 1L   
+   
+      MF = REFORM(MF, Yinfo[1], N_Time) 
+
+       numStart = where(abs((*AverStructure.Time)-Time[0]) LE Graph[GraphN].Step/2.0, IndexStart)
+       if (IndexStart le 0) then begin
+            AverStructure.Nres = -1L
+            return
+       endif
+
+      NStart = numStart[0]
+             
+      for i = NStart, NStart+NAve-1L  do begin 
+  
+         num = where(abs((*AverStructure.Time)[i]-Time) LE Graph[GraphN].Step/2.0, Index)
+         if (Index LE 0) AND (Graph[GraphN].Sampling[Item] LT Graph[GraphN].MaxSampling[Item]) then $
+              num = where(abs((*AverStructure.Time)[i]-Time) LE Graph[GraphN].MaxSampling[Item]/2.0, Index)
+
+         if (Index GT 0) then begin         
+                 for k = 0, Yinfo[1]-1 do begin
+                    if (Index GT 1) then begin
+                        temp = 0.0;
+                        numK = 0;
+                        for kk = 0, Index-1 do begin                         
+                            ;    if (finite(MF[k,num[kk]])) then begin
+                                if (MF[k,num[kk]] GT  Graph[GraphN].FillValue[Item]) then begin
+                                    temp += MF[k,num[kk]];
+                                    numK++;
+                                endif                        
+                         endfor                   
+                        (*AverStructure.Val)[i*Yinfo[1]+k]  = numK GT 0 ? temp/ numK : Graph[GraphN].FillValue[Item] ; 
+                    endif else begin 
+                         (*AverStructure.Val)[i*Yinfo[1]+k] = MF[k,num[0]]
+                     endelse
+               endfor                      
+         endif
+  
+      endfor  
+      AverStructure.Nres = Index GT 0 ? (num[Index-1]+1) < (N_Time-1) : -1L  
+      
+      if (Graph[GraphN].CurrentSec LT (Graph[GraphN].Sections-1)) AND (AverStructure.Nres ne -1) then begin
+            *(AverStructure.LastTime) = Time[AverStructure.NRes:*]; 
+            *(AverStructure.LastVal) = REFORM(MF[*,AverStructure.NRes:N_Time-1], Yinfo[1]*(N_Time-AverStructure.Nres))
+  endif
+      
+      if (Graph[GraphN].CurrentSec eq (Graph[GraphN].Sections-1)) then begin
+        if ptr_valid(AverStructure.LastTime) then ptr_free, AverStructure.LastTime
+        if ptr_valid(AverStructure.LastVal) then ptr_free, AverStructure.LastVal
+      endif
+      
+   endif else begin ;if (Step GT Sampling*2.0)
+
+     records = N_elements(*AverStructure.Time);
+
+     if (records eq 0) then begin
+         (*AverStructure.Time) = T; : [*AverStructure.Time,T]
+         (*AverStructure.Val) =  Val;
+     endif else begin 
+         numBefore = where((*AverStructure.Time)[records-1] GT T, IndexBefore);
+         if (IndexBefore eq 0) then begin
+                (*AverStructure.Time) =  [*AverStructure.Time,T];
+                (*AverStructure.Val) =   [*AverStructure.Val,Val];
+         endif else begin
+                (*AverStructure.Time) =  [*AverStructure.Time,T[IndexBefore-1:RetSize-1]];              
+                (*AverStructure.Val) =  [*AverStructure.Val,Val[Yinfo[1]*(IndexBefore-1):Yinfo[1]*RetSize-1]]  
+         endelse
+     endelse
+    endelse  
+
+return
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/ba2000_0_init.pro
@@ -0,0 +1,41 @@
+;----------------ENVELOPE-------------------------
+; 
+pro ba2000_0_init, CLUSTER1 = Cluster1, CLUSTER2 = Cluster2, CLUSTER3 = Cluster3, CLUSTER4 = Cluster4, $
+                     THEMISA = ThemisA, THEMISB = ThemisB, THEMISC = ThemisC, $
+                     THEMISD = ThemisD, THEMISE = ThemisE, $
+                     GEOTAIL = Geotail, $
+                     OMNI = OMNI, AMDA = AMDA, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+ 
+  cmd = 'ba2000_init, /X'
+  if KEYWORD_SET(FRAME) then cmd += ',FRAME = ' + frame
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(CLUSTER1) then cmd += ',/CLUSTER1' 
+  if KEYWORD_SET(CLUSTER2) then cmd += ',/CLUSTER2'  
+  if KEYWORD_SET(CLUSTER3) then cmd += ',/CLUSTER3'
+  if KEYWORD_SET(CLUSTER4) then cmd += ',/CLUSTER4'  
+  
+  if KEYWORD_SET(THEMISA) then cmd += ',/THEMISA' 
+  if KEYWORD_SET(THEMISB) then cmd += ',/THEMISB'  
+  if KEYWORD_SET(THEMISC) then cmd += ',/THEMISC'
+  if KEYWORD_SET(THEMISD) then cmd += ',/THEMISD'  
+  if KEYWORD_SET(THEMISE) then cmd += ',/THEMISE'
+
+  if KEYWORD_SET(GEOTAIL) then cmd += ',/GEOTAIL'
+
+  if KEYWORD_SET(OMNI) then cmd += ',/OMNI' 
+  if KEYWORD_SET(AMDA) then cmd += ',/AMDA' 
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/ba2000_1_init.pro
@@ -0,0 +1,41 @@
+;----------------ENVELOPE-------------------------
+; 
+pro ba2000_1_init, CLUSTER1 = Cluster1, CLUSTER2 = Cluster2, CLUSTER3 = Cluster3, CLUSTER4 = Cluster4, $
+                     THEMISA = ThemisA, THEMISB = ThemisB, THEMISC = ThemisC, $
+                     THEMISD = ThemisD, THEMISE = ThemisE, $
+                     GEOTAIL = Geotail, $
+                     OMNI = OMNI, AMDA = AMDA, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+ 
+  cmd = 'ba2000_init, /Y'
+  if KEYWORD_SET(FRAME) then cmd += ',FRAME = ' + frame
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(CLUSTER1) then cmd += ',/CLUSTER1' 
+  if KEYWORD_SET(CLUSTER2) then cmd += ',/CLUSTER2'  
+  if KEYWORD_SET(CLUSTER3) then cmd += ',/CLUSTER3'
+  if KEYWORD_SET(CLUSTER4) then cmd += ',/CLUSTER4'  
+  
+  if KEYWORD_SET(THEMISA) then cmd += ',/THEMISA' 
+  if KEYWORD_SET(THEMISB) then cmd += ',/THEMISB'  
+  if KEYWORD_SET(THEMISC) then cmd += ',/THEMISC'
+  if KEYWORD_SET(THEMISD) then cmd += ',/THEMISD'  
+  if KEYWORD_SET(THEMISE) then cmd += ',/THEMISE'
+
+  if KEYWORD_SET(GEOTAIL) then cmd += ',/GEOTAIL'
+
+  if KEYWORD_SET(OMNI) then cmd += ',/OMNI' 
+  if KEYWORD_SET(AMDA) then cmd += ',/AMDA' 
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/ba2000_2_init.pro
@@ -0,0 +1,41 @@
+;----------------ENVELOPE-------------------------
+; 
+pro ba2000_2_init, CLUSTER1 = Cluster1, CLUSTER2 = Cluster2, CLUSTER3 = Cluster3, CLUSTER4 = Cluster4, $
+                     THEMISA = ThemisA, THEMISB = ThemisB, THEMISC = ThemisC, $
+                     THEMISD = ThemisD, THEMISE = ThemisE, $
+                     GEOTAIL = Geotail, $
+                     OMNI = OMNI, AMDA = AMDA, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+ 
+  cmd = 'ba2000_init, /Z'
+  if KEYWORD_SET(FRAME) then cmd += ',FRAME = ' + frame
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(CLUSTER1) then cmd += ',/CLUSTER1' 
+  if KEYWORD_SET(CLUSTER2) then cmd += ',/CLUSTER2'  
+  if KEYWORD_SET(CLUSTER3) then cmd += ',/CLUSTER3'
+  if KEYWORD_SET(CLUSTER4) then cmd += ',/CLUSTER4'  
+  
+  if KEYWORD_SET(THEMISA) then cmd += ',/THEMISA' 
+  if KEYWORD_SET(THEMISB) then cmd += ',/THEMISB'  
+  if KEYWORD_SET(THEMISC) then cmd += ',/THEMISC'
+  if KEYWORD_SET(THEMISD) then cmd += ',/THEMISD'  
+  if KEYWORD_SET(THEMISE) then cmd += ',/THEMISE'
+
+  if KEYWORD_SET(GEOTAIL) then cmd += ',/GEOTAIL'
+
+  if KEYWORD_SET(OMNI) then cmd += ',/OMNI' 
+  if KEYWORD_SET(AMDA) then cmd += ',/AMDA' 
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/ba2000_init.pro
@@ -0,0 +1,110 @@
+; a2000 field along the orbit
+;  
+
+pro ba2000_init, CLUSTER1 = Cluster1, CLUSTER2 = Cluster2, CLUSTER3 = Cluster3, CLUSTER4 = Cluster4, $
+                     THEMISA = ThemisA, THEMISB = ThemisB, THEMISC = ThemisC, $
+                     THEMISD = ThemisD, THEMISE = ThemisE, $
+                     GEOTAIL = Geotail, $
+                     OMNI = OMNI, AMDA = AMDA, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+common GraphC, graph, GraphN, Item
+ 
+
+ names = ['B!DX', 'B!DY', 'B!DZ', 'B']
+ 
+ pen = !dnc-1 & comp = 100 & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;
+     
+     title = 'CL1';
+     graph[GraphN].SetCommonF[0] = 'c1_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getc1_xyz'; 
+ 
+  if KEYWORD_SET(CLUSTER2) then begin
+     title = 'CL2';
+     graph[GraphN].SetCommonF[0] = 'c2_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getc2_xyz'; 
+   endif
+  if KEYWORD_SET(CLUSTER3) then begin     
+     title = 'CL3';
+     graph[GraphN].SetCommonF[0] = 'c3_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getc3_xyz'; 
+   endif
+  if KEYWORD_SET(CLUSTER4) then begin     
+     title = 'CL4'
+     graph[GraphN].SetCommonF[0] = 'c4_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getc4_xyz'; 
+   endif
+  if KEYWORD_SET(THEMISA) then begin
+    title = 'THA'
+    graph[GraphN].SetCommonF[0] = 'tha_xyz_com'; 
+    graph[GraphN].GetDataF[0] ='gettha_xyz';
+   endif
+ if KEYWORD_SET(THEMISB) then begin
+    title = 'THB'
+    graph[GraphN].SetCommonF[0] = 'thb_xyz_com'; 
+    graph[GraphN].GetDataF[0] ='getthb_xyz';
+   endif
+ if KEYWORD_SET(THEMISC) then begin
+    title = 'THC'
+    graph[GraphN].SetCommonF[0] = 'thc_xyz_com'; 
+    graph[GraphN].GetDataF[0] ='getthc_xyz';
+   endif
+ if KEYWORD_SET(THEMISD) then begin
+    title = 'THD'
+    graph[GraphN].SetCommonF[0] = 'thd_xyz_com'; 
+    graph[GraphN].GetDataF[0] = 'getthd_xyz';
+   endif
+ if KEYWORD_SET(THEMISE) then begin
+    title = 'THE'
+    graph[GraphN].SetCommonF[0] = 'the_xyz_com'; 
+    graph[GraphN].GetDataF[0] = 'getthe_xyz';
+   endif
+ if KEYWORD_SET(GEOTAIL) then begin
+    title = 'GTL'
+    graph[GraphN].SetCommonF[0] = 'gtl_xyz_com'; 
+    graph[GraphN].GetDataF[0] = 'getgtl_xyz';
+   endif  
+
+
+a2000Function = KEYWORD_SET(OMNI) ? 'A2000OMNI' : 'A2000O'; 
+
+frame = "GSM";
+graph[GraphN].ProcessF= 'gse2gsm;' + a2000Function; 
+  
+if (KEYWORD_SET(FRAME)) then begin 
+  if (frame eq "GSE") then graph[GraphN].ProcessF= 'gse2gsm;' + a2000Function + ';gsm2gse';
+  if (frame eq "SM")  then graph[GraphN].ProcessF= 'gse2gsm;' + a2000Function + ';gse2sm';
+endif 
+
+yTitle = title + ' A2000 B '+frame+', nT';
+
+if (KEYWORD_SET(AY1TITLE)) then graph[GraphN].Ly.title = AY1Title $
+else graph[GraphN].Ly.title = yTitle
+
+graph[GraphN].Type = 'T96_T'
+graph[GraphN].DataN = 1 
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 120.0D0) < 36000.0D0 ; 
+graph[GraphN].AxisF = 'lineaxis,' + string(Graph[GraphN].TotalTime);
+
+if (KEYWORD_SET(X)) then comp = 0;
+if (KEYWORD_SET(Y)) then comp = 1;
+if (KEYWORD_SET(Z)) then comp = 2;
+if (KEYWORD_SET(TOTAL)) then comp = 3;
+
+graph[GraphN].PlotF = comp eq 100 ? 'PlotVec': 'PlotVec_Comp_New, ' + string(comp) +','+ string(pen ) +','+ string(symbol)
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) = (comp le 3) ? names : names[0:2];
+ 
+ 
+graph[GraphN].LastPlotF = 'nodatalast'  
+  
+return
+end
+++ a/amda_plus/amdalib/bcain_mso_0_init.pro
@@ -0,0 +1,28 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bcain_mso_0_init, MEX = Mex, MGS = mgs, MAVEN = maven, SAMPLING = sampling,     $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bcain_mso_init, /X, SAMPLING = ' + string(sampling)
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(MEX) then cmd += ',/MEX' 
+    
+  if KEYWORD_SET(MGS) then cmd += ',/MGS' 
+ 
+  if KEYWORD_SET(MAVEN) then cmd += ',/MAVEN'
+  
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bcain_mso_1_init.pro
@@ -0,0 +1,28 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bcain_mso_1_init, MEX = Mex, MGS = mgs, MAVEN = maven, SAMPLING = sampling,     $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bcain_mso_init, /Y, SAMPLING = ' + string(sampling)
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(MEX) then cmd += ',/MEX' 
+    
+  if KEYWORD_SET(MGS) then cmd += ',/MGS' 
+ 
+  if KEYWORD_SET(MAVEN) then cmd += ',/MAVEN'
+  
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bcain_mso_2_init.pro
@@ -0,0 +1,28 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bcain_mso_2_init, MEX = Mex, MGS = mgs, MAVEN = maven,   SAMPLING = sampling,   $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bcain_mso_init, /Z, SAMPLING = '+ string(sampling)
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(MEX) then cmd += ',/MEX' 
+    
+  if KEYWORD_SET(MGS) then cmd += ',/MGS' 
+ 
+  if KEYWORD_SET(MAVEN) then cmd += ',/MAVEN'
+  
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bcain_mso_init.pro
@@ -0,0 +1,73 @@
+; cain field along the orbit
+;  
+
+pro bcain_mso_init,   MEX = Mex, MGS = mgs, MAVEN = maven, SAMPLING = sampling, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e 
+
+common GraphC, graph, GraphN, Item
+ 
+ names = ['B!DX', 'B!DY', 'B!DZ', 'B']
+ 
+ pen = !dnc-1 & comp = 100 & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;
+     
+     title = 'MEX';
+     graph[GraphN].SetCommonF[0] = 'mex_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getmex_xyz'; 
+ 
+  if KEYWORD_SET(MGS) then begin
+     title = 'MGS';
+     graph[GraphN].SetCommonF[0] = 'xyz_mgs_mso_com';
+     graph[GraphN].GetDataF[0] = 'getxyz_mgs_mso'; 
+   endif
+
+  if KEYWORD_SET(MAVEN) then begin     
+     title = 'MAVEN';
+     graph[GraphN].SetCommonF[0] = 'mav_xyz_mso_com';
+     graph[GraphN].GetDataF[0] = 'getmav_xyz_mso'; 
+   endif
+    
+ 
+frame = "MSO";
+graph[GraphN].ProcessF =  "CAIN";
+;if (KEYWORD_SET(FRAME)) then begin 
+;  if (frame eq "AC") then graph[GraphN].ProcessF= 'mso2ac';
+;endif 
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+graph[GraphN].Type = 'T96_T'
+graph[GraphN].DataN = 1 
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 120.0D0) < 36000.0D0 ;
+
+yTitle = title + ' Cain B '+frame+', nT';
+
+if (KEYWORD_SET(AY1TITLE)) then graph[GraphN].Ly.title = AY1Title $
+else graph[GraphN].Ly.title = yTitle
+
+graph[GraphN].AxisF = 'lineaxis,' + string(Graph[GraphN].TotalTime)
+
+if (KEYWORD_SET(X)) then comp = 0;
+if (KEYWORD_SET(Y)) then comp = 1;
+if (KEYWORD_SET(Z)) then comp = 2;
+if (KEYWORD_SET(TOTAL)) then comp = 3;
+
+if (KEYWORD_SET(SAMPLING)) then begin 
+    if (Sampling eq 60) then  graph[GraphN].PlotF = comp eq 100 ? 'PlotVec' : 'PlotVec_Comp_New, ' + string(comp) +','+ string(pen ) +','+ string(symbol) $  
+    else graph[GraphN].PlotF = comp eq 100 ? 'PlotVecCain,' + string(sampling) : 'PlotVecCompCain, '  + string(sampling) +','+ string(comp) +','+ string(pen ) +','+ string(symbol)
+endif else begin
+  graph[GraphN].PlotF = comp eq 100 ? 'PlotVec' : 'PlotVec_Comp_New, ' + string(comp) +','+ string(pen ) +','+ string(symbol)
+endelse
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) = (comp le 3) ? names : names[0:2];
+ 
+graph[GraphN].LastPlotF = 'nodatalast' 
+return
+end
+++ a/amda_plus/amdalib/bcain_tot_init.pro
@@ -0,0 +1,28 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bcain_tot_init, MEX = Mex, MGS = mgs, MAVEN = maven,  SAMPLING = sampling,    $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bcain_mso_init, /TOT, SAMPLING = ' + string(sampling)
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(MEX) then cmd += ',/MEX' 
+    
+  if KEYWORD_SET(MGS) then cmd += ',/MGS' 
+ 
+  if KEYWORD_SET(MAVEN) then cmd += ',/MAVEN'
+  
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bmorsch_mso_0_init.pro
@@ -0,0 +1,28 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bmorsch_mso_0_init, MEX = Mex, MGS = mgs, MAVEN = maven, SAMPLING = sampling,     $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bmorsch_mso_init, /X, SAMPLING = ' + string(sampling)
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(MEX) then cmd += ',/MEX' 
+    
+  if KEYWORD_SET(MGS) then cmd += ',/MGS' 
+ 
+  if KEYWORD_SET(MAVEN) then cmd += ',/MAVEN'
+  
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bmorsch_mso_1_init.pro
@@ -0,0 +1,28 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bmorsch_mso_1_init, MEX = Mex, MGS = mgs, MAVEN = maven, SAMPLING = sampling,     $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bmorsch_mso_init, /Y, SAMPLING = ' + string(sampling)
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(MEX) then cmd += ',/MEX' 
+    
+  if KEYWORD_SET(MGS) then cmd += ',/MGS' 
+ 
+  if KEYWORD_SET(MAVEN) then cmd += ',/MAVEN'
+  
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bmorsch_mso_2_init.pro
@@ -0,0 +1,28 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bmorsch_mso_2_init, MEX = Mex, MGS = mgs, MAVEN = maven, SAMPLING = sampling,     $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bmorsch_mso_init, /Z, SAMPLING = ' + string(sampling)
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(MEX) then cmd += ',/MEX' 
+    
+  if KEYWORD_SET(MGS) then cmd += ',/MGS' 
+ 
+  if KEYWORD_SET(MAVEN) then cmd += ',/MAVEN'
+  
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bmorsch_mso_init.pro
@@ -0,0 +1,73 @@
+; morsch field along the orbit
+;  
+
+pro bmorsch_mso_init,   MEX = Mex, MGS = mgs, MAVEN = maven, SAMPLING = sampling, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e 
+
+common GraphC, graph, GraphN, Item
+ 
+ names = ['B!DX', 'B!DY', 'B!DZ', 'B']
+ 
+ pen = !dnc-1 & comp = 100 & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;
+     
+     title = 'MEX';
+     graph[GraphN].SetCommonF[0] = 'mex_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getmex_xyz'; 
+ 
+  if KEYWORD_SET(MGS) then begin
+     title = 'MGS';
+     graph[GraphN].SetCommonF[0] = 'xyz_mgs_mso_com';
+     graph[GraphN].GetDataF[0] = 'getxyz_mgs_mso'; 
+   endif
+
+  if KEYWORD_SET(MAVEN) then begin     
+     title = 'MAVEN';
+     graph[GraphN].SetCommonF[0] = 'mav_xyz_mso_com';
+     graph[GraphN].GetDataF[0] = 'getmav_xyz_mso'; 
+   endif
+    
+ 
+frame = "MSO";
+graph[GraphN].ProcessF =  "MORSCH";
+;if (KEYWORD_SET(FRAME)) then begin 
+;  if (frame eq "AC") then graph[GraphN].ProcessF= 'mso2ac';
+;endif 
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+graph[GraphN].Type = 'T96_T'
+graph[GraphN].DataN = 1 
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 120.0D0) < 36000.0D0 ;
+
+yTitle = title + ' Morschhauser B '+frame+', nT';
+
+if (KEYWORD_SET(AY1TITLE)) then graph[GraphN].Ly.title = AY1Title $
+else graph[GraphN].Ly.title = yTitle
+
+graph[GraphN].AxisF = 'lineaxis,' + string(Graph[GraphN].TotalTime)
+
+if (KEYWORD_SET(X)) then comp = 0;
+if (KEYWORD_SET(Y)) then comp = 1;
+if (KEYWORD_SET(Z)) then comp = 2;
+if (KEYWORD_SET(TOTAL)) then comp = 3;
+
+if (KEYWORD_SET(SAMPLING)) then begin 
+    if (Sampling eq 60) then  graph[GraphN].PlotF = comp eq 100 ? 'PlotVec' : 'PlotVec_Comp_New, ' + string(comp) +','+ string(pen ) +','+ string(symbol) $  
+    else graph[GraphN].PlotF = comp eq 100 ? 'PlotVecCain,' + string(sampling) : 'PlotVecCompCain, '  + string(sampling) +','+ string(comp) +','+ string(pen ) +','+ string(symbol)
+endif else begin
+  graph[GraphN].PlotF = comp eq 100 ? 'PlotVec' : 'PlotVec_Comp_New, ' + string(comp) +','+ string(pen ) +','+ string(symbol)
+endelse
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) = (comp le 3) ? names : names[0:2];
+ 
+graph[GraphN].LastPlotF = 'nodatalast' 
+return
+end
+++ a/amda_plus/amdalib/bmorsch_tot_init.pro
@@ -0,0 +1,28 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bmorsch_tot_init, MEX = Mex, MGS = mgs, MAVEN = maven,  SAMPLING = sampling,    $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bmorsch_mso_init, /TOT, SAMPLING = ' + string(sampling)
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(MEX) then cmd += ',/MEX' 
+    
+  if KEYWORD_SET(MGS) then cmd += ',/MGS' 
+ 
+  if KEYWORD_SET(MAVEN) then cmd += ',/MAVEN'
+  
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bt96_0_init.pro
@@ -0,0 +1,41 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bt96_0_init,     CLUSTER1 = Cluster1, CLUSTER2 = Cluster2, CLUSTER3 = Cluster3, CLUSTER4 = Cluster4, $
+                     THEMISA = ThemisA, THEMISB = ThemisB, THEMISC = ThemisC, $
+                     THEMISD = ThemisD, THEMISE = ThemisE, $
+                     GEOTAIL = Geotail, $
+                     OMNI = OMNI, AMDA = AMDA, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bt96_init, /X'
+  if KEYWORD_SET(FRAME) then cmd += ',FRAME = ' + frame
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(CLUSTER1) then cmd += ',/CLUSTER1' 
+  if KEYWORD_SET(CLUSTER2) then cmd += ',/CLUSTER2'  
+  if KEYWORD_SET(CLUSTER3) then cmd += ',/CLUSTER3'
+  if KEYWORD_SET(CLUSTER4) then cmd += ',/CLUSTER4'  
+  
+  if KEYWORD_SET(THEMISA) then cmd += ',/THEMISA' 
+  if KEYWORD_SET(THEMISB) then cmd += ',/THEMISB'  
+  if KEYWORD_SET(THEMISC) then cmd += ',/THEMISC'
+  if KEYWORD_SET(THEMISD) then cmd += ',/THEMISD'  
+  if KEYWORD_SET(THEMISE) then cmd += ',/THEMISE'
+
+  if KEYWORD_SET(GEOTAIL) then cmd += ',/GEOTAIL'
+
+  if KEYWORD_SET(OMNI) then cmd += ',/OMNI' 
+  if KEYWORD_SET(AMDA) then cmd += ',/AMDA' 
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bt96_1_init.pro
@@ -0,0 +1,41 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bt96_1_init,     CLUSTER1 = Cluster1, CLUSTER2 = Cluster2, CLUSTER3 = Cluster3, CLUSTER4 = Cluster4, $
+                     THEMISA = ThemisA, THEMISB = ThemisB, THEMISC = ThemisC, $
+                     THEMISD = ThemisD, THEMISE = ThemisE, $
+                     GEOTAIL = Geotail, $
+                     OMNI = OMNI, AMDA = AMDA, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bt96_init, /Y'
+  if KEYWORD_SET(FRAME) then cmd += ',FRAME = ' + frame
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(CLUSTER1) then cmd += ',/CLUSTER1' 
+  if KEYWORD_SET(CLUSTER2) then cmd += ',/CLUSTER2'  
+  if KEYWORD_SET(CLUSTER3) then cmd += ',/CLUSTER3'
+  if KEYWORD_SET(CLUSTER4) then cmd += ',/CLUSTER4'  
+  
+  if KEYWORD_SET(THEMISA) then cmd += ',/THEMISA' 
+  if KEYWORD_SET(THEMISB) then cmd += ',/THEMISB'  
+  if KEYWORD_SET(THEMISC) then cmd += ',/THEMISC'
+  if KEYWORD_SET(THEMISD) then cmd += ',/THEMISD'  
+  if KEYWORD_SET(THEMISE) then cmd += ',/THEMISE'
+
+  if KEYWORD_SET(GEOTAIL) then cmd += ',/GEOTAIL'
+
+  if KEYWORD_SET(OMNI) then cmd += ',/OMNI' 
+  if KEYWORD_SET(AMDA) then cmd += ',/AMDA' 
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bt96_2_init.pro
@@ -0,0 +1,41 @@
+;----------------ENVELOPE-------------------------
+; 
+pro bt96_2_init,     CLUSTER1 = Cluster1, CLUSTER2 = Cluster2, CLUSTER3 = Cluster3, CLUSTER4 = Cluster4, $
+                     THEMISA = ThemisA, THEMISB = ThemisB, THEMISC = ThemisC, $
+                     THEMISD = ThemisD, THEMISE = ThemisE, $
+                     GEOTAIL = Geotail, $
+                     OMNI = OMNI, AMDA = AMDA, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e  
+
+  cmd = 'bt96_init, /Z'
+  if KEYWORD_SET(FRAME) then cmd += ',FRAME = ' + frame
+
+  if KEYWORD_SET(COLOR) then cmd += ',COLOR = ' + color;
+  if KEYWORD_SET(AY1TITLE) then cmd += ',AY1TITLE = ' + AY1Title;
+
+  if KEYWORD_SET(CLUSTER1) then cmd += ',/CLUSTER1' 
+  if KEYWORD_SET(CLUSTER2) then cmd += ',/CLUSTER2'  
+  if KEYWORD_SET(CLUSTER3) then cmd += ',/CLUSTER3'
+  if KEYWORD_SET(CLUSTER4) then cmd += ',/CLUSTER4'  
+  
+  if KEYWORD_SET(THEMISA) then cmd += ',/THEMISA' 
+  if KEYWORD_SET(THEMISB) then cmd += ',/THEMISB'  
+  if KEYWORD_SET(THEMISC) then cmd += ',/THEMISC'
+  if KEYWORD_SET(THEMISD) then cmd += ',/THEMISD'  
+  if KEYWORD_SET(THEMISE) then cmd += ',/THEMISE'
+
+  if KEYWORD_SET(GEOTAIL) then cmd += ',/GEOTAIL'
+
+  if KEYWORD_SET(OMNI) then cmd += ',/OMNI' 
+  if KEYWORD_SET(AMDA) then cmd += ',/AMDA' 
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/bt96_init.pro
@@ -0,0 +1,108 @@
+; tsyganenko 96 field along the orbit
+;  
+
+pro bt96_init,   CLUSTER1 = Cluster1, CLUSTER2 = Cluster2, CLUSTER3 = Cluster3, CLUSTER4 = Cluster4, $
+                     THEMISA = ThemisA, THEMISB = ThemisB, THEMISC = ThemisC, $
+                     THEMISD = ThemisD, THEMISE = ThemisE, $
+                     GEOTAIL = Geotail, $
+                     OMNI = OMNI, AMDA = AMDA, $
+                     PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                    
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $                    
+                     ; all additional keywords
+                    _EXTRA = e 
+
+common GraphC, graph, GraphN, Item
+ 
+ names = ['B!DX', 'B!DY', 'B!DZ', 'B']
+ 
+ pen = !dnc-1 & comp = 100 & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;
+     
+     title = 'CL1';
+     graph[GraphN].SetCommonF[0] = 'c1_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getc1_xyz'; 
+ 
+  if KEYWORD_SET(CLUSTER2) then begin
+     title = 'CL2';
+     graph[GraphN].SetCommonF[0] = 'c2_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getc2_xyz'; 
+   endif
+  if KEYWORD_SET(CLUSTER3) then begin     
+     title = 'CL3';
+     graph[GraphN].SetCommonF[0] = 'c3_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getc3_xyz'; 
+   endif
+  if KEYWORD_SET(CLUSTER4) then begin     
+     title = 'CL4'
+     graph[GraphN].SetCommonF[0] = 'c4_xyz_com';
+     graph[GraphN].GetDataF[0] = 'getc4_xyz'; 
+   endif
+  if KEYWORD_SET(THEMISA) then begin
+    title = 'THA'
+    graph[GraphN].SetCommonF[0] = 'tha_xyz_com'; 
+    graph[GraphN].GetDataF[0] ='gettha_xyz';
+   endif
+ if KEYWORD_SET(THEMISB) then begin
+    title = 'THB'
+    graph[GraphN].SetCommonF[0] = 'thb_xyz_com'; 
+    graph[GraphN].GetDataF[0] ='getthb_xyz';
+   endif
+ if KEYWORD_SET(THEMISC) then begin
+    title = 'THC'
+    graph[GraphN].SetCommonF[0] = 'thc_xyz_com'; 
+    graph[GraphN].GetDataF[0] ='getthc_xyz';
+   endif
+ if KEYWORD_SET(THEMISD) then begin
+    title = 'THD'
+    graph[GraphN].SetCommonF[0] = 'thd_xyz_com'; 
+    graph[GraphN].GetDataF[0] = 'getthd_xyz';
+   endif
+ if KEYWORD_SET(THEMISE) then begin
+    title = 'THE'
+    graph[GraphN].SetCommonF[0] = 'the_xyz_com'; 
+    graph[GraphN].GetDataF[0] = 'getthe_xyz';
+   endif
+ if KEYWORD_SET(GEOTAIL) then begin
+    title = 'GTL'
+    graph[GraphN].SetCommonF[0] = 'gtl_xyz_com'; 
+    graph[GraphN].GetDataF[0] = 'getgtl_xyz';
+   endif
+
+t96Function = KEYWORD_SET(OMNI) ? 'T96OMNI' : 'T96';
+ 
+frame = "GSM";
+graph[GraphN].ProcessF= t96Function; 
+  
+if (KEYWORD_SET(FRAME)) then begin 
+  if (frame eq "GSE") then graph[GraphN].ProcessF= t96Function + ';gsm2gse';
+  if (frame eq "SM")  then graph[GraphN].ProcessF= t96Function + ';gse2sm';
+endif 
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+graph[GraphN].Type = 'T96_T'
+graph[GraphN].DataN = 1 
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 120.0D0) < 36000.0D0 ;
+
+yTitle = title + ' T96 B '+frame+', nT';
+
+if (KEYWORD_SET(AY1TITLE)) then graph[GraphN].Ly.title = AY1Title $
+else graph[GraphN].Ly.title = yTitle
+
+graph[GraphN].AxisF = 'lineaxis,' + string(Graph[GraphN].TotalTime)
+
+if (KEYWORD_SET(X)) then comp = 0;
+if (KEYWORD_SET(Y)) then comp = 1;
+if (KEYWORD_SET(Z)) then comp = 2;
+if (KEYWORD_SET(TOTAL)) then comp = 3;
+
+graph[GraphN].PlotF = comp eq 100 ? 'PlotVec': 'PlotVec_Comp_New, ' + string(comp) +','+ string(pen ) +','+ string(symbol)
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) = (comp le 3) ? names : names[0:2];
+ 
+graph[GraphN].LastPlotF = 'nodatalast' 
+return
+end
+++ a/amda_plus/amdalib/c1_hia_pad_init.pro
@@ -0,0 +1,33 @@
+
+pro c1_hia_pad_init,  PA = PA,  ENERGY = energy, LIN = lin, LOG = log
+
+Pitch = '0';
+PitchValue = [  '0-11.25','11.25-22.5','22.5-33.75','33.75-45', $
+		'45-56.25','56.25-67.5','67.5-78.75', $ 
+		'78.75-90','90-101.25','101.25-112.5','112.5-123.75', $
+		'123.75-135','135-146.25','146.25-157.5','157.5-168.75', $
+		'168.75-180'];
+ 
+if (KEYWORD_SET(PA)) then Pitch = string(PA-1);
+    
+En = 31;
+if (KEYWORD_SET(ENERGY)) then En = Energy;
+
+    
+       if (En LT 31) then begin
+          common GraphC, graph, GraphN, Item
+          graph[GraphN].Type = 'FGM_T' 
+	  graph[GraphN].DataN = 1
+	  graph[GraphN].SetCommonF[0] = 'pad_hia_com, 1' 
+	  graph[GraphN].GetDataF[0] ='get_pad_hia, 1'
+	  graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) <  18000.0D0 
+	  graph[GraphN].Ly.title = 'C1 hia Flux !CP_'+PitchValue[Pitch] + ' E eV '; 
+	  graph[GraphN].AxisF =  KEYWORD_SET(log) ?  'logaxis' : 'lineaxis, 0';
+	  graph[GraphN].PlotF = 'spec3D_comp_plot,' + string(Pitch) +','+string(En);
+   endif else begin
+	  functionName = 'P' + strtrim(Pitch,1) + '_hia_c1_init';
+	  call_procedure, functionName;
+   endelse
+
+return
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/c1_whi_ne_.pro
@@ -0,0 +1,50 @@
+
+; $Id: whi_ne_c1_.pro,v 1.1 2010/02/22 09:36:19 budnik Exp $
+;
+function c1_whi_ne_, StartTime, TimeInt, Time, RetSize, MaxSampling, close=close 
+ common GraphC, graph, graphN
+ common WHINEC1_COM, ID
+
+   FillValue = -1.0
+ 
+   if n_elements(ID) EQ 0 then ID = -1L
+   Sampling = 1.5
+   MaxSampling = 104.0
+   Size = long(graph[graphN].DeltaT/sampling)
+
+    VI = 'clust1:whi:whi'    
+    Val = fltarr(Size)
+    
+    Time = dblarr(Size)
+
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                             ID, VI, 'Density', StartTime, TimeInt, Size, Time, Val)
+ 
+     if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+     endif
+
+    if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, 0.0
+   endif
+   
+   Time=Time[0:Retsize-1]
+   Val=Val[0:Retsize-1]
+   
+   num_gap = where(Val ne FillValue, Cntr)
+  
+   if Cntr le 0 then begin
+       Retsize = -999
+       Time = 0.D0
+       return, [0.0, 0.0, 0.0]
+   endif
+   
+    Time=Time[num_gap]
+    Val=Val[num_gap]
+    Retsize=Cntr
+  
+return, Val
+end
+++ a/amda_plus/amdalib/c1_whi_ne_init.pro
@@ -0,0 +1,38 @@
+;-----------------WHISPER-PP_INIT--------------------------
+;
+pro c1_whi_ne_INIT,  LIN = LIN, LOG = LOG, PSYM = psym, COLOR = color, CONTRAST = CONTRAST, UNCERTAINTY =  UNCERTAINTY 
+
+common GraphC, graph, GraphN, Item
+common Request, List, StartTime, TimeInt
+   
+ names = ['Contrast', '0.7-1.0', '0.4-0.7', '0.0-0.4', 'other src']  
+ 
+ quality = '0';
+ error = '0';
+
+ pen = !dnc-1  & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;
+
+ if KEYWORD_SET(CONTRAST) then quality = '1'; 
+ if KEYWORD_SET(UNCERTAINTY) then error = '1'; 
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) =  names;
+graph[GraphN].Type = 'STAFF_T'
+graph[GraphN].DataN = 3
+graph[GraphN].SetCommonF[0] = 'whi_ne_com, 1'
+graph[GraphN].SetCommonF[1] = 'whi_ne_qual_com, 1'
+graph[GraphN].SetCommonF[2] = 'whi_ne_err_com, 1'
+graph[GraphN].GetDataF[0] ='get_whi_ne, 1'
+graph[GraphN].GetDataF[1] ='get_whi_ne_qual, 1'
+graph[GraphN].GetDataF[2] ='get_whi_ne_err, 1'
+graph[GraphN].DeltaT = (graph[GraphN].TotalTime + 10.0D0) <  18000.0D0 
+graph[GraphN].Ly.title = 'C1 WHI!CNe, cm!U-3 '
+graph[GraphN].AxisF =  KEYWORD_SET(LIN) ?  'lineaxis,0' : 'logaxis'
+graph[GraphN].PlotF = 'plot_whi_ne_new,' + quality + ','+ error + ','+ string(pen) + ',' + string(symbol)
+ 
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/c2_whi_ne_.pro
@@ -0,0 +1,50 @@
+
+; $Id: whi_ne_c2_.pro,v 1.1 2010/02/22 09:36:19 budnik Exp $
+;
+function c2_whi_ne_c, StartTime, TimeInt, Time, RetSize, MaxSampling, close=close 
+ common GraphC, graph, graphN
+ common WHINEC2_COM, ID
+
+   FillValue = -1.0
+ 
+   if n_elements(ID) EQ 0 then ID = -1L
+   Sampling = 1.5
+   MaxSampling = 104.0
+   Size = long(graph[graphN].DeltaT/sampling)
+
+    VI = 'clust2:whi:whi'    
+    Val = fltarr(Size)
+    
+    Time = dblarr(Size)
+
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                             ID, VI, 'Density', StartTime, TimeInt, Size, Time, Val)
+ 
+     if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+     endif
+
+    if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, 0.0
+   endif
+   
+   Time=Time[0:Retsize-1]
+   Val=Val[0:Retsize-1]
+   
+   num_gap = where(Val ne FillValue, Cntr)
+  
+   if Cntr le 0 then begin
+       Retsize = -999
+       Time = 0.D0
+       return, [0.0, 0.0, 0.0]
+   endif
+   
+    Time=Time[num_gap]
+    Val=Val[num_gap]
+    Retsize=Cntr
+  
+return, Val
+end
+++ a/amda_plus/amdalib/c2_whi_ne_init.pro
@@ -0,0 +1,38 @@
+;-----------------WHISPER-PP_INIT--------------------------
+;
+pro c2_whi_ne_INIT,  LIN = LIN, LOG = LOG, PSYM = psym, COLOR = color, CONTRAST = CONTRAST, UNCERTAINTY =  UNCERTAINTY 
+
+common GraphC, graph, GraphN, Item
+common Request, List, StartTime, TimeInt
+   
+ names = ['Contrast', '0.7-1.0', '0.4-0.7', '0.0-0.4', 'other src']  
+ 
+ quality = '0';
+ error = '0';
+
+ pen = !dnc-1  & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;
+
+ if KEYWORD_SET(CONTRAST) then quality = '1'; 
+ if KEYWORD_SET(UNCERTAINTY) then error = '1'; 
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) =  names;
+graph[GraphN].Type = 'STAFF_T'
+graph[GraphN].DataN = 3
+graph[GraphN].SetCommonF[0] = 'whi_ne_com, 2'
+graph[GraphN].SetCommonF[1] = 'whi_ne_qual_com, 2'
+graph[GraphN].SetCommonF[2] = 'whi_ne_err_com, 2'
+graph[GraphN].GetDataF[0] ='get_whi_ne, 2'
+graph[GraphN].GetDataF[1] ='get_whi_ne_qual, 2'
+graph[GraphN].GetDataF[2] ='get_whi_ne_err, 2'
+graph[GraphN].DeltaT = (graph[GraphN].TotalTime + 10.0D0) <  18000.0D0 
+graph[GraphN].Ly.title = 'C2 WHI!CNe, cm!U-3 '
+graph[GraphN].AxisF =  KEYWORD_SET(LIN) ?  'lineaxis,0' : 'logaxis'
+graph[GraphN].PlotF = 'plot_whi_ne_new,' + quality + ','+ error + ','+ string(pen) + ',' + string(symbol)
+ 
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/c3_hia_pad_init.pro
@@ -0,0 +1,33 @@
+
+pro c3_hia_pad_init,  PA = PA,  ENERGY = energy, LIN = lin, LOG = log
+
+Pitch = '0';
+PitchValue = [  '0-11.25','11.25-22.5','22.5-33.75','33.75-45', $
+		'45-56.25','56.25-67.5','67.5-78.75', $ 
+		'78.75-90','90-101.25','101.25-112.5','112.5-123.75', $
+		'123.75-135','135-146.25','146.25-157.5','157.5-168.75', $
+		'168.75-180'];
+ 
+if (KEYWORD_SET(PA)) then Pitch = string(PA-1);
+    
+En = 31;
+if (KEYWORD_SET(ENERGY)) then En = Energy;
+
+    
+       if (En LT 31) then begin
+          common GraphC, graph, GraphN, Item
+          graph[GraphN].Type = 'FGM_T' 
+	  graph[GraphN].DataN = 1
+	  graph[GraphN].SetCommonF[0] = 'pad_hia_com, 3' 
+	  graph[GraphN].GetDataF[0] ='get_pad_hia, 3'
+	  graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) <  18000.0D0 
+	  graph[GraphN].Ly.title = 'C3 hia Flux !CP_'+PitchValue[Pitch] + ' E eV '; 
+	  graph[GraphN].AxisF =  KEYWORD_SET(log) ?  'logaxis' : 'lineaxis, 0';
+	  graph[GraphN].PlotF = 'spec3D_comp_plot,' + string(Pitch) +','+string(En);
+   endif else begin
+	  functionName = 'P' + strtrim(Pitch,1) + '_hia_c3_init';
+	  call_procedure, functionName;
+   endelse
+
+return
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/c3_whi_ne_.pro
@@ -0,0 +1,50 @@
+
+; $Id: whi_ne_c3_.pro,v 1.1 2010/02/22 09:36:19 budnik Exp $
+;
+function c3_whi_ne_, StartTime, TimeInt, Time, RetSize, MaxSampling, close=close 
+ common GraphC, graph, graphN
+ common WHINEC3_COM, ID
+
+   FillValue = -1.0
+ 
+   if n_elements(ID) EQ 0 then ID = -1L
+   Sampling = 1.5
+   MaxSampling = 104.0
+   Size = long(graph[graphN].DeltaT/sampling)
+
+    VI = 'clust3:whi:whi'    
+    Val = fltarr(Size)
+    
+    Time = dblarr(Size)
+
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                             ID, VI, 'Density', StartTime, TimeInt, Size, Time, Val)
+ 
+     if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+     endif
+
+    if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, 0.0
+   endif
+   
+   Time=Time[0:Retsize-1]
+   Val=Val[0:Retsize-1]
+   
+   num_gap = where(Val ne FillValue, Cntr)
+  
+   if Cntr le 0 then begin
+       Retsize = -999
+       Time = 0.D0
+       return, [0.0, 0.0, 0.0]
+   endif
+   
+    Time=Time[num_gap]
+    Val=Val[num_gap]
+    Retsize=Cntr
+  
+return, Val
+end
+++ a/amda_plus/amdalib/c3_whi_ne_init.pro
@@ -0,0 +1,38 @@
+;-----------------WHISPER-PP_INIT--------------------------
+;
+pro c3_whi_ne_INIT,  LIN = LIN, LOG = LOG, PSYM = psym, COLOR = color, CONTRAST = CONTRAST, UNCERTAINTY =  UNCERTAINTY 
+
+common GraphC, graph, GraphN, Item
+common Request, List, StartTime, TimeInt
+   
+ names = ['Contrast', '0.7-1.0', '0.4-0.7', '0.0-0.4', 'other src']  
+ 
+ quality = '0';
+ error = '0';
+
+ pen = !dnc-1  & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;
+
+ if KEYWORD_SET(CONTRAST) then quality = '1'; 
+ if KEYWORD_SET(UNCERTAINTY) then error = '1'; 
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) =  names;
+graph[GraphN].Type = 'STAFF_T'
+graph[GraphN].DataN = 3
+graph[GraphN].SetCommonF[0] = 'whi_ne_com, 3'
+graph[GraphN].SetCommonF[1] = 'whi_ne_qual_com, 3'
+graph[GraphN].SetCommonF[2] = 'whi_ne_err_com, 3'
+graph[GraphN].GetDataF[0] ='get_whi_ne, 3'
+graph[GraphN].GetDataF[1] ='get_whi_ne_qual, 3'
+graph[GraphN].GetDataF[2] ='get_whi_ne_err, 3'
+graph[GraphN].DeltaT = (graph[GraphN].TotalTime + 10.0D0) <  18000.0D0 
+graph[GraphN].Ly.title = 'C3 WHI!CNe, cm!U-3 '
+graph[GraphN].AxisF =  KEYWORD_SET(LIN) ?  'lineaxis,0' : 'logaxis'
+graph[GraphN].PlotF = 'plot_whi_ne_new,' + quality + ','+ error + ','+ string(pen) + ',' + string(symbol)
+ 
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/c4_whi_ne_.pro
@@ -0,0 +1,50 @@
+
+; $Id: whi_ne_c4_.pro,v 1.1 2010/02/22 09:36:19 budnik Exp $
+;
+function c4_whi_ne_, StartTime, TimeInt, Time, RetSize, MaxSampling, close=close 
+ common GraphC, graph, graphN
+ common WHINEC4_COM, ID
+
+   FillValue = -1.0
+ 
+   if n_elements(ID) EQ 0 then ID = -1L
+   Sampling = 1.5
+   MaxSampling = 104.0
+   Size = long(graph[graphN].DeltaT/sampling)
+
+    VI = 'clust4:whi:whi'    
+    Val = fltarr(Size)
+    
+    Time = dblarr(Size)
+
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                             ID, VI, 'Density', StartTime, TimeInt, Size, Time, Val)
+ 
+     if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+     endif
+
+    if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, 0.0
+   endif
+   
+   Time=Time[0:Retsize-1]
+   Val=Val[0:Retsize-1]
+   
+   num_gap = where(Val ne FillValue, Cntr)
+  
+   if Cntr le 0 then begin
+       Retsize = -999
+       Time = 0.D0
+       return, [0.0, 0.0, 0.0]
+   endif
+   
+    Time=Time[num_gap]
+    Val=Val[num_gap]
+    Retsize=Cntr
+  
+return, Val
+end
+++ a/amda_plus/amdalib/c4_whi_ne_init.pro
@@ -0,0 +1,38 @@
+;-----------------WHISPER-PP_INIT--------------------------
+;
+pro c4_whi_ne_INIT,  LIN = LIN, LOG = LOG, PSYM = psym, COLOR = color, CONTRAST = CONTRAST, UNCERTAINTY =  UNCERTAINTY 
+
+common GraphC, graph, GraphN, Item
+common Request, List, StartTime, TimeInt
+   
+ names = ['Contrast', '0.7-1.0', '0.4-0.7', '0.0-0.4', 'other src']  
+ 
+ quality = '0';
+ error = '0';
+
+ pen = !dnc-1  & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;
+
+ if KEYWORD_SET(CONTRAST) then quality = '1'; 
+ if KEYWORD_SET(UNCERTAINTY) then error = '1'; 
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) =  names;
+graph[GraphN].Type = 'STAFF_T'
+graph[GraphN].DataN = 3
+graph[GraphN].SetCommonF[0] = 'whi_ne_com, 4'
+graph[GraphN].SetCommonF[1] = 'whi_ne_qual_com, 4'
+graph[GraphN].SetCommonF[2] = 'whi_ne_err_com, 4'
+graph[GraphN].GetDataF[0] ='get_whi_ne, 4'
+graph[GraphN].GetDataF[1] ='get_whi_ne_qual, 4'
+graph[GraphN].GetDataF[2] ='get_whi_ne_err, 4'
+graph[GraphN].DeltaT = (graph[GraphN].TotalTime + 10.0D0) <  18000.0D0 
+graph[GraphN].Ly.title = 'C4 WHI!CNe, cm!U-3 '
+graph[GraphN].AxisF =  KEYWORD_SET(LIN) ?  'lineaxis,0' : 'logaxis'
+graph[GraphN].PlotF = 'plot_whi_ne_new,' + quality + ','+ error + ','+ string(pen) + ',' + string(symbol)
+ 
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/cain.pro
@@ -0,0 +1,26 @@
+pro cain, Time, Orbit
+
+  FillValue = -1.e31
+  Number = N_elements(Time)
+  Bmso = fltarr(3)
+  Rm = 3393.0;
+
+  num31 = where(Orbit[0,*] ne FillValue, Index31)
+  if (Index31 le 0) then return
+
+  Time = Time[num31];
+  Orbit = Orbit[*,num31]*Rm;
+  err = call_external(!PROJLIB + '/idl_cain.so', 'CainSetUp');
+
+  for i = 0L, Index31-1L do begin
+          julday = doubl2jd(Time[i]);
+          err = call_external(!PROJLIB + '/idl_cain.so', 'CainField', julday, Orbit[*,i], Bmso);
+          if (Bmso[0] eq 0.0) then Bmso(*) = FillValue; 
+          Orbit[*,i] =  Bmso;                
+  endfor
+
+ return
+end
+
+  
+  
 \ No newline at end of file
+++ a/amda_plus/amdalib/doubl2jd.pro
@@ -0,0 +1,22 @@
+;----------------------------------------------------------
+;                  IMA COMMON LIB (IDL)
+;                     doubl2jd.pro
+;                         V.1.0
+;
+;    Function which Julian Day time from DD double time
+;
+;    Versions:
+;       Jan 25 2005, V.1.0
+;===========================================================
+ 
+function Doubl2JD, Time  
+     CTimeS =  call_external(!DDLIB + 'DD_idl.so','GetTime', Time, /S_VALUE)
+     Year =  fix(STRMID(CTimeS, 0, 4) )
+     Day =  double(STRMID(CTimeS, 4, 3))
+     Hour =  fix(STRMID(CTimeS, 7, 2))
+     Minute =  fix(STRMID(CTimeS, 9, 2))
+     Second =  fix(STRMID(CTimeS, 11, 2))
+     JDReal = JULDAY(01, 01, Year, Hour, Minute, Second)  +  Day
+   return,JDReal 
+end
+   
+++ a/amda_plus/amdalib/dst_init.pro
@@ -0,0 +1,23 @@
+;----------------DST_INIT--------------------------
+; 
+pro DST_INIT, LOG=LOG, LIN=LIN, RED=RED, BLUE=BLUE, GREEN=GREEN
+
+common GraphC, graph, GraphN, Item
+pen = 0
+
+ if KEYWORD_SET(RED) then pen = 3
+ if KEYWORD_SET(BLUE) then pen = 1  
+ if KEYWORD_SET(GREEN) then pen = 2
+   
+graph[GraphN].Type = 'DST_T'
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'dstcom' 
+graph[GraphN].GetDataF[0] ='getdst' 
+graph[GraphN].DeltaT = (graph[graphN].TotalTime+120.0D0) < 86400.0D0
+graph[GraphN].Ly.title = 'DST, nT'
+graph[GraphN].AxisF = 'lineaxis,' + string(graph[graphN].TotalTime)
+graph[GraphN].PlotF = 'dstplot,' + string(pen);
+ 
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/dstplot.pro
@@ -0,0 +1,33 @@
+;============================================================
+;
+;     DST Time Profile PLOT
+;
+;  
+;=============================================================
+
+  pro dstplot, pen
+; procedure to plot prepared data on prepared graph
+     
+  common GraphC, Graph, GraphN
+  colour = [!dnc - 1, !dnc*20/256, !dnc*90/256, !dnc - 2]  
+
+  Time =  Graph[GraphN].Time[0]
+  if (N_elements(*Time) eq 0) then return    
+  Val = (*Graph[GraphN].Val[0])
+  numMinMax = where(Val gt -10000.0, IndexMinMax)
+  if IndexMinMax LE 0 then return 
+  MinVal = min(Val[numMinMax])
+  MaxVal = max(Val[numMinMax]) 
+  if Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] then begin
+                  Graph[GraphN].Ly.range = [MinVal, MaxVal]
+                  Graph[GraphN].Ly.style = 0
+   endif else Graph[GraphN].Ly.style = 1 
+  
+    !y = Graph[GraphN].Ly
+    
+    if (execute(Graph[GraphN].AxisF) NE 1) then plot, *Time, Val, PSYM = 10, color = !DNC-1  else $ 
+    oplot, *Time, Val, PSYM = 10,  color=colour[pen];
+  
+return
+end
+;--------------------------------------------------------------------
+++ a/amda_plus/amdalib/e_mgs_omni_init.pro
@@ -0,0 +1,64 @@
+;---------------------------------------------------
+; $Id: e_mgs_omni_init.pro,v 1.3 2013/03/28 08:40:06 budnik Exp $
+;              MGS ER SPECTRA  OMNI init
+;                
+;----------------------------------------------------
+
+pro e_mgs_omni_init, LOG = log, LIN = lin, PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $
+                     ; spectra
+                     MINCOUNT = minCount, MAXCOUNT = maxCount, BACKGRDCHAN = backgrdChan, BACKGRDVAL = backgrdVal, $
+                     CHANNEL = channel, Energy = energy, ENERGYMIN = energymin, ENERGYMAX =  energymax, $
+                     
+                     ; vector                     
+                     DELAY = delay, $                      
+                     ; all additional keywords
+                    _EXTRA = e
+
+  common Request, List, StartTime, TimeInt
+  common GraphC, graph, GraphN, Item
+
+ pen = !dnc-1 & comp = 100 & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;
+   
+  graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+  *(graph[GraphN].names) = 'log Flux 1/(cm!U2!N*s*ster*eV)'
+  
+  graph[GraphN].Type = 'ER_T'
+  graph[GraphN].DataN = 1
+  graph[GraphN].SetCommonF[0] = 'e_mgs_omni_com'
+  graph[GraphN].GetDataF[0] = 'get_e_mgs_omni'
+
+  graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) < 36100.0D0
+
+EnMin = 11;
+EnMax = 16217;
+
+energyTitle = '';
+
+  if (KEYWORD_SET(energymin)) then begin
+   EnMin = energymin;
+   if (KEYWORD_SET(energymax)) then EnMax =  energymax else EnMax = energymin;
+   energyTitle = energymin eq energymax ? '!CE ' + strtrim(string(energymin),1) + 'eV' : '!CE ' + strtrim(string(energymin),1) + '-'+ strtrim(string(energymax),1) + ' eV'; 
+  endif else  energyTitle = '!Call energies';
+
+  options = string(EnMin) +','+string(EnMax);
+
+   
+  if (KEYWORD_SET(ENERGYMIN)) then begin
+	      graph[GraphN].AxisF =  KEYWORD_SET(log) ?  'logaxis' : 'lineaxis, 0';
+	      ytitle = 'MGS e-Flux (cm!U2!N*s*st*eV)!U-1' + energyTitle;
+	      PlotF = 'spec_comp_plot, '+ options;                
+  endif else begin
+     ytitle = 'MGS e-Flux!C Energy, eV' 
+     PlotF = 'e_mgs_omni_plot'
+  endelse
+
+ graph[GraphN].PlotF = PlotF;
+ if (KEYWORD_SET(AY1TITLE)) then graph[GraphN].Ly.title = AY1Title $
+ else graph[GraphN].Ly.title = yTitle
+ 
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/earthorbaxis.pro
@@ -0,0 +1,121 @@
+;==================================================
+;                 DD  
+;              ORBITAL OBJECTS
+;              mexorbaxis.pro
+;                   V.2.0
+; Modifications:
+;     19 May 2006: V.1.0, Fedorov
+;     09 Dec 2006: V.2.0, Fedorov, New approach for graphical environment
+;---------------------------------------------------
+; Implementation of all types of axis for orbital plot
+; Options:
+;   Key                  ToPlot
+; /CYL - cylindrical XR   0
+; /XY                     1
+; /XZ                     2
+; /YZ                     3
+;====================================================
+
+pro earthorbaxis, ToPlot
+common GraphC, graph, GraphN
+
+;------------ Set graphics Colors ------------------------
+ if (Graph[GraphN].Lx.range[0] eq Graph[GraphN].Lx.range[1]) then begin
+       Graph[GraphN].Lx.range[0] = -20;
+       Graph[GraphN].Lx.range[1] = 10; 
+ endif
+ if (Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1]) then begin
+       Graph[GraphN].Ly.range[0] = ToPlot eq 0 ? 0: -20;
+       Graph[GraphN].Ly.range[1] = 10;
+ endif
+ 
+ !x = Graph[GraphN].Lx
+ !y = Graph[GraphN].Ly
+
+  x = !x.range
+  y = !y.range
+  black = !DNC- 1
+  zero = fix(!DNC * 20.0/ 256.0 )
+  planet = fix(!DNC * 10.0/ 256.0)
+  bound =  fix(!DNC * 1.0/ 256.0)
+;------------------------------------------------------------
+
+  ThickPlanet = 2 
+  ThickBound =  2
+
+;------- main axis -----------------
+  case ToPlot of
+  0: begin
+      !x.title = 'X, R!DE'
+      !y.title = 'R, R!DE'
+     end
+  1: begin
+      !x.title = 'X, R!DE'
+      !y.title = 'Y, R!DE'
+     end
+  2: begin
+      !x.title = 'X, R!DE'
+      !y.title = 'Z, R!DE'
+     end
+  3: begin
+      !x.title = 'Y, R!DE'
+      !y.title = 'Z, R!DE'
+     end
+  endcase
+  
+  ;------ Plot main axis. Use common Graph system variables
+  plot,x,y,/nodata,/noerase,xstyle=1,ystyle=1,color=black,xminor=1, yminor=1, charsize=graph[GraphN].Lp.charsize
+     
+  ;------- plot Zero lines --------------------- 
+  oplot,x,[0,0],color=zero, thick = 1
+  oplot,[0,0],y,color=zero, thick = 1
+  
+  ;------- plot Earth ---------------------
+  dal = 2.0 * !PI / 60.
+  if ToPlot EQ 0 then $ 
+       al = dal * findgen(31) $
+  else al = dal * findgen(61) 
+    
+  xx = cos(al)
+  yy = sin(al)
+  oplot,xx,yy,color=black, thick = ThickPlanet
+    
+  if ToPlot EQ 0 then begin
+     ;---- Plot boundaries for cylindrical plot
+     ;----------------------- Bow Shock -----------------------------
+   RamPress = 1.5  
+   theta = (135.*!pi/180.)*findgen(100)/100.
+   
+   
+  ; rr = 2.04/(1.0+1.02*cos(theta))
+   xx = fltarr(100)
+   yy = fltarr(100)
+   
+   epsilon = 0.81
+   L0 = 24.8
+   ramP = 1.8 
+   L_BS = L0*(RamPress/ramP)^(-1./6.)
+
+    rr=L_BS/(1.+epsilon*COS(theta))  
+    xx = rr *cos(theta)
+    yy = rr*sin(theta)
+ 
+   oplot,xx,yy,color=bound,  thick = ThickBound
+   oplot,xx,-yy,color=bound, thick = ThickBound
+
+    ;----------------------- MP Shue-----------------------------
+    Bz = -1.0
+    
+    alpha=(.58-.007*Bz)*(1.+0.024*alog(RamPress))
+    r0=(10.22+1.29*tanh(0.184*(Bz+8.14)))*RamPress^(-1./6.6)
+    rr=r0*(2./(1.+COS(theta)))^alpha
+    
+     xx = rr*cos(theta)
+     yy = rr*sin(theta)
+   
+     oplot,xx,yy,color=bound,  thick = ThickBound
+     oplot,xx,-yy,color=bound, thick = ThickBound
+     
+  endif
+return
+end
+++ a/amda_plus/amdalib/elscntplot.pro
@@ -0,0 +1,133 @@
+;==========================================================
+;                DD DDMARS/DDVEX
+;                 elscntplot.pro
+;    
+;-----------------------------------------------------------
+
+pro elscntplot, Mission, anodStart, anodStop, enMin, enMax
+
+common VexELS, VexELS, VexElsPlus
+common MexELS, MexELS, MexElsPlus
+common GraphC, Graph,GraphN, Item
+common request, List, StartTime, TimeInt
+   
+   cmd = 'ElsPlus = ' + Mission + 'ElsPlus'
+   R = execute(cmd);  
+   
+ ; Expand DeltaT a little : Get Energy Calibration
+ 
+  StD = call_external(!DDLIB+'DD_idl.so','GetDTime', StartTime, /D_VALUE)
+  DtD = call_external(!DDLIB+'DD_idl.so','GetDTime', TimeInt, /D_VALUE)
+  StD = StD - 100.0D0
+  DtD = DtD + 200.0D0
+  stloc = call_external(!DDLIB+'DD_idl.so','GetTime',StD, /S_VALUE)
+  dtloc = call_external(!DDLIB+'DD_idl.so','GetTime',DtD, /S_VALUE)
+  Size = Long(DtD/4. + 2)
+  Nen = 128;
+  Energy = fltarr(Nen, Size)
+  Tm = dblarr(Size)
+  ID = -1L;
+  
+  err = call_external(!PROJLIB + 'libplanetcom.so','SetMission_idl', Mission);
+  RetSize = call_external(!PROJLIB + 'els_idl.so','GetRef',   $
+                                                 ID,    $
+                                                 stloc, $
+                                                 dtloc, $
+                                                 Size,  $
+                                                 Tm,    $
+                                                 Energy)
+  if (err GE 0) then Tm  -= Graph[GraphN].TimeBase
+  err = call_external(!DDLIB +'DD_idl.so','CloseID',ID)
+
+; data 
+ 
+   Time = (*Graph[GraphN].Time[0]); 
+   N_Time = N_elements(Time);
+
+   Val = (*Graph[GraphN].Val[0]);
+   Val = temporary(reform(Val, 16, Nen, N_Time)); 
+   
+   num = where(abs(Tm - Time[N_Time/2]) EQ min(abs(Tm - Time[N_Time/2])));
+
+   E = Energy[*,num[0]]*ElsPlus.Sens[7];
+ 
+   if ((enmin eq enmax) and (enmin eq 0)) then begin      
+            enmax = E[0];
+            enmin = E[126];       
+   endif
+       
+   Data =  anodStop GT anodStart ? total(Val[anodStart:anodStop,*,*],1) : reform(Val[anodStart,*,*], 128, N_Time);
+
+; if reduced Energy Table  
+
+   reducedEnergy = E[0] LT 1000.0 ? 1 : 0; 
+  ; if reduced Energy Table expand time
+     if (reducedEnergy) then begin
+	Nen = 32;
+	enMax = enMax < 158.0; 
+        newTime = dblarr(N_Time*4);
+        for i = 0, N_Time-1 do for j = 0, 3 do newTime[i*4+j] = Time[i] + double(j);
+        Time = newTime; 
+        Data = reform(Data, Nen, N_Time*4);
+        N_Time *= 4;    
+        E = E(0:Nen-1);                   
+     endif  
+     
+    E[Nen-1] = E[Nen-2]; 
+    EBounds = fltarr(Nen+1);
+    EBounds[0] = E[0];
+    EBounds[Nen] = E[Nen-1]; 
+    for ie = 1, Nen-1 do EBounds(ie) = sqrt(E(ie-1)*E(ie)); 
+   
+    num0 = where(EBounds ge enMin, Index0);
+    num1 = where(EBounds ge enMax, Index1);
+ 
+    if (Index1 eq Nen+1) then begin
+	  nodatafortimeinterval
+	  return
+    endif
+
+    if (Index1 eq 0)  then numStop = 0 else  numStop = num1[Index1-1];
+
+    if (Index0 eq Nen+1) then numStart = 31 else numStart = num0[Index0-1];
+
+    Data =  numStart ne numStop ? total(Data[numStop:numStart,*],1) : reform(Data[numStart,*], N_Time);
+ ;   Data = IndexE gt 1 ? total(Data[numE,*],1) : reform(Data[numE[0],*], N_Time);
+
+  ;  if (avEnergy) then Data /= IndexE;
+  ;  if (avAngle) then Data /= (anodStop-anodStart+1); 
+  
+    numMinMax = where(finite(Data), IndexMinMax)
+ 
+  ;--------------------- No DATA - just axes -----------------------
+ 
+  if (IndexMinMax LE 0) then begin
+        nodatafortimeinterval  
+        return
+  endif
+
+  if (Time[0] GT Graph[GraphN].TotalTime) then begin
+        nodatafortimeinterval
+        return
+  endif
+ ;-------------------------------------------------------------------
+  
+  ;TODO test mode   if E(1) EQ E(126) then begin ; The test mode, just plot everything???         
+    MinVal = min(Data[numMinMax]);   
+    MaxVal = max(Data[numMinMax]);
+   
+    if Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] then begin
+                Graph[GraphN].Ly.range = Graph[GraphN].AxisF eq 'logaxis' ? [MinVal, MaxVal] : [0.0, MaxVal]
+                Graph[GraphN].Ly.style = 0
+    endif else Graph[GraphN].Ly.style = 1
+    
+    !y = Graph[GraphN].Ly;
+    
+     
+    if (execute(Graph[GraphN].AxisF) NE 1) then plot, Time, Data, MAX_VALUE=MaxValue,  color = !dnc-1, /NODATA, /NOERASE 
+    oplot, Time, Data,  MAX_VALUE=MaxValue, color = !dnc-1 
+
+        
+return
+end  
+;--------------------------------------------------------------------
+++ a/amda_plus/amdalib/elsspecplot.pro
@@ -0,0 +1,160 @@
+;==========================================================
+;                DD DDMARS/DDVEX
+;                 elsetplot.pro
+;                     V.3.2
+;     Constructor of ELSET object
+;     Plot of the ET spectrogram of ELS count.
+; Mandatory arguments: FirstAnode, LastAnode
+; Keywords:            no keywords
+; Resource control the color scale
+;
+; Versions:
+;      28 Jan 2004: V.1.0
+;      18 han 2005: V.2.0  Fedorov, variable energy range
+;      30 May 2007: V.3.2, Fedorov. Anodes interval. Blue background
+;-----------------------------------------------------------
+
+pro elsspecplot, Mission, anodStart, anodStop 
+
+common VexELS, VexELS, VexElsPlus
+common MexELS, MexELS, MexElsPlus
+common GraphC, Graph,GraphN, Item
+common request, List, StartTime, TimeInt
+   
+   cmd = 'ElsPlus = ' + Mission + 'ElsPlus'
+   R = execute(cmd);  
+   
+ ; Expand DeltaT a little : Get Energy Calibration
+ 
+  StD = call_external(!DDLIB+'DD_idl.so','GetDTime', StartTime, /D_VALUE)
+  DtD = call_external(!DDLIB+'DD_idl.so','GetDTime', TimeInt, /D_VALUE)
+  StD = StD - 100.0D0
+  DtD = DtD + 200.0D0
+  stloc = call_external(!DDLIB+'DD_idl.so','GetTime',StD, /S_VALUE)
+  dtloc = call_external(!DDLIB+'DD_idl.so','GetTime',DtD, /S_VALUE)
+  Size = Long(DtD/4. + 2)
+  Energy = fltarr(128, Size)
+  Tm = dblarr(Size)
+  ID = -1L;
+  
+  err = call_external(!PROJLIB + 'libplanetcom.so','SetMission_idl', Mission);
+  RetSize = call_external(!PROJLIB + 'els_idl.so','GetRef',   $
+                                                 ID,    $
+                                                 stloc, $
+                                                 dtloc, $
+                                                 Size,  $
+                                                 Tm,    $
+                                                 Energy)
+
+ 
+  if (err GE 0) then Tm  -= Graph[GraphN].TimeBase
+  err = call_external(!DDLIB +'DD_idl.so','CloseID',ID)
+  
+   Time = (*Graph[GraphN].Time[0])
+   N_Time = N_elements(Time);
+       
+   EnumbFull = 128
+   EnumbReduced = 31
+
+   E = fltarr(EnumbFull)
+
+  ;-------- Cat in the Energy range ------------
+   E0 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ? 0.610171 : Graph[GraphN].Ly.range[0]
+   E1 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ? 20832.0: Graph[GraphN].Ly.range[1]
+   
+    Val = (*Graph[GraphN].Val[0])   
+    Val = temporary(reform(Val, 16, 128, N_Time)) 
+    LocVal =  anodStop GT anodStart ? total(Val[anodStart:anodStop,*,*],1) : reform(Val[anodStart,*,*], 128, N_Time);
+ 
+    numMinMax = where(finite(LocVal), IndexMinMax)
+ 
+   if (Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1]) then  $
+                                       Graph[GraphN].Ly.range = [E0, E1]
+         
+  Graph[GraphN].Ly.style = 1
+
+  !y = Graph[GraphN].Ly
+
+   plot_io, graph[GraphN].Lx.range, graph[GraphN].Ly.range, /nodata, /noerase, ystyle = 1, $
+                              ytickformat='yticks', charsize=graph[GraphN].Lp.charsize*0.8, color = !dnc-1
+                                 
+  ;--------------------- No DATA - just axes -----------------------
+ 
+  if (IndexMinMax LE 0) then begin
+        nodatafortimeinterval  
+         return
+  endif
+  if (Time[0] GT Graph[GraphN].TotalTime) then begin
+        nodatafortimeinterval
+         return
+  endif
+  MinVal = min(LocVal[numMinMax]) > 1.0
+  MaxVal = max(LocVal[numMinMax])                            
+
+    Step = Graph[GraphN].Sampling[0] > Graph[GraphN].Step;
+    StepMax = Graph[GraphN].MaxSampling[0] > Step
+ ;-----------   NOW PLOT SPECTRA
+    alogMin = alog(MinVal)
+    KF = float(!DNC - 2) /(alog(MaxVal) - alogMin)
+    xx = fltarr(4) & yy = fltarr(4)
+ 
+     for i = 0, N_Time - 2 do begin
+      
+        num = where(abs(Tm - Time[i]) EQ min(abs(Tm - Time[i])));
+        E =  anodStop GT anodStart ? Energy[*,num[0]]*ElsPlus.Sens[4] : Energy[*,num[0]]*ElsPlus.Sens[anodStart];
+        FullEnergy = 1;  
+        if (mission eq "VEX") then FullEnergy = E[0] GT 1000.0 ? 1 : 0; 
+ 
+         Enumb = FullEnergy ? EnumbFull :EnumbReduced
+
+        eindex = where((E[0:Enumb-1] GE E0) AND $
+			  (E[0:Enumb-1] LE E1) AND $
+			(E[0:Enumb-1] GT 0.0),counte)
+  
+        if (FullEnergy) then E[127] = E[126] 	    		    
+	    if (counte GT 0) AND (E[1] NE E[126]) then begin
+		EBounds = fltarr(ENumb+1)
+		EBounds(min(eindex)) = E(min(eindex))
+		EBounds(max(eindex+1)) = E(max(eindex))
+		for ie = min(eindex+1), max(eindex) do EBounds(ie) = sqrt(E(ie-1)*E(ie))
+	      endif else begin
+	      if E(1) EQ E(126) then begin ; The test mode, just plot everything
+		  eindex = indgen(128)
+		  step = exp(alog(E1/E0)/127.0)
+		  E = E0 * step^findgen(128)
+		  EBounds = fltarr(ENumb+1)
+		  EBounds(min(eindex)) = E(min(eindex))
+		  EBounds(max(eindex+1)) = E(max(eindex))
+		  for ie = min(eindex+1), max(eindex) do EBounds(ie) = sqrt(E(ie-1)*E(ie))
+	      endif else ERRFLAG = 1   
+	    endelse  
+        if (not FullEnergy) then for kk = 0, 30 do for kkk = 1, 3 do LocVal[kk,i] += LocVal[kk+kkk*32,i]/4.0
+
+        xx[0] = Time[i] > 0.0
+        xx[1] = xx[0]
+        xx[2] = Time[i+1] < Time[i] + StepMax
+        if (xx[2] GT graph[GraphN].TotalTime) then xx[2] = graph[GraphN].TotalTime
+        xx[3]= xx[2]
+      
+        if (xx[0] LT graph[GraphN].TotalTime AND xx[2] LE graph[GraphN].TotalTime) then begin
+           for k = min(eindex), max(eindex) do begin
+                  yy[0] = EBounds[k]
+                  yy[1] = EBounds[k+1]
+                  yy[2]=yy[1]
+                  yy[3]=yy[0]
+                  lcount = LocVal[k,i]
+                  cl = finite(lcount) ? fix(KF*(alog(lcount) - alogMin)) > 1 : 0;
+                  if (cl GT !DNC - 2) then  cl = !DNC - 2
+                  PolyFill, xx ,yy, color = cl, /data 
+            endfor          
+        endif
+     endfor
+
+;-------------------- Legend ----------------------------------------
+      spec_legend, MinVal, MaxVal
+;------------------ MIN-MAX------------------------------------------------
+      Graph[GraphN].Min  = MinVal
+ 
+return
+end  
+;--------------------------------------------------------------------
+++ a/amda_plus/amdalib/get_ica_spec.pro
@@ -0,0 +1,46 @@
+;===========================================================
+;=============================================================
+
+pro Get_Ica_Spec, StartTime, TimeInt 
+
+common IcaSpecCom, Spec
+common GraphC, Graph, GraphN, Item
+common Etab, Flag, Tab
+
+;------------Get data-----------------------------
+
+    ID = Long(graph[GraphN].DataID[Item])
+
+    Size = Spec.DefSize
+    TimeArr = dblarr(Size);
+    Spectra = fltarr(96, 16, 16, Size)
+;    PaccArr = lonarr(Size);
+    EnergyTab = intarr(Size);
+    VI = "ros:ica:all";
+    paramName = "no_counts_summed";
+
+    RetSize = call_external(!PROJLIB + '/libParam.so', 'getLocalData', $   
+                                 ID, VI, paramName, $
+                                 StartTime, TimeInt, Size, $
+                                 TimeArr, Spectra)
+  
+    RetSize = call_external(!PROJLIB + '/libParam.so', 'getLocalData', $   
+                                 ID, VI, "energyTab", $
+                                 StartTime, TimeInt, Size, $
+                                 TimeArr, EnergyTab)
+
+     if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+     if (RetSize GT 0) then begin
+         Spectra = temporary(reform(Spectra[*,*,*,0:RetSize-1], 96L*16L*16L, RetSize));
+         TimeArr -= graph[GraphN].TimeBase
+         graph[GraphN].DataFlags[Item] = 1
+         graph[GraphN].dataTime[Item] = TimeArr[0]
+         average, Spec, TimeArr[0:RetSize-1], Spectra
+         if (Flag eq GraphN) then immediate, Tab, TimeArr[0:RetSize-1], float(EnergyTab[0:RetSize-1]);
+     endif; else Graph[GraphN].ErrorFlag = 1
+   
+ 
+   
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/get_ima_o_spec.pro
@@ -0,0 +1,56 @@
+;===========================================================
+;       Get_Vex_H_Spec
+;
+;  18 Mars 2009
+;=============================================================
+
+pro Get_Ima_O_Spec, Mission, StartTime, TimeInt, RESTRICTED = restricted
+
+common VexOSpecCom, VexOSpec
+common MexOSpecCom, MexOSpec
+common GraphC, Graph, GraphN, Item
+common Etab, FlagVex, EtabVex, FlagMex, EtabMex 
+ 
+
+  cmd = 'Spec = ' + Mission + 'OSpec & Tab = Etab' + Mission;
+  R = execute(cmd);
+
+  Flag = Mission eq 'vex' ? FlagVex : FlagMex;
+  
+;------------Get data-----------------------------
+
+    ID = Long(graph[GraphN].DataID[Item])
+
+    Size = Spec.DefSize
+    TimeArr = dblarr(Size);
+    Spectra = fltarr(16, 96, 16, Size)
+    Product = 'HEAVYSPEC';
+    PaccArr = intarr(Size);
+    EnergyTab = intarr(Size);
+
+    if (ID eq -1L) then begin 
+           err = call_external(!PROJLIB + 'libplanetcom.so','SetMission_idl', Mission);
+           if (KEYWORD_SET(RESTRICTED)) then err = call_external(!PROJLIB + 'libplanetcom.so','SetRestriction');
+    endif
+    RetSize = call_external(!PROJLIB + 'libImaExtra.so','GetImaExtra_idl',$
+                                ID, StartTime, TimeInt,$
+                                Product, Size, $
+                                TimeArr, PaccArr, EnergyTab, Spectra)
+ 
+   
+     if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+     if (RetSize GT 1) then begin
+         Spectra = temporary(reform(Spectra[*,*,*,0:RetSize-1], 16L*96L*16L,RetSize));
+         TimeArr -= graph[GraphN].TimeBase
+         graph[GraphN].DataFlags[Item] = 1
+         graph[GraphN].dataTime[Item] = TimeArr[0]
+         average, Spec, TimeArr[0:RetSize-1], Spectra
+         if (Flag eq GraphN) then immediate, Tab, TimeArr[0:RetSize-1], float(EnergyTab[0:RetSize-1]);
+     endif; else Graph[GraphN].ErrorFlag = 1
+   
+   cmd = Mission + 'OSpec = Spec & Etab' + Mission +'= Tab'
+   R = execute(cmd);
+   
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/get_imam_h_spec.pro
@@ -0,0 +1,78 @@
+;===========================================================
+;       Get_Vex_H_Spec
+;
+;  18 Mars 2009
+;=============================================================
+
+pro Get_ImaM_H_Spec, StartTime, TimeInt, RESTRICTED = restricted
+
+common MexHSpecCom, Spec
+common InternalImaMex, Energy
+common GraphC, Graph, GraphN, Item
+common Etab, FlagVex, EtabVex, FlagMex, EtabMex
+
+FORWARD_FUNCTION gfl
+;------------Get data-----------------------------
+
+    ID = Long(graph[GraphN].DataID[Item])
+
+    Size =    Spec.DefSize
+    TimeArr = dblarr(Size);
+    SpectraTot = fltarr(16, 96, 16, 2, Size)
+    Product = 'HPSPEC';
+    PaccArr = intarr(Size);
+    EnergyTab = intarr(Size);
+
+    Mission = "mex";
+    
+    if (ID eq -1L) then begin
+         err = call_external(!PROJLIB + 'libplanetcom.so','SetMission_idl', Mission);
+         if (KEYWORD_SET(RESTRICTED)) then err = call_external(!PROJLIB + 'libplanetcom.so','SetRestriction');
+    endif
+
+    EnNum = 96
+    Energy = fltarr(EnNum)
+       
+    RetSize = call_external(!PROJLIB + 'libImaExtra.so','GetImaExtra_idl',$
+                                ID, StartTime, TimeInt,$
+                                Product, Size, $
+                                TimeArr, PaccArr, EnergyTab, SpectraTot)
+ 
+     if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+ 
+     err = call_external(!PROJLIB + 'imainfo.so','GetImaEner_idl', Energy, EnergyTab[0])
+     Energy = temporary(Energy(where(Energy GT 0))); 
+ 
+     if (RetSize GT 1) then begin
+ 
+         Spectra = fltarr(16, 96, 16, RetSize);
+   
+         for i = 0, RetSize - 1 do begin
+
+            if (i ne RetSize - 1) then $ 
+	      if (EnergyTab[i] ne EnergyTab[i+1]) then begin
+		err = call_external(!PROJLIB + 'imainfo.so','GetImaEner_idl', Energy, EnergyTab[i+1])
+		Energy = temporary(Energy(where(Energy GT 0))); 
+	    endif
+
+            for j = 0, N_elements(Energy) - 1 do begin
+               GFeff = gfl(Energy[j], PaccArr[i], 1.0);
+               GFeffG = gfl(Energy[j], PaccArr[i], 0.0); Ghost
+                 
+               coeffH = GFeff GT 0 ? 1.e-5/GFeff : 0.0;
+               coeffG = GFeffG GT 0 ? 1.e-5/GFeffG : 0.0;
+             
+               Spectra[*,j,*,i] = SpectraTot[*,j,*,0,i]*coeffH + SpectraTot[*,j,*,1,i] * coeffG;
+         endfor 
+      endfor
+
+         TimeArr -= graph[GraphN].TimeBase
+         graph[GraphN].DataFlags[Item] = 1
+         graph[GraphN].dataTime[Item] = TimeArr[0]
+         average, Spec, TimeArr[0:RetSize-1], reform(Spectra,16L*96L*16L,RetSize)
+	 if (FlagMex eq GraphN) then immediate, EtabMex, TimeArr[0:RetSize-1], float(EnergyTab[0:RetSize-1]);
+     endif else Graph[GraphN].ErrorFlag = 1
+ 
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/get_mode_peir.pro
@@ -0,0 +1,48 @@
+;-------------------- Get THEMIS Modeity ION -------------------------------
+;  
+;-----------------------------------------------------------------------------
+pro Get_Mode_Peir, THEMIS, st, dt   
+
+common GraphC, Graph, GraphN, Item
+common ThaModePeirCom, ThaMode
+common ThbModePeirCom, ThbMode
+common ThcModePeirCom, ThcMode
+common ThdModePeirCom, ThdMode
+common TheModePeirCom, TheMode 
+ 
+    cmd = 'ThMode = Th'+THEMIS+'Mode'   
+    R = execute(cmd) 
+      
+    sz = ThMode.DefSize
+    T=dblarr(sz)
+    mode = fltarr(3,sz);
+
+    ID = Long(graph[GraphN].DataID[Item])
+
+    VI = 'th'+THEMIS+':peir:mom'
+CONT:
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, 'Mode', st, dt, sz, T, mode)
+ 
+    if (Retsize eq -34) OR (Retsize eq -35) OR (Retsize eq -0) then begin  
+                   wait, 3               
+                   goto, CONT
+     endif 
+
+   
+  if (RetSize GT 1) then begin 
+ 
+      Val =  mode[2, 0:Retsize-1]
+      T -= graph[GraphN].TimeBase
+      graph[GraphN].DataID[Item] = ID
+      graph[GraphN].DataFlags[Item] = 1
+      graph[GraphN].dataTime[Item] = T[0] 
+      immediate, ThMode, T[0:RetSize-1], Val
+   endif
+
+  cmd = 'Th'+THEMIS+'Mode = ThMode' 
+  R = execute(cmd)
+       
+return
+end
+;-----------------------------------------------------------
+++ a/amda_plus/amdalib/get_multi.pro
@@ -0,0 +1,86 @@
+;;
+;---------------------------------------
+
+pro GET_MULTI, st, dt
+
+common Multi_Commom, getStart, N_VIs,  DataCommon
+common GraphC, graph, GraphN, Item 
+ 
+   FillValue = !Values.F_NAN 
+   
+   StartTime = call_external(!DDLIB +'DD_idl.so','GetDTime', st, /D_VALUE)	      
+   TimeCalc = dindgen(Graph[GraphN].ppp)*Graph[GraphN].Step + StartTime + Graph[GraphN].Step/2.0; 
+   k = 0;
+
+   for i = 0, N_VIs[graphN] - 1 do begin 
+
+     if(Graph[GraphN].CurrentSec EQ Graph[GraphN].Sections-1) then $
+         *(DataCommon[i, GraphN].data) = Call_Function(DataCommon[i, GraphN].name, st, dt, Time, RetSize, Sampling, /close) $
+     else *(DataCommon[i, GraphN].data) = Call_Function(DataCommon[i, GraphN].name, st, dt, Time, RetSize, Sampling)
+ 
+     if (RetSize GT 0) then begin
+    ;
+    ; Different processing if 1D array and RetSize = 1
+    ;
+	  Yinfo = size(*DataCommon[i, GraphN].data);
+	  if  (Yinfo[0] eq 1) then begin
+	    if (RetSize gt 1) then begin 
+		(*DataCommon[i, GraphN].data) = reform((*DataCommon[i, GraphN].data),1, RetSize);
+		  DataCommon[i, GraphN].size[1] = 1;
+	    endif else (*DataCommon[i, GraphN].data) = reform((*DataCommon[i, GraphN].data), Yinfo[1], RetSize);
+	  endif else DataCommon[i, GraphN].size[1] = Yinfo[1]
+ 
+	  *(DataCommon[i, GraphN].Time) = Time;
+ 
+; Several data sections 
+      if (Graph[GraphN].Sections gt 1) then begin 	 
+	  if (N_elements(*(DataCommon[i, GraphN].LastTime)) GT 0) then $
+	      if (*(DataCommon[i, GraphN].LastTime) LT Time[0]) then begin
+		*(DataCommon[i, GraphN].Time) = [*(DataCommon[i, GraphN].LastTime),Time]
+		TempArr = fltarr(DataCommon[i, GraphN].size[1],Retsize+1)
+		TempArr[*,0] = *(DataCommon[i, GraphN].LastData)
+		TempArr[*,1:*] = *(DataCommon[i, GraphN].data)
+		*(DataCommon[i, GraphN].data) = TempArr
+		*(DataCommon[i, GraphN].LastTime) = Time[RetSize-1];
+		*(DataCommon[i, GraphN].LastData) = (*DataCommon[i, GraphN].data)[*,RetSize-1];
+		  RetSize += 1;
+	    endif else begin
+	    *(DataCommon[i, GraphN].LastTime) = Time[RetSize-1];
+	    *(DataCommon[i, GraphN].LastData) = (*DataCommon[i, GraphN].data)[*,RetSize-1];
+	    endelse
+	  endif ;  if (Graph[GraphN].Sections gt 1)
+	  TempArr = 0.0  
+      endif ; if (RetSize GT 0)
+
+	DataCommon[i, GraphN].size[0] = RetSize;
+	DataCommon[i, GraphN].Sampling = Sampling;	           
+ 
+      if (DataCommon[i, GraphN].Sampling LT Graph[GraphN].Step/2.0) then begin		    
+	      Y = reform(replicate(FillValue, DataCommon[i, GraphN].size[1]*Graph[GraphN].ppp), DataCommon[i, GraphN].size[1], Graph[GraphN].ppp);
+	      Flag = intarr(N_elements(TimeCalc));   
+	      averaging_, *(DataCommon[i, GraphN].time), *(DataCommon[i, GraphN].data), TimeCalc, Y, Flag;
+	      *DataCommon[i, GraphN].data = Y;
+	      *(DataCommon[i, GraphN].Time) = TimeCalc;
+      endif  
+    
+      indices = *(DataCommon[i, GraphN].indices);
+ 
+      if (N_Elements(indices) gt 0) then begin
+	for j = 0, N_Elements(indices) - 1 do begin 	 
+		*(Graph[GraphN].Val)[k] =  getStart[GraphN] ? [*(Graph[GraphN].Val)[k],reform((*DataCommon[i, GraphN].data)[indices[j],*])] : reform((*DataCommon[i, GraphN].data)[indices[j],*]);
+                *(Graph[GraphN].Time)[k] = getStart[GraphN] ? [*(Graph[GraphN].Time)[k], *(DataCommon[i, GraphN].Time)-Graph[GraphN].TimeBase] : *(DataCommon[i, GraphN].Time)-Graph[GraphN].TimeBase;
+	         k += 1;
+      endfor
+      endif else begin
+	*(Graph[GraphN].Val)[k] =  getStart[GraphN] ? [*(Graph[GraphN].Val)[k],*(DataCommon[i, GraphN].data)] : *(DataCommon[i, GraphN].data);
+	*(Graph[GraphN].Time)[k] = getStart[GraphN] ? [*(Graph[GraphN].Time)[k], *(DataCommon[i, GraphN].Time)-Graph[GraphN].TimeBase] : *(DataCommon[i, GraphN].Time)-Graph[GraphN].TimeBase;
+        k += 1;
+      endelse
+ 
+ endfor
+
+  Graph[GraphN].DataFlags = 1
+  Graph[GraphN].ErrorFlag = 0
+  getStart[GraphN] = 1;
+return
+end
+++ a/amda_plus/amdalib/get_scatter.pro
@@ -0,0 +1,108 @@
+;;
+;---------------------------------------
+
+pro GET_SCATTER, st, dt
+
+common Scatter_Commom, getStart, N_param, indices, DataCommon
+common GraphC, graph, GraphN, Item 
+ 
+   FillValue = !Values.F_NAN 
+ 
+   for i = 0, N_param[graphN] - 1 do begin 
+     if(Graph[GraphN].CurrentSec EQ Graph[GraphN].Sections-1) then $
+         *(DataCommon[i, GraphN].data) = Call_Function(DataCommon[i, GraphN].name, st, dt, Time, RetSize, Sampling, /close) $
+     else *(DataCommon[i, GraphN].data) = Call_Function(DataCommon[i, GraphN].name, st, dt, Time, RetSize, Sampling)
+  
+    if (RetSize GT 0) then begin
+    ;
+    ; Different processing if 1D array and RetSize = 1
+    ;
+ 
+          Yinfo = size(*DataCommon[i, GraphN].data);
+          if  (Yinfo[0] eq 1) then begin
+            if (RetSize gt 1) then begin 
+                (*DataCommon[i, GraphN].data) = reform((*DataCommon[i, GraphN].data), 1, RetSize);
+                  DataCommon[i, GraphN].size[1] = 1;
+            endif else (*DataCommon[i, GraphN].data) = reform((*DataCommon[i, GraphN].data), Yinfo[1], RetSize);
+          endif else DataCommon[i, GraphN].size[1] = Yinfo[1]
+
+          *(DataCommon[i, GraphN].Time) = Time;
+ 
+; Several data sections 
+      if (Graph[GraphN].Sections gt 1) then begin       
+          if (N_elements(*(DataCommon[i, GraphN].LastTime)) GT 0) then $
+              if (*(DataCommon[i, GraphN].LastTime) LT Time[0]) then begin
+                *(DataCommon[i, GraphN].Time) = [*(DataCommon[i, GraphN].LastTime),Time]
+                TempArr = fltarr(DataCommon[i, GraphN].size[1],Retsize+1)
+                TempArr[*,0] = *(DataCommon[i, GraphN].LastData)
+                TempArr[*,1:*] = *(DataCommon[i, GraphN].data)
+                *(DataCommon[i, GraphN].data) = TempArr
+                *(DataCommon[i, GraphN].LastTime) = Time[RetSize-1];
+                *(DataCommon[i, GraphN].LastData) = (*DataCommon[i, GraphN].data)[*,RetSize-1];
+                  RetSize += 1;
+            endif else begin
+            *(DataCommon[i, GraphN].LastTime) = Time[RetSize-1];
+            *(DataCommon[i, GraphN].LastData) = (*DataCommon[i, GraphN].data)[*,RetSize-1];
+            endelse
+          endif else if (Yinfo[0] gt 1) then DataCommon[i, GraphN].size[1] = Yinfo[1]
+          TempArr = 0.0  
+      endif
+        DataCommon[i, GraphN].size[0] = RetSize;
+        DataCommon[i, GraphN].Sampling = Sampling;               
+    endfor
+ 
+; Processing Y 
+  if (DataCommon[0, GraphN].Sampling LT Graph[GraphN].Step/2.0) then begin
+             StartTime = call_external(!DDLIB+'DD_idl.so','GetDTime',st, /D_VALUE)
+             TimeInterval = call_external(!DDLIB+'DD_idl.so','GetDTime',dt, /D_VALUE)
+              numberPoints = long(TimeInterval/Graph[GraphN].Step + 0.5);
+	      toAverageY = 1;
+	      Time = dindgen(numberPoints)*Graph[GraphN].Step + StartTime+ Graph[GraphN].Step/2.0; 		    
+	      Y = reform(replicate(FillValue, DataCommon[0, GraphN].size[1]*numberPoints), DataCommon[0, GraphN].size[1],numberPoints);     
+              Flag = intarr(N_elements(Time));   
+	      averaging_, *(DataCommon[0, GraphN].time), *(DataCommon[0, GraphN].data), Time, Y, Flag;
+	      if (getStart[graphN]) then begin 
+                   *Graph[GraphN].Val[1] = [*Graph[GraphN].Val[1], reform(Y[indices[0,graphN],*],numberPoints)]                      
+              endif  else begin
+	                *Graph[GraphN].Val[1] =  reform(Y[indices[0,graphN],*],numberPoints);
+              endelse 
+  endif else begin 
+	      toAverageY = 0;
+	      Time = *DataCommon[0, GraphN].time;
+	      if (getStart[graphN]) then *Graph[GraphN].Val[1] = [*Graph[GraphN].Val[1], reform((*DataCommon[0, GraphN].data)[indices[0],*],DataCommon[0, GraphN].size[0])] $
+	      else *Graph[GraphN].Val[1] =  reform((*DataCommon[0, GraphN].data)[indices[0,graphN],*], DataCommon[0, GraphN].size[0]);
+  endelse
+
+; Processing X 
+; the same parameter - different component
+   if (N_param[graphN] eq 1) then begin
+	      if  (toAverageY) then X = reform(Y[indices[1,graphN],*],numberPoints) $
+	      else X = reform((*DataCommon[0, GraphN].data)[indices[1,graphN],*],DataCommon[0, GraphN].size[0]);
+	      if (getStart[graphN]) then *Graph[GraphN].Val[0] = [*Graph[GraphN].Val[0], X] $
+	      else *Graph[GraphN].Val[0] = X;
+   endif else begin
+	      if (toAverageY) then  step = Graph[GraphN].Step $
+	      else step = DataCommon[0, GraphN].Sampling; 
+ 
+	      X = reform(replicate(FillValue, DataCommon[1, GraphN].size[1]*N_Elements(Time)), DataCommon[1, GraphN].size[1], N_Elements(Time));
+	      Flag = intarr(N_elements(Time)); 
+	      if (DataCommon[1, GraphN].Sampling LT step/2.0) then $                 
+		 averaging_, *(DataCommon[1, GraphN].time), (*DataCommon[1, GraphN].data)[indices[1,graphN],*], Time, X, Flag $
+	      else $
+		  interpolation_, *(DataCommon[1, GraphN].time), (*DataCommon[1, GraphN].data)[indices[1,graphN],*], DataCommon[1, GraphN].Sampling*3, Time, X, Flag;
+	       if (getStart[graphN]) then *Graph[GraphN].Val[0] = [*Graph[GraphN].Val[0], reform(X)] $
+	       else *Graph[GraphN].Val[0] = reform(X);
+  endelse
+ 
+  if (toAverageY) then begin
+    if (Graph[GraphN].CurrentSec EQ 0) then *Graph[GraphN].Time[0] =  dindgen(Graph[GraphN].ppp)*Graph[GraphN].Step + StartTime + Graph[GraphN].Step/2.0;  
+  endif else begin 
+   if (getStart[graphN]) then  *Graph[GraphN].Time[0] = [*Graph[GraphN].Time[0], Time] $
+   else *Graph[GraphN].Time[0] = Time
+  endelse
+
+  Graph[GraphN].DataFlags = 1
+  Graph[GraphN].ErrorFlag = 0
+  getStart[GraphN] = 1;
+return
+end
+++ a/amda_plus/amdalib/get_th_e_sp.pro
@@ -0,0 +1,42 @@
+;-------------------- Get THEMIS SPECTRA ION -------------------------------
+;  
+;-----------------------------------------------------------------------------
+pro Get_Th_E_Sp, THEMIS, st, dt   
+
+common GraphC, Graph, GraphN, Item
+common ThaeSpecCom, ThaeSpec
+common ThbeSpecCom, ThbeSpec
+common ThceSpecCom, ThceSpec
+common ThdeSpecCom, ThdeSpec
+common TheeSpecCom, TheeSpec 
+ 
+    cmd = 'TheSpec = Th'+THEMIS+'eSpec'   
+    R = execute(cmd) 
+      
+    sz = TheSpec.DefSize
+    V=dblarr(32,sz)
+    T=dblarr(sz)
+ 
+    ID = Long(graph[GraphN].DataID[Item])
+
+    VI = 'th'+THEMIS+':esa:l2e'
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, 'Flux', st, dt, sz, T, V);
+ 
+  
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+   if (RetSize GT 1) then begin 
+      Val = float(V[*,0:RetSize-1])    
+      T -= graph[GraphN].TimeBase
+      graph[GraphN].DataFlags[Item] = 1
+      graph[GraphN].dataTime[Item] = T[0] 
+      average, TheSpec, T[0:RetSize-1], Val
+   endif else Graph[GraphN].ErrorFlag = 1
+ 
+  cmd = 'Th'+THEMIS+'eSpec = TheSpec' 
+  R = execute(cmd)
+       
+return
+end
+;-----------------------------------------------------------
+++ a/amda_plus/amdalib/get_vexmex_dens.pro
@@ -0,0 +1,58 @@
+;-------------------- Get THEMIS ION Velocity -------------------------------
+;  
+;-----------------------------------------------------------------------------
+pro Get_VexMex_Dens, Mission, Product, st, dt, RESTRICTED = restricted
+
+common GraphC, Graph, GraphN, Item
+common VexHDensCom, VexH
+common VexODensCom, VexO
+common MexHDensCom, MexH
+common MexODensCom, MexO
+
+    cmd = 'VEX = ' + Mission + Product; 
+    R = execute(cmd) 
+
+    sz = Vex.DefSize
+     
+    Quality = fltarr(sz)
+    Val = fltarr(sz)
+    Time = dblarr(sz)
+
+    ID = Long(graph[GraphN].DataID[Item])
+
+    VI = Mission + ':ima:param'
+
+    if (KEYWORD_SET(RESTRICTED)) then  VI = Mission + ':imaparam:restricted'
+    
+    ParamNameQ = Product eq "h" ? "Quality_P" : "Quality_O";
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamNameQ, st, dt, sz, Time, Quality)
+ 
+    ParamName = Product eq "h" ? "Density_P" : "Density_O";
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, st, dt, sz, Time, Val)
+ 
+ 
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+   if (RetSize GT 1) then begin 
+      stD = call_external(!DDLIB + 'DD_idl.so','GetDTime', st, /D_VALUE) 
+      numTime = where(Time GE stD AND Time LE (stD+graph[graphN].DeltaT), IndexTime)
+      if IndexTime EQ 0 then return
+      Time = temporary(Time[numTime])
+      Val = temporary(Val[numTime])
+      Quality = temporary(Quality[numTime])
+      Density = fltarr(2, IndexTime);
+      Density[0,*] =  Val
+      Density[1,*] =  Quality
+      Time -= graph[GraphN].TimeBase
+      graph[GraphN].DataFlags[Item] = 1
+      graph[GraphN].dataTime[Item] = Time[0] 
+      average, VEX, Time, Density
+   endif; else Graph[GraphN].ErrorFlag = 1
+ 
+  cmd = Mission + Product +' = VEX' 
+  R = execute(cmd)
+       
+return
+end
+;-----------------------------------------------------------
+++ a/amda_plus/amdalib/get_vexmex_qual.pro
@@ -0,0 +1,48 @@
+;-------------------- Get THEMIS ION Velocity -------------------------------
+;  
+;-----------------------------------------------------------------------------
+pro Get_VexMex_Qual, Mission, Product, st, dt, RESTRICTED = restricted
+
+common GraphC, Graph, GraphN, Item
+common VexHQualCom, VexH
+common VexOQualCom, VexO
+common MexHQualCom, MexH
+common MexOQualCom, MexO
+
+    cmd = 'VEX = ' + Mission + Product; 
+    R = execute(cmd) 
+
+    sz = Vex.DefSize
+
+    Val = fltarr(sz)
+    Time = dblarr(sz)
+
+    ID = Long(graph[GraphN].DataID[Item])
+ 
+    VI = Mission + ':ima:param';
+    if (KEYWORD_SET(RESTRICTED)) then  VI = Mission + ':imaparam:restricted';
+
+    ParamName = Product eq "h" ? "Quality_P" : "Quality_O";
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, st, dt, sz, Time, Val)
+ 
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+   if (RetSize GT 1) then begin 
+      stD = call_external(!DDLIB + 'DD_idl.so','GetDTime', st, /D_VALUE) 
+      numTime = where(Time GE stD AND Time LE (stD+graph[graphN].DeltaT), IndexTime)
+      if IndexTime EQ 0 then return
+      Time = temporary(Time[numTime])
+      Val = temporary(Val[numTime])
+      Time -= graph[GraphN].TimeBase
+      graph[GraphN].DataFlags[Item] = 1
+      graph[GraphN].dataTime[Item] = Time[0]
+      average, VEX, Time, Val
+   endif; else Graph[GraphN].ErrorFlag = 1
+ 
+  cmd = Mission + Product +' = VEX' 
+  R = execute(cmd)
+       
+return
+end
+;-----------------------------------------------------------
+++ a/amda_plus/amdalib/get_vexmex_scan.pro
@@ -0,0 +1,46 @@
+;-------------------- Get THEMIS ION Velocity -------------------------------
+;  
+;-----------------------------------------------------------------------------
+pro Get_VexMex_Scan, Mission,  st, dt, RESTRICTED = restricted
+
+common GraphC, Graph, GraphN, Item
+common VexScanCom, VexScan
+common MexScanCom, MexScan
+
+    cmd = 'VEX = ' + Mission + 'Scan';
+    R = execute(cmd) 
+
+    sz = Vex.DefSize
+
+    Val = intarr(sz)
+    Time = dblarr(sz)
+
+    ID = Long(graph[GraphN].DataID[Item])
+ 
+    VI = Mission + ':ima:param';
+    if (KEYWORD_SET(RESTRICTED)) then  VI = Mission + ':imaparam:restricted';
+
+    ParamName = "ScanFlag";
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, st, dt, sz, Time, Val)
+ 
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+   if (RetSize GT 0) then begin 
+      stD = call_external(!DDLIB + 'DD_idl.so','GetDTime', st, /D_VALUE) 
+      numTime = where(Time GE stD AND Time LE (stD+graph[graphN].DeltaT), IndexTime)
+      if IndexTime EQ 0 then return
+      Time = temporary(Time[numTime])
+      Val = temporary(Val[numTime])
+      Time -= graph[GraphN].TimeBase
+      graph[GraphN].DataFlags[Item] = 1
+      graph[GraphN].dataTime[Item] = Time[0]
+      immediate, VEX, Time, Val
+   endif; else Graph[GraphN].ErrorFlag = 1
+  
+  cmd = Mission +   'Scan  = VEX' 
+  R = execute(cmd)
+       
+return
+end
+;-----------------------------------------------------------
+++ a/amda_plus/amdalib/get_vexmex_temp.pro
@@ -0,0 +1,49 @@
+;-------------------- Get THEMIS ION Velocity -------------------------------
+;  
+;-----------------------------------------------------------------------------
+pro Get_VexMex_Temp, Mission, Product, st, dt, RESTRICTED = restricted
+
+common GraphC, Graph, GraphN, Item
+common VexHTempCom, VexH
+common VexOTempCom, VexO
+common MexHTempCom, MexH
+common MexOTempCom, MexO
+
+    cmd = 'VEX = ' + Mission + Product; 
+    R = execute(cmd) 
+
+    sz = Vex.DefSize
+
+    Val = fltarr(3,sz)
+    Time = dblarr(sz)
+
+    ID = Long(graph[GraphN].DataID[Item])
+ 
+    VI = Mission + ':ima:param'
+    if (KEYWORD_SET(RESTRICTED)) then  VI = Mission + ':imaparam:restricted';
+
+    ParamName = Product eq "h" ? "Temperature_P" : "Temperature_O";
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, st, dt, sz, Time, Val)
+
+
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+   if (RetSize GT 1) then begin 
+      stD = call_external(!DDLIB + 'DD_idl.so','GetDTime', st, /D_VALUE) 
+      numTime = where(Time GE stD AND Time LE (stD+graph[graphN].DeltaT), IndexTime)
+      if IndexTime EQ 0 then return
+      Time = temporary(Time[numTime])
+      Val = temporary(Val[*,numTime])
+      Time -= graph[GraphN].TimeBase
+      graph[GraphN].DataFlags[Item] = 1
+      graph[GraphN].dataTime[Item] = Time[0] 
+      average, VEX, Time, Val
+   endif; else Graph[GraphN].ErrorFlag = 1
+ 
+  cmd = Mission + Product +' = VEX' 
+  R = execute(cmd)
+       
+return
+end
+;-----------------------------------------------------------
+++ a/amda_plus/amdalib/get_vexmex_vel.pro
@@ -0,0 +1,50 @@
+;-------------------- Get THEMIS ION Velocity -------------------------------
+;  
+;-----------------------------------------------------------------------------
+pro Get_VexMex_Vel, Mission, Product, st, dt, RESTRICTED = restricted
+
+common GraphC, Graph, GraphN, Item
+common VexHVelCom, VexH
+common VexOVelCom, VexO
+common MexHVelCom, MexH
+common MexOVelCom, MexO
+
+    cmd = 'VEX = ' + Mission + Product; 
+    R = execute(cmd) 
+
+    sz = Vex.DefSize
+
+    Val = fltarr(3,sz)
+    Time = dblarr(sz)
+
+    ID = Long(graph[GraphN].DataID[Item])
+ 
+    VI = Mission + ':ima:param'
+
+    if (KEYWORD_SET(RESTRICTED)) then  VI = Mission + ':imaparam:restricted'
+
+    ParamName = Product eq "h" ? "Velocity_P" : "Velocity_O";
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, st, dt, sz, Time, Val)
+
+
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+   if (RetSize GT 1) then begin 
+      stD = call_external(!DDLIB + 'DD_idl.so','GetDTime', st, /D_VALUE) 
+      numTime = where(Time GE stD AND Time LE (stD+graph[graphN].DeltaT), IndexTime)
+      if IndexTime EQ 0 then return
+      Time = temporary(Time[numTime])
+      Val = temporary(Val[*,numTime])
+      Time -= graph[GraphN].TimeBase
+      graph[GraphN].DataFlags[Item] = 1
+      graph[GraphN].dataTime[Item] = Time[0] 
+      average, VEX, Time, Val
+   endif; else Graph[GraphN].ErrorFlag = 1
+ 
+  cmd = Mission + Product +' = VEX' 
+  R = execute(cmd)
+       
+return
+end
+;-----------------------------------------------------------
+++ a/amda_plus/amdalib/getcdfvar.pro
@@ -0,0 +1,63 @@
+function getCdfVar, files,  paramName, dims, StartDouble, IntDouble, Time, Data
+
+ Size = N_elements(Time);
+ TIMESHIFT = 62167219200.D0 
+ startIndex = 0L;
+ 
+ datapath = getenv("USER_DATA_PATH") eq '' ?  '../DATA/' : getenv("USER_DATA_PATH");
+
+ for i = 0, n_elements(files) - 1 do begin
+   
+     fileName = datapath + files[i];
+     id = CDF_OPEN(fileName);
+ 
+     CDF_CONTROL, id, SET_ZMODE=2
+
+    ; CDF_CONTROL, id, GET_VAR_INFO=info, VARIABLE=0, /ZVARIABLE 
+     dataType = "UNKNOWN"
+     timeId = 0L;
+     while (dataType ne "CDF_EPOCH" AND dataType ne "CDF_EPOCH16") do begin
+           rest = cdf_varinq(id, timeId, /ZVARIABLE) 
+           dataType = rest.datatype;
+           timeId += 1L;
+     endwhile
+ 
+     timeId -= 1L;
+     recs = getcdfnbrecs_(FileName, timeId);
+ 
+   ; CDF_CONTROL, id, GET_VAR_INFO=info, VARIABLE=0, /ZVARIABLE 
+    ;EXTENDRECS MAXALLOCREC MAXREC MAXRECS NINDEXENTRIES NINDEXRECORDS
+    
+    if (dataType eq "CDF_EPOCH16") then begin 
+     CDF_VARGET, id,  timeId, EPOCH16, REC_COUNT = recs+1, /ZVARIABLE
+     EPOCH = dblarr(recs);
+     for j = 0, recs -1 do begin
+        CDF_EPOCH16, EPOCH16[j], yr, mo, dy, hr, min, sec, milli, micro, pico, /BREAKDOWN_EPOCH 
+        CDF_EPOCH, EPOCH8, yr, mo, dy, hr, min, sec, milli, /COMPUTE_EPOCH 
+        EPOCH[j] = EPOCH8;
+     endfor
+    endif else CDF_VARGET, id,  timeId, EPOCH, REC_COUNT = recs+1, /ZVARIABLE
+       
+    Time1 =  reform(EPOCH/1000.D0 - TIMESHIFT) 
+ 
+      num = where(Time1 GE StartDouble AND Time1 LT StartDouble+IntDouble, Index);
+   
+      if ((Index eq 0 AND startIndex eq 0) OR (startIndex + Index gt Size))  then begin       
+             return, 0
+      endif
+ 
+      if (Index gt 0) then begin
+            Time[startIndex:startIndex+Index-1] = Time1[num];
+        
+            varId = CDF_VARNUM(id, paramName);
+            CDF_VARGET, id,  varId, data1, REC_COUNT = recs+1, /ZVARIABLE
+        
+            if (dims eq 1) then  Data[startIndex:startIndex+Index-1] = reform(Data1[num])  $
+            else Data[*, startIndex:startIndex+Index-1] = Data1[*, num]
+      endif
+      CDF_CLOSE, id
+      startIndex += Index;     
+  endfor  
+  
+ return,startIndex;
+end
+++ a/amda_plus/amdalib/getcefvar.pro
@@ -0,0 +1,48 @@
+function getCefVar, files,  paramName, dims, StartDouble, IntDouble, Time, Data
+ ; CEF_UNDEF, CEF_INT, CEF_FLOAT, CEF_DOUBLE, CEF_CHAR, CEF_ISO_TIME, CEF_ISO_TIME_RANGE,
+  
+ Size = N_elements(Time);
+ TIMESHIFT =  378691200.D0 ; secs 1970 - 1958
+ 
+ startIndex = 0L;
+
+ datapath = getenv("USER_DATA_PATH") eq '' ?  '../DATA/' : getenv("USER_DATA_PATH");
+
+ for i = 0, n_elements(files) - 1 do begin
+   
+     fileName = datapath + files[i];
+     nrec = CEF_READ(fileName);
+     if (nrec LE 0) then continue;
+
+     if (i EQ 0) then begin
+	dataType = "UNKNOWN"
+	timeId = 0L;
+        vars = CEF_VARNAMES();
+	while (dataType ne "ISO_TIME") do begin
+	      dataType = CEF_VATTR(vars[timeId],"VALUE_TYPE"); 
+	      timeId += 1L;
+	endwhile
+	timeId -= 1L;
+     endif
+
+      Time1 = CEF_VAR(vars[timeId])/1000.D0 -  TIMESHIFT;
+       
+      num = where(Time1 GE StartDouble AND Time1 LT StartDouble+IntDouble, Index);
+ 
+      if ((Index eq 0 AND startIndex eq 0) OR (startIndex + Index gt Size))  then return, 0
+        
+      if (Index gt 0) then begin
+            Time[startIndex:startIndex+Index-1] = Time1[num];
+        
+            Data1 = CEF_VAR(paramName);
+           ; dims = CEF_VATTR(paramName, "SIZES")
+            dims = n_elements(Data1)/n_elements(Time1);
+            if (dims eq 1) then  Data[startIndex:startIndex+Index-1] = Data1[num]  $
+            else Data[*, startIndex:startIndex+Index-1] = Data1[*, num]
+      endif
+      CEF_CLOSE
+      startIndex += Index;     
+  endfor 
+  
+ return,startIndex;
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/getmex_n_sw.pro
@@ -0,0 +1,46 @@
+;--------------------  -------------------------------
+;  
+;-----------------------------------------------------------------------------
+pro GetMex_N_SW, st, dt
+
+common GraphC, Graph, GraphN, Item
+common Mex_N_SW, Mex
+
+ 
+    sz = Mex.DefSize
+     
+    Quality = fltarr(sz)
+    Val = fltarr(sz)
+    Time = dblarr(sz)
+
+    ID = Long(graph[GraphN].DataID[Item]) 
+    VI = 'mex:ima:sw'
+      
+    RetSize = call_external(!PROJLIB + '/libParam.so', 'getLocalData', $
+                          ID, VI,  'Flag', st, dt, sz, Time, Quality)
+ 
+     
+    RetSize = call_external(!PROJLIB + '/libParam.so', 'getLocalData', $
+                          ID, VI, 'N', st, dt, sz, Time, Val)
+ 
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+   if (RetSize GE 1) then begin 
+      stD = call_external(!DDLIB + 'DD_idl.so','GetDTime', st, /D_VALUE) 
+      numTime = where(Time GE stD AND Time LE (stD + graph[graphN].DeltaT), IndexTime)
+      if IndexTime EQ 0 then return
+      Time = temporary(Time[numTime])
+      Val = temporary(Val[numTime])
+      Quality = temporary(Quality[numTime])
+      Density = fltarr(2, IndexTime);
+      Density[0,*] =  Val
+      Density[1,*] =  Quality
+      Time -= graph[GraphN].TimeBase
+      graph[GraphN].DataFlags[Item] = 1
+      graph[GraphN].dataTime[Item] = Time[0] 
+      (*Mex.Time) = N_elements(*Mex.Time) eq 0 ? Time : [*Mex.Time,Time]  
+      (*Mex.Val) =  N_elements(*Mex.Val)  eq 0 ?   reform(Density,2*IndexTime) : [*Mex.Val,reform(Density,2*IndexTime)]
+   endif ;else Graph[GraphN].ErrorFlag = 1
+         
+return
+end
+;-----------------------------------------------------------
+++ a/amda_plus/amdalib/getmex_v_sw.pro
@@ -0,0 +1,46 @@
+;--------------------  -------------------------------
+;  
+;-----------------------------------------------------------------------------
+pro GetMex_v_SW, st, dt
+
+common GraphC, Graph, GraphN, Item
+common Mex_v_SW, Mex
+
+ 
+    sz = Mex.DefSize
+     
+    Quality = fltarr(sz)
+    Val = fltarr(sz)
+    Time = dblarr(sz)
+
+    ID = Long(graph[GraphN].DataID[Item]) 
+    VI = 'mex:ima:sw'
+      
+    RetSize = call_external(!PROJLIB + '/libParam.so', 'getLocalData', $
+                          ID, VI,  'Flag', st, dt, sz, Time, Quality)
+ 
+     
+    RetSize = call_external(!PROJLIB + '/libParam.so', 'getLocalData', $
+                          ID, VI, 'V', st, dt, sz, Time, Val)
+ 
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID
+   if (RetSize GE 1) then begin 
+      stD = call_external(!DDLIB + 'DD_idl.so','GetDTime', st, /D_VALUE) 
+      numTime = where(Time GE stD AND Time LE (stD + graph[graphN].DeltaT), IndexTime)
+      if IndexTime EQ 0 then return
+      Time = temporary(Time[numTime])
+      Val = temporary(Val[numTime])
+      Quality = temporary(Quality[numTime])
+      Density = fltarr(2, IndexTime);
+      Density[0,*] =  Val
+      Density[1,*] =  Quality
+      Time -= graph[GraphN].TimeBase
+      graph[GraphN].DataFlags[Item] = 1
+      graph[GraphN].dataTime[Item] = Time[0] 
+      (*Mex.Time) = N_elements(*Mex.Time) eq 0 ? Time : [*Mex.Time,Time]  
+      (*Mex.Val) =  N_elements(*Mex.Val)  eq 0 ?   reform(Density,2*IndexTime) : [*Mex.Val,reform(Density,2*IndexTime)]
+   endif ;else Graph[GraphN].ErrorFlag = 1
+         
+return
+end
+;-----------------------------------------------------------
+++ a/amda_plus/amdalib/getncvar.pro
@@ -0,0 +1,53 @@
+function getNcVar, files,  paramName, dims, StartDouble, IntDouble, Time, Data
+ 
+ Size = N_elements(Time);
+
+ LastTime = -1;
+
+ startIndex = 0L;
+ datapath = getenv("USER_DATA_PATH") eq '' ?  '../DATA/' : getenv("USER_DATA_PATH");
+
+ for i = 0, n_elements(files) - 1 do begin
+   
+     fileName = datapath + files[i];
+     id = NCDF_OPEN(fileName);
+ 
+     
+     TimeID = NCDF_VARID(id, 'Time');
+     VarID =  NCDF_VARID(id, paramName);
+   
+     res = NCDF_VARINQ(id, VarID);     
+     dims = res.ndims;
+           
+      NCDF_VARGET, id, TimeId, Time1 
+ 
+      num = where(Time1 GE StartDouble AND Time1 LT StartDouble+IntDouble AND Time1 GT LastTime, Index);
+         
+      if (Index eq 0 AND startIndex eq 0 AND i eq (n_elements(files) - 1))  then return, 0
+
+      stopIndex = startIndex + Index;
+ 
+      if (stopIndex gt Size) then begin
+            Time1 = Time1[0:Index-(stopIndex-Size)]
+            num = where(Time1 GE StartDouble AND Time1 LT StartDouble+IntDouble, Index);
+            stopIndex =  Size ; 
+      endif    
+       
+      if (Index GT 0) then begin
+        Time[startIndex:stopIndex-1] = Time1[num];
+        
+        NCDF_VARGET, id, VarId, data1
+    
+        if (dims eq 1) then  Data[startIndex:stopIndex-1] = reform(Data1[num])  $
+        else Data[*, startIndex:stopIndex-1] = Data1[*, num]
+
+      endif 
+
+      NCDF_CLOSE, id
+      startIndex += Index;  
+      LastTime =  Time[stopIndex-1]  
+  endfor  
+ 
+ return, startIndex;
+
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/getros_lap_sweep.pro
@@ -0,0 +1,60 @@
+; $Id: getlocalparam.pro,v 1.5 2011/09/21 16:01:45 budnik Exp $
+;--------------------------------------
+;   GETNEWPARAM - TEMPLATE
+;
+; Fri Sep 21 11:57:01 CEST 2007
+;
+;---------------------------------------
+
+ pro GETros_lap_sweep, LAP, st, dt
+ 
+ 
+common ros_lap1_sweep_Common, Final1
+common ros_lap2_sweep_Common, Final2
+common GraphC, graph, GraphN, Item
+
+    cmd = 'Final = Final'+string(LAP,format='(i1)')
+    R = execute(cmd);
+
+    paramName = 'full_sweep';
+    
+    VI = 'ros:lap'+string(LAP,format='(i1)')+':sweep';
+    ID = Graph[GraphN].DataID[Item];
+
+    Size = Final.DefSize
+    flag = intarr(size); 
+    Time  = dblarr(size);
+
+  
+     RetSize = call_external(!PROJLIB + '/libParam.so', 'getLocalData', $         
+                          ID, VI, 'macros', st, dt, size, Time, flag);
+
+     if (RetSize GT 1) then begin 
+        flag = temporary(flag[0:RetSize-1]);
+           
+        dims = 241L;          
+        Data = fltarr(dims,retSize);       
+        Time  = temporary(Time[0:retSize-1]);    
+        Data1 =  replicate(Graph[GraphN].FillValue[Item], (dims+1L)*retSize);
+    
+        RetSize = call_external(!PROJLIB + '/libParam.so', 'getLocalData', $         
+                            ID, VI, paramName, st, dt, size, Time, Data);
+ 
+        
+         Data1 = temporary(reform(Data1,dims+1L,Retsize));
+         Data1[0:dims-1, *] =  Data[*, 0:Retsize - 1];                                  
+ 
+        Data1[dims,*] = float(flag);       
+        Time = Time - graph[GraphN].TimeBase;
+        graph[GraphN].DataFlags[Item] = 1;
+        graph[GraphN].dataTime[Item] = Time[0];
+         
+        averagelap, Final, Time[ 0:Retsize - 1],Data1  
+      
+    endif;  else Graph[GraphN].ErrorFlag = 1;
+
+    cmd = 'Final'+string(LAP,format='(i1)')+ ' = Final'
+    R = execute(cmd)
+ 
+return
+end
+++ a/amda_plus/amdalib/getros_mip_surv.pro
@@ -0,0 +1,55 @@
+; $Id: getlocalparam.pro,v 1.5 2011/09/21 16:01:45 budnik Exp $
+;--------------------------------------
+;   GETNEWPARAM - TEMPLATE
+;
+; Fri Sep 21 11:57:01 CEST 2007
+;
+;---------------------------------------
+
+ pro GETros_mip_surv, st, dt
+ 
+ 
+common ros_mip_surv_Common, Final 
+common GraphC, graph, GraphN, Item
+
+    paramName = 'SurveyFull';
+    dims = 92L;    
+    VI = 'ros:mip:surv';
+    Size = Final.DefSize
+
+    
+    Data = fltarr(dims,size);       
+    flag = lonarr(size); 
+    Time  =  dblarr(size)
+    Data1 = fltarr(dims+1,size);
+
+   ID = Long(graph[GraphN].DataID[Item])
+   Lag = Graph[GraphN].Lag[0];
+   Sampling =  Graph[GraphN].Sampling[Item];
+ 
+   if (Lag ne 0.0) then begin
+     rt = call_external(!DDLIB+'DD_idl.so', 'GetDTime', st, /D_VALUE)
+     rt -=  Lag
+     st = call_external(!DDLIB+'DD_idl.so', 'GetTime', double(rt), /S_VALUE)
+   endif
+ 
+   RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $         
+                          ID, VI, paramName, st, dt, size, Time, Data) 
+   
+   RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $         
+                          ID, VI, 'Frequency_table_idx', st, dt, size, Time, flag) 
+ 
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID    
+    
+   if (RetSize GT 1) then begin  
+         Data1[0:dims-1, 0:RetSize-1] = Data[*,0:RetSize-1];
+         Data1[dims,  0:RetSize-1] = float(flag[0:RetSize-1]);
+   
+         Time = Time - graph[GraphN].TimeBase + Lag;
+         graph[GraphN].DataFlags[Item] = 1;
+         graph[GraphN].dataTime[Item] = Time[0];
+         average, Final, Time[0:RetSize-1], Data1[*,0:RetSize-1]
+   endif;  else Graph[GraphN].ErrorFlag = 1;
+    
+return
+end
+++ a/amda_plus/amdalib/getros_mip_survphase.pro
@@ -0,0 +1,55 @@
+; $Id: getlocalparam.pro,v 1.5 2011/09/21 16:01:45 budnik Exp $
+;--------------------------------------
+;   GETNEWPARAM - TEMPLATE
+;
+; Fri Sep 21 11:57:01 CEST 2007
+;
+;---------------------------------------
+
+ pro GETros_mip_survphase, st, dt
+ 
+ 
+common ros_mip_survphase_Common, Final 
+common GraphC, graph, GraphN, Item
+
+    paramName = 'SurveyFull_Phase';
+    dims = 92L;    
+    VI = 'ros:mip:survphase';
+    Size = Final.DefSize
+
+    
+    Data = fltarr(dims,size);       
+    flag = lonarr(size); 
+    Time  =  dblarr(size)
+    Data1 = fltarr(dims+1,size);
+
+   ID = Long(graph[GraphN].DataID[Item])
+   Lag = Graph[GraphN].Lag[0];
+   Sampling =  Graph[GraphN].Sampling[Item];
+ 
+   if (Lag ne 0.0) then begin
+     rt = call_external(!DDLIB+'DD_idl.so', 'GetDTime', st, /D_VALUE)
+     rt -=  Lag
+     st = call_external(!DDLIB+'DD_idl.so', 'GetTime', double(rt), /S_VALUE)
+   endif
+ 
+   RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $         
+                          ID, VI, 'Frequency_table_idx', st, dt, size, Time, flag) 
+
+   RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $         
+                          ID, VI, paramName, st, dt, size, Time, Data) 
+ 
+   if (ID GE 0) then  Graph[GraphN].DataID[Item] = ID    
+    
+   if (RetSize GT 1) then begin  
+         Data1[0:dims-1, 0:RetSize-1] = Data[*,0:RetSize-1];
+         Data1[dims,  0:RetSize-1] = float(flag[0:RetSize-1]);
+   
+         Time = Time - graph[GraphN].TimeBase + Lag;
+         graph[GraphN].DataFlags[Item] = 1;
+         graph[GraphN].dataTime[Item] = Time[0];
+         averagemip, Final, Time[0:RetSize-1], Data1[*,0:RetSize-1]
+   endif ; else Graph[GraphN].ErrorFlag = 1;
+    
+return
+end
+++ a/amda_plus/amdalib/gettxtvar.pro
@@ -0,0 +1,40 @@
+function getTxtVar, files,  paramName, dims, StartDouble, IntDouble, Time, Data
+ 
+ Size = N_elements(Time);
+  
+ startIndex = 0L;
+ i_start = long(paramName);
+ i_stop = i_start + dims - 1;
+ 
+ datapath = getenv("USER_DATA_PATH") eq '' ?  '../DATA/' : getenv("USER_DATA_PATH");
+ for i = 0, n_elements(files) - 1 do begin
+ ; print,  i, "FILENEME", files[i]
+     openr, LUN, datapath+files[i], err = err, /GET_LUN ;
+     temp = " ";
+    
+     while not EOF(LUN) do begin
+        readf, LUN, temp
+        arr = strsplit(temp, " ", /EXTRACT);
+        Time1 = N_Elements(Time1) eq 0 ? double(arr[0]) : [Time1,double(arr[0])];   
+        Data1 = N_Elements(Data1) eq 0 ? float(arr[i_start:i_stop]) : [Data1, float(arr[i_start:i_stop])]        
+    endwhile
+
+    close,LUN
+    free_lun,LUN  
+  endfor 
+
+   num = where(Time1 GE StartDouble AND Time1 LT StartDouble+IntDouble, Index);
+ ;  print, "INDEX ",Index, " SIZE ", Size, " INDEX+ ", startIndex + Index 
+   if ((Index eq 0 AND startIndex eq 0) OR (startIndex + Index gt Size))  then return, 0;   
+
+   if (Index GT 0) then begin 
+        Time[startIndex:startIndex+Index-1] = Time1[num];
+        if (dims eq 1) then  Data[startIndex:startIndex+Index-1] = reform(Data1[num])  $
+        else begin 
+                Data1 = reform(Data1,dims,N_elements(Time1))
+                Data[*, startIndex:startIndex+Index-1] = Data1[*, num]
+        endelse      
+   endif
+  startIndex += Index; 
+ return, startIndex;
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/ica_spec_com.pro
@@ -0,0 +1,57 @@
+;===========================================================
+;
+;===================== ========================================
+pro ica_spec_com 
+
+common IcaSpecCom, Spec
+common GraphC, graph, GraphN, Item
+common Etab, Flag, Tab
+Common Record, Record
+
+ if (N_ELEMENTS(Flag) eq 0) then Flag = GraphN
+
+ FillValue =  !Values.F_NAN
+ Sampling = 192.0
+ Sz =  Long((Graph[GraphN].DeltaT)/Sampling + 1)
+
+  Spec = record
+  Spec.DefSize = long(Graph[GraphN].DeltaT/Sampling + 1)
+  Spec.Sampling = Sampling
+  Graph[GraphN].FillValue[Item] = FillValue
+  ; Graph[GraphN].Lag[Item] = double(Delay)
+
+  Spec.Time = PTR_NEW(/ALLOCATE_HEAP)
+  Spec.Val = PTR_NEW(/ALLOCATE_HEAP)
+  Spec.LastTime = PTR_NEW(/ALLOCATE_HEAP)
+  Spec.LastVal = PTR_NEW(/ALLOCATE_HEAP)
+
+  if (Flag eq GraphN) then begin
+      Tab = record;
+      Tab.DefSize = long(Graph[GraphN].DeltaT/Sampling + 1)
+      Tab.Sampling = Sampling
+      
+      Tab.Time = PTR_NEW(/ALLOCATE_HEAP)
+      Tab.Val = PTR_NEW(/ALLOCATE_HEAP)
+      Tab.LastTime = PTR_NEW(/ALLOCATE_HEAP)
+      Tab.LastVal = PTR_NEW(/ALLOCATE_HEAP)
+      Flag = GraphN;
+  endif
+
+  Graph[GraphN].Sampling[Item] =  Sampling
+  Graph[GraphN].Time[Item] =  Spec.Time
+  Graph[GraphN].Val[Item] = Spec.Val
+ 
+  if (Graph[GraphN].Step GT Sampling*2.0) then begin
+    *(Spec.Time) = indgen(Graph[GraphN].ppp)*Graph[GraphN].Step + Graph[GraphN].Step/2.0
+    *(Spec.Val)  =  replicate(FillValue, 16L*96L*16L*Graph[GraphN].ppp)
+    if (Flag eq GraphN) then begin
+	*(Tab.Time) = indgen(Graph[GraphN].ppp)*Graph[GraphN].Step + Graph[GraphN].Step/2.0
+	*(Tab.Val)  =  replicate(FillValue, Graph[GraphN].ppp)
+    endif
+  endif
+  
+
+return
+end
+;------------------------------------------------------------
+
+++ a/amda_plus/amdalib/icacntplot.pro
@@ -0,0 +1,104 @@
+pro icacntplot, azMin, azMax, enMin, enMax, elevation
+
+  common GraphC, Graph, GraphN, Item
+  common Etab, Flag, ETab
+    
+    ;------- Energy table -----------------    
+   
+    Tab = fix(*ETab[0].val); 
+    VI = "ros:ica:all";
+    ID = Long(graph[GraphN].DataID[Item]);
+    EnNum = 96
+    Energy = fltarr(EnNum)
+    numTab = uniq(Tab,sort(Tab));
+    uniqTab = Tab[numTab];
+    temp = fltarr(96L,2L); 
+
+
+    err = call_external(!PROJLIB + '/libParam.so', 'getInfo', $
+                             ID, VI, "Energy", 2L*96L, temp)
+
+    Energy = temp[*, uniqTab[0]];
+    Energy =  temporary(Energy(where(Energy GT 0)));
+    EnNum = N_elements(Energy)
+    
+  ;-------- Cat in the Energy range ------------
+   E0 =  min(Energy)
+   E1 =  max(Energy)
+   eindex = where((Energy GE E0) AND $
+                  (Energy LE E1) AND $
+                  (Energy GT 0.0),counte)
+   if counte GT 0 then begin
+      EBounds = fltarr(EnNum+1)
+      EBounds(min(eindex)) = Energy(min(eindex))
+      EBounds(max(eindex+1)) = Energy(max(eindex))
+      for ie = min(eindex+1), max(eindex) do EBounds(ie) = sqrt(Energy(ie-1)*Energy(ie))
+   endif else return
+
+;--------------------- No DATA - just axes -----------------------
+   if (N_Elements(*Graph[GraphN].Time[0]) eq 0) then begin
+         nodata
+         return
+   endif
+
+
+    Time = (*Graph[GraphN].Time[0])
+
+    Val = (*Graph[GraphN].Val[0])
+    N_Time = N_elements(Time);
+
+   
+   if (elevation) then begin 
+      Val = temporary(reform(Val, 96, 16, 16, N_Time)) 
+      LocVal =  azMax GT azMin ? total(total(Val[*, azMin:azMax, 0:15, *],2),2) : reform(total(Val[*,azMin,0:15,*],3), 96, N_time);
+      newTime = Time;  
+      deltaT = 96.0;
+      gapT = deltaT*10.0 ;
+  endif else begin
+; elevation as TIME
+      newTime = dblarr(N_Time*16);
+      for i = 0, N_Time - 1 do for j = 0, 15 do newTime[i*16 + j] = Time[i] + j*12.D0
+      N_Time *= 16;
+      Val = temporary(reform(Val, 96, 16, N_Time));
+      LocVal =  azMax GT azMin ? total(Val[*,azMin:azMax,*],2) : reform(Val[*,azMin,*], 96, N_time);  
+      deltaT = 6.0;
+      gapT = deltaT*10.0;
+    endelse
+      
+    num0 = where(EBounds ge enMin, Index0);
+    num1 = where(EBounds ge enMax, Index1);
+    
+    if (Index1 eq 0)  then numStop = 0 else  numStop = num1[Index1-1];
+    if (Index0 eq EnNum+1) then numStart = 95 else numStart = num0[Index0-1];
+    
+    Data =  numStart ne numStop ? total(LocVal[numStop:numStart,*],1) : reform(LocVal[numStart,*], N_Time);   
+    numMinMax = where(finite(Data), IndexMinMax)
+
+    MinVal = min(Data[numMinMax]) > 0.1;   
+    MaxVal = max(Data[numMinMax]);
+   
+     
+   if Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] then begin
+                Graph[GraphN].Ly.range = Graph[GraphN].AxisF eq 'logaxis' ? [MinVal, MaxVal] : [0.0, MaxVal]
+                Graph[GraphN].Ly.style = 0
+    endif else Graph[GraphN].Ly.style = 1
+    
+    !y = Graph[GraphN].Ly;  
+  
+    if (execute(Graph[GraphN].AxisF) NE 1) then plot, newTime, Data, MAX_VALUE=MaxValue,  color = !dnc-1, /NODATA, /NOERASE 
+        
+    plots, newTime[0]+deltaT, Data[0], /DATA
+ 
+    for i = 1, N_elements(newTime) - 1 do begin
+      if (newTime[i] LE Graph[GraphN].TotalTime) then begin
+          if (data[i] lt Graph[GraphN].Ly.range[0]) then data[i] = Graph[GraphN].Ly.range[0];
+          if (data[i] gt Graph[GraphN].Ly.range[1]) then data[i] = Graph[GraphN].Ly.range[1];  
+          if (newTime[i] - newTime[i-1] LE gapT) then $
+                    plots, newTime[i]+deltaT, data[i], color=col, /CONTINUE, /DATA $
+          else   plots, newTime[i]+deltaT, data[i], /DATA 
+        endif
+    endfor
+      
+return
+end
+;--------------------------------------------------------------------
+++ a/amda_plus/amdalib/icaspecplot.pro
@@ -0,0 +1,186 @@
+;============================================================
+;
+;=============================================================
+
+pro icaspecplot,  AzMin, AzMax, elevation, minCount, maxCount
+
+   common GraphC, Graph, GraphN, Item
+   common Etab, Flag, ETab
+
+   ;------- Energy table -----------------    
+    
+    VI = "ros:ica:all";
+    ID = Long(graph[GraphN].DataID[Item])
+ 
+    EnNum0 = 96
+    Energy = fltarr(EnNum0);
+    Energy1 = fltarr(EnNum0);
+; set min/max energy
+    EnergyMin = 40000.0;
+    EnergyMax = 0.0;
+ 
+    if (n_elements(*ETab[0].val) eq 0) then begin 
+     !y = Graph[GraphN].Ly
+      plot_io, graph[GraphN].Lx.range, [1., 40000.0], /nodata, /noerase, ystyle = 1, $
+                              ytickformat='yticks', charsize=graph[GraphN].Lp.charsize*0.8, color = !dnc-1
+      nodatafortimeinterval
+     return
+    endif
+
+    Tab = fix(*ETab[0].val); 
+    numTab = uniq(Tab,sort(Tab));
+    uniqTab = Tab[numTab];
+    temp = fltarr(96L,2L);
+
+
+    err = call_external(!PROJLIB + '/libParam.so', 'getInfo', $
+                             ID, VI, "Energy", 2L*96L, temp)
+
+    for i = 0, N_elements(uniqTab) - 1 do begin 
+        if uniqTab[i] lt 0 then continue;
+         Energy1 = temp[*, uniqTab[i]]
+         Energy1 =  temporary(Energy1(where(Energy1 GT 0)));
+        if (min(Energy1) LT EnergyMin) then EnergyMin = min(Energy1);
+        if (max(Energy1) GT EnergyMax) then EnergyMax = max(Energy1);
+    endfor 
+ 
+    Energy = temp[*, uniqTab[0]]
+    Energy =  temporary(Energy(where(Energy GT 0)));
+    EnNum = N_elements(Energy)
+  
+  ;-------- Cat in the Energy range ------------
+   E0 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ?  min(Energy) : Graph[GraphN].Ly.range[0]
+   E1 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ?  max(Energy) : Graph[GraphN].Ly.range[1]
+
+   eindex = where((Energy GE E0) AND $
+                  (Energy LE E1) AND $
+                  (Energy GT 0.0),counte)
+
+   if counte GT 0 then begin
+      EBounds = fltarr(EnNum+1)
+      EBounds(min(eindex)) = Energy(min(eindex))
+      EBounds(max(eindex+1)) = Energy(max(eindex))
+      for ie = min(eindex+1), max(eindex) do EBounds[ie] = sqrt(Energy[ie-1]*Energy[ie])
+   endif else begin
+         !y = Graph[GraphN].Ly
+         plot_io, graph[GraphN].Lx.range, [1., 40000.0], /nodata, /noerase, ystyle = 1, $
+                              ytickformat='yticks', charsize=graph[GraphN].Lp.charsize*0.8, color = !dnc-1
+        nodatafortimeinterval
+        return
+   endelse
+ 
+;--------------------- No DATA - just axes -----------------------
+   if (N_Elements(*Graph[GraphN].Time[0]) eq 0) then begin
+         nodata
+         return
+   endif
+
+    Time = (*Graph[GraphN].Time[0])
+
+    Val = (*Graph[GraphN].Val[0])
+    N_Time = N_elements(Time);
+ 
+   if (elevation) then begin 
+      Val = temporary(reform(Val, 96, 16, 16, N_Time)) 
+      LocVal =  azMax GT azMin ? total(total(Val[*, azMin:azMax, 0:15, *],2),2) : reform(total(Val[*,azMin,0:15,*],3), 96, N_time);
+      newTime = Time;      
+  endif else begin
+; elevation as TIME
+      newTime = dblarr(N_Time*16);
+      for i = 0, N_Time - 1 do for j = 0, 15 do newTime[i*16 + j] = Time[i] + j*12.D0
+      N_Time *= 16;
+      Val = temporary(reform(Val, 96, 16, N_Time));
+      LocVal =  azMax GT azMin ? total(Val[*,azMin:azMax,*],2) : reform(Val[*,azMin,*], 96, N_time);  
+    endelse
+
+  
+   numMinMax = where(LocVal gt 0.0, IndexMinMax)   
+ 
+
+  if (Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1]) then  $
+                                      Graph[GraphN].Ly.range = [EnergyMin, EnergyMax]
+       
+   Graph[GraphN].Ly.style = 1
+
+  !y = Graph[GraphN].Ly
+ 
+   plot_io, graph[GraphN].Lx.range, graph[GraphN].Ly.range, /nodata, /noerase, ystyle = 1, $
+                              ytickformat='yticks', charsize=graph[GraphN].Lp.charsize*0.8, color = !dnc-1
+
+   Step = Graph[GraphN].Sampling[0] > Graph[GraphN].Step;  
+   if (not elevation) then Step /= 16.0;    
+
+ ;--------------------- No DATA - just axes -----------------------
+  if (IndexMinMax LE 0) then begin
+         nodatafortimeinterval
+         return
+  endif
+  if (Time[0] GT Graph[GraphN].TotalTime) then begin
+        nodatafortimeinterval
+         return
+  endif
+
+   MinVal = (MinCount ge 0) ? MinCount > 0.1: min(Val[numMinMax]) > 0.1
+   MaxVal = (MaxCount ge 0) ? MaxCount : max(Val[numMinMax]) 
+  
+ ; MinVal = min(LocVal[numMinMax]) > 0.1
+ ; MaxVal = max(LocVal[numMinMax])   
+ 
+ ;-----------   NOW PLOT SPECTRA
+    alogMin = alog(MinVal)
+    KF = float(!DNC - 2) /(alog(MaxVal) - alogMin)
+    xx = fltarr(4) & yy = fltarr(4)
+ 
+     for i = 0, N_Time - 2 do begin
+
+	ii = elevation ? i : i/16 < (N_elements(Tab) - 2);
+ 
+	if (Tab[ii+1] ne Tab[ii] AND Tab[ii+1] GE 0) then begin                
+	      Energy = temp[*,Tab[ii+1]]
+	      Energy =  temporary(Energy(where(Energy GT 0)));
+	      EnNum = N_elements(Energy)
+	      
+	    ;-------- Cat in the Energy range ------------
+	    E0 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ? min(Energy) : Graph[GraphN].Ly.range[0]
+	    E1 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ? max(Energy) : Graph[GraphN].Ly.range[1]
+	    eindex = where((Energy GE E0) AND $
+			    (Energy LE E1) AND $
+			    (Energy GT 0.0),counte)
+	    if counte GT 0 then begin
+		EBounds = fltarr(EnNum+1)
+		EBounds(min(eindex)) = Energy(min(eindex))
+		EBounds(max(eindex+1)) = Energy(max(eindex))
+		for ie = min(eindex+1), max(eindex) do EBounds(ie) = sqrt(Energy(ie-1)*Energy(ie))
+	    endif else return
+	endif
+        xx[0] = newTime[i]  > 0.0
+        xx[1] = xx[0]
+        xx[2] = newTime[i+1] < newTime[i] + Step
+        xx[3] = xx[2]
+
+        if (xx[2] LE graph[GraphN].TotalTime AND xx[0] LT graph[GraphN].TotalTime) then begin
+           for k = min(eindex), max(eindex) do begin
+                  yy[0] = EBounds[k]
+                  yy[1] = EBounds[k+1]
+                  yy[2] = yy[1]
+                  yy[3] = yy[0]
+                  lcount = LocVal[k,i] ;< 1.e5
+                  if finite(lcount) then begin
+                   cl = lcount GT MinVal ? fix(KF*(alog(lcount) - alogMin)) : 1;
+                  endif else cl = 0;
+                  
+                  if (cl GT !DNC - 2) then  cl = !DNC - 2
+                  PolyFill, xx ,yy, color = cl, /data 
+            endfor
+        endif
+     endfor
+  
+;-------------------- Legend ----------------------------------------
+       spec_legend, MinVal, MaxVal
+;------------------ MIN-MAX------------------------------------------------
+      Graph[GraphN].Min  = MinVal
+      Graph[GraphN].Max  = MaxVal
+
+return
+end
+;--------------------------------------------------------------------
+++ a/amda_plus/amdalib/ima_h_spec_com.pro
@@ -0,0 +1,72 @@
+;===========================================================
+;
+;                   VEX_H_SPEC COM
+;  18 Mars 2009
+;===================== ========================================
+pro ima_h_spec_com, Mission
+
+common VexHSpecCom, VexHSpec
+common MexHSpecCom, MexHSpec
+
+common GraphC, graph, GraphN, Item
+Common Record, Record
+common Etab, FlagVex, EtabVex, FlagMex, EtabMex
+
+
+ if (Mission eq 'vex') then begin
+      if (N_ELEMENTS(FlagVex) eq 0) then Flag = GraphN $
+      else Flag = FlagVex;
+ endif else begin
+      if (N_ELEMENTS(FlagMex) eq 0) then Flag = GraphN $
+      else Flag = FlagMex;
+ endelse 
+   
+ FillValue =  !Values.F_NAN
+ Sampling = 192.0
+ Sz =  Long((Graph[GraphN].DeltaT)/Sampling + 1)
+
+  Spec = record
+  Spec.DefSize = long(Graph[GraphN].DeltaT/Sampling + 1)
+  Spec.Sampling = Sampling
+  Graph[GraphN].FillValue[Item] = FillValue
+  ; Graph[GraphN].Lag[Item] = double(Delay)
+
+  Spec.Time = PTR_NEW(/ALLOCATE_HEAP)
+  Spec.Val = PTR_NEW(/ALLOCATE_HEAP)
+  Spec.LastTime = PTR_NEW(/ALLOCATE_HEAP)
+  Spec.LastVal = PTR_NEW(/ALLOCATE_HEAP)
+
+  Graph[GraphN].Sampling[Item] =  Sampling
+  Graph[GraphN].Time[Item] =  Spec.Time
+  Graph[GraphN].Val[Item] = Spec.Val
+
+  if (Flag eq GraphN) then begin
+	Tab = record;
+	Tab.DefSize = long(Graph[GraphN].DeltaT/Sampling + 1)
+	Tab.Sampling = Sampling
+	
+	Tab.Time = PTR_NEW(/ALLOCATE_HEAP)
+	Tab.Val = PTR_NEW(/ALLOCATE_HEAP)
+	Tab.LastTime = PTR_NEW(/ALLOCATE_HEAP)
+	Tab.LastVal = PTR_NEW(/ALLOCATE_HEAP)
+	Flag = GraphN;
+    endif
+ 
+  if (Graph[GraphN].Step GT Sampling*2.0) then begin
+    *(Spec.Time) = indgen(Graph[GraphN].ppp)*Graph[GraphN].Step + Graph[GraphN].Step/2.0
+    *(Spec.Val)  = replicate(FillValue, 16L*96L*16L*Graph[GraphN].ppp)
+    if (Flag eq GraphN) then begin
+	*(Tab.Time) = indgen(Graph[GraphN].ppp)*Graph[GraphN].Step + Graph[GraphN].Step/2.0
+	*(Tab.Val)  =  replicate(FillValue, Graph[GraphN].ppp)
+    endif
+  endif
+  
+  if (Flag eq GraphN) then cmd = Mission + 'HSpec = Spec & Etab' + Mission +'= Tab & Flag'+ Mission +'= Flag' $
+  else cmd = Mission + 'HSpec = Spec';
+
+  R = execute(cmd);
+  
+return
+end
+;------------------------------------------------------------
+
+++ a/amda_plus/amdalib/ima_o_spec_com.pro
@@ -0,0 +1,70 @@
+;===========================================================
+;
+;                   VEX_o_SPEC COM
+;  18 Mars 2009
+;===================== ========================================
+pro ima_o_spec_com, Mission
+
+common VexOSpecCom, VexOSpec
+common MexOSpecCom, MexOSpec
+common GraphC, graph, GraphN, Item
+common Etab, FlagVex, EtabVex, FlagMex, EtabMex
+Common Record, Record
+
+ if (Mission eq 'vex') then begin
+      if (N_ELEMENTS(FlagVex) eq 0) then Flag = GraphN $
+      else Flag = FlagVex;
+ endif else begin
+      if (N_ELEMENTS(FlagMex) eq 0) then Flag = GraphN $
+      else Flag = FlagMex;
+ endelse  
+
+ FillValue =  !Values.F_NAN
+ Sampling = 192.0
+ Sz =  Long((Graph[GraphN].DeltaT)/Sampling + 1)
+
+  Spec = record
+  Spec.DefSize = long(Graph[GraphN].DeltaT/Sampling + 1)
+  Spec.Sampling = Sampling
+  Graph[GraphN].FillValue[Item] = FillValue
+  ; Graph[GraphN].Lag[Item] = double(Delay)
+
+  Spec.Time = PTR_NEW(/ALLOCATE_HEAP)
+  Spec.Val = PTR_NEW(/ALLOCATE_HEAP)
+  Spec.LastTime = PTR_NEW(/ALLOCATE_HEAP)
+  Spec.LastVal = PTR_NEW(/ALLOCATE_HEAP)
+
+  if (Flag eq GraphN) then begin
+      Tab = record;
+      Tab.DefSize = long(Graph[GraphN].DeltaT/Sampling + 1)
+      Tab.Sampling = Sampling
+      
+      Tab.Time = PTR_NEW(/ALLOCATE_HEAP)
+      Tab.Val = PTR_NEW(/ALLOCATE_HEAP)
+      Tab.LastTime = PTR_NEW(/ALLOCATE_HEAP)
+      Tab.LastVal = PTR_NEW(/ALLOCATE_HEAP)
+      Flag = GraphN;
+  endif
+
+  Graph[GraphN].Sampling[Item] =  Sampling
+  Graph[GraphN].Time[Item] =  Spec.Time
+  Graph[GraphN].Val[Item] = Spec.Val
+ 
+  if (Graph[GraphN].Step GT Sampling*2.0) then begin
+    *(Spec.Time) = indgen(Graph[GraphN].ppp)*Graph[GraphN].Step + Graph[GraphN].Step/2.0
+    *(Spec.Val)  =  replicate(FillValue, 16L*96L*16L*Graph[GraphN].ppp)
+    if (Flag eq GraphN) then begin
+	*(Tab.Time) = indgen(Graph[GraphN].ppp)*Graph[GraphN].Step + Graph[GraphN].Step/2.0
+	*(Tab.Val)  =  replicate(FillValue, Graph[GraphN].ppp)
+    endif
+  endif
+  
+  if (Flag eq GraphN) then cmd = Mission + 'OSpec = Spec & Etab' + Mission +'= Tab & Flag'+ Mission +'= Flag' $
+  else cmd = Mission + 'OSpec = Spec';
+    
+  R = execute(cmd);
+
+return
+end
+;------------------------------------------------------------
+
+++ a/amda_plus/amdalib/imacntplot.pro
@@ -0,0 +1,93 @@
+pro imacntplot, Mission, azMin, azMax, enMin, enMax, elevation
+
+  common GraphC, Graph, GraphN, Item
+  common Etab, FlagVex, EtabVex, FlagMex, EtabMex
+
+   ;------- Energy table -----------------
+    if (Mission eq 'vex') then Tab = fix(*EtabVex[0].val) else Tab = fix(*EtabMex[0].val);
+    EnNum = 96
+    Energy = fltarr(EnNum)
+     
+    RetSize = call_external(!PROJLIB + 'imainfo.so','GetImaEner_idl', Energy, Tab[0])
+    Energy = temporary(Energy(where(Energy GT 0)));
+    EnNum = N_elements(Energy)
+    
+  ;-------- Cat in the Energy range ------------
+   E0 =  min(Energy)
+   E1 =  max(Energy)
+   eindex = where((Energy GE E0) AND $
+                  (Energy LE E1) AND $
+                  (Energy GT 0.0),counte)
+   if counte GT 0 then begin
+      EBounds = fltarr(EnNum+1)
+      EBounds(min(eindex)) = Energy(min(eindex))
+      EBounds(max(eindex+1)) = Energy(max(eindex))
+      for ie = min(eindex+1), max(eindex) do EBounds(ie) = sqrt(Energy(ie-1)*Energy(ie))
+   endif else return
+
+
+  ;  FillValue = Graph[GraphN].FillValue[0]
+;--------------------- No DATA - just axes -----------------------
+   if (N_Elements(*Graph[GraphN].Time[0]) eq 0) then begin
+         nodata
+         return
+   endif
+
+    Time = (*Graph[GraphN].Time[0])
+
+    Val = (*Graph[GraphN].Val[0])
+    N_Time = N_elements(Time);
+
+    if (elevation) then begin 
+      Val = temporary(reform(Val, 16, 96, 16, N_Time))  
+      LocVal =  azMax GT azMin ? total(total(Val[azMin:azMax,*,0:15, *],1),2) : reform(total(Val[azMin,*,0:15,*],3), 96, N_time);   
+      newTime = Time;  
+      deltaT = 96.0;
+      gapT = deltaT*10.0 ;
+  endif else begin
+; elevation as TIME
+      newTime = dblarr(N_Time*16);
+      for i = 0, N_Time - 1 do for j = 0, 15 do newTime[i*16 + j] = Time[i] + j*12.D0
+      N_Time *= 16;
+      Val = temporary(reform(Val, 16, 96, N_Time));
+      LocVal =  azMax GT azMin ? total(Val[azMin:azMax,*,*],1) : reform(Val[azMin,*,*], 96, N_time);  
+      deltaT = 6.0;
+      gapT = deltaT*10.0;
+    endelse
+  
+    num0 = where(EBounds ge enMin, Index0);
+    num1 = where(EBounds ge enMax, Index1);
+    
+    if (Index1 eq 0)  then numStop = 0 else  numStop = num1[Index1-1];
+    if (Index0 eq EnNum+1) then numStart = 95 else numStart = num0[Index0-1];
+    
+    Data =  numStart ne numStop ? total(LocVal[numStop:numStart,*],1) : reform(LocVal[numStart,*], N_Time);   
+    numMinMax = where(finite(Data), IndexMinMax)
+    MinVal = min(Data[numMinMax]) > 0.1;   
+    MaxVal = max(Data[numMinMax]);
+   
+     
+   if Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] then begin
+                Graph[GraphN].Ly.range = Graph[GraphN].AxisF eq 'logaxis' ? [MinVal, MaxVal] : [0.0, MaxVal]
+                Graph[GraphN].Ly.style = 0
+    endif else Graph[GraphN].Ly.style = 1
+    
+    !y = Graph[GraphN].Ly;  
+   
+    if (execute(Graph[GraphN].AxisF) NE 1) then plot, newTime, Data, MAX_VALUE=MaxValue,  color = !dnc-1, /NODATA, /NOERASE 
+
+    plots, newTime[0]+deltaT, Data[0], /DATA
+ 
+    for i = 1, N_elements(newTime) - 1 do begin
+      if (newTime[i] LE Graph[GraphN].TotalTime) then begin
+          if (data[i] lt Graph[GraphN].Ly.range[0]) then data[i] = Graph[GraphN].Ly.range[0];
+          if (data[i] gt Graph[GraphN].Ly.range[1]) then data[i] = Graph[GraphN].Ly.range[1];  
+          if (newTime[i] - newTime[i-1] LE gapT) then $
+                    plots, newTime[i]+deltaT, data[i], color=col, /CONTINUE, /DATA $
+          else   plots, newTime[i]+deltaT, data[i], /DATA 
+        endif
+    endfor
+      
+return
+end
+;--------------------------------------------------------------------
+++ a/amda_plus/amdalib/imaspecplot.pro
@@ -0,0 +1,172 @@
+;============================================================
+;
+;=============================================================
+
+pro imaspecplot, Mission, AzMin, AzMax, elevation
+
+   common GraphC, Graph, GraphN, Item
+   common Etab, FlagVex, EtabVex, FlagMex, EtabMex
+
+   ;------- Energy table -----------------
+    if (Mission eq 'vex') then begin
+        numTab = where(finite(*EtabVex[0].val) eq 0, indexTab);
+        if (indexTab GT 0) then (*EtabVex[0].val)[numTab] = -1;
+        Tab = fix(*EtabVex[0].val); 
+    endif else begin
+        numTab = where(finite(*EtabMex[0].val) eq 0, indexTab);
+        if (indexTab GT 0) then (*EtabMex[0].val)[numTab] = -1;
+        Tab = fix(*EtabMex[0].val);  
+    endelse
+ 
+    EnNum0 = 96
+    Energy = fltarr(EnNum0);
+    Energy1 = fltarr(EnNum0);
+
+    numTab = uniq(Tab,sort(Tab));
+    uniqTab = Tab[numTab];
+; set min/max energy
+    EnergyMin = 10000.0;
+    EnergyMax = 0.0;
+    for i = 0, N_elements(uniqTab) - 1 do begin 
+        if uniqTab[i] lt 0 then continue;
+        RetSize = call_external(!PROJLIB + 'imainfo.so','GetImaEner_idl', Energy1,  uniqTab[i]); 
+        Energy1 =  temporary(Energy1(where(Energy1 GT 0)));
+        if (min(Energy1) LT EnergyMin) then EnergyMin = min(Energy1);
+        if (max(Energy1) GT EnergyMax) then EnergyMax = max(Energy1);
+    endfor 
+     
+    RetSize = call_external(!PROJLIB + 'imainfo.so','GetImaEner_idl', Energy, Tab[0])
+  
+    Energy =  temporary(Energy(where(Energy GT 0)));
+    EnNum = N_elements(Energy)
+   
+  ;-------- Cat in the Energy range ------------
+   E0 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ?  min(Energy) : Graph[GraphN].Ly.range[0]
+   E1 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ?  max(Energy) : Graph[GraphN].Ly.range[1]
+
+   eindex = where((Energy GE E0) AND $
+                  (Energy LE E1) AND $
+                  (Energy GT 0.0),counte)
+
+   if counte GT 0 then begin
+      EBounds = fltarr(EnNum+1)
+      EBounds(min(eindex)) = Energy(min(eindex))
+      EBounds(max(eindex+1)) = Energy(max(eindex))
+      for ie = min(eindex+1), max(eindex) do EBounds[ie] = sqrt(Energy[ie-1]*Energy[ie])
+   endif else begin
+        nodata
+        return
+   endelse
+
+;--------------------- No DATA - just axes -----------------------
+   if (N_Elements(*Graph[GraphN].Time[0]) eq 0) then begin
+         nodata
+         return
+   endif
+
+    Time = (*Graph[GraphN].Time[0])
+
+    Val = (*Graph[GraphN].Val[0])
+    N_Time = N_elements(Time);
+ 
+   if (elevation) then begin
+      Val = temporary(reform(Val, 16, 96, 16, N_Time)) 
+      LocVal =  azMax GT azMin ? total(total(Val[azMin:azMax,*,0:15, *],1),2) : reform(total(Val[azMin,*,0:15,*],3), 96, N_time);
+      newTime = Time;       
+  endif else begin
+; elevation as TIME
+      newTime = dblarr(N_Time*16);
+      for i = 0, N_Time - 1 do for j = 0, 15 do newTime[i*16 + j] = Time[i] + j*12.D0
+      N_Time *= 16;
+      Val = temporary(reform(Val,16, 96, N_Time));
+      LocVal =  azMax GT azMin ? total(Val[azMin:azMax,*,*],1) : reform(Val[azMin,*,*], 96, N_time);  
+    endelse
+      
+   numMinMax = where(LocVal gt 0.5, IndexMinMax)   
+ 
+  if (Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1]) then  $
+                                      Graph[GraphN].Ly.range = [EnergyMin, EnergyMax]
+       
+   Graph[GraphN].Ly.style = 1
+
+  !y = Graph[GraphN].Ly
+
+   plot_io, graph[GraphN].Lx.range, graph[GraphN].Ly.range, /nodata, /noerase, ystyle = 1, $
+                              ytickformat='yticks', charsize=graph[GraphN].Lp.charsize*0.8, color = !dnc-1
+
+   Step = Graph[GraphN].Sampling[0] > Graph[GraphN].Step;  
+   if (not elevation) then Step /= 16.0;    
+
+ ;--------------------- No DATA - just axes -----------------------
+  if (IndexMinMax LE 0) then begin
+         nodatafortimeinterval
+         return
+  endif
+  if (Time[0] GT Graph[GraphN].TotalTime) then begin
+        nodatafortimeinterval
+         return
+  endif
+  
+  
+  MinVal = min(LocVal[numMinMax]) > 0.5
+  MaxVal = max(LocVal[numMinMax])   
+ 
+ ;-----------   NOW PLOT SPECTRA
+    alogMin = alog(MinVal)
+    KF = float(!DNC - 2) /(alog(MaxVal) - alogMin)
+    xx = fltarr(4) & yy = fltarr(4)
+ 
+     for i = 0, N_Time - 2 do begin
+
+	ii = elevation ? i : i/16 < (N_elements(Tab) - 2);
+ 
+	if (Tab[ii+1] ne Tab[ii] AND Tab[ii+1] GE 0) then begin
+              Energy = fltarr(EnNum0)  
+	      RetSize = call_external(!PROJLIB + 'imainfo.so','GetImaEner_idl', Energy, Tab[ii+1])
+	      Energy =  temporary(Energy(where(Energy GT 0)));
+	      EnNum = N_elements(Energy)
+	      
+	    ;-------- Cat in the Energy range ------------
+	    E0 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ? min(Energy) : Graph[GraphN].Ly.range[0]
+	    E1 = Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1] ? max(Energy) : Graph[GraphN].Ly.range[1]
+	    eindex = where((Energy GE E0) AND $
+			    (Energy LE E1) AND $
+			    (Energy GT 0.0),counte)
+	    if counte GT 0 then begin
+		EBounds = fltarr(EnNum+1)
+		EBounds(min(eindex)) = Energy(min(eindex))
+		EBounds(max(eindex+1)) = Energy(max(eindex))
+		for ie = min(eindex+1), max(eindex) do EBounds(ie) = sqrt(Energy(ie-1)*Energy(ie))
+	    endif else return
+	endif
+        xx[0] = newTime[i]  > 0.0
+        xx[1] = xx[0]
+        xx[2] = newTime[i+1] < newTime[i] + Step
+        xx[3] = xx[2]
+
+        if (xx[2] LE graph[GraphN].TotalTime AND xx[0] LT graph[GraphN].TotalTime) then begin
+           for k = min(eindex), max(eindex) do begin
+                  yy[0] = EBounds[k]
+                  yy[1] = EBounds[k+1]
+                  yy[2] = yy[1]
+                  yy[3] = yy[0]
+                  lcount = LocVal[k,i] < 1.e5
+                  if finite(lcount) then begin
+                   cl = lcount GT MinVal ? fix(KF*(alog(lcount) - alogMin)) : 1;
+                  endif else cl = 0;
+                  
+                  if (cl GT !DNC - 2) then  cl = !DNC - 2
+                  PolyFill, xx ,yy, color = cl, /data 
+            endfor
+        endif
+     endfor
+
+;-------------------- Legend ----------------------------------------
+      spec_legend, MinVal, MaxVal
+;------------------ MIN-MAX------------------------------------------------
+      Graph[GraphN].Min  = MinVal
+      Graph[GraphN].Max  = MaxVal
+
+return
+end
+;--------------------------------------------------------------------
+++ a/amda_plus/amdalib/legend_multi.pro
@@ -0,0 +1,11 @@
+pro legend_multi, Names, Colors
+common GraphC, Graph, GraphN, Item
+
+    xy = Graph[GraphN].lp.region
+    y0 = xy[1] + (xy[3] - xy[1])/2. + 0.02*(N_elements(Names)-1)/2
+ 
+   for i = 0, N_elements(Names)-1 do $
+      xyouts, xy[2]+0.02, y0-0.03*i ,Names[i], /NORMAL, charsize=graph[GraphN].Lp.charsize*0.8, color=Colors[i]    
+ 
+return
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/lep_g_0_init.pro
@@ -0,0 +1,24 @@
+;-----------------GTLN_INIT--------------------------
+; 
+pro lep_g_0_INIT, LIN = LIN, LOG = LOG,  RED = RED, BLUE = BLUE, GREEN = GREEN
+
+common GraphC, graph, GraphN, Item
+ 
+ pen = 0
+
+ if KEYWORD_SET(RED) then pen = 3
+ if KEYWORD_SET(BLUE) then pen = 1  
+ if KEYWORD_SET(GREEN) then pen = 2   
+ 
+graph[GraphN].Type = 'HIAM_T' 
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'gtlcom' 
+graph[GraphN].GetDataF[0] ='getgtl' 
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime+30.0D0) <  18000.0D0
+graph[GraphN].Ly.title = 'GTL !C N, cm!U-3'
+graph[GraphN].AxisF =  KEYWORD_SET(LIN) ?  'lineaxis,0' : 'logaxis'
+graph[GraphN].PlotF = 'plotscal, 3,'+ string(pen)
+graph[GraphN].LastPlotF = 'nodatalast'
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/lep_g_1_3_init.pro
@@ -0,0 +1,42 @@
+;-----------------GTLV_INIT--------------------------
+; ---------------------------------------------------------
+pro lep_g_1_3_INIT, X = X, Y = Y, Z = Z, TOT = TOT, GSE = GSE, GSM = GSM, SM = SM, RED = RED, GREEN = GREEN, BLUE = BLUE
+
+common GraphC, graph, GraphN, Item
+  
+   names = ['V!DX', 'V!DY', 'V!DZ', 'V']
+   comp = 4 & pen = 0
+   
+   if KEYWORD_SET(X) then comp = 0
+   if KEYWORD_SET(Y) then comp = 1  
+   if KEYWORD_SET(Z) then comp = 2  
+   if KEYWORD_SET(TOT) then comp = 3
+    if KEYWORD_SET(BLUE) then pen = 1
+    if KEYWORD_SET(GREEN) then pen = 2
+    if KEYWORD_SET(RED) then pen = 3
+
+frame = 'GSE'
+ 
+if KEYWORD_SET(GSM) then begin
+  frame = 'GSM'
+  graph[GraphN].ProcessF='gse2gsm'
+endif
+if KEYWORD_SET(SM) then begin
+  frame = 'SM'
+  graph[GraphN].ProcessF='gse2sm'
+endif 
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP) 
+graph[GraphN].Type = 'HIAM_T'
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'gtlcom' 
+graph[GraphN].GetDataF[0] ='getgtl'
+graph[GraphN].DeltaT =  (Graph[GraphN].TotalTime+30.0D0) <  18000.0D0
+graph[GraphN].Ly.title =  'GTL !CV '+frame+', km/s'
+graph[GraphN].AxisF = 'lineaxis,' + string(Graph[GraphN].TotalTime)
+*(graph[GraphN].names) = (comp le 3) ? names : names[0:2]
+graph[GraphN].PlotF = (comp le 3) ? 'plotvec_comp_, '+string(comp)+',' +string(pen) : 'plotvec_'
+graph[GraphN].LastPlotF = 'nodatalast'
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/lep_g_1_init.pro
@@ -0,0 +1,17 @@
+;----------------ENVELOPE-------------------------
+; 
+pro lep_g_1_init, GSE = GSE, GSM = GSM, SM = SM, RED = RED, GREEN = GREEN, BLUE = BLUE
+
+  cmd = 'lep_g_1_3_init, /X'
+  if KEYWORD_SET(GSM) then cmd += ',/GSM'
+  if KEYWORD_SET(SM) then cmd += ',/SM'
+
+  if KEYWORD_SET(BLUE) then cmd += ',/BLUE'
+  if KEYWORD_SET(RED) then cmd += ',/RED'
+  if KEYWORD_SET(GREEN) then cmd += ',/GREEN'
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/lep_g_2_init.pro
@@ -0,0 +1,17 @@
+;----------------ENVELOPE-------------------------
+; 
+pro lep_g_2_init, GSE = GSE, GSM = GSM, SM = SM, RED = RED, GREEN = GREEN, BLUE = BLUE
+
+  cmd = 'lep_g_1_3_init, /Y'
+  if KEYWORD_SET(GSM) then cmd += ',/GSM'
+  if KEYWORD_SET(SM) then cmd += ',/SM'
+
+  if KEYWORD_SET(BLUE) then cmd += ',/BLUE'
+  if KEYWORD_SET(RED) then cmd += ',/RED'
+  if KEYWORD_SET(GREEN) then cmd += ',/GREEN'
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/lep_g_3_init.pro
@@ -0,0 +1,17 @@
+;----------------ENVELOPE-------------------------
+; 
+pro lep_g_3_init, GSE = GSE, GSM = GSM, SM = SM, RED = RED, GREEN = GREEN, BLUE = BLUE
+
+  cmd = 'lep_g_1_3_init, /Z'
+  if KEYWORD_SET(GSM) then cmd += ',/GSM'
+  if KEYWORD_SET(SM) then cmd += ',/SM'
+
+  if KEYWORD_SET(BLUE) then cmd += ',/BLUE'
+  if KEYWORD_SET(RED) then cmd += ',/RED'
+  if KEYWORD_SET(GREEN) then cmd += ',/GREEN'
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/lep_g_4_init.pro
@@ -0,0 +1,15 @@
+;----------------ENVELOPE-------------------------
+; 
+pro lep_g_4_init, GSE = GSE, GSM = GSM, SM = SM, RED = RED, GREEN = GREEN, BLUE = BLUE
+
+  cmd = 'lep_g_1_3_init, /TOT'
+  
+  if KEYWORD_SET(BLUE) then cmd += ',/BLUE'
+  if KEYWORD_SET(RED) then cmd += ',/RED'
+  if KEYWORD_SET(GREEN) then cmd += ',/GREEN'
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/lep_g_5_6_init.pro
@@ -0,0 +1,29 @@
+;-----------------GTLT_INIT--------------------------
+;         14 Apr, 2004 V2.5
+;-------------------------------------------------------
+pro lep_g_5_6_INIT, PAR=PAR, PERP=PERP, RED = RED, GREEN = GREEN, BLUE = BLUE, LIN = LIN, LOG = LOG
+
+common GraphC, graph, GraphN, Item
+names = ['T!DYY', 'T!DZZ']
+comp = 2 & pen = 0
+ 
+ if KEYWORD_SET(BLUE) then pen = 1
+ if KEYWORD_SET(GREEN) then pen = 2
+ if KEYWORD_SET(RED) then pen = 3
+ if KEYWORD_SET(PAR) then comp = 0
+ if KEYWORD_SET(PERP) then comp = 1
+
+graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP) 
+graph[GraphN].Type = 'HIAM_T' 
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'gtlcom' 
+graph[GraphN].GetDataF[0] = 'getgtl'
+graph[GraphN].DeltaT =  (Graph[GraphN].TotalTime+30.0D0)  <  18000.0D0
+graph[GraphN].Ly.title = 'GTL !C T, eV'
+graph[GraphN].AxisF = 'logaxis'
+*(graph[GraphN].names) = names 
+graph[GraphN].PlotF = (comp le 1) ? 'plotten_comp, 4, 2,'+string(comp)+',' +string(pen) : 'plotten, 4, 2'
+graph(GraphN).LastPlotF = 'nodatalast'
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/lep_g_5_init.pro
@@ -0,0 +1,15 @@
+;----------------ENVELOPE-------------------------
+; 
+pro lep_g_5_init,  RED = RED, GREEN = GREEN, BLUE = BLUE, LIN = LIN, LOG = LOG
+
+  cmd = 'lep_g_5_6_init, /PAR'
+ 
+  if KEYWORD_SET(BLUE) then cmd += ',/BLUE'
+  if KEYWORD_SET(RED) then cmd += ',/RED'
+  if KEYWORD_SET(GREEN) then cmd += ',/GREEN'
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/lep_g_6_init.pro
@@ -0,0 +1,15 @@
+;----------------ENVELOPE-------------------------
+; 
+pro lep_g_6_init,  RED = RED, GREEN = GREEN, BLUE = BLUE, LIN = LIN, LOG = LOG
+
+  cmd = 'lep_g_5_6_init, /PERP'
+ 
+  if KEYWORD_SET(BLUE) then cmd += ',/BLUE'
+  if KEYWORD_SET(RED) then cmd += ',/RED'
+  if KEYWORD_SET(GREEN) then cmd += ',/GREEN'
+
+  R = execute(cmd)
+
+return
+end
+;------------------------------------------------------------
 \ No newline at end of file
+++ a/amda_plus/amdalib/median_.pro
@@ -0,0 +1,18 @@
+function median_, Data, T_AVER
+
+common  InternalCommon, FinalTime, num
+
+     Npoints = N_elements(FinalTime);
+     Nreal = n_elements(num);
+     Timy = FinalTime(num);
+ 
+     FillValue = !Values.F_NAN
+     outPut = replicate(FillValue, Nreal); outPut = fltarr(Nreal);
+ 
+     for  tt = FinalTime[0], FinalTime[Npoints-1]+T_AVER, T_AVER do begin
+         nummi = where(abs(tt - Timy) LE T_AVER/2.0, Index)  
+         if (Index GT 0) then OutPut[nummi] = median(Data[nummi])
+      endfor
+
+return, OutPut
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/mercorbaxis.pro
@@ -0,0 +1,110 @@
+;==================================================
+;                 DD MERCURY
+;              ORBITAL OBJECTS
+;              mercorbaxis.pro
+;                    
+;---------------------------------------------------
+; Implementation of all types of axis for orbital plot
+; Options:
+;   Key                  ToPlot
+; /CYL - cylindrical XR   0
+; /XY                     1
+; /XZ                     2
+; /YZ                     3
+;====================================================
+
+pro mercorbaxis, ToPlot
+common GraphC, graph, GraphN
+
+;------------ Set graphics Colors ------------------------
+ if (Graph[GraphN].Lx.range[0] eq Graph[GraphN].Lx.range[1]) then begin
+       Graph[GraphN].Lx.range[0] = -5;
+       Graph[GraphN].Lx.range[1] = 5; 
+ endif
+ if (Graph[GraphN].Ly.range[0] eq Graph[GraphN].Ly.range[1]) then begin
+       Graph[GraphN].Ly.range[0] = ToPlot eq 0 ? 0: -5;
+       Graph[GraphN].Ly.range[1] = 5;
+ endif
+ 
+ !x = Graph[GraphN].Lx
+ !y = Graph[GraphN].Ly
+
+  x = !x.range
+  y = !y.range
+  black = !DNC- 1
+  zero = fix(!DNC * 20.0/ 256.0 )
+  planet = fix(!DNC * 10.0/ 256.0)
+  bound =  fix(!DNC * 1.0/ 256.0)
+;------------------------------------------------------------
+
+  ThickPlanet = 2 
+  ThickBound =  2
+
+;------- main axis -----------------
+  case ToPlot of
+  0: begin
+      !x.title = 'X, R!DM'
+      !y.title = 'R, R!DM'
+     end
+  1: begin
+      !x.title = 'X, R!DM'
+      !y.title = 'Y, R!DM'
+     end
+  2: begin
+      !x.title = 'X, R!DM'
+      !y.title = 'Z, R!DM'
+     end
+  3: begin
+      !x.title = 'Y, R!DM'
+      !y.title = 'Z, R!DM'
+     end
+  endcase
+  
+  ;------ Plot main axis. Use common Graph system variables
+  plot,x,y,/nodata,/noerase,xstyle=1,ystyle=1,color=black,xminor=1, yminor=1, charsize=graph[GraphN].Lp.charsize
+     
+  ;------- plot Zero lines --------------------- 
+  oplot,x,[0,0],color=zero, thick = 1
+  oplot,[0,0],y,color=zero, thick = 1
+  
+  ;------- plot Mars ---------------------
+  dal = 2.0 * !PI / 60.
+  if ToPlot EQ 0 then $ 
+       al = dal * findgen(31) $
+  else al = dal * findgen(61) 
+    
+  xx = cos(al)
+  yy = sin(al)
+  oplot,xx,yy,color=black, thick = ThickPlanet
+    
+  if ToPlot EQ 0 then begin
+     ;---- Plot boundaries for cylindrical plot
+     ;----------------------- Bow Shock -----------------------------
+   theta = (135.*!pi/180.)*findgen(100)/100.
+
+   rr = 2.96/(1.0+1.02*cos(theta))
+   xx = fltarr(100)
+   yy = fltarr(100)
+
+   for i = 0, 99 do begin
+     xx(i) =  0.5 + rr(i)*cos(theta(i))
+     yy(i) =  rr(i)*sin(theta(i))
+   endfor
+
+   oplot,xx,yy,color=bound,  thick = ThickBound
+   oplot,xx,-yy,color=bound, thick = ThickBound
+
+    ;----------------------- MP Shue-----------------------------
+    alpha = 0.5;
+    r0 = 1.45;
+
+    rr=r0*(2./(1.+COS(theta)))^alpha
+    
+     xx = rr*cos(theta)
+     yy = rr*sin(theta)  
+       
+     oplot,xx,rr,color=bound,  thick = ThickBound
+     oplot,xx,-rr,color=bound, thick = ThickBound
+  endif
+return
+end
+++ a/amda_plus/amdalib/mex_els_spec_init.pro
@@ -0,0 +1,61 @@
+;==========================================================
+;                DD DDMARS/DDVEX
+;                elspec_init.pro
+;-----------------------------------------------------------
+
+pro mex_els_spec_init, ANODE = Anode, Stop_Anode = Stop_Anode, $
+                      LOG = log, LIN = lin, PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $
+                     ; spectra
+                     MINCOUNT = minCount, MAXCOUNT = maxCount, BACKGRDCHAN = backgrdChan, BACKGRDVAL = backgrdVal, $
+                     CHANNEL = channel, Energy = energy, ENERGYMIN = energymin, ENERGYMAX =  energymax, $
+                     
+                     DELAY = delay, $                      
+                     ; all additional keywords
+                    _EXTRA = e 
+
+common Request, List, StartTime, TimeInt
+common GraphC, graph, GraphN, Item
+
+AnMin = 0
+AnMax = 15
+ 
+EnMin = 1.0;
+EnMax = 2800.0;
+
+if (KEYWORD_SET(ANODE)) then begin
+   AnMin = Anode
+   if (KEYWORD_SET(Stop_Anode)) then AnMax = Stop_Anode $
+   else AnMax = -1;
+endif
+      
+graph[GraphN].Type = 'ELS_T'
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'elscom, "mex"'
+graph[GraphN].GetDataF[0] = 'getels, "mex"'
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) < 36100.0D0
+energyTitle = '';
+ 
+if (KEYWORD_SET(energymin)) then begin
+   EnMin = energymin;
+   EnMax =  energymax;  
+   energyTitle = energymin eq energymax ? '!C E ' + string(energymin) + 'eV' : '!C E ' + string(energymin) + '-'+  string(energymax) + 'eV'; 
+      
+      options = '"mex", ' + string(AnMin) +','+string(AnMax) +','  +string(EnMin) +','+string(EnMax);
+      angleTitle =  AnMax GT AnMin ? 'MEX ELS An'+ string(AnMin,format='(i2.2)') +'-'+ string(AnMax,format='(i2.2)') : 'MEX ELS An'+ string(AnMin,format='(i2.2)');
+ 
+      graph[GraphN].AxisF =  KEYWORD_SET(lin) ?  'lineaxis, 0' : 'logaxis';
+      graph[GraphN].Ly.title =  'MEX ELS Cnts ' +  angleTitle + energyTitle					    
+      graph[GraphN].PlotF = 'elscntplot,' + options
+endif else begin 
+      options = '"mex", ' + string(AnMin) +','+string(AnMax);
+
+      graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+      *(graph[GraphN].names) = 'log Counts'
+ 
+      graph[GraphN].Ly.title = AnMax GT AnMin ? 'MEX ELS An'+ string(AnMin,format='(i2.2)') +'-'+ string(AnMax,format='(i2.2)') +'!C E, eV' : 'MEX ELS An'+ string(AnMin,format='(i2.2)') + '!C E, eV'
+      graph[GraphN].PlotF = 'elsspecplot,' + options
+ 
+endelse
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/mex_h_dens_.pro
@@ -0,0 +1,54 @@
+function MEX_H_DENS_, StartTime, TimeInt, Time, RetSize, Sampling, close=close
+common GraphC, graph, GraphN
+common MEXHDENS_COM, ID
+
+   FillValue = !Values.F_NAN 
+   if n_elements(ID) EQ 0 then ID = -1L
+   sampling = 192.0
+   Size = long(graph[GraphN].DeltaT/sampling)
+     
+   Bdata = fltarr(Size)
+   Time = dblarr(Size)
+   Quality = fltarr(Size)
+   Time_Q = dblarr(Size)
+   
+   Mission ="mex";
+ 
+    VI = Mission + ':ima:param'
+    ParamName = "Density_P";
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, StartTime, TimeInt, Size, Time, Bdata)
+   
+    ParamName = "Quality_P";
+
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, StartTime, TimeInt, Size, Time_Q, Quality)
+                          
+   if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+   endif                        
+  
+   if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, 0.0
+   endif
+   
+   Time=Time[0:Retsize-1]
+   Bdata=Bdata[0:Retsize-1]
+   Quality = Quality[0:Retsize-1]
+ ;  num_gap = where(Quality GE 0.7, Cntr)
+  
+ ;  if Cntr le 0 then begin
+ ;      Retsize = -999
+ ;      Time = 0.D0
+  ;     return, 0.0
+  ; endif
+   
+ ;   Time=Time[num_gap]
+  ;  Bdata=Bdata[num_gap];
+  ;  Retsize=Cntr
+
+return, Bdata
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/mex_h_dens_init.pro
@@ -0,0 +1,24 @@
+;--------------------------------------
+;
+;---------------------------------------
+
+pro mex_h_dens_INIT,  LIN = LIN, LOG = LOG, COLOR = color, PSYM = psym 
+
+common GraphC, graph, GraphN, Item
+
+names = ['Poor', 'Satisfactory', 'Good', 'Excellent']
+
+Graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) =  names;
+graph[GraphN].Type = 'FGM_T' 
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'vexmex_dens_com, "mex", "h"' 
+
+graph[GraphN].GetDataF[0] ='get_vexmex_dens, "mex", "h"' 
+
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) <  36100.0D0
+graph[GraphN].Ly.title = 'MEX ima-H N, cm!U-3'
+graph[GraphN].AxisF =  KEYWORD_SET(LIN) ?  'lineaxis, 0' : 'logaxis'
+graph[GraphN].PlotF = 'PlotImaDens'
+return
+end
+++ a/amda_plus/amdalib/mex_h_dens_restr_.pro
@@ -0,0 +1,52 @@
+function MEX_H_DENS_RESTR_, StartTime, TimeInt, Time, RetSize, Sampling, close=close
+common GraphC, graph, GraphN
+common MEXHDENS_COM, ID
+
+   FillValue = !Values.F_NAN 
+   if n_elements(ID) EQ 0 then ID = -1L
+   sampling = 192.0
+   Size = long(graph[GraphN].DeltaT/sampling)
+     
+   Bdata = fltarr(Size)
+   Time = dblarr(Size)
+   Quality = fltarr(Size)
+   Time_Q = dblarr(Size)
+    
+    VI = 'mex:imaparam:restricted'
+    ParamName = "Density_P";
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, StartTime, TimeInt, Size, Time, Bdata)
+   
+    ParamName = "Quality_P";
+
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, StartTime, TimeInt, Size, Time_Q, Quality)
+                          
+   if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+   endif                        
+  
+   if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, 0.0
+   endif
+   
+   Time=Time[0:Retsize-1]
+   Bdata=Bdata[0:Retsize-1]
+   Quality = Quality[0:Retsize-1]
+ ;  num_gap = where(Quality GE 0.7, Cntr)
+  
+ ;  if Cntr le 0 then begin
+ ;      Retsize = -999
+ ;      Time = 0.D0
+  ;     return, 0.0
+  ; endif
+   
+ ;   Time=Time[num_gap]
+  ;  Bdata=Bdata[num_gap];
+  ;  Retsize=Cntr
+
+return, Bdata
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/mex_h_dens_restr_init.pro
@@ -0,0 +1,24 @@
+;--------------------------------------
+;
+;---------------------------------------
+
+pro mex_h_dens_restr_INIT,  LIN = LIN, LOG = LOG, COLOR = color, PSYM = psym 
+
+common GraphC, graph, GraphN, Item
+
+names = ['Poor', 'Satisfactory', 'Good', 'Excellent']
+
+Graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+*(graph[GraphN].names) =  names;
+graph[GraphN].Type = 'FGM_T' 
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'vexmex_dens_com, "mex", "h"' 
+
+graph[GraphN].GetDataF[0] = 'get_vexmex_dens, "mex", "h", /RESTRICTED' 
+
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) <  36100.0D0
+graph[GraphN].Ly.title = 'MEX ima-H N, cm!U-3'
+graph[GraphN].AxisF =  KEYWORD_SET(LIN) ?  'lineaxis, 0' : 'logaxis'
+graph[GraphN].PlotF = 'PlotImaDens'
+return
+end
+++ a/amda_plus/amdalib/mex_h_qual_.pro
@@ -0,0 +1,49 @@
+function MEX_H_QUAL_, StartTime, TimeInt, Time, RetSize, Sampling, close=close
+common GraphC, graph, GraphN
+common MEXHQUAL_COM, ID
+
+   FillValue = !Values.F_NAN 
+   if n_elements(ID) EQ 0 then ID = -1L
+   sampling = 192.0
+   Size = long(graph[GraphN].DeltaT/sampling)
+     
+   Bdata=fltarr(Size)
+   Time=dblarr(Size)
+   
+   Mission ="mex";
+ 
+    VI = Mission + ':ima:param'
+    ParamName = "Quality_P";
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, StartTime, TimeInt, Size, Time, Bdata)
+                          
+   if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+   endif                       
+  
+   if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, 0.0
+   endif
+   
+   Time=Time[0:Retsize-1]
+   Bdata=Bdata[0:Retsize-1]
+   
+   num_gap = where(finite(Bdata), Cntr)
+  
+   if Cntr le 0 then begin
+       Retsize = -999
+       Time = 0.D0
+       return, 0.0
+   endif
+   
+    Time=Time[num_gap]
+    Bdata=Bdata[num_gap];
+    Retsize=Cntr
+
+   
+   
+return, Bdata
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/mex_h_qual_init.pro
@@ -0,0 +1,26 @@
+;--------------------------------------
+;
+;---------------------------------------
+
+pro mex_h_qual_INIT, $
+                     LOG = log, LIN = lin, PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                   
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $
+                     
+                     ; all additional keywords
+                    _EXTRA = e   
+
+common GraphC, graph, GraphN, Item
+
+
+graph[GraphN].Type = 'FGM_T' 
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'vexmex_qual_com, "mex", "h"' 
+graph[GraphN].GetDataF[0] ='get_vexmex_qual, "mex", "h"'
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) <  18000.0D0 
+graph[GraphN].Ly.title = 'MEX ima-H Quality'
+graph[GraphN].AxisF = 'lineaxis, 0'
+graph[GraphN].PlotF = 'PlotImaQual'
+return
+end
+++ a/amda_plus/amdalib/mex_h_qual_restr_.pro
@@ -0,0 +1,47 @@
+function MEX_H_QUAL_RESTR_, StartTime, TimeInt, Time, RetSize, Sampling, close=close
+common GraphC, graph, GraphN
+common MEXHQUAL_COM, ID
+
+   FillValue = !Values.F_NAN 
+   if n_elements(ID) EQ 0 then ID = -1L
+   sampling = 192.0
+   Size = long(graph[GraphN].DeltaT/sampling)
+     
+   Bdata=fltarr(Size)
+   Time=dblarr(Size)
+   
+    VI = 'mex:imaparam:restricted'
+    ParamName = "Quality_P";
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, StartTime, TimeInt, Size, Time, Bdata)
+                          
+   if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+   endif                       
+  
+   if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, 0.0
+   endif
+   
+   Time=Time[0:Retsize-1]
+   Bdata=Bdata[0:Retsize-1]
+   
+   num_gap = where(finite(Bdata), Cntr)
+  
+   if Cntr le 0 then begin
+       Retsize = -999
+       Time = 0.D0
+       return, 0.0
+   endif
+   
+    Time=Time[num_gap]
+    Bdata=Bdata[num_gap];
+    Retsize=Cntr
+
+   
+   
+return, Bdata
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/mex_h_qual_restr_init.pro
@@ -0,0 +1,26 @@
+;--------------------------------------
+;
+;---------------------------------------
+
+pro mex_h_qual_restr_INIT, $
+                     LOG = log, LIN = lin, PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                   
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $
+                     
+                     ; all additional keywords
+                    _EXTRA = e   
+
+common GraphC, graph, GraphN, Item
+
+
+graph[GraphN].Type = 'FGM_T' 
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'vexmex_qual_com, "mex", "h"' 
+graph[GraphN].GetDataF[0] ='get_vexmex_qual, "mex", "h", /RESTRICTED'
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) <  18000.0D0 
+graph[GraphN].Ly.title = 'MEX ima-H Quality'
+graph[GraphN].AxisF = 'lineaxis, 0'
+graph[GraphN].PlotF = 'PlotImaQual'
+return
+end
+++ a/amda_plus/amdalib/mex_h_spec_.pro
@@ -0,0 +1,63 @@
+function mex_h_spec_, StartTime, TimeInt, Time, RetSize, Sampling, close=close
+common GraphC, graph, GraphN
+common MEX_H_COM, ID, Energy
+
+FORWARD_FUNCTION gfl
+
+   FillValue =  !Values.F_NAN
+   
+   if n_elements(ID) EQ 0 then  begin
+     ID = -1L
+     err = call_external(!PROJLIB + 'libplanetcom.so','SetMission_idl', 'mex');
+     
+     EnNum = 96
+     Energy = fltarr(EnNum)
+     err = call_external(!PROJLIB + 'imainfo.so','GetImaEner_idl', Energy, 0)
+     Energy =  temporary(Energy(where(Energy GT 0)));
+   endif
+   
+   sampling = 192.0;
+  
+   Size = long(graph[GraphN].DeltaT/sampling)
+
+   Time=dblarr(Size)
+   SpectraTot = fltarr(16, 96, 16, 2, Size)
+   Product = 'HPSPEC';
+   PaccArr = intarr(Size);
+   EnergyTab = intarr(Size);
+
+    RetSize = call_external(!PROJLIB + 'libImaExtra.so','GetImaExtra_idl',$
+                                ID, StartTime, TimeInt,$
+                                Product, Size, $
+                                Time, PaccArr, EnergyTab, SpectraTot)
+                                
+  
+   if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+   endif                              
+   
+   if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, replicate(0.0, 96)
+   endif
+ 
+   Spectra = fltarr(16, 96, 16, RetSize);
+   
+         for i = 0, RetSize - 1 do $
+            for j = 0, N_elements(Energy) - 1 do begin
+               GFeff = gfl(Energy[j], PaccArr[i], 1.0);
+               GFeffG = gfl(Energy[j], PaccArr[i], 0.0); Ghost
+                 
+               coeffH = GFeff GT 0 ? 1.e-5/GFeff : 0.0;
+               coeffG = GFeffG GT 0 ? 1.e-5/GFeffG : 0.0;
+             
+               Spectra[*,j,*,i] = SpectraTot[*,j,*,0,i]*coeffH + SpectraTot[*,j,*,1,i] * coeffG;
+         endfor 
+
+   
+   SpectraOut = total(total(Spectra,1),2);
+   Time=Time[0:Retsize-1]
+   
+return, SpectraOut
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/mex_h_spec_init.pro
@@ -0,0 +1,71 @@
+;====================================================
+;              VEX_h_spec_init
+;                 1.0
+;   Init module of "h_spec" object
+;   HPSPEC from IMAEXTRA product VEX
+;
+;              18 Mars 2009
+;----------------------------------------------------
+
+pro mex_h_spec_init, ANODE = Anode, STOP_ANODE = Stop_Anode, elevationAsTime = elevationAsTime, $
+                     LOG = log, LIN = lin, PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $
+                     ; spectra
+                     MINCOUNT = minCount, MAXCOUNT = maxCount, BACKGRDCHAN = backgrdChan, BACKGRDVAL = backgrdVal, $
+                     CHANNEL = channel, Energy = energy, ENERGYMIN = energymin, ENERGYMAX =  energymax, $
+                     
+                     DELAY = delay, $                      
+                     ; all additional keywords
+                    _EXTRA = e 
+
+  common Request, List, StartTime, TimeInt
+  common GraphC, graph, GraphN, Item
+
+  AzMin = 0
+  AzMax = 15
+
+   elevation = 1;
+   if (KEYWORD_SET(elevationAsTime)) then elevation = 0;
+
+  if (KEYWORD_SET(ANODE)) then begin
+    AzMin = Anode
+   if (KEYWORD_SET(STOP_ANODE)) then AzMax = Stop_Anode $
+   else AzMax = -1;
+  endif
+    
+  graph(GraphN).Type = 'IMA_T'
+  graph(GraphN).DataN = 1
+  graph(GraphN).SetCommonF(0) = 'ima_h_spec_com, "mex"'
+  graph(GraphN).GetDataF(0) = 'get_imam_h_spec'
+  graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) < 36100.0D0
+
+EnMin = 1.0;
+EnMax = 25000.0;
+
+energyTitle = '';
+
+if (KEYWORD_SET(ENERGYMIN)) then begin
+   EnMin = energymin;
+   EnMax =  energymax;
+   energyTitle = energymin eq energymax ? '!C E ' + string(energymin,format='(f7.1)') + 'eV' : '!C E ' + string(energymin,format='(f7.1)') + '-'+  string(energymax,format='(f7.1)') + 'eV'; 
+    
+    options = '"mex", ' + string(AzMin) +','+string(AzMax) +','+string(EnMin) +','+string(EnMax)+ ',' + string(elevation);;
+
+    angleTitle =   AzMax GT AzMin ?  'MEX IMA H+ An'+ string(AzMin,format='(i2.2)') +'-'+ string(AzMax,format='(i2.2)') : 'MEX IMA H+ An'+ string(AzMin,format='(i2.2)');
+    
+    graph[GraphN].AxisF =  KEYWORD_SET(lin) ?  'lineaxis, 0' : 'logaxis';
+    graph[GraphN].Ly.title =  'MEX H Cnts ' +  angleTitle + energyTitle
+					  
+    graph[GraphN].PlotF = 'imacntplot,' + options
+endif else begin
+    options = '"mex",'+string(AzMin) +','+string(AzMax) + ',' + string(elevation);
+    
+    graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+    *(graph[GraphN].names) = 'log Counts'
+    
+    
+    graph[GraphN].Ly.title = AzMax GT AzMin ?  'MEX IMA H+ An'+ string(AzMin,format='(i2.2)') +'-'+ string(AzMax,format='(i2.2)') +'!CE/Q, eV' : 'MEX IMA H+ An'+ string(AzMin,format='(i2.2)') + '!C E/Q, eV'
+    graph[GraphN].PlotF =  'imaspecplot,' + options
+endelse
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/mex_h_spec_restr_.pro
@@ -0,0 +1,63 @@
+function mex_h_spec_restr_, StartTime, TimeInt, Time, RetSize, Sampling, close=close
+common GraphC, graph, GraphN
+common MEX_H_COM, ID, Energy
+
+FORWARD_FUNCTION gfl
+
+   FillValue =  !Values.F_NAN
+   
+   if n_elements(ID) EQ 0 then  begin
+     ID = -1L
+     err = call_external(!PROJLIB + 'libplanetcom.so','SetMission_idl', 'mex');
+     err = call_external(!PROJLIB + 'libplanetcom.so','SetRestriction');
+     EnNum = 96
+     Energy = fltarr(EnNum)
+     err = call_external(!PROJLIB + 'imainfo.so','GetImaEner_idl', Energy, 0)
+     Energy =  temporary(Energy(where(Energy GT 0)));
+   endif
+   
+   sampling = 192.0;
+  
+   Size = long(graph[GraphN].DeltaT/sampling)
+
+   Time=dblarr(Size)
+   SpectraTot = fltarr(16, 96, 16, 2, Size)
+   Product = 'HPSPEC';
+   PaccArr = intarr(Size);
+   EnergyTab = intarr(Size);
+
+    RetSize = call_external(!PROJLIB + 'libImaExtra.so','GetImaExtra_idl',$
+                                ID, StartTime, TimeInt,$
+                                Product, Size, $
+                                Time, PaccArr, EnergyTab, SpectraTot)
+                                
+  
+   if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+   endif                              
+   
+   if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, replicate(0.0, 96)
+   endif
+ 
+   Spectra = fltarr(16, 96, 16, RetSize);
+   
+         for i = 0, RetSize - 1 do $
+            for j = 0, N_elements(Energy) - 1 do begin
+               GFeff = gfl(Energy[j], PaccArr[i], 1.0);
+               GFeffG = gfl(Energy[j], PaccArr[i], 0.0); Ghost
+                 
+               coeffH = GFeff GT 0 ? 1.e-5/GFeff : 0.0;
+               coeffG = GFeffG GT 0 ? 1.e-5/GFeffG : 0.0;
+             
+               Spectra[*,j,*,i] = SpectraTot[*,j,*,0,i]*coeffH + SpectraTot[*,j,*,1,i] * coeffG;
+         endfor 
+
+   
+   SpectraOut = total(total(Spectra,1),2);
+   Time=Time[0:Retsize-1]
+   
+return, SpectraOut
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/mex_h_spec_restr_init.pro
@@ -0,0 +1,74 @@
+;====================================================
+;              VEX_h_spec_init
+;                 1.0
+;   Init module of "h_spec" object
+;   HPSPEC from IMAEXTRA product VEX
+;
+;              18 Mars 2009
+;----------------------------------------------------
+
+pro mex_h_spec_restr_init, ANODE = Anode, STOP_ANODE = Stop_Anode, elevationAsTime = elevationAsTime, $
+                     LOG = log, LIN = lin, PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $
+                     ; spectra
+                     MINCOUNT = minCount, MAXCOUNT = maxCount, BACKGRDCHAN = backgrdChan, BACKGRDVAL = backgrdVal, $
+                     CHANNEL = channel, Energy = energy, ENERGYMIN = energymin, ENERGYMAX =  energymax, $
+                     
+                     DELAY = delay, $                      
+                     ; all additional keywords
+                    _EXTRA = e 
+
+  common Request, List, StartTime, TimeInt
+  common GraphC, graph, GraphN, Item
+
+  AzMin = 0
+  AzMax = 15
+
+   elevation = 1;
+   if (KEYWORD_SET(elevationAsTime)) then elevation = 0;
+
+  if (KEYWORD_SET(ANODE)) then begin
+    AzMin = Anode
+   if (KEYWORD_SET(STOP_ANODE)) then AzMax = Stop_Anode $
+   else AzMax = -1;
+  endif
+  
+  
+  graph(GraphN).Type = 'IMA_T'
+  graph(GraphN).DataN = 1
+  graph(GraphN).SetCommonF(0) = 'ima_h_spec_com, "mex"'
+  graph(GraphN).GetDataF(0) = 'get_imam_h_spec, /RESTRICTED'
+  graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) < 36100.0D0
+
+
+EnMin = 1.0;
+EnMax = 25000.0;
+
+energyTitle = '';
+
+
+if (KEYWORD_SET(ENERGYMIN)) then begin
+   EnMin = energymin;
+   EnMax =  energymax;
+   energyTitle = energymin eq energymax ? '!C E ' + string(energymin,format='(f7.1)') + 'eV' : '!C E ' + string(energymin,format='(f7.1)') + '-'+  string(energymax,format='(f7.1)') + 'eV'; 
+    
+    options = '"mex", ' + string(AzMin) +','+string(AzMax) +','+string(EnMin) +','+string(EnMax) + ',' + string(elevation);
+
+    angleTitle =   AzMax GT AzMin ?  'MEX IMA H+ An'+ string(AzMin,format='(i2.2)') +'-'+ string(AzMax,format='(i2.2)') : 'MEX IMA H+ An'+ string(AzMin,format='(i2.2)');
+    
+    graph[GraphN].AxisF =  KEYWORD_SET(lin) ?  'lineaxis, 0' : 'logaxis';
+    graph[GraphN].Ly.title =  'MEX H Cnts ' +  angleTitle + energyTitle
+					  
+    graph[GraphN].PlotF = 'imacntplot,' + options
+endif else begin
+    options = '"mex",'+string(AzMin) +','+string(AzMax) + ',' + string(elevation);
+    
+    graph[GraphN].names = PTR_NEW(/ALLOCATE_HEAP)
+    *(graph[GraphN].names) = 'log Counts'
+    
+    
+    graph[GraphN].Ly.title = AzMax GT AzMin ?  'MEX IMA H+ An'+ string(AzMin,format='(i2.2)') +'-'+ string(AzMax,format='(i2.2)') +'!CE/Q, eV' : 'MEX IMA H+ An'+ string(AzMin,format='(i2.2)') + '!C E/Q, eV'
+    graph[GraphN].PlotF =  'imaspecplot,' + options
+endelse
+return
+end
+;------------------------------------------------------------
+++ a/amda_plus/amdalib/mex_h_temp_init.pro
@@ -0,0 +1,33 @@
+;--------------------------------------
+;
+;---------------------------------------
+
+pro mex_h_temp_INIT,  $
+                     LOG = log, LIN = lin, PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                   
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $
+                     
+                     ; all additional keywords
+                    _EXTRA = e
+
+common GraphC, graph, GraphN, Item
+ 
+   pen = !dnc-1  & symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;    
+ 
+graph[GraphN].Type = 'FGM_T' 
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'vexmex_temp_com, "mex", "h"' 
+graph[GraphN].GetDataF[0] ='get_vexmex_temp, "mex", "h"'
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) <  18000.0D0 
+graph[GraphN].Ly.title = 'MEX ima-H T VSO, eV'
+graph[GraphN].PlotF = 'PlotScal_New, 0,' + string(pen)  + ',' + string(symbol);  
+if (KEYWORD_SET(LOG)) then AxisF = 'logaxis' $
+else AxisF = 'lineaxis, ' + string(Graph[GraphN].TotalTime);
+graph[GraphN].AxisF = AxisF; 
+
+return
+end
+++ a/amda_plus/amdalib/mex_h_temp_restr_.pro
@@ -0,0 +1,48 @@
+function MEX_H_TEMP_RESTR_, StartTime, TimeInt, Time, RetSize, Sampling, close=close
+common GraphC, graph, GraphN
+common MEXHDTEMP_COM, ID
+
+   FillValue = !Values.F_NAN 
+   if n_elements(ID) EQ 0 then ID = -1L
+   sampling = 192.0
+   Size = long(graph[GraphN].DeltaT/sampling)
+     
+   Bdata=fltarr(3,Size)
+   Time=dblarr(Size)
+   
+    VI = 'mex:imaparam:restricted'
+    ParamName = "Temperature_P";
+ 
+    RetSize = call_external(!PROJLIB + 'libParam.so', 'getLocalData', $
+                          ID, VI, ParamName, StartTime, TimeInt, Size, Time, Bdata)
+                          
+  
+   if KEYWORD_SET(CLOSE) then begin
+      error = call_external(!DDLIB+'DD_idl.so','CloseID',ID)
+      ID = -1L
+   endif                        
+  
+   if (RetSize LE 0) then begin
+      Time = 0.D0
+      return, 0.0
+   endif
+   
+   Time=Time[0:Retsize-1]
+   Bdata=Bdata[*,0:Retsize-1]
+   
+   num_gap = where(finite(Bdata[0,*]), Cntr)
+  
+   if Cntr le 0 then begin
+       Retsize = -999
+       Time = 0.D0
+       return, 0.0
+   endif
+   
+    Time=Time[num_gap]
+    Bdata=Bdata[*,num_gap];
+    Retsize=Cntr
+    B = fltarr(Retsize)
+    
+    for i = 0, Retsize - 1 do B[i] = min(Bdata[*,i])
+return, B
+end
 \ No newline at end of file
+++ a/amda_plus/amdalib/mex_h_temp_restr_init.pro
@@ -0,0 +1,31 @@
+;--------------------------------------
+;
+;---------------------------------------
+
+pro mex_h_temp_RESTR_INIT,  $
+                     LOG = log, LIN = lin, PSYM = psym, COLOR = color, AY1TITLE = AY1Title, $                   
+                     ; vector
+                     X = x, Y = y, Z = z, TOTAL = total, FRAME = frame, $
+                     DELAY = delay, $
+                     
+                     ; all additional keywords
+                    _EXTRA = e
+
+common GraphC, graph, GraphN, Item
+ 
+   pen = !dnc-1 &  symbol = -3;
+ 
+ if (KEYWORD_SET(COLOR)) then pen = color; 
+ if (KEYWORD_SET(PSYM)) then symbol = psym;    
+ 
+graph[GraphN].Type = 'FGM_T' 
+graph[GraphN].DataN = 1
+graph[GraphN].SetCommonF[0] = 'vexmex_temp_com, "mex", "h"' 
+graph[GraphN].GetDataF[0] ='get_vexmex_temp, "mex", "h", /RESTRICTED'
+graph[GraphN].DeltaT = (Graph[GraphN].TotalTime + 60.0D0) <  18000.0D0 
+graph[GraphN].Ly.title = 'MEX ima-H T VSO, eV'
+graph[GraphN].PlotF = 'PlotScal_New, 0,' + string(pen)  + ',' + string(symbol);  
+graph[GraphN].AxisF = 'lineaxis,' + string(graph[graphN].TotalTime)
+
+return
+end