#text(100, 0, "185", horizontalalignment='center', verticalalignment='top', fontsize=colorbar_fontsize,rotation=0)
#text(10, 10, "104", horizontalalignment='center', verticalalignment='top', fontsize=colorbar_fontsize,rotation=0)
#savefig(sys.argv[1]+'-moll.png',dpi=dpires, facecolor=bgcolor, edgecolor=bgcolor); #savefig(sys.argv[1]+'-moll.eps')
#savefig(sys.argv[1]+'-moll.gif',dpi=dpires, facecolor=bgcolor, edgecolor=bgcolor); #savefig(sys.argv[1]+'-moll.eps')
#savefig(sys.argv[1]+'-moll.eps',dpi=dpires, facecolor=bgcolor, edgecolor=bgcolor); #savefig(sys.argv[1]+'-moll.eps')
#show()
#exit
if plotNS == '1' :
# north polar projection.
m2 = Basemap(llcrnrlon=-45.,llcrnrlat=10.,urcrnrlon=135.,urcrnrlat=10., projection='stere', lat_0=90.,lon_0=0.,lat_ts=90.)
# resolution='c',area_thresh=10000.,
fig2 = m2.createfigure()
x, y = m2(lons, lats)
cs = m2.contourf(x,y,mapM,color_num,cmap=cm.jet)
if mymeridians == '1' :
delat = DeltaLatitude
circles = arange(0.,90.+delat,delat).tolist()+arange(-delat,-90.-delat,-delat).tolist()
m2.drawparallels(circles,labels=[1,1,1,1])
# draw meridians
delon = DeltaLongitude
meridians = arange(-180,180,delon)
m2.drawmeridians(meridians,labels=[1,1,1,1])
if mytitle != 'notitle' :
title(mytitle,color=fgcolor)
