[ncl-talk] Aerosol optical depth from WRF-chem

[Adam Phillips]

Hi Setareh,
The array that you are plotting, AOD550_2D,  is a scalar value. I do not
believe that gsn_csm_contour_map can plot a single value.
The end of your script:
AOD550_3D = a2* ((400/550)^angstrom_exponent)
 printVarSummary(AOD550_3D)
AOD550_2D = sum(AOD550_3D)
 printVarSummary(AOD550_2D)
res at tiMainString = "WRF-CHEM (DUST_5) " + times(nt)
res at gsnLeftString = x at description
plot = gsn_csm_contour_map(wks,AOD550_2D,res)

Your AO550_3D array is dimensioned 35 x 40 x 40, but there is no coordinate
variable information. I would suggest adding
copy_VarMeta(a2,AOD550_3D)
to copy metadata from your a2 array to AOD550_3D. To plot the spatial
array, you will also need to read in the XLONG/XLAT arrays (which are 3D)
and set them to lat1d/lon1d arrays attached as attributes:
AOD550_3D at lat2d = f->XLAT(nt,:,:)
AOD550_3D at lon2d = f->XLONG(nt,:,:)

Hope that gets you going down the right path. If you have any further
questions, please respond to the ncl-talk email list.
Adam




On Wed, Aug 10, 2022 at 3:58 PM Setareh Rahimi via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Dear NCL users,
>
> I need to plot aerosol optical depth(AOD) from WRF-chem outputs. I follow
> the structure( script attached), however,
>  wrfout_d01_2010-07-14_00:00:00
> <https://drive.google.com/file/d/1-9zPQ-P6cmvrYs86CCeiNB-VU66Nuluw/view?usp=drive_web>
> the NCL returned this error:
> Error: scalar_field: If the input data is 1-dimensional, you must set
> sfXArray and sfYArray to 1-dimensional arrays of the same length.
> My question is how can I set sfXArray and sfYArray to 1-dimensional
> arrays of the same length, please?
> Many thanks in advance,
> Best wishes,
> --
> S.Rahimi
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220811/50b78b85/attachment.htm>