Changeset 3664

Timestamp:

Feb 28, 2025, 11:16:41 AM (4 months ago)

Author:

afalco

Message:

z2sig: pa and preff are input parameters.
AF

File:

• trunk/UTIL/z2sig/build_sig.py (modified) (5 diffs)

trunk/UTIL/z2sig/build_sig.py

@@ -9 +9 @@
 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 +39 @@
     data=np.zeros(nbline-1,dtype='f')
     i=0
+    nbl = 0
     for line in mylines:
         s=str.split(line)
@@ -28 +46 @@
         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 +114 @@
 # 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 +166 @@
 
 if plot:
+    clev=np.linspace(1,nbl,nbl)
     ### Plot pressure levels
     mpl.figure(figsize=(50/2.54, 30/2.54),facecolor='w')