if 0: #Alboran
xs = 160.
xh = xh[xs:]
bot = bot[:, xs:]
if 1: #Alboran run90+
xs = 140.
xh = xh[xs:]
bot = bot[:, xs:]
bot[bot == 0.] = 0.1
bot[0, :] = 0.1
bot[-1, :] = 0.1 # dirty
md.create_grd(grdname, yh, xh, 'x')
grid = NF(grdname, 'a')
grid.variables['xh'][:] = xh
grid.variables['yh'][:] = yh
grid.variables['D'][:] = bot
grid.close()
#
print ""
print "lenlat: " + str(np.diff(yh)[0] * yh.shape[0])
print "lenlon: " + str(np.diff(xh)[0] * xh.shape[0])
print ""
print "lowlat: " + str(min(yh))
print "westlon: " + str(min(xh))
#
print "size: " + str(bot.shape)
if 0: #Alboran
xs=160.;
xh=xh[xs:]
bot=bot[:,xs:]
if 1: #Alboran run90+
xs=140.;
xh=xh[xs:]
bot=bot[:,xs:]
bot[bot==0.]=0.1
bot[0,:]=0.1
bot[-1,:]=0.1 # dirty
md.create_grd(grdname,yh,xh,'x')
grid=NF(grdname,'a')
grid.variables['xh'][:]=xh
grid.variables['yh'][:]=yh
grid.variables['D'][:]=bot
grid.close()
#
print ""
print "lenlat: " + str(np.diff(yh)[0]*yh.shape[0])
print "lenlon: " + str(np.diff(xh)[0]*xh.shape[0])
print ""
print "lowlat: " + str(min(yh))
print "westlon: " + str(min(xh))
#
print "size: " + str(bot.shape)
x2=np.arange(np.size(yi))
for j in np.arange(0,np.size(yi)):
x1[j]=np.size(x)*(np.size(y)-yi[j])+xi[0];
x2[j]=np.size(x)*(np.size(y)-yi[j])+xi[-1];
tz=np.nan*np.ones([np.size(yi),np.size(xi)]);
if 1:
for j in np.arange(0,np.size(yi)):
tz[j]=ff.variables['z'][x1[j]:x2[j]+1];
f = sci_interp.RectBivariateSpline(x[xi],y[yi], tz.transpose())
h = -f(lon, lat).transpose()
md.create_grd(grdname,lat,lon,'l')
grid=NF(grdname,'a')
grid.variables['lonh'][:]=lon
grid.variables['lath'][:]=lat
grid.variables['D'][:]=h
grid.close()
print "lenlat: " + str(lat[-1]-lat[0])
print "lenlon: " + str(lon[-1]-lon[0])
print ""
print "lowlat: " + str(min(lat))
print "westlon: " + str(min(lon))
print "size: " + str(h.shape)
# for ii in range(13,30): # de=bot[ii,w0-len(sl)] # nu=(ii-12.)/(30-13) # print nu # de2=de*nu+0.1*(1.-nu) # de2=bot[29,w0]*nu+0.1*(1.-nu) # bot[ii,w0-len(sl):w1+len(sl)]=de2*gu+de*(1.-gu) # de2=bot[0:30,-1] # gr=de2.shape[0] # for ii in range(gr): # de=bot[ii+13,w0-len(sl)] # bot[ii+13,w0-len(sl):w1+len(sl)]=de2[ii]*gu+de*(1.-gu) md.create_grd(grdname, yh, xh, "x") grid = NF(grdname, "a") grid.variables["xh"][:] = xh grid.variables["yh"][:] = yh grid.variables["D"][:] = bot grid.close() # print "" print "lenlat: " + str(np.diff(yh)[0] * yh.shape[0]) print "lenlon: " + str(np.diff(xh)[0] * xh.shape[0]) print "" print "lowlat: " + str(min(yh)) print "westlon: " + str(min(xh)) # print "size: " + str(bot.shape)
