Changeset 827 for palm/trunk/SCRIPTS/NCL

Timestamp:

Feb 20, 2012 1:09:00 PM (13 years ago)

Author:

heinze

Message:

allow plotting of data with very small time increments

Location:

palm/trunk/SCRIPTS/NCL

Files:

• cross_sections.ncl (modified) (1 diff)
• profiles.ncl (modified) (2 diffs)
• spectra.ncl (modified) (1 diff)
• timeseries.ncl (modified) (1 diff)

palm/trunk/SCRIPTS/NCL/cross_sections.ncl

@@ -419 +419 @@
    ; *********************************************
 
-   t_all = f[:]->time
-   nt = dimsizes(t_all)
-   delta_t = t_all(nt-1)/nt
+   t_all         = f[:]->time
+   nt            = dimsizes(t_all)
+
+   if (nt .EQ. 1)then
+      delta_t=t_all(nt-1)/nt
+   else
+      delta_t_array = new(nt-1,double)
+
+      do i=0,nt-2
+         delta_t_array(i) = t_all(i+1)-t_all(i)
+      end do
+
+      delta_t = min(delta_t_array)
+      delete(delta_t_array)
+   end if
+
 
    ; ****************************************************       

palm/trunk/SCRIPTS/NCL/profiles.ncl

@@ -1010 +1010 @@
    nt    = dimsizes(t_all)
 
-   
    if (nt .EQ. 1)then
       delta_t=t_all(nt-1)/nt
    else
-      delta_t=(t_all(nt-1)-t_all(0))/(nt-1)
-   end if
+      delta_t_array = new(nt-1,double)
+
+      do i=0,nt-2
+         delta_t_array(i) = t_all(i+1)-t_all(i)
+      end do
+
+      delta_t = min(delta_t_array)
+      delete(delta_t_array)
+   end if
+
 
    ; ****************************************************       
@@ -1205 +1212 @@
    resP                        = True
    resP@txFont                 = "helvetica"
-   resP@txString               = f_att@title
+   ;resP@txString               = f_att@title
    resP@txFuncCode             = "~"
    resP@txFontHeightF          = 0.015

palm/trunk/SCRIPTS/NCL/spectra.ncl

@@ -220 +220 @@
    t_all = f[:]->time
    nt    = dimsizes(t_all)
-   delta_t=t_all(nt-1)/nt
+
+   if (nt .EQ. 1)then
+      delta_t=t_all(nt-1)/nt
+   else
+      delta_t_array = new(nt-1,double)
+
+      do i=0,nt-2
+         delta_t_array(i) = t_all(i+1)-t_all(i)
+      end do
+
+      delta_t = min(delta_t_array)
+      delete(delta_t_array)
+   end if
    
    k_x=f_att->k_x

palm/trunk/SCRIPTS/NCL/timeseries.ncl

@@ -193 +193 @@
    t_all   = f[:]->time
    nt      = dimsizes(t_all)
-   delta_t = t_all(nt-1)/nt
+
+   if (nt .EQ. 1)then
+      delta_t=t_all(nt-1)/nt
+   else
+      delta_t_array = new(nt-1,double)
+
+      do i=0,nt-2
+         delta_t_array(i) = t_all(i+1)-t_all(i)
+      end do
+
+      delta_t = min(delta_t_array)
+      delete(delta_t_array)
+   end if
 
    if (nt .LE. 1) then