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

Dennis Shea shea at ucar.edu
Fri Nov 1 13:01:25 MDT 2019


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/a5a15c0d/attachment.html>


More information about the ncl-talk mailing list