eta_ini=-eta; u=np.squeeze(0.*eta[0:2,...]); lay=[1027.6,1029.1] md.create_init(initname,grdname,'x') init=NF(initname,'a') init.variables['LAYER'][:]=lay; init.variables['ETA'][:]=eta_ini init.variables['u'][:]=u init.variables['v'][:]=u init.close() eta_bdy=-eta; md.create_bdy(bdyname,grdname,'x') bdy=NF(bdyname,'a') bdy.variables['LAYER'][:]=lay; bdy.variables['ETA'][:]=eta_ini bdy.variables['u'][:]=u bdy.variables['v'][:]=u bdy.close() plt.close() plt.close() plt.close() plt.close() plt.close() plt.close() plt.close() plt.figure(1)
eta_ini = -eta u = np.squeeze(0. * eta[0:2, ...]) lay = [1027.6, 1030.1] md.create_init(initname, grdname, 'x') init = NF(initname, 'a') init.variables['LAYER'][:] = lay init.variables['ETA'][:] = eta_ini init.variables['u'][:] = u init.variables['v'][:] = u init.close() eta_bdy = -eta md.create_bdy(bdyname, grdname, 'x') bdy = NF(bdyname, 'a') bdy.variables['LAYER'][:] = lay bdy.variables['ETA'][:] = eta_ini bdy.variables['u'][:] = u bdy.variables['v'][:] = u bdy.close() plt.close() plt.close() plt.close() plt.close() plt.close() plt.close() plt.close() plt.figure(1)
eta_ini = -eta u = np.squeeze(0.0 * eta[0:2, ...]) lay = [1027.6, 1030.3] md.create_init(initname, grdname, "x") init = NF(initname, "a") init.variables["LAYER"][:] = lay init.variables["ETA"][:] = eta_ini init.variables["u"][:] = u init.variables["v"][:] = u init.close() eta_bdy = -eta md.create_bdy(bdyname, grdname, "x") bdy = NF(bdyname, "a") bdy.variables["LAYER"][:] = lay bdy.variables["ETA"][:] = eta_ini bdy.close() plt.close() plt.close() plt.close() plt.close() plt.close() plt.close() plt.close() plt.figure(1) ax = plt.axes() # p=ax.contourf(lon,lat,eta_ini[1,...],20)