Transform the lambert conformal to cylindrical equitidistant

From: gaby colorado <gaby.colorado_at_nyahnyahspammersnyahnyah>
Date: Tue, 21 Nov 2006 22:25:55 -0600

hello, i want to transform the lambert conformal to cylindrical
equitidistant, i saw some examples, but my cuestion is how do it considering
the time, one part the example is:

> NLAT = 46
> MLON = 180

> dlat = 90./(NLAT-1) ; output lat
> lat = ispan ( 0,NLAT-1,1 )*dlat
> lat!0 = "lat"
> lat_at_units = "degrees_north"
> lat&lat = lat

> dlon = 360./MLON ; output lon
> lon = ispan ( 0,MLON-1,1 )*dlon - 180.
> lon!0 = "lon"
> lon_at_units = "degrees_east"
> lon&lon = lon

> kl = 10 ; arbitrary level, in this case 400hPa
> X2D = triple2grid(ndtooned(lon2d),ndtooned(lat2d),ndtooned(x(kl,:,:)), lon,lat, False)

and i try the follow:
   ntime = dimsizes(tim)
   NLAT = 90
   MLON = 180

   dlat = 90./(NLAT-1) ; output lat
   lat = ispan ( 0,NLAT-1,1 )*dlat
   lat!0 = "lat"
   lat_at_units = "degrees_north"
   lat&lat = lat

   dlon = 360./MLON ; output lon
   lon = ispan ( 0,MLON-1,1 )*dlon - 180.
   lon!0 = "lon"
   lon_at_units = "degrees_east"
   lon&lon = lon

  X2D = new ((/ntime,MLON,NLAT/), float,-999)
  do i = 0,ntime-1,1
      X2D(i,:,:) =
triple2grid(ndtooned(lon2d),ndtooned(lat2d),ndtooned(x(i,:,:)), lon,lat,
False)
      X2D!0 = "lat" ; triple2grid regrids based on a
nearest neighbor routine.
      X2D&lat = lat
      X2D!1 = "lon"
      X2D&lon = lon
      X2D_at_units = x_at_units
   end do

I know that is wrong but i don't have other idea, the error tha it send me
es:

fatal:Dimension size mismatch on subscript #1, left-hand and right-hand side
dimensions do not match
fatal:Execute: Error occurred at or near line 49 in file conversion.ncl

thanks for your help

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Nov 21 2006 - 21:25:55 MST

This archive was generated by hypermail 2.2.0 : Tue Nov 28 2006 - 15:20:59 MST