Re: How to use rip_cape_2d with radiosonde dataset?

Thank you for your help Denis,Sending my files through ftp link failed.
So I attached a zipped file containing data which I am using.Louis--- En date de : Sam 1.11.08, Dennis Shea <shea_at_ucar.edu> a écrit :
De: Dennis Shea <shea_at_ucar.edu>
Objet: Re: [ncl-talk] How to use rip_cape_2d with radiosonde dataset?
À: appopson_at_yahoo.fr
Cc: "Mary Haley" <haley_at_ucar.edu>
Date: Samedi 1 Novembre 2008, 15h54

Louis

offline

Can you ftp the input data files to:

ftp ftp.cgd.ucar.edu
anonymous
email
cd incoming
put ...
put ...

quit

Then send me an email telling me when you have done this.

THX
D

louis Vonder wrote:
>
> Thank you Mary,
>
> Here errors appearing when I am trying to use rip_cape_2d
>
> Copyright (C) 1995-2007 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 5.0.0
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> *** glibc detected *** ncl: double free or corruption (out): 0x0aeb87d8
***
> ======= Backtrace: =========
> /lib/libc.so.6[0xb7d25c37]
> /lib/libc.so.6(cfree+0x89)[0xb7d278f9]
> ncl(NclFree+0x21)[0x810922f]
> ncl(rip_cape_2d_W+0xf0a)[0x8200b5a]
> ncl(CallINTRINSIC_FUNC_CALL+0x47)[0x817f2c0]
> ncl(_NclExecute+0x377)[0x8182305]
> ncl(yyparse+0x76b)[0x80bdbed]
> ncl(main+0xbc4)[0x80b9464]
> /lib/libc.so.6(__libc_start_main+0xe0)[0xb7cd4450]
> ncl(vfprintf+0x9fd)[0x80b8811]
> ======= Memory map: ========
> 08048000-08d6b000 r-xp 00000000 08:05 4162260 /usr/local/ncl/bin/ncl
> 08d6b000-08e9a000 rwxp 00d22000 08:05 4162260 /usr/local/ncl/bin/ncl
> 08e9a000-0aee8000 rwxp 08e9a000 00:00 0 [heap]
> b7b00000-b7b21000 rwxp b7b00000 00:00 0
> b7b21000-b7c00000 ---p b7b21000 00:00 0
> b7c6b000-b7c75000 r-xp 00000000 08:05 987649 /lib/libgcc_s.so.1
> b7c75000-b7c76000 rwxp 0000a000 08:05 987649 /lib/libgcc_s.so.1
> b7c9a000-b7c9b000 rwxp b7c9a000 00:00 0
> b7c9b000-b7c9f000 r-xp 00000000 08:05 3772997 /usr/lib/libXdmcp.so.6.0.0
> b7c9f000-b7ca0000 rwxp 00003000 08:05 3772997 /usr/lib/libXdmcp.so.6.0.0
> b7ca0000-b7ca1000 rwxp b7ca0000 00:00 0
> b7ca1000-b7ca3000 r-xp 00000000 08:05 3772995 /usr/lib/libXau.so.6.0.0
> b7ca3000-b7ca4000 rwxp 00001000 08:05 3772995 /usr/lib/libXau.so.6.0.0
> b7ca4000-b7cbb000 r-xp 00000000 08:05 3772904 /usr/lib/libxcb.so.1.0.0
> b7cbb000-b7cbc000 rwxp 00016000 08:05 3772904 /usr/lib/libxcb.so.1.0.0
> b7cbc000-b7cbd000 r-xp 00000000 08:05 3772908
/usr/lib/libxcb-xlib.so.0.0.0
> b7cbd000-b7cbe000 rwxp 00000000 08:05 3772908
/usr/lib/libxcb-xlib.so.0.0.0
> b7cbe000-b7df1000 r-xp 00000000 08:05 987776 /lib/libc-2.7.so
> b7df1000-b7df2000 r-xp 00133000 08:05 987776 /lib/libc-2.7.so
> b7df2000-b7df4000 rwxp 00134000 08:05 987776 /lib/libc-2.7.so
> b7df4000-b7df7000 rwxp b7df4000 00:00 0
> b7df7000-b7df9000 r-xp 00000000 08:05 987782 /lib/libdl-2.7.so
> b7df9000-b7dfb000 rwxp 00001000 08:05 987782 /lib/libdl-2.7.so
> b7dfb000-b7e1e000 r-xp 00000000 08:05 987783 /lib/libm-2.7.so
> b7e1e000-b7e20000 rwxp 00022000 08:05 987783 /lib/libm-2.7.so
> b7e20000-b7e21000 rwxp b7e20000 00:00 0
> b7e21000-b7e2e000 r-xp 00000000 08:05 3773011 /usr/lib/libXext.so.6.4.0
> b7e2e000-b7e2f000 rwxp 0000d000 08:05 3773011 /usr/lib/libXext.so.6.4.0
> b7e2f000-b7f13000 r-xp 00000000 08:05 3772912 /usr/lib/libX11.so.6.2.0
> b7f13000-b7f16000 rwxp 000e4000 08:05 3772912 /usr/lib/libX11.so.6.2.0
> b7f34000-b7f3c000 rwxp b7f34000 00:00 0
> b7f3c000-b7f3d000 r-xp b7f3c000 00:00 0 [vdso]
> b7f3d000-b7f57000 r-xp 00000000 08:05 3817509 /lib/ld-2.7.so
> b7f57000-b7f59000 rwxp 00019000 08:05 3817509 /lib/ld-2.7.so
> bfa2b000-bfa65000 rwxp bffc6000 00:00 0 [stack]
> Aborted (core dumped)
>
>
> --- En date de : *Ven 31.10.08, Mary Haley /<haley_at_ucar.edu>/* a
écrit :
>
> De: Mary Haley <haley_at_ucar.edu>
> Objet: Re: [ncl-talk] How to use rip_cape_2d with radiosonde dataset?
> À: appopson_at_yahoo.fr
> Cc: "Dennis Shea" <shea_at_ucar.edu>
> Date: Vendredi 31 Octobre 2008, 18h24
>
> Hi Louis,
>
> I'm not sure what you are asking about. Is this script giving you
an error
> message?
> It's hard to just look at a script and determine what's wrong
with it
> without
> more information. It would be helpful to email the exact error message
you are
> getting, or if not, to email a description of what you mean by
> "there's something wrong with it".
>
> --Mary
>
> On Oct 31, 2008, at 9:28 AM, louis Vonder wrote:
>
> > > Sorry Denis,
> > > In a previous message I send the script.
> > > > Here the script I am talking about
> > > > Thanks
> > > > > ;
***********************************************
> > 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
> > nrow = numAsciiRow("pres_04_05_07")
> > ncol = numAsciiCol("pres_04_05_07")
> > > pres = asciiread("pres_04_05_07", (/nrow,
ncol/),
> "float") ; pressure values
> > PP = new ( (/nrow, 1, 1/), typeof(pres), getFillValue(pres))
> > > z = asciiread("z_04_05_07", (/nrow, ncol/),
"float")
> ; geopotential heights
> > ZZ = new ( (/nrow, 1, 1/), typeof(z), getFillValue(z))
> > > t = asciiread("t_04_05_07", (/nrow, ncol/),
"float")
> ; temperature values
> > TT = new ( (/nrow, 1, 1/), typeof(t), getFillValue(t))
> > > rhum = asciiread("rhum_04_05_07", (/nrow,
ncol/),
> "float") ; relative humidity
> > td = asciiread("td_04_05_07", (/nrow, ncol/),
> "float") ; dew point temperature
> > > psfc = asciiread("psfc_04_05_07", (/1/),
"float") ;
> pressure at surface
> > PS = new ( (/1, 1/), typeof(psfc), getFillValue(psfc))
> > > zsfc = asciiread("zsfc_04_05_07", (/1/),
"float") ;
> geopotential at surface
> > ZS = new ( (/1, 1/), typeof(zsfc), getFillValue(zsfc))
> > > > q = mixhum_ptd (pres, td, 2) ; specific
humidities
> > QQ = new ( (/nrow, 1, 1/), typeof(q), getFillValue(q))
> > > > cinfo = rip_cape_2d(PP, TT, QQ, ZZ, ZS, PS,
True)
> > > cape = cinfo(0)
> > cin = cinfo(1)
> > lcl = cinfo(2)
> > lfc = cinfo(3)
> > > print(cape)
> > > end
> > > > --- En date de : Ven 31.10.08, Dennis Shea
<shea_at_ucar.edu> a écrit
> :
> > De: Dennis Shea <shea_at_ucar.edu>
> > Objet: Re: [ncl-talk] How to use rip_cape_2d with radiosonde
dataset?
> > À: appopson_at_yahoo.fr
> > Date: Vendredi 31 Octobre 2008, 14h02
> > > Louis You did not send the script. Most of us are not
paid by the NCL
> project. We want to help people but we must do our own jobs also.
Please try
> print statements. That is what we would do. Then, sfter you have tried
different
> thing, send us a *clean* script. Often, people send us scripts with
mnay
> commented lines and the code is poor. Our time is valuable too.
Regards D louis
> Vonder wrote: > > Hi every body, > > Someone found what is
wrong
> with my script? > > Thanks > > --- En date de : *Mer
29.10.08, Mary
> Haley /<haley_at_ucar.edu>/* a écrit : > > De: Mary Haley
> <haley_at_ucar.edu> > Objet: Re: [ncl-talk] How to use
rip_cape_2d with
> radiosonde dataset? > À: "Dennis Shea"
<shea_at_ucar.edu> >
> Cc: appopson_at_yahoo.fr > Date: Mercredi 29 Octobre 2008, 19h31 >
> Minor
> correction, conform_dims is in V5.0.1, not 5.0.0. We hope to >
release 5.0.1
> in late November. > > --Mary > > > On Wed, 29 Oct 2008,
Dennis
> Shea wrote: > > > TT = new ( (/1,nlevel,1,1/), typeof(T),
> getFillValue(T)) > > > or, with the yet to be released, 5.0.0
> >
> > TT = conform_dims( (/1,nlevel,1,1/), T, 1) > > > Good
luck >
> > > louis Vonder wrote: > >> >> Many Thanks
Dennis, >
> >> >> Data I am using is already ordered top to bottom as
indicated
> in the > document. > >> >> But, how to create
nominal
> multi-dimensional arrays ? > >> >> Thanks > >>
>>
> --- En date de : *Mer 29.10.08, Dennis Shea /<shea_at_ucar.edu>/* a
>
> écrit : > >> >> De: Dennis Shea <shea_at_ucar.edu>
>
> >> Objet: Re: [ncl-talk] How to use rip_cape_2d with radiosonde
>
> dataset? > >> À: appopson_at_yahoo.fr > >> Cc:
ncl-talk_at_ucar.edu
> > >> Date: Mercredi 29 Octobre 2008, 15h21 > >>
>> As
> noted in the documentation, there is very limited support >
>> for this
> routine. > >> >> [1] > >> You could create
nominal
> multi-dimensional arrays. > >> >> T[nlevel] ==>
TT(nlevel,1,1)
> or TT(1,nlevel,1,1) > >> etc > >> >> [2] >
>>
> >> As noted in the documentation, the data must be top to bottom
>
> [sfc] > >> >> if not, usee the ::-1 syntax to reverse
the order.
> > >> >> TT(:,0,0) = T(::-1) > >> >> Good
luck >
> >> >> louis Vonder wrote: > >> > Dear ncl
users, >
> >> > > I want to use the NCL "rip_cape_2d" to
>
> calculate CAPE and CIN. > >> > > For this purpose I am
using a
> radio sonde data. > >> > > But I found the following
errors >
> >> > > fatal:rip_cape_2d: The p, t, q, and z arrays must
be > 3
> or 4-dimensional > >> > fatal:Execute: Error occurred at
or near
> line 27 in file > cape_crash.ncl > >> > > This
function
> cannot be used for 1D dataset? > >> > > Here my script
>
> >> > > ; ***********************************************
>
> >> > 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 > >> > nrow =
> numAsciiRow("pres_04_05_07") > >> > ncol =
> numAsciiCol("pres_04_05_07") > >> > > pres =
> asciiread("pres_04_05_07", > (/nrow, ncol/), >
>>
> "float") ; pressure values > >> > z =
> asciiread("z_04_05_07", (/nrow, ncol/), >
"float") >
> >> ; geopotential heights > >> > t =
> asciiread("t_04_05_07", (/nrow, ncol/), >
"float") >
> >> ; temperature values > >> > rhum =
> asciiread("rhum_04_05_07", (/nrow, ncol/), > >>
> "float") ; relative humidity > >> > td =
> asciiread("td_04_05_07", (/nrow, ncol/), > >>
> "float") ; dew point temperature > >> > >
psfc =
> asciiread("psfc_04_05_07", (/1/), > "float") ;
>
> >> pressure at surface > >> > zsfc =
> asciiread("zsfc_04_05_07", (/1/), > "float") ;
>
> >> geopotential at surface > >> > > q =
mixhum_ptd
> (pres, td, 2) ; specific humidities > >> > > cinfo =
> rip_cape_2d(pres, t, q, z, zsfc, psfc, > True) > >> >
> cape =
> cinfo(0) > >> > cin = cinfo(1) > >> > lcl =
cinfo(2)
> > >> > lfc = cinfo(3) > >> > > end >
>> >
> > Thanks > >> > > > > >
>
------------------------------------------------------------------------ >
> >> > > _______________________________________________
> >>
> > ncl-talk mailing list > >> > List instructions,
subscriber
> options, unsubscribe: > >> >
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk > >>
>> >
> > > >
>
------------------------------------------------------------------------ >
> > _______________________________________________ > ncl-talk
mailing list
> > List instructions, subscriber options, unsubscribe: >
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> > > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk