Re: (no subject)

Assuming line 84 is:
       alpha(:,:,:,:) = (RGAS*t(:,:,:,:))/t&level

Use conform to 'broadcast' the one-dimensional
t&level to all dimensions of t t(0,1,2,3)
       http://www.ncl.ucar.edu/Document/Functions/Built-in/conform.shtml

       alpha(:,:,:,:) = (RGAS*t(:,:,:,:))/conform(t,t&level,1)

or, better

        aplha = (RGAS*t)/conform(t,t&level,1)

The latter is more efficient.

On 7/22/12 3:02 AM, gibies george wrote:
> Dear friends,
>
> I tried to convert omega to w using the following part of script.
>
>
> RGAS = 287. ; J/(kg-K) => m2/(s2 K)
> GRAV = 9.8 ; m/s2
> alpha(:,:,:,:) = (RGAS*t(:,:,:,:))/t&level ;;;;; alpha=1/rho unit =
> m3kg-1
>
> w = -omega*alpha/(GRAV)
>
> But as expected it shows error as
>
> fatal:Div: Number of dimensions do not match, can't continue
> fatal:Execute: Error occurred at or near line 84 in file /work/wind.ncl
>
> Can anybody suggest an efficient way to solve this problem ?
>
>
> --
> *Gibies George,
> CSIR-Senior Research Fellow,*
> *Monsoon Group,
> Climate and Global Modelling Division,
> Indian Institute of Tropical Meteorology,
> Dr. Homi Bhabha Road,
> NCL (P. O.), Pashan,
> Pune 411008, India.
> Tel. Off. +912025904306*
> *Mobile : +919028523910
> +919881707796
> *
>
> *http://sites.google.com/site/gibiesge/*
>
>
>
>
>
>
>
>
> _______________________________________________
> 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 Sun Jul 22 06:33:32 2012