Re: area average across 0deg Meridian

Use: http://www.ncl.ucar.edu/Document/Functions/Contributed/lonFlip.shtml

This is located in
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

Note: if your data are sst(time,lat,lon) then you do not need to use
         the dimension names. Dimension names should be used when
         you are reshaping [permuting the array order.]

         Also, if you are using NCL's coordimnate subscripting you must
         use the {...} syntax.

         sst = lonflip( sst )
        
         sstAvg = wgt_areaave( sst(:,{-30:30},{-60:15}), clat, 1.0, 0)

====
Also, you may want to use

   http://www.ncl.ucar.edu/Document/Functions/Built-in/wgt_areaave.shtml

====
If you are learning NCL, I suggest that you ndownload the mini-language
manual.

http://www.ncl.ucar.edu/Document/index.shtml

This looks best in colr but black and white is fine. Also,
print on both sides of paper.

Good luck.
Kabumbwe Hansingo wrote:
> Hi,
> am trying to get an area average for an area 60degWest to 15degEast ,
> 30degSouth to 30degNorth. My data has longitudes from 0 to 359. I used
> sst(time|:,latitude|-30:30,longitude|-60:15) to select the data I need
> for the average. But the timeseries I get is different from that in,
> say GrADS. It appears that the problem is from crossing the 0deg
> Meridian. How do I get round this problem?
>
> thanks
> Hansingo
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Apr 21 2008 - 06:58:02 MDT