Skip to content

Commit 7b2f5f7

Browse files
committed
added test case for ascontiguous bug
1 parent ddc3a6d commit 7b2f5f7

1 file changed

Lines changed: 34 additions & 0 deletions

File tree

Lines changed: 34 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,34 @@
1+
import cdms2,MV2,regrid2,cdutil
2+
3+
4+
a=MV2.reshape(MV2.sin(MV2.arange(20000)),(2,1,100,100))
5+
6+
lon=cdms2.createAxis(MV2.arange(100)*3.6)
7+
lon.designateLongitude()
8+
lon.units="degrees_east"
9+
lon.id="longitude"
10+
11+
lat = cdms2.createAxis(MV2.arange(100)*1.8-90.)
12+
lat.id="latitude"
13+
lat.designateLatitude()
14+
lat.units="degrees_north"
15+
16+
lev = cdms2.createAxis([1000.])
17+
lev.id="plev"
18+
lev.designateLevel()
19+
lev.units="hPa"
20+
21+
t=cdms2.createAxis([0,31.])
22+
t.id="time"
23+
t.designateTime()
24+
t.units="days since 2014"
25+
26+
cdutil.setTimeBoundsMonthly(t)
27+
a.setAxisList((t,lev,lat,lon))
28+
a=MV2.masked_less(a,.5)
29+
grd=cdms2.createGaussianGrid(64)
30+
31+
a=a.ascontiguous()
32+
a=a.regrid(grd,regridTool="regrid2")
33+
a=cdutil.averager(a,axis='txy')
34+
assert a[0]==0.7921019540305255

0 commit comments

Comments
 (0)