schain-jc Commit - r877:f62d4806c545

4

5

from figure import Figure, isRealtime, isTimeInHourRange

5

from figure import Figure, isRealtime, isTimeInHourRange

6

from plotting_codes import *

6

from plotting_codes import *

7

from numpy import NaN

7

8

9

class MomentsPlot(Figure):

9

class MomentsPlot(Figure):

10

446

# tmin = None

446

# tmin = None

447

# tmax = None

447

# tmax = None

448

449

450

x = dataOut.getTimeRange()

450

x = dataOut.getTimeRange1(dataOut.outputInterval)

451

y = dataOut.~~getHeiRange~~()

451

y = dataOut.heightList

452

#x = dataOut.getTimeRange1(dataOut.outputInterval)

452

z = dataOut.data_output.copy()

453

#y = dataOut.heightList

454

455

z = dataOut.data_output.copy()

456

#print 'dataOut_JI',z

457

nplots = z.shape[0] #Number of wind dimensions estimated

453

nplots = z.shape[0] #Number of wind dimensions estimated

458

nplotsw = nplots

454

nplotsw = nplots

459

455

472

468

473

for i in range(nplotsw):

469

for i in range(nplotsw):

474

z[i,ind] = numpy.nan

470

z[i,ind] = numpy.nan

475

471

476

477

# showprofile = False

478

# thisDatetime = dataOut.datatime

479

thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.ltctime)

472

thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.ltctime)

473

#thisDatetime = datetime.datetime.now()

480

title = wintitle + "Wind"

474

title = wintitle + "Wind"

481

xlabel = ""

475

xlabel = ""

482

ylabel = "Height (km)"

476

ylabel = "Height (km)"

539

533

540

z1 = z[i,:].reshape((1,-1))*windFactor[i]

534

z1 = z[i,:].reshape((1,-1))*windFactor[i]

541

#z1=numpy.ma.masked_where(z1==0.,z1)

535

#z1=numpy.ma.masked_where(z1==0.,z1)

542

536

543

axes.pcolorbuffer(x, y, z1,

537

axes.pcolorbuffer(x, y, z1,

544

xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=zminVector[i], zmax=zmaxVector[i],

538

xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=zminVector[i], zmax=zmaxVector[i],

545

xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,

539

xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,

549

i += 1

543

i += 1

550

title = "Signal Noise Ratio (SNR): %s" %(thisDatetime.strftime("%Y/%m/%d %H:%M:%S"))

544

title = "Signal Noise Ratio (SNR): %s" %(thisDatetime.strftime("%Y/%m/%d %H:%M:%S"))

551

axes = self.axesList[i*self.__nsubplots]

545

axes = self.axesList[i*self.__nsubplots]

552

546

SNRavgdB = SNRavgdB.reshape((1,-1))

553

SNRavgdB = SNRavgdB.reshape((1,-1))

554

555

axes.pcolorbuffer(x, y, SNRavgdB,

547

axes.pcolorbuffer(x, y, SNRavgdB,

556

xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=SNRmin, zmax=SNRmax,

548

xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=SNRmin, zmax=SNRmax,

557

xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,

549

xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages

             from figure import Figure, isRealtime, isTimeInHourRange
             from plotting_codes import *
-            from numpy import NaN
             class MomentsPlot(Figure):
             #         tmin = None
             #         tmax = None
-                    x = dataOut.getTimeRange()
+                    x = dataOut.getTimeRange1(dataOut.outputInterval)
-                    y = dataOut.getHeiRange()
+                    y = dataOut.heightList
-                    #x = dataOut.getTimeRange1(dataOut.outputInterval)
+                    z = dataOut.data_output.copy()
-                    #y = dataOut.heightList
-                    z = dataOut.data_output.copy()
-                    #print 'dataOut_JI',z
                     nplots = z.shape[0]    #Number of wind dimensions estimated
                     nplotsw = nplots
                         for i in range(nplotsw):
                             z[i,ind] = numpy.nan
-            #         showprofile = False
-            #        thisDatetime = dataOut.datatime
                     thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.ltctime)
+                    #thisDatetime = datetime.datetime.now()
                     title = wintitle + "Wind"
                     xlabel = ""
                     ylabel = "Height (km)"
                         z1 = z[i,:].reshape((1,-1))*windFactor[i]
                         #z1=numpy.ma.masked_where(z1==0.,z1)
                         axes.pcolorbuffer(x, y, z1,
                                     xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=zminVector[i], zmax=zmaxVector[i],
                                     xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,
                         i += 1
                         title = "Signal Noise Ratio (SNR): %s" %(thisDatetime.strftime("%Y/%m/%d %H:%M:%S"))
                         axes = self.axesList[i*self.__nsubplots]
+                        SNRavgdB = SNRavgdB.reshape((1,-1))
-                        SNRavgdB = SNRavgdB.reshape((1,-1))
                         axes.pcolorbuffer(x, y, SNRavgdB,
                                     xmin=self.xmin, xmax=self.xmax, ymin=ymin, ymax=ymax, zmin=SNRmin, zmax=SNRmax,
                                     xlabel=xlabel, ylabel=ylabel, title=title, rti=True, XAxisAsTime=True,