Changeset 827


Ignore:
Timestamp:
Feb 20, 2012 1:09:00 PM (11 years ago)
Author:
heinze
Message:

allow plotting of data with very small time increments

Location:
palm/trunk/SCRIPTS/NCL
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/NCL/cross_sections.ncl

    r770 r827  
    419419   ; *********************************************
    420420
    421    t_all = f[:]->time
    422    nt = dimsizes(t_all)
    423    delta_t = t_all(nt-1)/nt
     421   t_all         = f[:]->time
     422   nt            = dimsizes(t_all)
     423
     424   if (nt .EQ. 1)then
     425      delta_t=t_all(nt-1)/nt
     426   else
     427      delta_t_array = new(nt-1,double)
     428
     429      do i=0,nt-2
     430         delta_t_array(i) = t_all(i+1)-t_all(i)
     431      end do
     432
     433      delta_t = min(delta_t_array)
     434      delete(delta_t_array)
     435   end if
     436
    424437
    425438   ; ****************************************************       
  • palm/trunk/SCRIPTS/NCL/profiles.ncl

    r783 r827  
    10101010   nt    = dimsizes(t_all)
    10111011
    1012    
    10131012   if (nt .EQ. 1)then
    10141013      delta_t=t_all(nt-1)/nt
    10151014   else
    1016       delta_t=(t_all(nt-1)-t_all(0))/(nt-1)
    1017    end if
     1015      delta_t_array = new(nt-1,double)
     1016
     1017      do i=0,nt-2
     1018         delta_t_array(i) = t_all(i+1)-t_all(i)
     1019      end do
     1020
     1021      delta_t = min(delta_t_array)
     1022      delete(delta_t_array)
     1023   end if
     1024
    10181025
    10191026   ; ****************************************************       
     
    12051212   resP                        = True
    12061213   resP@txFont                 = "helvetica"
    1207    resP@txString               = f_att@title
     1214   ;resP@txString               = f_att@title
    12081215   resP@txFuncCode             = "~"
    12091216   resP@txFontHeightF          = 0.015
  • palm/trunk/SCRIPTS/NCL/spectra.ncl

    r769 r827  
    220220   t_all = f[:]->time
    221221   nt    = dimsizes(t_all)
    222    delta_t=t_all(nt-1)/nt
     222
     223   if (nt .EQ. 1)then
     224      delta_t=t_all(nt-1)/nt
     225   else
     226      delta_t_array = new(nt-1,double)
     227
     228      do i=0,nt-2
     229         delta_t_array(i) = t_all(i+1)-t_all(i)
     230      end do
     231
     232      delta_t = min(delta_t_array)
     233      delete(delta_t_array)
     234   end if
    223235   
    224236   k_x=f_att->k_x
  • palm/trunk/SCRIPTS/NCL/timeseries.ncl

    r769 r827  
    193193   t_all   = f[:]->time
    194194   nt      = dimsizes(t_all)
    195    delta_t = t_all(nt-1)/nt
     195
     196   if (nt .EQ. 1)then
     197      delta_t=t_all(nt-1)/nt
     198   else
     199      delta_t_array = new(nt-1,double)
     200
     201      do i=0,nt-2
     202         delta_t_array(i) = t_all(i+1)-t_all(i)
     203      end do
     204
     205      delta_t = min(delta_t_array)
     206      delete(delta_t_array)
     207   end if
    196208
    197209   if (nt .LE. 1) then
Note: See TracChangeset for help on using the changeset viewer.