schain-jc Commit - r870:bd8def453202

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

8

import matplotlib.pyplot as plt

9

8

10

class MomentsPlot(Figure):

9

class MomentsPlot(Figure):

11

10

447

# tmin = None

446

# tmin = None

448

# tmax = None

447

# tmax = None

449

448

450

x = dataOut.getTimeRange1(dataOut.outputInterval)

449

451

y = dataOut.~~heightList~~

450

x = dataOut.getTimeRange()

451

y = dataOut.getHeiRange()

452

#x = dataOut.getTimeRange1(dataOut.outputInterval)

453

#y = dataOut.heightList

454

452

z = dataOut.data_output.copy()

455

z = dataOut.data_output.copy()

456

#print 'dataOut_JI',z

453

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

457

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

454

455

nplotsw = nplots

458

nplotsw = nplots

456

459

460

457

#If there is a SNR function defined

461

#If there is a SNR function defined

458

if dataOut.data_SNR is not None:

462

if dataOut.data_SNR is not None:

459

nplots += 1

463

nplots += 1

472

476

473

# showprofile = False

477

# showprofile = False

474

# thisDatetime = dataOut.datatime

478

# thisDatetime = dataOut.datatime

475

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

479

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

476

thisDatetime = datetime.datetime.now()

477

title = wintitle + "Wind"

480

title = wintitle + "Wind"

478

xlabel = ""

481

xlabel = ""

479

ylabel = "Height (km)"

482

ylabel = "Height (km)"

492

495

493

self.xmin, self.xmax = self.getTimeLim(x, xmin, xmax, timerange)

496

self.xmin, self.xmax = self.getTimeLim(x, xmin, xmax, timerange)

494

497

495

#if ymin == None: ymin = numpy.nanmin(y)

498

if ymin == None: ymin = numpy.nanmin(y)

496

#if ymax == None: ymax = numpy.nanmax(y)

499

if ymax == None: ymax = numpy.nanmax(y)

497

500

498

if zmax == None: zmax = numpy.nanmax(abs(z[range(2),:]))

501

if zmax == None: zmax = numpy.nanmax(abs(z[range(2),:]))

499

#if numpy.isnan(zmax): zmax = 50

502

#if numpy.isnan(zmax): zmax = 50

503

if zmax_ver == None: zmax_ver = numpy.nanmax(abs(z[2,:]))

506

if zmax_ver == None: zmax_ver = numpy.nanmax(abs(z[2,:]))

504

if zmin_ver == None: zmin_ver = -zmax_ver

507

if zmin_ver == None: zmin_ver = -zmax_ver

505

508

506

# if dataOut.data_SNR is not None:

509

if dataOut.data_SNR is not None:

507

# if SNRmin == None: SNRmin = numpy.nanmin(SNRavgdB)

510

if SNRmin == None: SNRmin = numpy.nanmin(SNRavgdB)

508

# if SNRmax == None: SNRmax = numpy.nanmax(SNRavgdB)

511

if SNRmax == None: SNRmax = numpy.nanmax(SNRavgdB)

509

512

510

513

511

self.FTP_WEI = ftp_wei

514

self.FTP_WEI = ftp_wei

520

523

521

self.setWinTitle(title)

524

self.setWinTitle(title)

522

525

523

#if ((self.xmax - x[1]) < (x[1]-x[0])):

526

if ((self.xmax - x[1]) < (x[1]-x[0])):

524

# x[1] = self.xmax

527

x[1] = self.xmax

525

528

526

strWind = ['Zonal', 'Meridional', 'Vertical']

529

strWind = ['Zonal', 'Meridional', 'Vertical']

527

strCb = ['Velocity (m/s)','Velocity (m/s)','Velocity (cm/s)']

530

strCb = ['Velocity (m/s)','Velocity (m/s)','Velocity (cm/s)']

535

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

538

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

536

539

537

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

540

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

541

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

538

542

539

axes.pcolorbuffer(x, y, z1,

543

axes.pcolorbuffer(x, y, z1,

540

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

544

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

541

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

545

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

542

ticksize=9, cblabel=strCb[i], cbsize="1%", colormap="~~RdBu_r~~" )

546

ticksize=9, cblabel=strCb[i], cbsize="1%", colormap="seismic" )

543

547

544

if dataOut.data_SNR is not None:

548

if dataOut.data_SNR is not None:

545

i += 1

549

i += 1

563

thisDatetime=thisDatetime,

567

thisDatetime=thisDatetime,

564

update_figfile=update_figfile)

568

update_figfile=update_figfile)

565

569

566

if ~~False~~ ~~and~~ dataOut.ltctime + dataOut.outputInterval >= self.xmax:

570

if dataOut.ltctime + dataOut.outputInterval >= self.xmax:

567

self.counter_imagwr = wr_period

571

self.counter_imagwr = wr_period

568

self.isConfig = False

572

self.isConfig = False

569

update_figfile = True

573

update_figfile = True

780

784

781

785

782

786

783

class ParametersPlot(Figure):

787

class Parameters1Plot(Figure):

784

788

785

__isConfig = None

789

__isConfig = None

786

__nsubplots = None

790

__nsubplots = None

	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
-            import matplotlib.pyplot as plt
             class MomentsPlot(Figure):
             #         tmin = None
             #         tmax = None
-                    x = dataOut.getTimeRange1(dataOut.outputInterval)
-                    y = dataOut.heightList
+                    x = dataOut.getTimeRange()
+                    y = dataOut.getHeiRange()
+                    #x = dataOut.getTimeRange1(dataOut.outputInterval)
+                    #y = dataOut.heightList
                     z = dataOut.data_output.copy()
+                    #print 'dataOut_JI',z
                     nplots = z.shape[0]    #Number of wind dimensions estimated
                     nplotsw = nplots
                     #If there is a SNR function defined
                     if dataOut.data_SNR is not None:
                         nplots += 1
             #         showprofile = False
             #        thisDatetime = dataOut.datatime
-                    #thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.ltctime)
+                    thisDatetime = datetime.datetime.utcfromtimestamp(dataOut.ltctime)
-                    thisDatetime = datetime.datetime.now()
                     title = wintitle + "Wind"
                     xlabel = ""
                     ylabel = "Height (km)"
                         self.xmin, self.xmax = self.getTimeLim(x, xmin, xmax, timerange)
-                        #if ymin == None: ymin = numpy.nanmin(y)
+                        if ymin == None: ymin = numpy.nanmin(y)
-                        #if ymax == None: ymax = numpy.nanmax(y)
+                        if ymax == None: ymax = numpy.nanmax(y)
                         if zmax == None: zmax = numpy.nanmax(abs(z[range(2),:]))
                         #if numpy.isnan(zmax): zmax = 50
                             if zmax_ver == None: zmax_ver = numpy.nanmax(abs(z[2,:]))
                             if zmin_ver == None: zmin_ver = -zmax_ver
-            #             if dataOut.data_SNR is not None:
+                        if dataOut.data_SNR is not None:
-            #                 if SNRmin == None:  SNRmin = numpy.nanmin(SNRavgdB)
+                            if SNRmin == None:  SNRmin = numpy.nanmin(SNRavgdB)
-            #                 if SNRmax == None:  SNRmax = numpy.nanmax(SNRavgdB)
+                            if SNRmax == None:  SNRmax = numpy.nanmax(SNRavgdB)
                         self.FTP_WEI = ftp_wei
                     self.setWinTitle(title)
-                    #if ((self.xmax - x[1]) < (x[1]-x[0])):
+                    if ((self.xmax - x[1]) < (x[1]-x[0])):
-                    #    x[1] = self.xmax
+                        x[1] = self.xmax
                     strWind = ['Zonal', 'Meridional', 'Vertical']
                     strCb = ['Velocity (m/s)','Velocity (m/s)','Velocity (cm/s)']
                         axes = self.axesList[i*self.__nsubplots]
                         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,
-                                    ticksize=9, cblabel=strCb[i], cbsize="1%", colormap="RdBu_r" )
+                                    ticksize=9, cblabel=strCb[i], cbsize="1%", colormap="seismic" )
                     if dataOut.data_SNR is not None:
                         i += 1
                               thisDatetime=thisDatetime,
                               update_figfile=update_figfile)
-                    if False and dataOut.ltctime + dataOut.outputInterval >= self.xmax:
+                    if dataOut.ltctime + dataOut.outputInterval >= self.xmax:
                         self.counter_imagwr = wr_period
                         self.isConfig = False
                         update_figfile = True
-            class ParametersPlot(Figure):
+            class Parameters1Plot(Figure):
                 __isConfig = None
                 __nsubplots = None