Re: hello! about f2fsh_Wrap with missing values

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Tue, 23 May 2006 10:59:32 -0600 (MDT)

>
> I have a HadISST1 data :
> lat = 180 ;
> lon = 360 ;
> lat:units = "degrees_north" ;
> lon:units = "degrees_east" ;
> lat = -89.5S to 89.5N ;
> lon = -179.5W to 179.5E ; all increase by 1
>degree.
>
>I want to change the SST(180*360) to T31(48*96), but
>with missing value inside, the function f2fsh_Wrap
>(SST, (/48,96/) ) will not work, do there have any way
>can do the job with missing value inside?
>
===============================
Use bilinear interpolation.

SST([ntim,]180,360) ==> sst([ntim,]48,96)

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

   nlat = 48
   mlon = 96
   
   lat = latGau (nlat, "lat", "latitude" , "degrees_north")
   lon = lonGlobeF (mlon, "lon", "longitude", "degrees_east")
   
   sst = linint2_Wrap(SST&lon, SST&lat, SST, True, lon, lat, 0)
   printVarSummary( sst )
   printMinMax(sst, True)
   

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue May 23 2006 - 10:59:32 MDT

This archive was generated by hypermail 2.2.0 : Tue May 23 2006 - 14:11:01 MDT