Skip to content

Comments

Add get_climatology_stdev()#16

Open
tramsauer wants to merge 6 commits intopygeo:devfrom
tramsauer:dev
Open

Add get_climatology_stdev()#16
tramsauer wants to merge 6 commits intopygeo:devfrom
tramsauer:dev

Conversation

@tramsauer
Copy link

Added get_climatology_stdev() as discrete method to also extract standard deviation in climatological values.

mask=(np.isnan(clim) | (n < nmin) |
np.logical_not(
(np.logical_not(self.data.mask)).
mean(axis=0))))

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E126 continuation line over-indented for hanging indent

Copy link
Contributor

@bulli92 bulli92 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for the PR.

Can you please
a) take care of the stickler-bot message and correct the line mentioned

b) what I am missing is a unittest, testing the routine. Please implement one and then redo the PR

#~ for i in xrange(n):
#~ try:
#~ yy = stats.mstats.linregress(x, dat[:, i])
#~ try:

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

#~ try:
#~ yy = stats.mstats.linregress(x, dat[:, i])
#~ try:
#~ yy = stats.mstats.linregress(x, dat[:, i])

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

#~ print x
#~ print dat[:,i]
#~ stop
#res[i] = stats.mstats.linregress(x, dat[:, i])

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

#~ print dat[:,i]
#~ stop
#res[i] = stats.mstats.linregress(x, dat[:, i])
#~ except:

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

#~ stop
#res[i] = stats.mstats.linregress(x, dat[:, i])
#~ except:
#~ print x

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

#res[i] = stats.mstats.linregress(x, dat[:, i])
#~ except:
#~ print x
#~ print dat[:,i]

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

#~ except:
#~ print x
#~ print dat[:,i]
#~ stop

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

#~ for i in xrange(len(self.time)):
#~ Y,M,D,h,m,s = self._split_time_float(self.time[i])
#~ newdate.append(basedate + relativedelta.relativedelta(years=Y,months=M,days=D,hours=h,minutes=m,seconds=s))
#~ Y,M,D,h,m,s = self._split_time_float(self.time[i])

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

#~ Y,M,D,h,m,s = self._split_time_float(self.time[i])
#~ newdate.append(basedate + relativedelta.relativedelta(years=Y,months=M,days=D,hours=h,minutes=m,seconds=s))
#~ Y,M,D,h,m,s = self._split_time_float(self.time[i])
#~ newdate.append(basedate + relativedelta.relativedelta(years=Y,months=M,days=D,hours=h,minutes=m,seconds=s))

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

#~ ----------
#~ t : float
#~ scalar time indicator
#~ scalar time indicator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

E265 block comment should start with '# '

information on areas from a classic ESRI shapefile
"""
assert isinstance(shape,shp.Reader)
assert isinstance(shape, shp.Reader)

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

F821 undefined name 'shp'

"""
assert isinstance(shape,shp.Reader)

assert isinstance(shape, shp.Reader)

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

F821 undefined name 'shp'

@tramsauer
Copy link
Author

  • don't know how to deal with shp.Reader --> add import?
  • errors in Travis related to older tests already implemented.
  • my autopep formatter changed code formatting.. if not desired i can use not formatted version and resubmit

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants