Hi Caroline,
You don't need to regrid the data at all in order to plot it in a
different projection. If you have the lat2d & lon2d arrays set,
NCL can plot the data as-is using any projection you like.
If you do something like this, I believe it will work:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
a = addfile(v_file,"r")
t = a->TMP_GDS3_HTGL_10(:,:)
t@lat2d= a->g3_lat_0
t@lon2d= a->g3_lon_1
wks =gsn_open_wks("ps","tempe")
res = True
res@gsnAddCyclic = False
res@mpProjection = "mercator"
plot=gsn_csm_contour_map(wks,tt,res)
end
Cheers,
--Seth
---- Seth McGinnis mcginnis@ucar.edu Associate Scientist ISSE / IMAGe / NCAR ---- >Date: Fri, 29 Jan 2010 10:22:06 +0100 >From: Caroline Serraud <cs@meteoconsult.fr> >Subject: Data from a Lambert native grid reprojected on a > Mercator grid. >To: ncl-talk@ucar.edu >Message-ID: <4B62A8BE.2090706@meteoconsult.fr> >Content-Type: text/plain; charset="iso-8859-1" > >Hi everyone, > >I am trying to plot data from a .nc file with a Lambert native grid, but >in a Mercator projection. >I can't find any working solution for several days: I try first with >rcrm2rgrid function, but when I launch my ncl script, nothing happens >and ncl seems to run for hours.. > > >Am I doing something wrong? >Any remark would be very helpful!! > >Thanks in advance for your time >Caroline > _______________________________________________ ncl-talk mailing list List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Fri Jan 29 12:48:31 2010
This archive was generated by hypermail 2.1.8 : Mon Feb 01 2010 - 08:05:34 MST