Re: linint2_Wrap works for half domain

Hi Ahmed,

You're almost there. You did the correct thing of calling lonFlip to convert your longitude values from a 0 to 360 range to a -180 to 180 range. But, you then passed "lon_era" to linint2_Wrap, which is still in the 0 to 360 range.

Since "lonFlip" also takes care of flipping your longitude coordinate array, you should use this, instead of "lon_era".

Change this line:

era = linint2_Wrap(lon_era,lat_era,avg_tem_era, True,lon_model,lat_model,0)

to be:

era = linint2_Wrap(avg_tem_era&longitude,avg_tem_era&latitude,avg_tem_era, True,lon_model,lat_model,0)

--Mary

On Nov 8, 2013, at 3:44 PM, Ahmed Lasheen <ahmed4kernel@gmail.com> wrote:

> Hello All
> I am trying to interpolate era interim to coarse domain model output. the lat, lon points are 73,144 for era data and 24,36 for model output. Half of interpolated era data are missing values.
> I have attached the script and data.
> Thanks in advance
> ---------------
> Ahmed L. Shaaban | Graduate Student
> Department of Atmospheric and Environmental science(DEAS).
> State University of New York, University (SUNY) at Albany
> 1400 Washington Avenue, Albany, NY 12222
>
>
> <era_40_control_model_bias_forum.ncl><Era_40_1958.nc><ControlRunATM551.1958-1958ij.nc>_______________________________________________
> 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 Mon Nov 11 15:28:15 2013