[ncl-talk] Vertical Feature Mask (VFM) diagram using CALIPSO

Sun Nov 27 15:52:07 MST 2022

You are  looping over multiple files.and each file has a different lat/lon
arrangement?

SInce you are using:

lat = hdf_file->Latitude(1650:1500,0)

?? Each file has two-dimensional lat/lon arrays with different values??

*rcm2points*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2points.shtml>:
where the output (target) lat/lon include you desired latitude?

Maybe others have a clearer idea of what is desired?

On Sun, Nov 27, 2022 at 3:30 AM Setareh Rahimi <setareh.rahimi at gmail.com>
wrote:

> Dear Dennis,
> Thanks for your email. I mentioned "exactly for latitude between 32 N
> and  38 N, please", because when I add one of the CALIPSO files, say file
> number 1, to the script and set the latitude as follows:
> lat = hdf_file->Latitude(1650:1500,0)
> profile = fcf(1650:1500,1165:5514)
>
> I gave a plot for, say X to Y. After adding another CALIPSO file, say file
> number 2, and making no changes to the above command lines I got the plot
> for another latitude say M-N.
> That is why I mentioned "exactly for latitude between 32 N and  38 N,
> please"
>
>
> On Sat, Nov 26, 2022 at 1:32 AM Dennis Shea <shea at ucar.edu> wrote:
>
>> re: "exactly for latitude between 32 N and  38 N, please?"
>>
>> The words "exactly" and "between" makes this unclear to me.
>>
>> Do you mean "exactly" 36N?  (32N + 36N)/2 = 36N
>>
>> Do you mean all latitudes "between" 32N and 36N?
>>
>> I am not very familiar with CALIPSO data grids.
>>
>> You may have to interpolate the data to some standard grid which includes
>> the latitude:
>>
>> One dimensional coordinate variable:
>>
>> NCL's coordinate subscripting:   {...}
>>   x = data(*{*36*}*,:)   ;  all data at or hear 36N
>>
>> To find indices (subscripts) corresponding to specific 'targets
>>
>> One-dimensional coordinates
>>
>> *ind* <https://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml>
>>
>>    il = ind(lat.ge.32 .and. lat.le.36)    ; lat(lat)
>>
>> Two-dimensional coordinates
>> *region_ind*
>> <https://www.ncl.ucar.edu/Document/Functions/Contributed/region_ind.shtml>
>> *getind_latlon2d*
>> <https://www.ncl.ucar.edu/Document/Functions/Contributed/region_ind.shtml>
>>
>> =============
>> THe HDF EOS group has many examples using Matlab, NCL, Python and IDL
>> *    http://hdfeos.org/zoo/* <http://hdfeos.org/zoo/>
>> perhaps, one of the NCL examples could help.
>>
>>
>> On Wed, Nov 23, 2022 at 2:12 AM Setareh Rahimi <setareh.rahimi at gmail.com>
>> wrote:
>>
>>> Dear Dennis,
>>> I could find the script I was looking for (
>>> http://hdfeos.org/zoo/MORE/LaRC/CALIPSO/CAL_LID_L2_VFM-Standard-V4-10.2009-07-13T08-27-14ZD.hdf.v.ncl
>>> ) under this page: http://hdfeos.org/zoo/LaRC_CALIPSO_ncl.php.
>>>
>>> I need to plot the following latitude: 32 N - 38 N.
>>>
>>> I should plots CALIPSO data for different months for my study area.
>>> However, as CALIPSO overpasses are not the same for a certain area, after
>>> setting the below command lines in the script I got plots for another
>>> latitude.
>>>
>>> lat = hdf_file->Latitude(1650:1500,0)
>>> profile = fcf(1650:1500,1165:5514)
>>>
>>> So how can I have plots exactly for latitude between 32 N and  38 N,
>>> please?
>>>
>>> Many thanks in advance,
>>> Best wishes,
>>>
>>>
>>>
>>>
>>>
>>> On Fri, Nov 18, 2022 at 3:42 AM Dennis Shea <shea at ucar.edu> wrote:
>>>
>>>> re: "... need to know how to indicate the names of available particles
>>>> in the atmosphere (like dust, Black carbon,...) using NCL"
>>>> ===
>>>> I must admit, I doin't understand "names of available particles"
>>>>
>>>> D
>>>>
>>>>
>>>> On Wed, Nov 16, 2022 at 11:04 PM Setareh Rahimi <
>>>> setareh.rahimi at gmail.com> wrote:
>>>>
>>>>> Dear Dennis,
>>>>> Thank you for your advice. I already visited this section of NCL page,
>>>>> but I need to know how to indicate the names of available particles in the
>>>>> atmosphere (like dust, Black carbon,...) using NCL.
>>>>> I appreciate any suggestions.
>>>>> Best wishes,
>>>>>
>>>>> On Wed, Nov 16, 2022 at 9:44 PM Dennis Shea <shea at ucar.edu> wrote:
>>>>>
>>>>>> To me, it looks like NCL was used to generate the plot.
>>>>>>
>>>>>> *https://www.ncl.ucar.edu/Applications/*
>>>>>> <https://www.ncl.ucar.edu/Applications/>
>>>>>>
>>>>>> Under Datasets: click* Calipso*
>>>>>>
>>>>>> These types of plots can be 'complicated'
>>>>>>
>>>>>> Maybe you can contact the person/group that created the VFM.png plot
>>>>>>
>>>>>>
>>>>>>  ..... Good Luck
>>>>>>
>>>>>> On Wed, Nov 16, 2022 at 9:37 AM Setareh Rahimi via ncl-talk <
>>>>>> ncl-talk at mailman.ucar.edu> wrote:
>>>>>>
>>>>>>>
>>>>>>> Dear all NCL users,
>>>>>>>
>>>>>>> I want to plot a Vertical Feature Mask (VFM) diagram using CALIPSO
>>>>>>> data, like the attached one. But I am not familiar with the type of
>>>>>>> data. So please advise me in this regard.
>>>>>>>
>>>>>>> Thank you so much for your attention and participation.
>>>>>>> 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
>>>>>>>
>>>>>>
>>>>>
>>>>> --
>>>>> S.Rahimi
>>>>>
>>>>>
>>>
>>> --
>>> S.Rahimi
>>>
>>>
>
> --
> S.Rahimi
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20221127/5d15adf1/attachment.htm>