[ncl-talk] 3d field from satellite L2 data

Debasish Hazra debasish.hazra5 at gmail.com
Fri Nov 1 14:53:22 MDT 2019


Thanks Dennis



On Fri, Nov 1, 2019 at 3:01 PM Dennis Shea <shea at ucar.edu> wrote:

> To what do the data dimensions refer? This is your responsibility.
> *The golden rule of data processing is: Know your data!*
> =====
> It is the user's responsibility to know dimension information.... not
> ncl-talk's responsibility.
> =====
> The ext variable [*Column_Optical_Depth_Aerosols_532*] look more like a *time
> series* than a vertical cross section.
> I am not sure to what fakeDim44=*1 *may refer. Perhaps. the '*nadir*'
> location.
>
> *ext* = hdf4_file->*Column_Optical_Depth_Aerosols_532*
> printVarSummary(ext)   ; [fakeDim43 | *3856*] x [fakeDim44 | *1*]     ;
> <== fakeDim44=*1 *may refer to the *'nadir' *of the satellite
> printMinMax(ext,0)
> ; ****This is NOT a vertical cross section****
> print("===")
>
> Latitude = hdf4_file->Latitude
> printVarSummary(Latitude)  ; [fakeDim2 |* 3856]* x [fakeDim3 |* 3*]  ;
> <== fakeDim3=*3 *may refer to the left/nadir/right scan locations
> printMinMax(Latitude,0) ;
> print("===")
>
> Longitude = hdf4_file->Longitude
> printVarSummary(Longitude)   [fakeDim2 |* 3856]* x [fakeDim3 |* 3*]
> printMinMax(Longitude,0)
>
> time    = hdf4_file->Profile_Time
> time at units = "seconds since 1993-01-01 00:00"
>
> ptinVarSummary(time)   ; [fakeDim6 | *3856]* x [fakeDim3 |* 3*]; ; <==
> fakeDim6=*3 *may refer to the left/nadir/right scan times
> printMinMax(time,0) ; Info:
> print("===")
>
>
>
>
> On Fri, Nov 1, 2019 at 10:56 AM Debasish Hazra via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
>> Thanks Rick for pointing that out. But unlike the examples in (
>> http://www.ncl.ucar.edu/Applications/Scripts/calipso_3.ncl), variable I
>> want to see has different dimensions. Example page dimensions are
>>
>> [fakeDim70 | 3001] x [fakeDim71 | 399]
>>
>> which is why altitude reading from text file can be assigned to one of
>> the dimension. Here it is not.
>>
>> Variable: ext
>> Type: float
>> Total Size: 15424 bytes
>>             3856 values
>> Number of Dimensions: 2
>> Dimensions and sizes:   [fakeDim43 | 3856] x [fakeDim44 | 1]
>> Coordinates:
>> Number Of Attributes: 4
>>   units :       NoUnits
>>   format :      Float_32
>>   valid_range : 0...3
>>   hdf_name :    Column_Optical_Depth_Aerosols_532
>>
>> How Do I modify that.
>>
>> Thanks
>> Debasish
>>
>>
>>
>>
>> On Fri, Nov 1, 2019 at 11:44 AM Rick Brownrigg <brownrig at ucar.edu> wrote:
>>
>>> From the docs on fspan:
>>>
>>> *    npts*
>>>
>>>     Number of equally-spaced points desired between *start* and *finish*.
>>> This value must be >= 2.
>>>
>>> http://ncl.ucar.edu/Document/Functions/Built-in/fspan.shtml
>>>
>>> You saying you want one -point between 0 and the length of your arrays.
>>>
>>> Rick
>>>
>>> On Fri, Nov 1, 2019 at 9:11 AM Debasish Hazra via ncl-talk <
>>> ncl-talk at ucar.edu> wrote:
>>>
>>>> Hi,
>>>> I am following NCL examples of Calipso to plot vertical distribution of
>>>> a variable. But I am getting error in fspan as (but I checked the variable
>>>> dimension) :
>>>> fatal:fspan: number of elements parameter is less-than-or-equal-to one,
>>>> can't continue
>>>> Attached is the code and also sent the file through ftp.
>>>>
>>>> Thanks in advance
>>>> Debasish
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191101/d39c9d76/attachment.html>


More information about the ncl-talk mailing list