source: trunk/src/vnx1.py @ 54

vnx1=np.array([60,22.5,14.4,10,8,7.2,6,5,4.5,4,3.6,3,3,3,2.5,2.5,2.5,2,2,2,2,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.2,1.2,1.2,1.2,1.2,1.2,1.2,1,1,1,1])
vnx1=vnx1*10

[1,2,3,4,5,6,8,9,10,12,15,16,18,20,24,25,30,36,
40,45,48,50,60,72,75,80,90,100,120,144,150,180,200,225,240,300,360,
400,450,600,720,900,1200,1800,3600]



nd_mean16=zeros([ny], float)
nd_max16=zeros([ny], float)
nd_min16=zeros([ny], float)
for i in range(0,ny):
    nd_mean16[i]=mean(tbnngrid_t16[i,:])
    nd_max16[i]=max(tbnngrid_t16[i,:])
    nd_min16[i]=min(tbnngrid_t16[i,:])

plot(nd_min16)
plt.show()

nd_lat=zeros([ny], float)
nd_lat=sum(tbngrid_t60)

##########################
inlat = (y >= -90) & (y < -88)
inz = (z >= 100) & (z <= 300)
innn = inlat & inz

zz = z[innn]
llat=y[innn]
llon=x[innn]

plot(llon,llat)
plt.show()

tt1='point_tb_SSMIS_CH16'
t1='scatter_point_tb_SSMIS_CH16_centre.png'



# on recupere les coordonnes xy de la projection orth standard
m1 = Basemap(projection='ortho', lat_0 = -90, lon_0 = 0,
resolution = 'l')

xii, yii = m1(*np.meshgrid(xvec,yvec))
width = m1.urcrnrx - m1.llcrnrx
height = m1.urcrnry - m1.llcrnry

coef = 0.6
width = width*coef
height = height*coef

# on fait une nouvelle projection en zoomant sur l'antarctique

map = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='l',\
llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height)
x1,y1=map(llon,llat)
clevs=arange(100,300,1)
cs=map.scatter(x1,y1,16,zz,marker="o",cmap=cm.s3pcpn_l_r,vmin=100, vmax=300,alpha=0.7)
cbar =colorbar(cs)
plt.title(tt1)
# draw coastlines, country boundaries, fill continents.
map.drawcoastlines(linewidth=1)
# draw the edge of the map projection region (the projection limb)
map.drawmapboundary()
# draw lat/lon grid lines every 30 degrees.
map.drawmeridians(np.arange(0, 360, 1))
map.drawparallels(np.arange(-90, 90, 1))
#plt.show()
plt.savefig(t1)
close()