mat[:,44:54]= x1; mat[:,54:]=0.; eta2=eta2b*mat+eta2a*(1-mat); #eta2=eta1+10. ### rest #eta2=75*np.ones(eta2.shape) ### eta3=h; eta=np.array([eta1,eta2,eta3]) 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()
mat[:, 44:54] = x1 mat[:, 54:] = 0. eta2 = eta2b * mat + eta2a * (1 - mat) #eta2=eta1+10. ### rest #eta2=75*np.ones(eta2.shape) ### eta3 = h eta = np.array([eta1, eta2, eta3]) 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()
mat[:, 24:34] = x1 mat[:, 34:] = 0.0 eta2 = eta2b * mat + eta2a * (1 - mat) # eta2=eta1+10. ### rest # eta2=75*np.ones(eta2.shape) ### eta3 = h eta = np.array([eta1, eta2, eta3]) 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()