[ncl-talk] ncl-talk Digest, Vol 147, Issue 8
Benjamin Alabi
alabibenjamin2 at gmail.com
Mon Feb 8 04:03:22 MST 2016
Goodday sir
I think i've added 75 to ' ir_u at _FillValue= 255 '. pls check this script
for any correction:
; A script to convert binary data to netcdf
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
file_name_read = "merg_2012090102_4km-pixel"
ir_u = byte2flt( fbindirread (file_name_read ,0, (/2,3298,9896/),"ubyte"))
+ 75.
lon_g = fspan(0.0182,360,9896)
lat_g = fspan(59.982,59.982-0.036383683*3298,3298)
lon_g at units = "degrees_east"
lon_g at long_name = "longitude"
lon_g at short_name = "longitude"
lon_g!0 = "lon"
lon_g&lon = lon_g
lat_g at units = "degrees_north"
lat_g at long_name = "latitude"
lat_g at short_name = "latitude"
lat_g!0 = "lat"
lat_g&lat = lat_g
ir_u!0 = "time"
ir_u!1 = "lat"
ir_u!2 = "lon"
ir_u&lat = lat_g
ir_u&lon = lon_g
ir_u at _FillValue= 255+75.
wks = gsn_open_wks("X11","")
res = True
res at cnFillOn = True
res at cnFillMode = "RasterFill"
res at cnLinesOn = False
plot = gsn_csm_contour_map(wks, ir_u(0,:,:), res)
ncout = addfile(file_name_read+".nc","c")
ncout->ir_data = ir_u
end
On Sat, Feb 6, 2016 at 8:00 PM, <ncl-talk-request at ucar.edu> wrote:
> Send ncl-talk mailing list submissions to
> ncl-talk at ucar.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> or, via email, send a message with subject or body 'help' to
> ncl-talk-request at ucar.edu
>
> You can reach the person managing the list at
> ncl-talk-owner at ucar.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of ncl-talk digest..."
>
>
> Today's Topics:
>
> 1. Re: Data conversion (Alan Brammer)
> 2. Re: gsn_csm_vector_scalar_map (Xi Chang)
> 3. Re: gsn_csm_vector_scalar_map (alesall)
> 4. Re: read and write netcdf-4 compressed files (David Brown)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 5 Feb 2016 14:10:24 -0500
> From: Alan Brammer <abrammer at albany.edu>
> Subject: Re: [ncl-talk] Data conversion
> To: Benjamin Alabi <alabibenjamin2 at gmail.com>
> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> Message-ID: <564569E1-65BD-4130-B4EA-872812D72277 at albany.edu>
> Content-Type: text/plain; charset="utf-8"
>
> > And you don't want to try to add 75.0 to the fill value in the script.
>
> Just to respond to this and clarify, the missing value in the
> documentation that was linked to is 330. (
> http://www.cpc.ncep.noaa.gov/products/global_precip/html/README <
> http://www.cpc.ncep.noaa.gov/products/global_precip/html/README> )
>
> When I put the example script together I first posted, I just saw what was
> obviously a missing value at 255, having forgotten the offset. 255 is a
> very reasonable value in the world of brightness temps, so apologies for
> that, in my haste I made a mistake there in the first script.
>
> So the real _FillValue should be 330 (i.e. 255+75 )
>
>
>
>
> > On 5 Feb 2016, at 09:15, HAVENER, KEVIN F GS-12 USAF ACC 14 WS/WXED <
> kevin.havener at us.af.mil> wrote:
> >
> > And you don't want to try to add 75.0 to the fill value in the script.
> >
> > Kevin Havener
> >
> > -----Original Message-----
> > From: ncl-talk-bounces at ucar.edu [mailto:ncl-talk-bounces at ucar.edu] On
> Behalf Of alesall
> > Sent: Friday, February 05, 2016 9:03 AM
> > To: Benjamin Alabi
> > Cc: ncl-talk at ucar.edu
> > Subject: Re: [ncl-talk] Data conversion
> >
> >
> > the error message is:
> >
> > fatal:Undefined identifier: (byte2flt) is undefined, can't continue
> >
> > it is not finding the function ?byte2flt"
> > on the documentation page, this one: <
> http://ncl.ucar.edu/Document/Functions/Contributed/byte2flt.shtml>
> > it says that such function is in the ?contributed? file.
> >
> > you need to add statements such as:
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl <
> http://ncl.ucar.edu/Document/Functions/Contributed/contrib.shtml> "
> > at the beginning of your script, actually before the ?begin?
> >
> > I have all of these in my scripts:
> >
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> >
> > hope this helps!
> >
> >
> >
> > --
> > Alessandra Giannini
> > IRI for Climate and Society - The Earth Institute at Columbia University
> P.O. Box 1000, Palisades NY 10964-8000
> > phone/fax: +1 845 680-4473/4864 - email: alesall at iri.columbia.edu
> >
> >
> > On Feb 5, 2016, at 8:57 AM, Benjamin Alabi <
> alabibenjamin2 at gmail.com> wrote:
> >
> > Goodday sir,
> >
> > I appreciate your effort in guiding me. I'm new in this platform,
> please dont be offended by sending the message three times, I'm just
> learning to get help from here.
> >
> > I used the script,
> > begin
> > file_name_read = "merg_2012090101_4km-pixel"
> >
> > ir_u = byte2flt( fbindirread (file_name_read ,0,
> (/2,3298,9896/),"ubyte")) +75.0
> >
> >
> > lon_g = fspan(0.0182,360,9896)
> > lat_g = fspan(59.982,59.982-0.036383683*3298,3298)
> >
> > lon_g at units = "degrees_east"
> > lon_g at long_name = "longitude"
> > lon_g at short_name = "longitude"
> > lon_g!0 = "lon"
> > lon_g&lon = lon_g
> >
> > lat_g at units = "degrees_north"
> > lat_g at long_name = "latitude"
> > lat_g at short_name = "latitude"
> > lat_g!0 = "lat"
> > lat_g&lat = lat_g
> >
> >
> > ir_u!0 = "time"
> > ir_u!1 = "lat"
> > ir_u!2 = "lon"
> >
> > ir_u&lat = lat_g
> > ir_u&lon = lon_g
> >
> > ir_u at _FillValue= 255+75.0
> >
> >
> > wks = gsn_open_wks("X11","")
> > res = True
> > res at cnFillOn = True
> > res at cnFillMode = "RasterFill"
> > res at cnLinesOn = False
> > plot = gsn_csm_contour_map(wks, ir_u(0,:,:), res)
> >
> > ncout = addfile(file_name_read+".nc","c")
> > ncout->ir_data = ir_u
> >
> > end
> >
> >
> > But it was generating this error:
> > fatal:Undefined identifier: (byte2flt) is undefined, can't continue
> > fatal:["Execute.c":7741]:Execute: Error occurred at or near line 41
> >
> >
> > Waiting for guide, Thanks.
> >
> >
> >
> > On Thu, Feb 4, 2016 at 8:07 PM, Alan Brammer <abrammer at albany.edu>
> wrote:
> >
> >
> > Forgot the offset, which is mentioned in the link Kevin
> sent.
> >
> >
> > ir_u = byte2flt( fbindirread (file_name_read ,0,
> (/2,3298,9896/),"ubyte")) +75.
> > ir_u at _FillValue= 255+75.
> >
> > "Each record is a 9896 x 3298 Fortran array of IR
> brightness temperatures that have been scaled to fit into 1-byte by
> subtracting "75" from each datum. Therefore it is necessary for the user to
> add a value of "75" to each data value when using the data."
> >
> >
> >
> > On 4 Feb 2016, at 13:30, Alan Brammer <
> abrammer at albany.edu> wrote:
> >
> > ok, wow. Don?t send the same message 3 times, if
> the first answers weren?t complete at least expand and attempt something
> not just copy paste the same lines again.
> >
> > as said before .Z means it?s compressed. gunzip
> merg_2016020123_4km-pixel.Z will uncompress it. Then you have a binary
> file as a small amount of googling would have told you.
> > The same search would have linked you to various
> messages on here and various scripts where other people have included
> simple examples of this file type (myself included). Here is an example of
> reading and making a simple plot. There are plenty of examples on how to
> write netcdf on the NCL website so you should develop code from there.
> Most simple example at bottom of script below.
> >
> https://www.ncl.ucar.edu/Applications/o-netcdf.shtml
> >
> > Read the documentation on the below functions and
> understand what is occurring .
> >
> >
> > begin
> > file_name_read = "merg_2016020123_4km-pixel"
> >
> > ir_u = byte2flt( fbindirread (file_name_read ,0,
> (/2,3298,9896/),"ubyte"))
> >
> >
> > lon_g = fspan(0.0182,360,9896)
> > lat_g = fspan(59.982,59.982-0.036383683*3298,3298)
> >
> > lon_g at units = "degrees_east"
> > lon_g at long_name = "longitude"
> > lon_g at short_name = "longitude"
> > lon_g!0 = "lon"
> > lon_g&lon = lon_g
> >
> > lat_g at units = "degrees_north"
> > lat_g at long_name = "latitude"
> > lat_g at short_name = "latitude"
> > lat_g!0 = "lat"
> > lat_g&lat = lat_g
> >
> >
> > ir_u!0 = "time"
> > ir_u!1 = "lat"
> > ir_u!2 = "lon"
> >
> > ir_u&lat = lat_g
> > ir_u&lon = lon_g
> >
> > ir_u at _FillValue= 255
> >
> >
> > wks = gsn_open_wks("X11","")
> > res = True
> > res at cnFillOn = True
> > res at cnFillMode = "RasterFill"
> > res at cnLinesOn = False
> > plot = gsn_csm_contour_map(wks, ir_u(0,:,:), res)
> >
> > ncout = addfile(file_name_read+?.nc?,?c?)
> > ncout->ir_data = ir_u
> >
> > end
> >
> >
> >
> >
> >
> >
> > On 4 Feb 2016, at 13:04, Benjamin Alabi <
> alabibenjamin2 at gmail.com> wrote:
> >
> > Thanks for your assistance sir. I've tried
> the cdo to get general information on the data, but it didn't work
> > I need help from someone that can guide me
> to convert this data format - 'merg_2012091000_4km-pixel.Z' to netcdf
> format.
> > I download Gobally-merged (60N-60S)
> pixel-resolution IR brightness temperature data (equivalent blackbody
> temps), merged from all available geostationary satellites (GOES-8/10,
> METEOSAT-7/5)
> > This is the link:
> >
> http://mirador.gsfc.nasa.gov/cgi-bin/mirador/presentNavigation.pl?tree=project&dataset=Global-merged%20IR%20Brightness%20Temperature%20Data&project=TRMM&dataGroup=Ancillary&version=001&CGISESSID=25c6337b80163b3743e8bea85d015558
> .
> >
> > I really want it to be in netcdf format.
> >
> >
> > Regard
> >
> >
> >
> >
> >
> _______________________________________________
> > 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
> >
> >
> >
> > _______________________________________________
> > 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/20160205/2fdc6d64/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Fri, 5 Feb 2016 21:19:10 +0100
> From: Xi Chang <xi.chang01 at gmail.com>
> Subject: Re: [ncl-talk] gsn_csm_vector_scalar_map
> To: Karin Meier-Fleischer <meier-fleischer at dkrz.de>
> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> Message-ID:
> <CAGfRhteOJRK4=Rvrnoua8a+JDA2HQwfe8Q2jMhN=
> cmtmrj7gVg at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi Karin,,
>
> I did, but it doesnt work for gsn_csm_vector_scalar_map, no idea why..
>
> Chang.
>
> On Fri, Feb 5, 2016 at 7:53 PM, Karin Meier-Fleischer <
> meier-fleischer at dkrz.de> wrote:
>
> > Hi Xi,
> >
> > have you taken a look at the first example of the vector examples page
> > http://ncl.ucar.edu/Applications/Scripts/vector_1.ncl
> >
> > When you want to control the minimum/maximum and interval of the contour
> > lines use
> >
> > res at cnLevelSelectionMode = "ManualLevels" ; set manual contour
> levels
> > res at cnMinLevelValF = 24.0 ; set min contour level
> > res at cnMaxLevelValF = 29 ; set max contour level
> > res at cnLevelSpacingF = 0.10 ; set contour spacing
> >
> > Bye,
> > Karin
> >
> > Am 05.02.16 um 19:00 schrieb Xi Chang:
> >
> > Hello...
> >
> > Could you please tell me how to control the contour interval for color
> > shading once im using this function *gsn_csm_vector_scalar_map*.
> >
> > ;===========PLOT=================
> >
> > wks = gsn_open_wks("png","overlay")
> >
> > res = True ; plot mods desired
> > plot = new(16,graphic)
> >
> > res at gsnDraw = False ; don't draw
> > res at gsnFrame = False ; don't advance frame
> >
> > ; vector rescoures
> > res at vcRefLengthF = 0.045
> > res at vcRefMagnitudeF = 20.0
> > res at vcRefAnnoOrthogonalPosF = -0.12
> > res at vcRefAnnoParallelPosF = 0.997
> > res at vcRefAnnoFontHeightF = 0.015
> > res at vcMinDistanceF = 0.03
> > res at lbLabelBarOn = False
> >
> > res at mpProjection = "AzimuthalEquidistant" ; choose map projection
> > res at mpLimitMode = "LatLon" ; use lat/lon coordinates to
> > limit area
> > res at mpMinLatF = 25.
> > res at mpMaxLatF = 75.
> > res at mpMinLonF = -80.
> > res at mpMaxLonF = 20.
> > res at gsnMaximize = True
> >
> >
> >
> >
> >
> > * res at gsnScalarContour = True res at cnFillOn =
> True
> > res at cnFillPalette = "gui_default" ; set color map
> > res at cnLinesOn = False ; turn off contour
> lines *
> >
> >
> > do i=0,15
> > plot(i) = gsn_csm_vector_scalar_map(wks,u(i,:,:),v(i,:,:), \
> > p(i,:,:),res)
> > end do
> >
> > resP = True ; modify the panel plot
> > resP at gsnMaximize = True
> > gsn_panel(wks,plot,(/4,4/),resP) ; now draw as one plot
> >
> >
> > end
> >
> >
> >
> >
> > _______________________________________________
> > ncl-talk mailing listncl-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/20160205/4e881daa/attachment-0001.html
>
> ------------------------------
>
> Message: 3
> Date: Fri, 5 Feb 2016 15:21:44 -0500
> From: alesall <alesall at iri.columbia.edu>
> Subject: Re: [ncl-talk] gsn_csm_vector_scalar_map
> To: Xi Chang <xi.chang01 at gmail.com>
> Cc: "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
> Message-ID: <257ABA37-BFD7-49C5-B6CB-193B6F81F39F at iri.columbia.edu>
> Content-Type: text/plain; charset="us-ascii"
>
> Hi, is it complaining that the spacing is too fine compared to the wide
> range between min and max values, and reverting to doing what it thinks
> best?
>
>
>
> --
> Alessandra Giannini
> IRI for Climate and Society - The Earth Institute at Columbia University
> P.O. Box 1000, Palisades NY 10964-8000
> phone/fax: +1 845 680-4473/4864 - email: alesall at iri.columbia.edu <mailto:
> alesall at iri.columbia.edu>
> > On Feb 5, 2016, at 3:19 PM, Xi Chang <xi.chang01 at gmail.com> wrote:
> >
> > Hi Karin,,
> >
> > I did, but it doesnt work for gsn_csm_vector_scalar_map, no idea why..
> >
> > Chang.
> >
> > On Fri, Feb 5, 2016 at 7:53 PM, Karin Meier-Fleischer <
> meier-fleischer at dkrz.de <mailto:meier-fleischer at dkrz.de>> wrote:
> > Hi Xi,
> >
> > have you taken a look at the first example of the vector examples page
> > http://ncl.ucar.edu/Applications/Scripts/vector_1.ncl <
> http://ncl.ucar.edu/Applications/Scripts/vector_1.ncl>
> >
> > When you want to control the minimum/maximum and interval of the contour
> lines use
> > res at cnLevelSelectionMode = "ManualLevels" ; set manual contour
> levels
> > res at cnMinLevelValF = 24.0 ; set min contour level
> > res at cnMaxLevelValF = 29 ; set max contour level
> > res at cnLevelSpacingF = 0.10 ; set contour spacing
> > Bye,
> > Karin
> >
> > Am 05.02.16 um 19:00 schrieb Xi Chang:
> >> Hello...
> >>
> >> Could you please tell me how to control the contour interval for color
> shading once im using this function gsn_csm_vector_scalar_map.
> >>
> >> ;===========PLOT=================
> >>
> >> wks = gsn_open_wks("png","overlay")
> >>
> >> res = True ; plot mods desired
> >> plot = new(16,graphic)
> >>
> >> res at gsnDraw = False ; don't draw
> >> res at gsnFrame = False ; don't advance frame
> >>
> >> ; vector rescoures
> >> res at vcRefLengthF = 0.045
> >> res at vcRefMagnitudeF = 20.0
> >> res at vcRefAnnoOrthogonalPosF = -0.12
> >> res at vcRefAnnoParallelPosF = 0.997
> >> res at vcRefAnnoFontHeightF = 0.015
> >> res at vcMinDistanceF = 0.03
> >> res at lbLabelBarOn = False
> >>
> >> res at mpProjection = "AzimuthalEquidistant" ; choose map projection
> >> res at mpLimitMode = "LatLon" ; use lat/lon coordinates to
> limit area
> >> res at mpMinLatF = 25.
> >> res at mpMaxLatF = 75.
> >> res at mpMinLonF = -80.
> >> res at mpMaxLonF = 20.
> >> res at gsnMaximize = True
> >>
> >> res at gsnScalarContour = True
> >> res at cnFillOn = True
> >> res at cnFillPalette = "gui_default" ; set color map
> >> res at cnLinesOn = False ; turn off contour
> lines
> >>
> >>
> >>
> >> do i=0,15
> >> plot(i) = gsn_csm_vector_scalar_map(wks,u(i,:,:),v(i,:,:), \
> >> p(i,:,:),res)
> >> end do
> >>
> >> resP = True ; modify the panel plot
> >> resP at gsnMaximize = True
> >> gsn_panel(wks,plot,(/4,4/),resP) ; now draw as one plot
> >>
> >>
> >> end
> >>
> >>
> >>
> >>
> >> _______________________________________________
> >> 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>
> >
> >
> > _______________________________________________
> > 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>
> >
> >
> > _______________________________________________
> > 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/20160205/8c35e975/attachment-0001.html
>
> ------------------------------
>
> Message: 4
> Date: Fri, 5 Feb 2016 14:37:26 -0700
> From: David Brown <dbrown at ucar.edu>
> Subject: Re: [ncl-talk] read and write netcdf-4 compressed files
> To: Dan Chen <dchen at ucar.edu>, ncl-talk <ncl-talk at ucar.edu>
> Message-ID:
> <
> CAC92F7fDd6iTx4N401jbbN-G6ieUzzgFEuW5TG9QhMHLaRvx3A at mail.gmail.com>
> Content-Type: text/plain; charset=UTF-8
>
> Hi Dan,
> The problem is a dimension mismatch. I am not sure why it does not get
> reported as such, but here is the issue:
> Your call to wrf_user_getvar asks for time step 0, granted the only
> time step since the Time dimension has only 1 element. But the result
> is that the return variable 'bb' is reduced from 4 dimensions (t,
> z,y,x) in the file variable to 3 (z,y,x) in the local variable.
>
> Then you are trying to assign this 3d variable to a 4d variable in the
> file. This causes the error.
> There are two ways to fix it.
> 1) Substitute -1 for 0 in the 3rd argument to wrf_user_getvar. You
> will then get Time as the leftmost dimension of a 4d variable, and the
> write can proceed without error.
> 2) Or leave the wrf_user_getvar call alone. You can write the 3d local
> variable into the 4d file variable by subscripting the left hand side
> of the assignment:
>
> in->PM2_5_DRY(0,:,:,:) = (/bb/)
>
> We will look into why this is not caught earlier with a more apropos
> "dimension mismatch" error to avoid the trip into the NetCDF library
> code and the resulting low-level error messages.
> -dave
>
>
> On Thu, Feb 4, 2016 at 8:16 PM, Dan Chen <dchen at ucar.edu> wrote:
> > Hi,
> >
> > My name is Dan Chen working at MMM, NCAR. I'm using ncl to read a
> variable
> > and write into the netcdf-4 compressed file on Yellowstone, but failed.
> >
> > The script is as blow:
> > ========================================
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> > load "$NCARG_ROOT//lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
> >
> > begin
> >
> > fils = "/glade/scratch/dchen/ncl_netcdf4/analysis.ensmean"
> > in = addfile(fils,"w")
> > bb = wrf_user_getvar(in,"PM2_5_DRY",0)
> > bb = bb*10
> > fils->PM2_5_DRY = (/bb/)
> >
> > end
> > ==========================================
> >
> >
> > The error is as below
> > =========================
> > Copyright (C) 1995-2015 - All Rights Reserved
> > University Corporation for Atmospheric Research
> > NCAR Command Language Version 6.3.0
> > The use of this software is governed by a License Agreement.
> > See http://www.ncl.ucar.edu/ for more details.
> > fatal:["NclNetCDF4.c":3073]:NclNetCDF4.c: Error in nc_def_var_chunking in
> > file (/glade/scratch/dchen/ncl_netcdf4/analysis.ensmean) for writing, at
> > line: 3073
> >
> > line 3074 of NclNetCDF4.c: NetCDF: Attempt to define var properties, like
> > deflate, after enddef.
> > ============================
> >
> >
> > If I change the input file to a normal netcdf file without the
> compression
> > option, it works fine. But I do need to process those netcdf-4 compressed
> > files. Can you let me know if there is anything I can try?
> >
> > Many thanks,
> > Dan
> >
> >
> >
> >
> >
> >
> > _______________________________________________
> > 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
>
>
> End of ncl-talk Digest, Vol 147, Issue 8
> ****************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160208/95b6eae0/attachment.html
More information about the ncl-talk
mailing list