Index: trunk/UTIL/interpolate_1d_tprofile_z2sig.py
===================================================================
--- trunk/UTIL/interpolate_1d_tprofile_z2sig.py	(revision 3730)
+++ trunk/UTIL/interpolate_1d_tprofile_z2sig.py	(revision 3731)
@@ -4,10 +4,19 @@
 import numpy as np
 import pandas as pd
+import sys
+
+if sys.argv[1] in ("-h", "--help"):
+    print("Must be provided:\n"
+          "profile.in   : vertical profile to interpolate\n"
+          "z2sig.def.in : input altitude profile\n"
+          "z2sig.def    : output altitude profile\n"
+    )
+    exit(0)
 
 temperature_in = pd.read_csv("profile.in", header=None).iloc[:,0]
 print(len(temperature_in), temperature_in[0])
-z2sig_in = pd.read_csv("z2sig.def.in", names=range(40), header=None, skiprows=1, delimiter=" ").iloc[:,0]
+z2sig_in = pd.read_csv("z2sig.def.in", names=range(40), header=None, skiprows=1, delimiter=" +").iloc[:,0]
 print(len(z2sig_in), z2sig_in[0])
-z2sig_out = pd.read_csv("z2sig.def", names=range(40), header=None, skiprows=1, delimiter=" ").iloc[:,0]
+z2sig_out = pd.read_csv("z2sig.def", names=range(40), header=None, skiprows=1, delimiter=" +").iloc[:,0]
 print(len(z2sig_out), z2sig_out[0])
 
@@ -16,4 +25,12 @@
 z2sig_out = np.insert(z2sig_out,0,0)
 
+print("Input temperature")
+print(temperature_in)
+print("Input altitude")
+print(z2sig_in)
+print("Output altitude")
+print(z2sig_out)
 temperature_out = np.interp(z2sig_out, z2sig_in, temperature_in)
+print("Output temperature")
+print(temperature_out)
 np.savetxt("profile", temperature_out, fmt="%f")