Python create_init示例

示例#1

文件： initial.py 项目： poidl/article_gib1

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()

示例#2

文件： initial.py 项目： zatuch/article_gib1

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()

示例#3

文件： initial.py 项目： poidl/poster_egu2012a

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()