Segmentation fault (core dumped)

From: Sabeerali(sebi) <sabeerl_at_nyahnyahspammersnyahnyah>
Date: Mon Jun 07 2010 - 23:25:35 MDT

Hi Mary and Dennis,

 The size of each variable is giving here. I commented all the lines after
rf_clim=clmDayTLL(rf,yyyyddd). Still the problem is there..If we comment
rf_clim=clmDayTLL(rf,yyyyddd) then no problem..So I think problem will be in
clmDayTLL(rf,yyyyddd)..Also I am not getting the size of rf_clim...When I
include printVarSummary(rf_clim), it shows Segmentation fault.
Thanks once again..

Variable: rf
Type: float
Total Size: -1902903296 bytes
            598016000 values

Variable: time
Type: double
Total Size: 146000 bytes
            18250 values

Variable: TIME
Type: float
Total Size: 438000 bytes
            109500 values
Variable: year
Type: integer
Total Size: 73000 bytes
            18250 values

Variable: month
Type: integer
Total Size: 73000 bytes
            18250 values
Variable: day
Type: integer
Total Size: 73000 bytes
            18250 values
Variable: ddd
Type: integer
Total Size: 73000 bytes
            18250 values

Variable: yyyyddd
Type: integer
Total Size: 73000 bytes
            18250 values

On 7 June 2010 20:33, <ncl-talk-request@ucar.edu> wrote:

> Send ncl-talk mailing list submissions to
> ncl-talk@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@ucar.edu
>
> You can reach the person managing the list at
> ncl-talk-owner@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: ncl-talk Digest, Vol 79, Issue 10 (Sabeerali(sebi))
> 2. Segmentation fault (core dumped). (Sabeerali(sebi))
> 3. How to calculate POP (Principal Oscillation Pattern) with NCL
> (???)
> 4. ncl fatal:ContourPlotDraw: CTTMTL/ICAEDG - LOGIC ERROR
> (Ekin Akoglu)
> 5. ncl fatal:ContourPlotDraw: CTTMTL/ICAEDG - LOGIC ERROR
> (Ekin Akoglu)
> 6. Re: length of full path (Fan Fang)
> 7. Re: How to calculate POP (Principal Oscillation Pattern) with
> NCL (Dennis Shea)
> 8. Re: Segmentation fault (core dumped). (Dennis Shea)
> 9. Re: Segmentation fault (core dumped). (Mary Haley)
> 10. Contouring temperature with GFS dataset. (Caroline Serraud)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 7 Jun 2010 09:44:17 +0530
> From: "Sabeerali(sebi)" <sabeerl@gmail.com>
> Subject: Re: ncl-talk Digest, Vol 79, Issue 10
> To: ncl-talk@ucar.edu
> Message-ID:
> <AANLkTinKDoVxW74uNOrIA_QJkYJirtuUANTT-flX1HZY@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear Dennis,
>
> Thanks for your reply..I don't think my variable size is more than 2GB.
> Also
> I am not using any reordering and more recent *_n_Wrap function.. The
> script
> is given below..You can go through...
> Thanks again...
>
> 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/diagnostics_cam.ncl"
> begin
>
> in2=addfile("/home/sabeer/AR4_20c3m/ncar_csm/pr_A2.nc","r")
> rf=in2->pr ; daily data from 1950 to 1999
>
> time=in2->time
> TIME = ut_calendar(time, 0)
> year = floattointeger( TIME(:,0) )
> month = floattointeger( TIME(:,1) )
> day = floattointeger( TIME(:,2) )
> ddd = day_of_year(year, month, day)
> yyyyddd = year*1000 + ddd
>
> rf_clim=clmDayTLL(rf,yyyyddd)
> smth_rf_clim=smthClmDayTLL(rf_clim,3)
> rf_anom=calcDayAnomTLL(rf,yyyyddd,smth_rf_clim)
>
> spd=1
> bpf=(/20,100,201/)
> rfa_f=band_pass_latlon_time (rf_anom,spd, bpf,False)
> ;================saving to netcdf==============================
> diro ="./"
> filo = "dly.rfa.filt.nc"
> system("/bin/rm -f "+diro+filo)
> f = addfile(diro+filo, "c")
> f->rfa_filt =rfa_f
> end
>
>
> On 6 June 2010 23:30, <ncl-talk-request@ucar.edu> wrote:
>
> > Send ncl-talk mailing list submissions to
> > ncl-talk@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@ucar.edu
> >
> > You can reach the person managing the list at
> > ncl-talk-owner@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. Segmentation fault (core dumped). (Sabeerali(sebi))
> > 2. Re: Segmentation fault (core dumped). (Dennis Shea)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Sun, 6 Jun 2010 13:51:30 +0530
> > From: "Sabeerali(sebi)" <sabeerl@gmail.com>
> > Subject: Segmentation fault (core dumped).
> > To: ncl-talk@ucar.edu
> > Message-ID:
> > <AANLkTil5IaeeJ7kX-8cTVGmse970Y89-elbqxZDPwA5j@mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > I am having the daily data of 60 years...I am not able to do any
> > calculation
> > by using this large data set in NCL. ...It shows segmentation fault (core
> > dumped).. I think it is because that I encountered memory limit..Actually
> I
> > want to the band pass filter..It works nicely if I reduce size of data
> set
> > (means for 30 years daily) . Is there any solution if I encounter the
> > memory
> > limit?
> >
> > Thanks for your time..
> >
> > --
> > **********************************
> > Sabeerali
> > Climate and Global Modeling Division
> > Indian Institute of Tropical Meteorology
> > Pashan, Pune, 411 008
> > ****************************************
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> >
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100606/6addf6e9/attachment.html
> >
> > ------------------------------
> >
> > Message: 2
> > Date: Sun, 06 Jun 2010 08:33:44 -0600
> > From: Dennis Shea <shea@ucar.edu>
> > Subject: Re: Segmentation fault (core dumped).
> > To: "Sabeerali(sebi)" <sabeerl@gmail.com>
> > Cc: ncl-talk@ucar.edu
> > Message-ID: <4C0BB1C8.6010900@ucar.edu>
> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed
> >
> > If you are using more memory than is available on your machine,
> > the answer is "no".
> >
> > [1] Is the variable size > 2GB? Currently, NCL does not allow
> > variables to be > 2GB.
> >
> > [2] Are you using functions that require "dimension reordering"
> > or, the more recent *_n [ *_n_Wrap] function that do not.
> >
> > D
> >
> >
> >
> > On 6/6/10 2:21 AM, Sabeerali(sebi) wrote:
> > > I am having the daily data of 60 years...I am not able to do any
> > > calculation by using this large data set in NCL. ...It shows
> > > segmentation fault (core dumped).. I think it is because that I
> > > encountered memory limit..Actually I want to the band pass filter..It
> > > works nicely if I reduce size of data set (means for 30 years daily) .
> > > Is there any solution if I encounter the memory limit?
> > >
> > > Thanks for your time..
> > >
> > > --
> > > **********************************
> > > Sabeerali
> > > Climate and Global Modeling Division
> > > Indian Institute of Tropical Meteorology
> > > Pashan, Pune, 411 008
> > > ****************************************
> > >
> > >
> > >
> > > _______________________________________________
> > > ncl-talk mailing list
> > > List instructions, subscriber options, unsubscribe:
> > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> >
> > ------------------------------
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk@ucar.edu
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> >
> > End of ncl-talk Digest, Vol 79, Issue 10
> > ****************************************
> >
>
>
>
> --
> **********************************
> Sabeerali
> Climate and Global Modeling Division
> Indian Institute of Tropical Meteorology
> Pashan, Pune, 411 008
> ****************************************
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100607/b1b38c7a/attachment.html
>
> ------------------------------
>
> Message: 2
> Date: Mon, 7 Jun 2010 09:45:26 +0530
> From: "Sabeerali(sebi)" <sabeerl@gmail.com>
> Subject: Segmentation fault (core dumped).
> To: ncl-talk@ucar.edu
> Message-ID:
> <AANLkTik77g4-UFpUofYOCbxslSEBRtc4iK0PpJ2xw0Ml@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> >
> >
> > Dear Dennis,
> >
> > Thanks for your reply..I don't think my variable size is more than 2GB.
> > Also I am not using any reordering and more recent *_n_Wrap function..
> The
> > script is given below..You can go through...
> > Thanks again...
> >
> > 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/diagnostics_cam.ncl"
> > begin
> >
> > in2=addfile("/home/sabeer/AR4_20c3m/ncar_csm/pr_A2.nc","r")
> > rf=in2->pr ; daily data from 1950 to 1999
> >
> > time=in2->time
> > TIME = ut_calendar(time, 0)
> > year = floattointeger( TIME(:,0) )
> > month = floattointeger( TIME(:,1) )
> > day = floattointeger( TIME(:,2) )
> > ddd = day_of_year(year, month, day)
> > yyyyddd = year*1000 + ddd
> >
> > rf_clim=clmDayTLL(rf,yyyyddd)
> > smth_rf_clim=smthClmDayTLL(rf_clim,3)
> > rf_anom=calcDayAnomTLL(rf,yyyyddd,smth_rf_clim)
> >
> > spd=1
> > bpf=(/20,100,201/)
> > rfa_f=band_pass_latlon_time (rf_anom,spd, bpf,False)
> > ;================saving to netcdf==============================
> > diro ="./"
> > filo = "dly.rfa.filt.nc"
> > system("/bin/rm -f "+diro+filo)
> > f = addfile(diro+filo, "c")
> > f->rfa_filt =rfa_f
> > end
> >
> >
> >
> > On 6 June 2010 23:30, <ncl-talk-request@ucar.edu> wrote:
> >
> >> Send ncl-talk mailing list submissions to
> >> ncl-talk@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@ucar.edu
> >>
> >> You can reach the person managing the list at
> >> ncl-talk-owner@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. Segmentation fault (core dumped). (Sabeerali(sebi))
> >> 2. Re: Segmentation fault (core dumped). (Dennis Shea)
> >>
> >>
> >> ----------------------------------------------------------------------
> >>
> >> Message: 1
> >> Date: Sun, 6 Jun 2010 13:51:30 +0530
> >> From: "Sabeerali(sebi)" <sabeerl@gmail.com>
> >> Subject: Segmentation fault (core dumped).
> >> To: ncl-talk@ucar.edu
> >> Message-ID:
> >> <AANLkTil5IaeeJ7kX-8cTVGmse970Y89-elbqxZDPwA5j@mail.gmail.com>
> >> Content-Type: text/plain; charset="iso-8859-1"
> >>
> >> I am having the daily data of 60 years...I am not able to do any
> >> calculation
> >> by using this large data set in NCL. ...It shows segmentation fault
> (core
> >> dumped).. I think it is because that I encountered memory
> limit..Actually
> >> I
> >> want to the band pass filter..It works nicely if I reduce size of data
> set
> >> (means for 30 years daily) . Is there any solution if I encounter the
> >> memory
> >> limit?
> >>
> >> Thanks for your time..
> >>
> >> --
> >> **********************************
> >> Sabeerali
> >> Climate and Global Modeling Division
> >> Indian Institute of Tropical Meteorology
> >> Pashan, Pune, 411 008
> >> ****************************************
> >> -------------- next part --------------
> >> An HTML attachment was scrubbed...
> >> URL:
> >>
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100606/6addf6e9/attachment.html
> >>
> >> ------------------------------
> >>
> >> Message: 2
> >> Date: Sun, 06 Jun 2010 08:33:44 -0600
> >> From: Dennis Shea <shea@ucar.edu>
> >> Subject: Re: Segmentation fault (core dumped).
> >> To: "Sabeerali(sebi)" <sabeerl@gmail.com>
> >> Cc: ncl-talk@ucar.edu
> >> Message-ID: <4C0BB1C8.6010900@ucar.edu>
> >> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
> >>
> >> If you are using more memory than is available on your machine,
> >> the answer is "no".
> >>
> >> [1] Is the variable size > 2GB? Currently, NCL does not allow
> >> variables to be > 2GB.
> >>
> >> [2] Are you using functions that require "dimension reordering"
> >> or, the more recent *_n [ *_n_Wrap] function that do not.
> >>
> >> D
> >>
> >>
> >>
> >> On 6/6/10 2:21 AM, Sabeerali(sebi) wrote:
> >> > I am having the daily data of 60 years...I am not able to do any
> >> > calculation by using this large data set in NCL. ...It shows
> >> > segmentation fault (core dumped).. I think it is because that I
> >> > encountered memory limit..Actually I want to the band pass filter..It
> >> > works nicely if I reduce size of data set (means for 30 years daily) .
> >> > Is there any solution if I encounter the memory limit?
> >> >
> >> > Thanks for your time..
> >> >
> >> > --
> >> > **********************************
> >> > Sabeerali
> >> > Climate and Global Modeling Division
> >> > Indian Institute of Tropical Meteorology
> >> > Pashan, Pune, 411 008
> >> > ****************************************
> >> >
> >> >
> >> >
> >> > _______________________________________________
> >> > ncl-talk mailing list
> >> > List instructions, subscriber options, unsubscribe:
> >> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
> >>
> >> ------------------------------
> >>
> >> _______________________________________________
> >> ncl-talk mailing list
> >> ncl-talk@ucar.edu
> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
> >>
> >> End of ncl-talk Digest, Vol 79, Issue 10
> >> ****************************************
> >>
> >
> >
> >
> > --
> > **********************************
> > Sabeerali
> > Climate and Global Modeling Division
> > Indian Institute of Tropical Meteorology
> > Pashan, Pune, 411 008
> > ****************************************
> >
>
>
>
> --
> **********************************
> Sabeerali
> Climate and Global Modeling Division
> Indian Institute of Tropical Meteorology
> Pashan, Pune, 411 008
> ****************************************
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100607/954d3ce2/attachment.html
>
> ------------------------------
>
> Message: 3
> Date: Mon, 7 Jun 2010 17:29:38 +0800
> From: ??? <annddyshow@gmail.com>
> Subject: How to calculate POP (Principal Oscillation
> Pattern) with NCL
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <AANLkTinsed9w5OmxvoCFxj8aLuiA7Z7Om6R0Ucqqgyoy@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear Madam/Sir,
> How to use NCL to do POP (Principal Oscillation Pattern) analysis? Is there
> any example?
>
> In the mailing list, I've got a similar thread
> http://www.ncl.ucar.edu/Support/talk_archives/2009/0265.html
>
> Thanks a lot!
> Xiao
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100607/03d46efc/attachment.html
>
> ------------------------------
>
> Message: 4
> Date: Mon, 07 Jun 2010 13:45:25 +0300
> From: Ekin Akoglu <ekin@ims.metu.edu.tr>
> Subject: ncl fatal:ContourPlotDraw: CTTMTL/ICAEDG - LOGIC
> ERROR
> To: ncl-talk@ucar.edu
> Message-ID: <4C0CCDC5.50502@ims.metu.edu.tr>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear Users,
>
> What does this "ncl fatal:ContourPlotDraw: CTTMTL/ICAEDG - LOGIC ERROR"
> error message mean?
>
> Regards,
> --
> *Ekin Akoglu*
>
> Institute of Marine Sciences
> Middle East Technical University
> P.O. Box 28, 33731
> Erdemli, Mersin
> Turkey
>
> Web: www.ims.metu.edu.tr
> Email: ekin@ims.metu.edu.tr <mailto:ekin@ims.metu.edu.tr>
> Phone: +90 324 521 34 34
> GSM: +90 506 554 03 90
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100607/e9695035/attachment.html
>
> ------------------------------
>
> Message: 5
> Date: Mon, 07 Jun 2010 14:37:29 +0300
> From: Ekin Akoglu <ekin@ims.metu.edu.tr>
> Subject: ncl fatal:ContourPlotDraw: CTTMTL/ICAEDG - LOGIC
> ERROR
> To: ncl-talk@ucar.edu
> Message-ID: <4C0CD9F9.6040000@ims.metu.edu.tr>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Sorry for cross-posting but I think the full error message will be more
> meaningful. Here the error I get:
>
> fatal:ContourPlotDraw: CTTMTL/ICAEDG - LOGIC ERROR
> ERROR 5 IN CTTMTL/ICAEDG - LOGIC ERROR
> fatal:ContourPlotDraw: CTMESH - UNCLEARED PRIOR ERROR
> fatal:ContourPlotDraw: draw error
> warning:WorkstationDeactivate: workstation not active or not opened
>
> And the code is:
>
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> >
> > begin
> >
> > inf = addfile("marmara.nc","r")
> >
> > scale_factor = inf->z@scale_factor
> > add_offset = inf->z@add_offset
> > topo = scale_factor * (inf->z) + add_offset
> >
> > lons = inf->x_range
> > lats = inf->y_range
> > latarray = fspan(lats(0),lats(1),108781)
> > lonarray = fspan(lons(0),lons(1),108781)
> >
> > wks = gsn_open_wks("x11","Black_Sea")
> > gsn_define_colormap (wks,"gui_default")
> >
> > res = True
> > res@cnFillOn = True ; turn on color
> > res@cnLinesOn = False ; no contour lines
> > res@gsnSpreadColors = True ; use full color map
> > res@lbLabelAutoStride = True ; every other label
> > res@cnRasterModeOn = True
> > res@gsnAddCyclic = False ; regional data
> >
> > res@cnLevelSelectionMode = "AutomaticLevels" ; automatic levels
> > ;res@cnMinLevelValF = 100.
> > ;res@cnMaxLevelValF = 3500.
> > ;res@cnLevelSpacingF = 100.
> >
> >
> > res@tiMainString = "Black Sea"
> > res@sfXArray = lonarray ; uses longitude as
> > plot x-axis
> > res@sfYArray = latarray ; uses latitude as
> > plot y-axis
> >
> > res@tfDoNDCOverlay = True ; do not transform data
> > res@pmTickMarkDisplayMode = "Always" ; get a box around the
> > plot
> >
> > plot = gsn_csm_contour(wks,topo,res)
> >
> > end
>
>
> Regards,
> --
> *Ekin Akoglu*
>
> Institute of Marine Sciences
> Middle East Technical University
> P.O. Box 28, 33731
> Erdemli, Mersin
> Turkey
>
> Web: www.ims.metu.edu.tr
> Email: ekin@ims.metu.edu.tr <mailto:ekin@ims.metu.edu.tr>
> Phone: +90 324 521 34 34
> GSM: +90 506 554 03 90
>
> --
> *Ekin Akoglu*
>
> Institute of Marine Sciences
> Middle East Technical University
> P.O. Box 28, 33731
> Erdemli, Mersin
> Turkey
>
> Web: www.ims.metu.edu.tr
> Email: ekin@ims.metu.edu.tr <mailto:ekin@ims.metu.edu.tr>
> Phone: +90 324 521 34 34
> GSM: +90 506 554 03 90
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100607/187d8371/attachment.html
>
> ------------------------------
>
> Message: 6
> Date: Mon, 7 Jun 2010 08:06:04 -0400
> From: Fan Fang <Fan.Fang@nasa.gov>
> Subject: Re: length of full path
> To: Dennis Shea <shea@ucar.edu>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID: <4C0CE0AC.9040500@nasa.gov>
> Content-Type: text/plain; charset="ISO-8859-1"; format=flowed
>
> I found that the limit actually has nothing to do with looping several
> plots. It seems to exist and be much smaller for a simple xy-plot, when
> loading the additional "gsn_csm.ncl". It is consistent with what I
> found from the more complicated time series script - 85 characters in my
> Linux system including the ".eps" extension.
>
> Is there any way in NCL to enlarge this limit?
>
> -Fan
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>
>
> nchar = 82 ; FAIL
> nchar = 81
>
> N=100
> x = new(N, integer)
> y = new(N, integer)
> do i=0,N-1
> x(i) = rand()
> y(i) = rand()
> end do
>
> char = new ( nchar, "character")
> ca = stringtochar("a")
> char = ca(0)
>
> pltName = chartostring( char )
> pltType = "eps"
>
> res = True
>
> wks = gsn_open_wks(pltType, pltName)
> plot = gsn_csm_xy(wks,x,y,res)
>
> Dennis Shea wrote:
> > It appears to be 220 characters on my MAC.
> >
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> >
> >
> > ;nchar = 221 ; FAIL
> > nchar = 220
> >
> > char = new ( nchar, "character")
> > ca = stringtochar("a")
> > char = ca(0)
> >
> > pltName = chartostring( char )
> > pltType = "eps"
> >
> > wks = gsn_open_wks(pltType, pltName)
> > gsn_draw_colormap(wks)
> >
> >
> > On 6/4/10 8:51 AM, Fang, Fan (GSFC-610.2)[ADNET SYSTEMS INC] wrote:
> >
> >> What is the character limit of a full path name of a graphic output file
> (e.g. .eps file) that can be specified in NCL?
> >>
> >> We had "segmentation fault" if the file name with path is too long.
> >>
> >> -Fan
> >> _______________________________________________
> >> ncl-talk mailing list
> >> List instructions, subscriber options, unsubscribe:
> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
>
>
>
> ------------------------------
>
> Message: 7
> Date: Mon, 07 Jun 2010 07:44:39 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: How to calculate POP (Principal Oscillation
> Pattern) with NCL
> To: ??? <annddyshow@gmail.com>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID: <4C0CF7C7.30301@ucar.edu>
> Content-Type: text/plain; charset=GB2312
>
> To my knowledge, no one has develeoped an NCL based
> Principal Oscillation Pattern function for NCL. I am sure
> a WWW search may help you find code in (say) fortran.
>
> Good luck
>
>
>
> On 6/7/10 3:29 AM, ??? wrote:
> > Dear Madam/Sir,
> > How to use NCL to do POP (Principal Oscillation Pattern) analysis? Is
> > there any example?
> > In the mailing list, I've got a similar thread
> > http://www.ncl.ucar.edu/Support/talk_archives/2009/0265.html
> > Thanks a lot!
> > Xiao
> >
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> ------------------------------
>
> Message: 8
> Date: Mon, 07 Jun 2010 07:55:01 -0600
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: Segmentation fault (core dumped).
> To: "Sabeerali(sebi)" <sabeerl@gmail.com>
> Cc: ncl-talk@ucar.edu
> Message-ID: <4C0CFA35.7090004@ucar.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Actually, you are using "dimension" reordering".
>
> EG: the "contributed.ncl" function: smthClmDay
>
> eg: z = clmDay($dn(1)$|:,$dn(2)$|:,$dn(0)$|:)
>
> ===
> I speculate that you are exceeding the memory of your system.
>
> ===
> The *_n and *_n_Wrap functions were only recently available.
> Hence, many functions in contributed still use the older
> 'dimension reordering'.
>
>
>
> On 6/6/10 10:15 PM, Sabeerali(sebi) wrote:
> >
> > Dear Dennis,
> >
> > Thanks for your reply..I don't think my variable size is more than
> > 2GB. Also I am not using any reordering and more recent *_n_Wrap
> > function.. The script is given below..You can go through...
> > Thanks again...
> >
> > 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/diagnostics_cam.ncl"
> > begin
> >
> > in2=addfile("/home/sabeer/AR4_20c3m/ncar_csm/pr_A2.nc","r")
> > rf=in2->pr ; daily data from 1950 to 1999
> >
> > time=in2->time
> > TIME = ut_calendar(time, 0)
> > year = floattointeger( TIME(:,0) )
> > month = floattointeger( TIME(:,1) )
> > day = floattointeger( TIME(:,2) )
> > ddd = day_of_year(year, month, day)
> > yyyyddd = year*1000 + ddd
> >
> > rf_clim=clmDayTLL(rf,yyyyddd)
> > smth_rf_clim=smthClmDayTLL(rf_clim,3)
> > rf_anom=calcDayAnomTLL(rf,yyyyddd,smth_rf_clim)
> >
> > spd=1
> > bpf=(/20,100,201/)
> > rfa_f=band_pass_latlon_time (rf_anom,spd, bpf,False)
> > ;================saving to netcdf==============================
> > diro ="./"
> > filo = "dly.rfa.filt.nc <http://dly.rfa.filt.nc>"
> > system("/bin/rm -f "+diro+filo)
> > f = addfile(diro+filo, "c")
> > f->rfa_filt =rfa_f
> > end
> >
> >
> >
> > On 6 June 2010 23:30, <ncl-talk-request@ucar.edu
> > <mailto:ncl-talk-request@ucar.edu>> wrote:
> >
> > Send ncl-talk mailing list submissions to
> > ncl-talk@ucar.edu <mailto:ncl-talk@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@ucar.edu <mailto:ncl-talk-request@ucar.edu>
> >
> > You can reach the person managing the list at
> > ncl-talk-owner@ucar.edu <mailto:ncl-talk-owner@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. Segmentation fault (core dumped). (Sabeerali(sebi))
> > 2. Re: Segmentation fault (core dumped). (Dennis Shea)
> >
> >
> >
> ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Sun, 6 Jun 2010 13:51:30 +0530
> > From: "Sabeerali(sebi)" <sabeerl@gmail.com
> > <mailto:sabeerl@gmail.com>>
> > Subject: Segmentation fault (core dumped).
> > To: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>
> > Message-ID:
> > <AANLkTil5IaeeJ7kX-8cTVGmse970Y89-elbqxZDPwA5j@mail.gmail.com
> > <mailto:
> AANLkTil5IaeeJ7kX-8cTVGmse970Y89-elbqxZDPwA5j@mail.gmail.com>>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > I am having the daily data of 60 years...I am not able to do any
> > calculation
> > by using this large data set in NCL. ...It shows segmentation
> > fault (core
> > dumped).. I think it is because that I encountered memory
> > limit..Actually I
> > want to the band pass filter..It works nicely if I reduce size
> > of data set
> > (means for 30 years daily) . Is there any solution if I
> > encounter the memory
> > limit?
> >
> > Thanks for your time..
> >
> > --
> > **********************************
> > Sabeerali
> > Climate and Global Modeling Division
> > Indian Institute of Tropical Meteorology
> > Pashan, Pune, 411 008
> > ****************************************
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> >
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100606/6addf6e9/attachment.html
> >
> > ------------------------------
> >
> > Message: 2
> > Date: Sun, 06 Jun 2010 08:33:44 -0600
> > From: Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>>
> > Subject: Re: Segmentation fault (core dumped).
> > To: "Sabeerali(sebi)" <sabeerl@gmail.com <mailto:
> sabeerl@gmail.com>>
> > Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>
> > Message-ID: <4C0BB1C8.6010900@ucar.edu
> > <mailto:4C0BB1C8.6010900@ucar.edu>>
> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed
> >
> > If you are using more memory than is available on your machine,
> > the answer is "no".
> >
> > [1] Is the variable size > 2GB? Currently, NCL does not allow
> > variables to be > 2GB.
> >
> > [2] Are you using functions that require "dimension reordering"
> > or, the more recent *_n [ *_n_Wrap] function that do not.
> >
> > D
> >
> >
> >
> > On 6/6/10 2:21 AM, Sabeerali(sebi) wrote:
> > > I am having the daily data of 60 years...I am not able to do
> any
> > > calculation by using this large data set in NCL. ...It shows
> > > segmentation fault (core dumped).. I think it is because that
> I
> > > encountered memory limit..Actually I want to the band pass
> > filter..It
> > > works nicely if I reduce size of data set (means for 30 years
> > daily) .
> > > Is there any solution if I encounter the memory limit?
> > >
> > > Thanks for your time..
> > >
> > > --
> > > **********************************
> > > Sabeerali
> > > Climate and Global Modeling Division
> > > Indian Institute of Tropical Meteorology
> > > Pashan, Pune, 411 008
> > > ****************************************
> > >
> > >
> > >
> > > _______________________________________________
> > > ncl-talk mailing list
> > > List instructions, subscriber options, unsubscribe:
> > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> >
> > ------------------------------
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> >
> > End of ncl-talk Digest, Vol 79, Issue 10
> > ****************************************
> >
> >
> >
> >
> > --
> > **********************************
> > Sabeerali
> > Climate and Global Modeling Division
> > Indian Institute of Tropical Meteorology
> > Pashan, Pune, 411 008
> > ****************************************
> >
> >
> >
> >
> > --
> > **********************************
> > Sabeerali
> > Climate and Global Modeling Division
> > Indian Institute of Tropical Meteorology
> > Pashan, Pune, 411 008
> > ****************************************
> >
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> ------------------------------
>
> Message: 9
> Date: Mon, 7 Jun 2010 08:43:59 -0600
> From: Mary Haley <haley@ucar.edu>
> Subject: Re: Segmentation fault (core dumped).
> To: "Sabeerali(sebi)" <sabeerl@gmail.com>
> Cc: ncl-talk@ucar.edu
> Message-ID: <A7E707D1-3E6E-426F-A124-1701DA5ADC3F@ucar.edu>
> Content-Type: text/plain; charset="us-ascii"
>
> Please add:
>
> printVarSummary(rf)
>
> after you read in pr and let us know what the size of the variable is.
>
> Also, do the same for other variables, like "time".
>
> Thanks,
>
> --Mary
>
> On Jun 6, 2010, at 10:15 PM, Sabeerali(sebi) wrote:
>
> >
> > Dear Dennis,
> >
> > Thanks for your reply..I don't think my variable size is more than 2GB.
> Also I am not using any reordering and more recent *_n_Wrap function.. The
> script is given below..You can go through...
> > Thanks again...
> >
> > 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/diagnostics_cam.ncl"
> > begin
> >
> > in2=addfile("/home/sabeer/AR4_20c3m/ncar_csm/pr_A2.nc","r")
> > rf=in2->pr ; daily data from 1950 to 1999
> >
> > time=in2->time
> > TIME = ut_calendar(time, 0)
> > year = floattointeger( TIME(:,0) )
> > month = floattointeger( TIME(:,1) )
> > day = floattointeger( TIME(:,2) )
> > ddd = day_of_year(year, month, day)
> > yyyyddd = year*1000 + ddd
> >
> > rf_clim=clmDayTLL(rf,yyyyddd)
> > smth_rf_clim=smthClmDayTLL(rf_clim,3)
> > rf_anom=calcDayAnomTLL(rf,yyyyddd,smth_rf_clim)
> >
> > spd=1
> > bpf=(/20,100,201/)
> > rfa_f=band_pass_latlon_time (rf_anom,spd, bpf,False)
> > ;================saving to netcdf==============================
> > diro ="./"
> > filo = "dly.rfa.filt.nc"
> > system("/bin/rm -f "+diro+filo)
> > f = addfile(diro+filo, "c")
> > f->rfa_filt =rfa_f
> > end
> >
> >
> >
> > On 6 June 2010 23:30, <ncl-talk-request@ucar.edu> wrote:
> > Send ncl-talk mailing list submissions to
> > ncl-talk@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@ucar.edu
> >
> > You can reach the person managing the list at
> > ncl-talk-owner@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. Segmentation fault (core dumped). (Sabeerali(sebi))
> > 2. Re: Segmentation fault (core dumped). (Dennis Shea)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Sun, 6 Jun 2010 13:51:30 +0530
> > From: "Sabeerali(sebi)" <sabeerl@gmail.com>
> > Subject: Segmentation fault (core dumped).
> > To: ncl-talk@ucar.edu
> > Message-ID:
> > <AANLkTil5IaeeJ7kX-8cTVGmse970Y89-elbqxZDPwA5j@mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > I am having the daily data of 60 years...I am not able to do any
> calculation
> > by using this large data set in NCL. ...It shows segmentation fault (core
> > dumped).. I think it is because that I encountered memory limit..Actually
> I
> > want to the band pass filter..It works nicely if I reduce size of data
> set
> > (means for 30 years daily) . Is there any solution if I encounter the
> memory
> > limit?
> >
> > Thanks for your time..
> >
> > --
> > **********************************
> > Sabeerali
> > Climate and Global Modeling Division
> > Indian Institute of Tropical Meteorology
> > Pashan, Pune, 411 008
> > ****************************************
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100606/6addf6e9/attachment.html
> >
> > ------------------------------
> >
> > Message: 2
> > Date: Sun, 06 Jun 2010 08:33:44 -0600
> > From: Dennis Shea <shea@ucar.edu>
> > Subject: Re: Segmentation fault (core dumped).
> > To: "Sabeerali(sebi)" <sabeerl@gmail.com>
> > Cc: ncl-talk@ucar.edu
> > Message-ID: <4C0BB1C8.6010900@ucar.edu>
> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed
> >
> > If you are using more memory than is available on your machine,
> > the answer is "no".
> >
> > [1] Is the variable size > 2GB? Currently, NCL does not allow
> > variables to be > 2GB.
> >
> > [2] Are you using functions that require "dimension reordering"
> > or, the more recent *_n [ *_n_Wrap] function that do not.
> >
> > D
> >
> >
> >
> > On 6/6/10 2:21 AM, Sabeerali(sebi) wrote:
> > > I am having the daily data of 60 years...I am not able to do any
> > > calculation by using this large data set in NCL. ...It shows
> > > segmentation fault (core dumped).. I think it is because that I
> > > encountered memory limit..Actually I want to the band pass filter..It
> > > works nicely if I reduce size of data set (means for 30 years daily) .
> > > Is there any solution if I encounter the memory limit?
> > >
> > > Thanks for your time..
> > >
> > > --
> > > **********************************
> > > Sabeerali
> > > Climate and Global Modeling Division
> > > Indian Institute of Tropical Meteorology
> > > Pashan, Pune, 411 008
> > > ****************************************
> > >
> > >
> > >
> > > _______________________________________________
> > > ncl-talk mailing list
> > > List instructions, subscriber options, unsubscribe:
> > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> >
> > ------------------------------
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk@ucar.edu
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> >
> > End of ncl-talk Digest, Vol 79, Issue 10
> > ****************************************
> >
> >
> >
> > --
> > **********************************
> > Sabeerali
> > Climate and Global Modeling Division
> > Indian Institute of Tropical Meteorology
> > Pashan, Pune, 411 008
> > ****************************************
> >
> >
> >
> > --
> > **********************************
> > Sabeerali
> > Climate and Global Modeling Division
> > Indian Institute of Tropical Meteorology
> > Pashan, Pune, 411 008
> > ****************************************
> > _______________________________________________
> > ncl-talk mailing list
> > 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/20100607/0d6b18e7/attachment.html
>
> ------------------------------
>
> Message: 10
> Date: Mon, 07 Jun 2010 17:03:06 +0200
> From: Caroline Serraud <cs@meteoconsult.fr>
> Subject: Contouring temperature with GFS dataset.
> To: ncl-talk@ucar.edu
> Message-ID: <4C0D0A2A.6010201@meteoconsult.fr>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hello everyone,
>
> I am facing some trouble to plot contour of 2m-temperature from GFS data
> downloaded from NCEP FTP site (for example, gfs.t18z.pgrb2f06)
> NCL crashes when I run my script with default contour fill mode (area
> fill mode) but is ok with RasterFill mode.
>
> But with a 0,5? resolution I would really appreciate area mode because
> raster mode is not that pretty (see example attached)
> Did someone already experience plotting filled contour for
> 2m-temperature from GFS?
>
> Caroline
>
> --
> Caroline Serraud - Ing?nieur Etudes & D?veloppements
> ------------------------------------------------------
> METEO CONSULT / La Chaine M?t?o - Groupe Figaro
> D?partement Informatique
> Domaine de Marsinval F-78540 Vernouillet, FRANCE
> T?l : 01 39 28 1990 - Fax : 01 39 71 85 31
> e-mail : cs@meteoconsult.fr
> ------------------------------------------------------
> Toute la m?t?o sur le Web : http://www.meteoconsult.fr
> Toute la m?t?o par t?l?phone : 3201
>
> -------------- next part --------------
> An embedded and charset-unspecified text was scrubbed...
> Name: tempe_gfs.ncl
> Url:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100607/b25be7b8/attachment.pl
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: 1_tmp2m_2010060709.png
> Type: image/png
> Size: 46690 bytes
> Desc: not available
> Url :
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100607/b25be7b8/attachment.png
>
> ------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk@ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> End of ncl-talk Digest, Vol 79, Issue 11
> ****************************************
>

-- 
**********************************
Sabeerali
Climate and Global Modeling Division
Indian Institute of Tropical Meteorology
Pashan, Pune, 411 008
****************************************

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 7 23:25:47 2010

This archive was generated by hypermail 2.1.8 : Fri Jun 11 2010 - 09:33:24 MDT