[751] | 1 | from dynamico.meshes import zeros |
---|
[758] | 2 | from dynamico import meshes |
---|
[751] | 3 | import numpy as np |
---|
| 4 | import netCDF4 as cdf |
---|
| 5 | import argparse |
---|
| 6 | |
---|
| 7 | parser = argparse.ArgumentParser() |
---|
| 8 | |
---|
| 9 | parser.add_argument("-nx", type=int, |
---|
[758] | 10 | default=64, choices=None, |
---|
[751] | 11 | help="number of x points") |
---|
| 12 | parser.add_argument("-ny", type=int, |
---|
[758] | 13 | default=64, choices=None, |
---|
[751] | 14 | help="number of y points") |
---|
| 15 | parser.add_argument("-Lx", type=float, |
---|
| 16 | default=8., choices=None, |
---|
| 17 | help="Lx") |
---|
| 18 | parser.add_argument("-Ly", type=float, |
---|
| 19 | default=8., choices=None, |
---|
| 20 | help="Ly") |
---|
| 21 | parser.add_argument("-llm", type=int, |
---|
| 22 | default=1, choices=[1], |
---|
| 23 | help="number of vertical levels") |
---|
| 24 | args = parser.parse_args() |
---|
| 25 | nx, ny, Lx, Ly, llm, nqdyn = args.nx, args.ny,args.Lx, args.Ly, args.llm, 1 |
---|
| 26 | |
---|
| 27 | dx,dy=Lx/nx,Ly/ny |
---|
[758] | 28 | mesh = meshes.Cartesian_mesh(nx,ny,llm,nqdyn,Lx,Ly,0.) |
---|
[751] | 29 | print('Successfully initialized Cartesian mesh') |
---|
[758] | 30 | mesh.ncwrite('cart_%03d_%03d.nc'%(nx,ny)) |
---|
[751] | 31 | print('Successfully written Cartesian mesh to NetCDF File') |
---|