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)