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

[Lyndz]

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*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200323/7c4b793c/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2020-03-23 at 11.54.00 AM.png
Type: image/png
Size: 170608 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200323/7c4b793c/attachment.png>