remove mean for intermagnet plots (temp fix)

Author: rweigel

hapiplot/hapiplot.py

@@ -491,6 +491,11 @@ def hapiplot(*args, **kwargs):
             else:
                 y = np.asarray(data[name])
 
+            remove_mean = False
+            if 'uk/GIN_' in meta['x_server']:
+                remove_mean = True
+                y_mean = np.nanmean(y, axis=0)
+
 
             if 'fill' in meta["parameters"][i] and meta["parameters"][i]['fill']:
                 if  ptype == 'isotime' or ptype == 'string':
@@ -586,13 +591,19 @@ def hapiplot(*args, **kwargs):
                         if nodata:
                             col_name = col_name + " [no data in interval]"
 
+                        if remove_mean:
+                            col_name = "{0:s} - {1:.2f}".format(col_name, y_mean[l])
+
                         if 'label' in meta['parameters'][i] and \
                             type(meta['parameters'][i]['label']) == list and \
                             len(meta['parameters'][i]['label']) > l and \
                             meta['parameters'][i]['label'][l].strip() != '':
                                 col_name = meta['parameters'][i]['label'][l]
                                 if nodata:
                                     col_name = col_name + " [no data in interval]"
+                                else:
+                                    if remove_mean:
+                                        col_name = "{0:s} - {1:.2f}".format(col_name, y_mean[l])
 
                         if type(units) == list:
                             if len(units) == 1:
@@ -633,8 +644,12 @@ def hapiplot(*args, **kwargs):
             if nodata == True:
                 tsopts['nodata'] = True
 
-            with rc_context(rc=opts['rcParams']):
-                fig = timeseries(Time, y, **tsopts)
+            if remove_mean:
+                with rc_context(rc=opts['rcParams']):
+                    fig = timeseries(Time, y-y_mean, **tsopts)
+            else:
+                with rc_context(rc=opts['rcParams']):
+                    fig = timeseries(Time, y, **tsopts)
 
             meta["parameters"][i]['hapiplot']['figure'] = fig