[ncl-talk] Plotting multiple .nc4 into one plot

[Rashed Mahmood]

Hi Lyndz,
I do not exactly what do you mean by calling each file and combining them
in one plot? do you mean a panel plot of all time-steps which could be done
using a simple loop, or do you intend to avg/sum all time steps? Both
should be straightforward when using addfiles.

Cheers,
Rashed




On Sun, Mar 22, 2020 at 9:08 PM Lyndz via ncl-talk <ncl-talk at ucar.edu>
wrote:

> Dear NCL-Experts,
>
> *[Description of the Problem]*
>
> I have a folder containing* interpolated* level 2 TRMM data (swath
> data).  Attached is the screenshot of the directory.
>
> The file names have the following format.
>
> *********************************************
> For the interpolated data, the convention is:
>
> *[version]_[orbit].[yyyymmdd].[hhmmss]_region.nc <http://region.nc> where*
>
> *[version]* = currently 'TPR7_uw1' which means this is based on TRMM-PR
> version 7 data and the UW's version 1 product
> *[orbit]* = five-digit orbit number
> *[yyyymmdd*] = year, month and day
> *[hhmmss]* = hour, minute and second
> *r* = three letter abbreviation for the region
> <http://trmm.atmos.washington.edu/regions.html> of interest
> *************************************
>
> Ncdump of one file gives the following:
> ***********************************************
>
> dimensions:
>
> time = UNLIMITED ; // (1 currently)
>
> longitude = 274 ;
>
> latitude = 235 ;
>
> altitude = 80 ;
>
> altitude_lh = 19 ;
>
> variables:
>
> int time(time) ;
>
> time:units = "seconds since 1970-01-01 00:00:00 +0000" ;
>
> float corr_Zfactor(time, altitude, latitude, longitude) ;
>
> corr_Zfactor:units = "dBZ" ;
>
> corr_Zfactor:long_name = "2A25 corrected Z factor" ;
>
> corr_Zfactor:range = "0.f, 80.f" ;
>
> corr_Zfactor:missing_value = -999.f ;
>
> float longitude(longitude) ;
>
> longitude:units = "degrees" ;
>
> longitude:long_name = "longitude" ;
>
> float latitude(latitude) ;
>
> latitude:units = "degrees" ;
>
> latitude:long_name = "latitude" ;
>
> float altitude(altitude) ;
>
> altitude:units = "km" ;
>
> altitude:long_name = "altitude for corr_Zfactor MSL" ;
>
> float rain_type(time, latitude, longitude) ;
>
> rain_type:units = "none" ;
>
> rain_type:long_name = "simple 2A23 Rain Type" ;
>
> rain_type:stratiform = 1 ;
>
> rain_type:convective = 2 ;
>
> rain_type:others = 3 ;
>
> rain_type:missing_value = -999 ;
>
> int rain_type_original(time, latitude, longitude) ;
>
> rain_type_original:units = "none" ;
>
> rain_type_original:long_name = "copy of 2A23 Rain Type field" ;
>
> rain_type_original:stratiform_range = "100, 170" ;
>
> rain_type_original:convective_range = "200, 291" ;
>
> rain_type_original:other_range = "300, 313" ;
>
> rain_type_original:no_rain_value = -88 ;
>
> rain_type_original:missing_value = -999 ;
>
> float surf_rain(time, latitude, longitude) ;
>
> surf_rain:units = "mm/hr" ;
>
> surf_rain:long_name = "2A25 rainfall rate near the surface" ;
>
> surf_rain:missing_value = -999.f ;
> ***********************************************************
>
> *[What I want to do]*
> [1] I would like to plot all the files into a single shaded map. Is this
> possible in NCL?
> [2] I know I should use the *addfiles *when reading the files, I just
> don't know how to call each file and combine them into one plot.
>
>
> I'll appreciate any help on this.
>
> Sincerely,
>
> *Lyndz*
>
> _______________________________________________
> 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/20200323/a2aad008/attachment.html>