[ncl-talk] plot the wrfoutput

[Dennis Shea]

I speculate:

replace

lat1=lat(61:148,:)
lon1=lon(:,36:168)

with

lat1=lat(61:148,36:168)
lon1=lon(61:148,36:168)

===========
Also, a good programming practice it to use variables rather than hard code
numbers

sLat = 61
nLat = 148
wLon = 36
eLon = 168


no3_1_w=ana->nh4_1(:,:,sLat:nLat,wLon:eLon)
no3_2_w=ana->nh4_2(:,:,sLat:nLat,wLon:eLon)
no3_3_w=ana->nh4_3(:,:,sLat:nLat,wLon:eLon)
no3_4_w=ana->nh4_4(:,:,sLat:nLat,wLon:eLon)

ana=addfile("nh4_w_wo_co3.nc","r")
no3_1_wo=ana->nh4_1(:,:,sLat:nLat,wLon:eLon)
no3_2_wo=ana->nh4_2(:,:,sLat:nLat,wLon:eLon)
no3_3_wo=ana->nh4_3(:,:,sLat:nLat,wLon:eLon)
no3_4_wo=ana->nh4_4(:,:,sLat:nLat,wLon:eLon)

lat1=lat(sLat:nLat,wLon:eLon)
lon1=lon(sLat:nLat,wLon:eLon)



On Tue, Apr 13, 2021 at 10:33 AM Anahita Amiri Farahani via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hello,
>
> I tried to plot wrfout but for a smaller domain than the original one (sub
> domain) but it does not work. I have attached the code and the plot, also
> the desired plot.
>
> Thank you,
> Ana
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210414/d1360676/attachment.html>