Re: How to use rip_cape_2d with radiosonde dataset?

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
Received on Fri Oct 31 2008 - 06:50:30 MDT