Looked into the projection. We are using he Lambert Equal Area projection. We can simply change the projection in the plotting script to use Lamber Conformal Conic. simply by changing the line:
m = Basemap(projection='laea', resolution='h', lat_1=lat1, lat_2=lat2, lat_0=lat0, lon_0=lon0, lon_1=lon1,
to
m = Basemap(projection='lcc', resolution='h', lat_1=lat1, lat_2=lat2, lat_0=lat0, lon_0=lon0, lon_1=lon1,
@dtonggmu
Looked into the projection. We are using he Lambert Equal Area projection. We can simply change the projection in the plotting script to use Lamber Conformal Conic. simply by changing the line:
to
This will give a plot something like this