salem

https://salem.readthedocs.io/en/latest/wrf.html

In [12]: fpath = get_demo_file('namelist_mercator.wps')

In [13]: with open(fpath, 'r') as f: # this is just to show the file
....: print(f.read())
....:
&geogrid
parent_id = 1, 1, 2, 3,
parent_grid_ratio = 1, 3, 5, 5,
i_parent_start = 1, 28, 54, 40,
j_parent_start = 1, 24, 97, 50,
e_we = 211, 178, 101, 136,
e_sn = 131, 220, 141, 136,
dx = 30000,
dy = 30000,
map_proj = 'mercator',
ref_lat = -8.0,
ref_lon = -62.0,
truelat1 = -8.0,
/

In [14]: from salem import geogrid_simulator

In [15]: g, maps = geogrid_simulator(fpath)
I am densified (external_values, 115 elements)
I am densified (external_values, 8 elements)
I am densified (external_values, 1 elements)
I am densified (external_values, 1 elements)

In [16]: maps[0].set_rgb(natural_earth='lr') # add a background image

In [17]: maps[0].visualize(title='Domains 1 to 4')

domain