linear fit

Hi all:

 I want to plot a curve with its linear fit. Which function can be work?

  Thanks.

Xiangjun


>
>
> -----ԭʼÓʼþ-----
> ´Ó: "AdamPhillips"
> µ½: "MateusdaSilvaTeixeira"
> ³­ËÍ: ncl-talk_at_ucar.edu
> ÈÕÆÚ: 2006/10/05(ËÄ) 01:39am
> Ö÷Ìâ: Re: [ncl-talk] meta data
>
> Hi Mateus,
>
> I think your coordinate variables are getting lost when you divide by 10.
> Try this:
> geo2 = geo(0, 22,{-55:15},{270:330})
> geo2 = geo2/10
> plot = gsn_csm_contour_map_ce(xwks,geo2,res)
>
> Adam
>
> Mateus da Silva Teixeira wrote:
> > Dear NCL users,
> >
> > I'm trying to plot height fields instead of geopotential ones, so I do
> > the following
> >
> > plot = gsn_csm_contour_map_ce(xwks, *geo( 0, 22,
> > {-55:15},{270:330})/10*, res)
> >
> > however, ncl gives warning messages:
> >
> > (0) check_for_y_lat_coord: Warning: Data either does not contain a
> > valid latitude coordinate array or doesn't contain one at all.
> > (0) A valid latitude coordinate array should have a 'units'
> > attribute equal to one of the following values:
> > (0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees
> > north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN'
> > (0) check_for_lon_coord: Warning: Data either does not contain a
> > valid longitude coordinate array or doesn't contain one at all.
> > (0) A valid longitude coordinate array should have a 'units'
> > attribute equal to one of the following values:
> > (0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
> > 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE'
> >
> > and the fields are wrong. Is there a way to put all meta data
> > (coordinates, units, etc) into a new variable in a mathematical
> > operation ( X = Y/10 ) or I should give them manually?
> >
> > Thanks
> > Mateus
> >
> >
> > ------------------------------------------------------------------------
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk_at_ucar.edu
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> --
> --------------------------------------------------------------
> Adam Phillips asphilli_at_ucar.edu
> National Center for Atmospheric Research tel: (303) 497-1726
> ESSL/CGD/CAS fax: (303) 497-1333
> P.O. Box 3000
> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>


****************************************
* Xiangjun Tian *
* Institute of Atmospheric Physics *
* Chinese Academy of Sciences *
* P.O.Box 9804, *
* Beijing 100029, China *
* E-mail: tianxj_at_mail.iap.ac.cn *
* Tel: 86-10-62063912 (Office) *
****************************************

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Oct 04 2006 - 13:13:22 MDT