Re: Meridonal Temperature Gradient

Hi Soumik,

The error message is indicating that the dimension number you're using in the call to center_finite_diff_n is invalid. You are using index=3, but your input variable is being subscripted such that it is only three dimensions, so the index must be either 0, 1, or 2.

Dimension numbering starts at 0, so if you want the leftmost dimension, use index=0. If you want the second from the left dimension, use index=1, and so on.

I'll take a look at example 6 on the function doc page. Thanks for pointing this out.

--Mary

On May 2, 2011, at 11:58 AM, sbasu@iarc.uaf.edu wrote:

> Hi,
>
> I am trying to use canter_finite _difference_n function to calculate
> meridional temperature gradient.
> But I am getting an error saying
>
> "Invalid dimension index for calculating the center finite difference"
>
> Here is the portion of my script where I am doing this:
>
> nlat = 128
> nlon =256
> dlon = (lon(2)-lon(1))*0.0174533 ; convert to radians
> ; pre-allocate space
> dLon = dble2flt(dlon)
> Lon = dble2flt(lon)
>
> dTdY = new((/124,15,nlat,nlon/),typeof(tempm))
>
> do nl =0,nlon-1 ; loop over each longitude
> dY = 6378388.*cos(0.0174533*Lon(nl))*dLon ; constant at this latitude
> printVarSummary(dY)
> dTdY(:,:,:,nl)= center_finite_diff_n (tempm(:,:,:,nl),dY,True,0,3)
> end do
>
> dtdy = flt2dble(dTdY)
> printVarSummary(dtdy)
> printMinMax(dtdy, True)
> printMinMax(dTdY, True)
>
> Also in the example 6 of documentation of "center_finite_difference_n"
> there are typeo in this line:
>
> dTdX(:,:,nl,:) = center_finite_diff_n (T(:,:,nl;,:), dX , True,0,0)
> it should be I think
> dTdX(:,:,nl,:) = center_finite_diff_n (T(:,:,nl,:), dX , True,0,2)
>
> Thanks,
>
> --
> *************************************
> Soumik Basu
> Graduate Student, Research Assistant
> International Arctic Research Center
> University of Alaska Fairbanks
> Fairbanks,Alaska,USA
> **************************************
>
>
> -----------------------------------------
> This email was sent using SquirrelMail.
> "Webmail for nuts!"
> http://squirrelmail.org/
>
> _______________________________________________
> 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 Tue May 3 13:59:57 2011