[ncl-talk] Read multiple nc files

Wang Xinyue wang.xinyue.q6 at dc.tohoku.ac.jp
Sat Dec 23 20:36:36 MST 2017


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 <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/ <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 <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 <mailto: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 list
>>>> ncl-talk at ucar.edu <mailto:ncl-talk at ucar.edu>
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk <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
>>> Fax:      +49 (0)40 460094 270
>>> E-Mail:   meier-fleischer at dkrz.de <mailto:meier-fleischer at dkrz.de>
>>> URL:      www.dkrz.de <http://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 <mailto:ncl-talk at ucar.edu>
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk <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/d251ad61/attachment.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/d251ad61/attachment.png>


More information about the ncl-talk mailing list