Changeset 585 for palm/trunk/SCRIPTS/NCL/timeseries.ncl

Timestamp:

Oct 14, 2010 11:15:06 AM (14 years ago)

Author:

heinze

Message:

Bugfix: enable plot of data if it is of kind double instead of kind float

File:

• palm/trunk/SCRIPTS/NCL/timeseries.ncl (modified) (2 diffs)

palm/trunk/SCRIPTS/NCL/timeseries.ncl

@@ -170 +170 @@
    ListSetType(f,"cat")
 
-   vNam=getfilevarnames(f_att)
-   vNam_static=vNam
+   vNam  =getfilevarnames(f_att)
+   vType = getfilevartypes(f_att,vNam) 
+   vNam_static = vNam
+
+   if ((all(vType .eq. "double"))) then ;distinction if data is double or float
+      check_vType = True
+   else
+      check_vType = False
+   end if
 
    print(" ")
@@ -356 +363 @@
                                  
    n=0 
-   minE=1.E27
-   maxE=-1.E27
-   minus=1.E27
-   maxus=-1.E27
-   minu=1.E27
-   maxu=-1.E27
-   minz=1.E27
-   maxz=-1.E27
-   minw=1.E27
-   maxw=-1.E27
-   minp=1.E27
-   maxp=-1.E27
-   mins=1.E27
-   maxs=-1.E27
-
-   data   = new((/dim,(end_time_step-start_time_step)+1/),float)
-   unit   = new(dim,string)
-   data_0 = new((end_time_step-start_time_step)+1,float)
+   if (check_vType) then
+      minE  =  100000.d
+      maxE  = -100000.d
+      minus =  100000.d
+      maxus = -100000.d
+      minu  =  100000.d
+      maxu  = -100000.d
+      minz  =  100000.d
+      maxz  = -100000.d
+      minw  =  100000.d
+      maxw  = -100000.d
+      minp  =  100000.d
+      maxp  = -100000.d
+      mins  =  100000.d
+      maxs  = -100000.d
+   else
+      minE  =  100000.
+      maxE  = -100000.
+      minus =  100000.
+      maxus = -100000.
+      minu  =  100000.
+      maxu  = -100000.
+      minz  =  100000.
+      maxz  = -100000.
+      minw  =  100000.
+      maxw  = -100000.
+      minp  =  100000.
+      maxp  = -100000.
+      mins  =  100000.
+      maxs  = -100000.
+   end if
+
+   if (check_vType) then
+      data   = new((/dim,(end_time_step-start_time_step)+1/),double)
+      data_0 = new((end_time_step-start_time_step)+1,double) 
+      mini   = new(dim,double)
+      maxi   = new(dim,double)
+   else
+      data   = new((/dim,(end_time_step-start_time_step)+1/),float) 
+      data_0 = new((end_time_step-start_time_step)+1,float) 
+      mini   = new(dim,float)
+      maxi   = new(dim,float)
+   end if 
+   unit   = new(dim,string)   
    data_0 = 0.0
-   mini   = new(dim,float)
-   maxi   = new(dim,float)
+  
    
    if (over .EQ. 1) then