Re: linint2_Wrap works for half domain

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Mon Nov 11 2013 - 15:28:06 MST

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)


On Nov 8, 2013, at 3:44 PM, Ahmed Lasheen <> 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><><>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:

ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
Received on Mon Nov 11 15:28:15 2013

This archive was generated by hypermail 2.1.8 : Fri Nov 22 2013 - 09:36:32 MST