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

[Setareh Rahimi]

Dear Adam,
So many thanks for your help and advice. And also sorry for the slow reply.
Best wishes,

On Fri, Aug 12, 2022 at 1:08 AM Adam Phillips <asphilli at ucar.edu> wrote:

> 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>
>


-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220820/6de848a7/attachment.htm>