Index: trunk/UTIL/z2sig/build_sig.py
===================================================================
--- trunk/UTIL/z2sig/build_sig.py	(revision 3663)
+++ trunk/UTIL/z2sig/build_sig.py	(revision 3664)
@@ -9,4 +9,21 @@
 import datetime
 from matplotlib import pyplot
+
+
+######################################################################
+# Input Parameters
+######################################################################
+input_file = 'z2sig.def_47'
+ps=1.1 # in Pa
+pa=0.5
+preff=2
+
+# New levels
+mynbl=201
+n=float(mynbl)
+
+plot = True
+# plot = False
+
 
 ######################################################################
@@ -22,4 +39,5 @@
     data=np.zeros(nbline-1,dtype='f')
     i=0
+    nbl = 0
     for line in mylines:
         s=str.split(line)
@@ -28,20 +46,9 @@
         else:
           data[i-1]=float(s[0])
+          nbl = nbl+1
         i=i+1
-    return data,hh
+    return data,hh,nbl
 
-######################################################################
-# Input Parameters
-######################################################################
-nbl=47
-data,hh=readfile('z2sig.def_47')
-# nb of levels
-clev=np.linspace(1,nbl,nbl)
-ps=1.1
-# New levels
-mynbl=201
-n=float(mynbl)
-plot = True
-plot = False
+data,hh,nbl = readfile(input_file)
 
 #####################
@@ -107,6 +114,4 @@
 # hybrid pressure:
 nsig=np.zeros(nbl,dtype='f')
-pa=0.5 #ps
-preff=2 #ps
 for i in range(nbl):
   nsig[i]=sig_hybrid(sig[i],pa,preff)
@@ -161,4 +166,5 @@
 
 if plot:
+    clev=np.linspace(1,nbl,nbl)
     ### Plot pressure levels
     mpl.figure(figsize=(50/2.54, 30/2.54),facecolor='w')