[ncl-talk] Read multiple nc files

Krishnamoorthy Chandramouli chandrakrishna.90 at gmail.com
Sat Dec 23 20:49:42 MST 2017


Hi

>From your error message it looks like you want to read ERA data

I think ERA interim data is in grib format.

I used to convert it to netcdf using ncl_convert2nc

then i read the file.

If you want to read nc4 file. you can do it directly using addfile.

I have read the nc4 IMERG data recently, i never encountered a problem.



*Regards,*
*Krishnamoorthy Chandramouli*

On Sun, Dec 24, 2017 at 9:06 AM, Wang Xinyue via ncl-talk <ncl-talk at ucar.edu
> wrote:

> Hi Karin,
>
> Thanks for the reply.
> However, it did not work, I got the same warning:
>
> On Dec 22, 2017, at 23:59, Karin Meier-Fleischer <meier-fleischer at dkrz.de>
> wrote:
>
> Try:
>
> fils = systemfunc("ls *.nc")
> setfileoption("nc", "FileStructure", "Advanced")
> f = addfiles(fils,"r")
>
> -Karin
>
> Am 22.12.17 um 12:20 schrieb Wang Xinyue:
>
> Hi Karin,
>
> Thanks for your prompt reply!
> The 199902 dataset is just one of the datasets I want to read, I selected
> that one randomly and every single dataset got the same warning as:
>
> warning:_NclOpenFile: Can not open file </.../.../merra2/MERRA2_200.
> tavgM_3d_cld_Np.199902.nc4.nc>; file format not supported or file is
> corrupted
>
>
>
> The ncl_filedump output is like this:
>
> Variable: f
> Type: file
> filename: MERRA2_400.tavgM_3d_cld_Np.201112.nc4
> path: /.../MERRA2_400.tavgM_3d_cld_Np.201112.nc4.nc
>    file global attributes:
>       History : Original file generated: Sat Jul  4 21:58:25 2015 GMT
>       Filename : MERRA2_400.tavgM_3d_cld_Np.201112.nc4
>       Comment : GMAO filename: d5124_m2_jan10.tavg3_3d_cld_
> Np.monthly.201112.nc4
>       Conventions : CF-1
>       Institution : NASA Global Modeling and Assimilation Office
>       References : http://gmao.gsfc.nasa.gov
>       Format : NetCDF-4/HDF-5
>       SpatialCoverage : global
>       VersionID : 5.12.4
>       TemporalRange : 1980-01-01 -> 2016-12-31
>       identifier_product_doi_authority : http://dx.doi.org/
>       ShortName : M2TMNPCLD
>       RangeBeginningDate : 2011-12-01
>       RangeEndingDate : 2011-12-31
>       GranuleID : MERRA2_400.tavgM_3d_cld_Np.201112.nc4
>       ProductionDateTime : Original file generated: Sat Jul  4 21:58:25
> 2015 GMT
>       LongName : MERRA2 tavg3_3d_cld_Np: 3d,3-Hourly,Time-Averaged,
> Pressure-Level,Assimilation,Cloud Diagnostics Monthly Mean
>       Title : MERRA2 tavg3_3d_cld_Np: 3d,3-Hourly,Time-Averaged,
> Pressure-Level,Assimilation,Cloud Diagnostics Monthly Mean
>       SouthernmostLatitude : -90.0
>       NorthernmostLatitude : 90.0
>       WesternmostLongitude : -180.0
>       EasternmostLongitude : 179.375
>       LatitudeResolution : 0.5
>       LongitudeResolution : 0.625
>       DataResolution : 0.5 x 0.625 (42 pressure levels)
>       Source : CVS tag: GEOSadas-5_12_4
>       Contact : http://gmao.gsfc.nasa.gov
>       identifier_product_doi : 10.5067/J9R0LXGH48JR
>       RangeBeginningTime : 00:00:00.000000
>       RangeEndingTime : 23:59:59.000000
>       DODS_EXTRA.Unlimited_Dimension : time
>       history : 2017-12-10 06:40:28 GMT Hyrax-1.13.4
> https://goldsmr5.gesdisc.eosdis.nasa.gov:443/opendap/
> MERRA2_MONTHLY/M2TMNPCLD.5.12.4/2011/MERRA2_400.tavgM_3d_
> cld_Np.201112.nc4.nc?CLOUD[0:0][0:41][0:360][0:575],time,
> lev,lat[0:360],lon[0:575
> <https://goldsmr5.gesdisc.eosdis.nasa.gov/opendap/MERRA2_MONTHLY/M2TMNPCLD.5.12.4/2011/MERRA2_400.tavgM_3d_cld_Np.201112.nc4.nc?CLOUD[0:0][0:41][0:360][0:575],time,lev,lat[0:360],lon[0:575>
> ]
>    dimensions:
>       time = 1
>       lev = 42
>       lat = 361
>       lon = 576
>    variables:
>       float CLOUD ( time, lev, lat, lon )
>          long_name : cloud_fraction_for_radiation
>          units : 1
>          _FillValue : 1e+15
>          missing_value : 1e+15
>          fmissing_value : 1e+15
>          vmax : 1e+15
>          vmin : -1e+15
>          valid_range : ( -1e+15, 1e+15 )
>          origname : CLOUD
>          fullnamepath : /CLOUD
>
>       double lat ( lat )
>          CLASS : DIMENSION_SCALE
>          NAME : lat
>          long_name : latitude
>          units : degrees_north
>          vmax : 1e+15
>          vmin : -1e+15
>          valid_range : ( -1e+15, 1e+15 )
>          origname : lat
>          fullnamepath : /lat
>
>       double lev ( lev )
>          CLASS : DIMENSION_SCALE
>          NAME : lev
>          long_name : vertical level
>          units : hPa
>          positive : down
>          vmax : 1e+15
>          vmin : -1e+15
>          valid_range : ( -1e+15, 1e+15 )
>          origname : lev
>          fullnamepath : /lev
>
>       double lon ( lon )
>          CLASS : DIMENSION_SCALE
>          NAME : lon
>          long_name : longitude
>          units : degrees_east
>          vmax : 1e+15
>          vmin : -1e+15
>          valid_range : ( -1e+15, 1e+15 )
>          origname : lon
>          fullnamepath : /lon
>
>       integer time ( time )
>          CLASS : DIMENSION_SCALE
>          NAME : time
>          long_name : time
>          units : minutes since 2011-12-01 01:30:00
>          time_increment : 60000
>          begin_date : 20111201
>          begin_time : 13000
>          vmax : 1e+15
>          vmin : -1e+15
>          valid_range : ( -1e+15, 1e+15 )
>          origname : time
>          fullnamepath : /time
> Thanks,
> Xinyue
>
> On Dec 22, 2017, at 19:16, Karin Meier-Fleischer via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
> Hi Xinyue,
>
> please, send the output of *ncl_filedump* from MERRA2_200.tavgM_3d_cld_Np.
> 199902.nc4.nc
>
> -Karin
>
>
> Am 22.12.17 um 10:16 schrieb Wang Xinyue via ncl-talk:
>
> Hi there,
>
> I met with a problem in reading nc files.
> I can use the following commands to read NCEP datasets:
>
> fils = systemfunc("ls /.../.../ncep/omega/*.nc”)
> f = addfiles(fils, "r”)
>
> But when I use the same commands to read ERA-I datasets and MERRA2
> datasets, it warned like this:
>
> warning:_NclOpenFile: Can not open file </.../.../merra2/MERRA2_200.
> tavgM_3d_cld_Np.199902.nc4.nc>; file format not supported or file is
> corrupted
>
> Is anyone there had similar experience before? And how did you fix it?
>
> Looking forward to your reply.
>
> Thanks,
> Xinyue
>
>
>
> _______________________________________________
> ncl-talk mailing listncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> --
> Dipl. Geophys. Karin Meier-Fleischer
> Visualization, NCL, CDO
> Application Support
>
> Deutsches Klimarechenzentrum GmbH (DKRZ)
> Bundesstrasse 45a - D20146 Hamburg - Germany
>
> Phone:    +49 (0)40 460094 126 <+49%2040%20460094126>
> Fax:      +49 (0)40 460094 270 <+49%2040%20460094270>
> E-Mail:   meier-fleischer at dkrz.de
> URL:      www.dkrz.de
>
> Geschäftsführer: Prof. Dr. Thomas Ludwig
> Sitz der Gesellschaft: Hamburg
> Amtsgericht Hamburg HRB 39784
>
> _______________________________________________
> 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/20171224/59db17e8/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PastedGraphic-3.png
Type: image/png
Size: 181076 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171224/59db17e8/attachment-0001.png>


More information about the ncl-talk mailing list