Re: Plotting RSM output

Hi Mary

I really appreciate your help.

Mary Haley escribió:
> Angel,
>
> After I got your data offline and took at look at it, I realized that
> your data already had coordinate arrays, so there's no need to create
> them again.
> In order to preserve the metadata, I read your variable with two lines
> instead of one:
>
> T = f->TMP_GDS1_SFC_10 ; Get variable and metadata
> T = T-273.15 ; Convert to different units
>
> I've attached a new version of your script.
>
> However, the problem still persists. My first inclination is to say
> there is something wrong with the lat/lon values on your file, as I
> tried other things to get the data to line up, with no luck.
>
> I'm going to ask one of our developers about this.
>
> --Mary
>
>
> ------------------------------------------------------------------------
>
>
> On Jul 28, 2010, at 11:02 AM, Ángel G. Muñoz wrote:
>
>> Thanks Mary. I understand.
>>
>> Nonetheless, I can´t solve the problem yet....
>>
>> Mary Haley escribió:
>>> Hi Angel,
>>>
>>> Since you have lat/lon coordinate arrays, you don't want to set this:
>>>
>>>
>>>> res@tfDoNDCOverlay = True ; do not transform data
>>>>
>>>
>>> This causes NCL to assume your data has already been transformed to some map projection, and hence it ignores your coordinate arrays.
>>>
>>> Also, since you've gone through the trouble of setting up coordinate arrays, you don't need these two lines either:
>>>
>>>
>>>> res@sfXArray=lon
>>>> res@sfYArray=lat
>>>>
>>>
>>> These lines are redundant.
>>>
>>> --Mary
>>>
>>>
>>> On Jul 28, 2010, at 7:43 AM, Ángel G. Muñoz wrote:
>>>
>>>
>>>> Dear all,
>>>>
>>>> Has anyone used NCL for RSM plots?
>>>>
>>>> I have tryied different options, but I always obtain the field in a different "projection" than the "map" borders: there is no coincidence between the physical field and the continent drawn by gsn_csm_contour_map.
>>>>
>>>> This is my script:
>>>>
>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>> ;*****************************************************
>>>> begin
>>>> f = addfile("r_pgb.ft00.nc","r") ; I have transformed the file to nc
>>>> T = f->TMP_GDS1_SFC_10-273.15 lat = f->g1_lat_0
>>>> lon = f->g1_lon_1
>>>> T@long_name = "Temperatura" T!0="lat"
>>>> T!1="lon"
>>>> T&lat=lat
>>>> T&lon=lon
>>>> T&lat@units="degrees_north"
>>>> T&lon@units="degrees_east"
>>>> ;***************************************
>>>> ; create individual plots
>>>> ;***************************************
>>>> wks = gsn_open_wks("ps","RSM") ; open ps file
>>>> gsn_define_colormap(wks,"gui_default") ; choose colormap
>>>>
>>>> res = True ; plot mods desired
>>>> res@gsnMaximize = True ; uncomment to maximize size
>>>> res@cnLineThicknessF = 2.0 ; make twice as thick as default
>>>> res@cnLineLabelBackgroundColor = "background" ; no lines through labels
>>>> res@cnFillOn = True ; color contours
>>>> res@cnLinesOn = False ; default is True
>>>> res@cnLineLabelsOn = False ; default is True
>>>> res@gsnSpreadColors = True ; use full range of colormap
>>>> res@gsnAddCyclic = False
>>>> res@sfXArray=lon
>>>> res@sfYArray=lat
>>>> res@mpProjection = "mercator" ; projection
>>>> res@mpLimitMode = "Corners" ; method to zoom
>>>> res@mpLeftCornerLatF = min(lat)
>>>> res@mpLeftCornerLonF = min(lon)
>>>> res@mpRightCornerLatF = max(lat)
>>>> res@mpRightCornerLonF = max(lon)
>>>> res@mpOutlineDrawOrder = "PostDraw" ; draw continental outline last
>>>> res@mpFillDrawOrder = "PreDraw"
>>>> res@mpOutlineBoundarySets = "National" ; state boundaries
>>>> plot = gsn_csm_contour_map(wks,T(:,:),res)
>>>> end
>>>>
>>>>
>>>> The output can be seen here:
>>>> http://cmc.org.ve/descargas/OA/RSM1.png
>>>>
>>>> Thanks
>>>> <agmunoz.vcf>_______________________________________________
>>>> ncl-talk mailing list
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>
>>>
>>>
>>>
>> <agmunoz.vcf>
>