Re: Defined lat and long

As Debashish has noted, the WRF_map_c [same as WRF_mapres_c]
sets map resources. It sets them to plot the entire region.
They are overwriting your settings.

My advice would be to ad your resources after WRF_map_c

Debasish Pai Mazumder wrote:
> Dear Moti,
>
> Your script looks fine except line 49 where you called the function
> WRF_map_c (f, res, 0). I guess the function WRF_map_c (WRF_mapres_c ??) is
> caused the problem. You may try to plot without using this function
>
> Please check the webpage
>
>
> http://www.ncl.ucar.edu/Document/Functions/WRF_contributed/wrf_mapres_c.shtm
> l
>
> Thanks
>
> Debasish
>
>
>
> -----Original Message-----
> From: ncl-talk-bounces_at_ucar.edu [mailto:ncl-talk-bounces_at_ucar.edu] On Behalf
> Of Mittal, Moti
> Sent: Thursday, April 23, 2009 1:39 PM
> To: ncl-talk_at_ucar.edu
> Subject: Defined lat and long
>
> In the following script, I want to plot O3 between 12.5 to 13.5 lat and 79.5
> to 80.5 longitude.
> I included res_at_mpMinLatF = 12.5 etc. (line 41-44) but still
> the plot gives me 11.5 to 14.5 latitude and 78 to 82 longitude which is the
> output from wrfout file.
> How should I change my script?
>
> Thank you.
>
> Moti Mittal
>
>
> Script;
> 1 ;*************************************************
> 2 ; WRF: color over map with lat/lon labels
> 3 ;************************************************
> 4 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> 5 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> 6 load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
> 7 begin
> 8 ;************************************************
> 9 ; open file and read in data
> 10 ;************************************************
> 11 f = addfile
> ("/home/mmittal/CHENNAI/CHEM/WRFOUT/MAY/wrfout_d01_2007-05-03_06:
> 00:00.nc", "r")
> 12 ;************************************************
> 13 ; Read character variable Times; Convert to string for plots
> 14 ; Read vertical coordinate for plot labels
> 15 ;************************************************
> 16 times = chartostring(f->Times) ; built-in function
> 17 znu = f->ZNU(0,:) ; (Time, bottom_top)
> 18 lat2d = f->XLAT(0,:,:)
> 19 lon2d = f->XLONG(0,:,:)
> 20 ;************************************************
> 21 ; Read ozone at all times and levels
> 22 ;************************************************
> 23 x = f->o3 ; (Time, bottom_top,
> south_north, w est_east)
> 24
> 25 ;************************************************
> 26 ; create plots
> 27 ;************************************************
> 28 wks = gsn_open_wks("x11" ,"WRF_po") ; ps,pdf,x11,ncgm,eps
> 29 gsn_define_colormap(wks,"BlAqGrYeOrReVi200") ; select color map
> 30
> 31 res = True ; plot mods desired
> 32 res_at_gsnMaximize = True ; uncomment to maximize
> size
> 33 res_at_gsnSpreadColors = True ; use full range of
> colormap
> 34 res_at_cnFillOn = True ; color plot desired
> 35 res_at_cnLinesOn = True ; turn off contour lines
> 36 res_at_cnLineLabelsOn = True ; turn off contour labels
> 37 res_at_cnLevelSelectionMode = "ManualLevels"
> 38 res_at_cnMinLevelValF = .025
> 39 res_at_cnMaxLevelValF = .035
> 40 res_at_cnLevelSpacingF = .002
> 41 res_at_mpMinLatF = 12.5
> 42 res_at_mpMaxLatF = 13.5
> 43 res_at_mpMinLonF = 79.5
> 44 res_at_mpMaxLonF = 80.5
> 45
> 46 ;************************************************
> 47 ; Use WRF_contributed procedure to set map resources
> 48 ;************************************************
> 49 WRF_map_c (f, res, 0)
> 50 ;************************************************
> 51 ; set True for native projection (faster)
> 52 ;************************************************
> 53 res_at_tfDoNDCOverlay = False
> 54 x_at_lat2d = lat2d
> 55 x_at_lon2d = lon2d
> 56
> 57
> 58 ;************************************************
> 59 ; Turn on lat / lon labeling
> 60 ;************************************************
> 61 res_at_pmTickMarkDisplayMode = "Always" ; turn on tickmarks
> 62
> 63 ;************************************************
> 64 ; Loop over all times and levels ( uncomment )
> 65 ; Demo: one arbitrarily closen time and level
> 66 ;************************************************
> 67 dimx = dimsizes(x) ; dimensions of x
> 68 ntim = dimx(0) ; number of time steps
> 69 klev = dimx(1) ; number of "bottom_top"
> levels
> 70
> 71 ; nt = ntim/2 ; arbitrary time
> 72 ; kl = 1 ; " level
> 73 do nt=0,ntim-1 ; uncomment for loop
> 74 do kl=0,klev-1
> 75 res_at_tiMainString = times(nt)
> 76 res_at_gsnLeftString = x_at_description+" z="+znu(kl)
> 77 plot = gsn_csm_contour_map(wks,x(nt,kl,:,:),res)
> 78 ; now change the size of the tickmark labels
> 79 res_at_tmXBLabelFontHeightF = 0.02 ; resize tick labels
> 80 res_at_tmYLLabelFontHeightF = 0.02
> 81 plot = gsn_csm_contour_map(wks,x(nt,kl,:,:),res)
> 82 end do
> 83 end do
> 84 end
>
>
> Professor, Department of Environment and Occupational Health
> College of Public Health
> University of South Florida
> 13201 Bruce B. Downs Blvd., MDC56
> Tampa, Fl 33612-3805
>
> Tel: (813) 974-9571
> Fax: (813) 974-4986
> Email: mmittal_at_health.usf.edu
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk