From wrfhelp@ucar.edu Mon Jan 07 16:04:50 2013
Received: from [128.117.88.60] (account wrfhelp@mail.mmm.ucar.edu)
  by mail.ucar.edu (CommuniGate Pro WEBUSER 5.4.8)
  with HTTP id 44051892 for milmilch@hotmail.com; Mon, 07 Jan 2013 16:04:50 -0700
From: "WRF Help" <wrfhelp@ucar.edu>
Subject: Re: NCL - plotgrids.ncl
To: "Camilla Cheng" <milmilch@hotmail.com>
X-Mailer: CommuniGate Pro WebUser v5.4.8
Date: Mon, 07 Jan 2013 16:04:50 -0700
Message-ID: <web-44051893@mail.ucar.edu>
In-Reply-To: <SNT107-W2390B449AB293B4F096E48C4340@phx.gbl>
References: <SNT107-W2390B449AB293B4F096E48C4340@phx.gbl>
MIME-Version: 1.0
Content-Type: text/plain;charset=utf-8; format="flowed"
Content-Transfer-Encoding: 8bit

Camilla,

Your plotgrids.ncl script will work if you remove the stand_lon line 
from your namelist file.  The mercator projection does not use 
stand_lon and having it in your namelist is confusing the program.  We 
will fix this in future releases of NCL to avoid this happening again. 
 Another note, you might want to change your truelat1 to 20 or even 10 
degrees because with the projection you are using your map factors 
will be quite large.  Also, your choice in domains is quite confusing. 
 Do you realize that they are all overlapping and have the same mother 
grid?

wrfhelp

On Sun, 23 Dec 2012 21:28:26 +0800
  Camilla Cheng <milmilch@hotmail.com> wrote:
> 
> Hi wrfhelp,
> 
> I use ncl to preview the simulation domain, the plotgrids.ncl script 
>can't plot all 4 domains set in my namelist.wps. Would mind giving me 
>some advice? Attached the namelist.wps_1 and plotgrids.ncl. I didn't 
>change anything in plotgrids.ncl download from the WRF/NCL website.
> 
> Thank you very much!
> 
> 
> 		 	   		  

From wrfhelp@ucar.edu Mon Jan 07 11:04:17 2013
Received: from [128.117.88.60] (account wrfhelp@mail.mmm.ucar.edu)
  by mail.ucar.edu (CommuniGate Pro WEBUSER 5.4.8)
  with HTTP id 44046306 for ee07a2jr@leeds.ac.uk; Mon, 07 Jan 2013 11:04:17 -0700
From: "WRF Help" <wrfhelp@ucar.edu>
Subject: Re: wrf_user_intrp3d problem
To: "Alexander Roberts" <ee07a2jr@leeds.ac.uk>
X-Mailer: CommuniGate Pro WebUser v5.4.8
Date: Mon, 07 Jan 2013 11:04:17 -0700
Message-ID: <web-44046307@mail.ucar.edu>
In-Reply-To: <20121217130824.151p2rnf3swsko4o@webmail.leeds.ac.uk>
References: <20121217130824.151p2rnf3swsko4o@webmail.leeds.ac.uk>
MIME-Version: 1.0
Content-Type: text/plain;charset=utf-8; format="flowed"
Content-Transfer-Encoding: 8bit

Alex,

The issue comes when you interpolate the data.  Because the function 
is interpolating, not extrapolating, there will be no data below your 
lowest half level.  We are considering putting in code that would 
extrapolate, but it would need to be used with extreme caution as it 
is obviously not real model data.  As of now, the lowest half level 
will have no data.

On Mon, 17 Dec 2012 13:08:24 +0000
  Alexander Roberts <ee07a2jr@leeds.ac.uk> wrote:
> Hi
> I am trying to plot some cross sections from a WRF-ARW case study 
>simulation using NCL. I have been using the wrf_user_intrp3d function 
>and have been able to produce some cross sections that I am almost 
>happy with. The only issue I have been having is that the 
>interpolation onto height levels seems to produce missing values 
>close to the surface. For the work that I'm doing I need to look into 
>nocturnal boundary layer formation, therefore this is quite a 
>problem. I have seen several other people mention this issue on the 
>NCL support pages but the advice given seems to be to contact WRF 
>help. Hopefully you will be able to help me.
> 
> Regards
> 
> Alex Roberts

From wrfhelp@ucar.edu Fri Dec 21 13:18:14 2012
Received: from [70.89.195.30] (account wrfhelp@mail.mmm.ucar.edu)
  by mail.ucar.edu (CommuniGate Pro WEBUSER 5.4.8)
  with HTTP id 43876422; Fri, 21 Dec 2012 13:18:14 -0700
From: "WRF Help" <wrfhelp@ucar.edu>
Subject: Re: [ncl-talk] Trying to get correct map over rotated lat-lon
 data
To: "David Jones" <jonesd647@yahoo.ca>, "David Brown" <dbrown@ucar.edu>
Cc: "wrfhelp@ucar.edu" <wrfhelp@ucar.edu>
X-Mailer: CommuniGate Pro WebUser v5.4.8
Date: Fri, 21 Dec 2012 13:18:14 -0700
Message-ID: <web-43876423@mail.ucar.edu>
In-Reply-To: <1354288361.41195.YahooMailNeo@web160102.mail.bf1.yahoo.com>
References: <1351091273.31536.YahooMailNeo@web163901.mail.gq1.yahoo.com>
  <1354128615.83337.YahooMailNeo@web160105.mail.bf1.yahoo.com>
 <5FC68FBB-7544-445A-8C2E-E1CF12CBCF9C@ucar.edu>
 <1354288361.41195.YahooMailNeo@web160102.mail.bf1.yahoo.com>
MIME-Version: 1.0
Content-Type: text/plain;charset=iso-8859-1; format="flowed"
Content-Transfer-Encoding: 8bit

David,

I would like to give you a heads up about plotting rotated lat-lon 
grids in NCL.  We are aware there is a problem and we are working 
putting a fix together.  We will post a note when it is fixed.

Thank you

On Fri, 30 Nov 2012 07:12:41 -0800 (PST)
  David Jones <jonesd647@yahoo.ca> wrote:
> Dear Dave,
> 
> thank you very much for your response. I really appreciate your 
>help. 
> 
> 
> Thanks for the file you sent. I made some slight changes (i.e. exact 
>corners, etc.) (attached) and was able to get it to work reasonably 
>well, but it's not really completely representing the borders in 
>places. For example, if I plot the variable LANDMASK instead of PSFC, 
>there are gaps near the edges of, for example, the Arabian Peninsula. 
>(Ignore the strange looking landmask in the polar regions, Greenland 
>etc., as that is an artifact of the reanalysis data i'm using). I 
>think that the method you've provided could well be the best that it 
>can work using the non-wrf routines in NCL. 
> 
> 
> I've uploaded the file I'm using on the ftp (it's called 
>wrfout_d01_2009-12-01_00:00:00) so you can see this yourself if you 
>like. I've reduced it to just one time step to make the file smaller. 
> 
> 
> When I try to use the 2D arrays XLAT and XLONG to define where the 
>map borders should be drawn I get an error message like 
>"warning:ScalarFieldSetValues: 2d coordinate array sfXArray has an 
>incorrect dimension size: defaulting sfXArray". Note that those 2D 
>arrays give the geographical lat and lon for each grid point in the 
>domain. My hope was that there would be some way to do this using the 
>WRF routines, as WRF output is not fully NetCDF compliant - which 
>could cause problems for some of the routines maybe. I thought that 
>maybe this could be done through attributes like LatLonOverlay, 
>sfXArray and sfYArray. 
> 
> 
>FYI, ref_lat and ref_lon define the geographical latitude/longitude 
>of a specific grid point in the domain. By default it is the middle 
>of the domain, however in my case I have shifted it northwards to the 
>grid point defined by (ref_x,ref_y) as listed in the comments at the 
>bottom of my original script (re-attached).
> 
> Many thanks for your assistance,
> 
> Dave. 
> 
> 
> 
> ________________________________
> From: David Brown <dbrown@ucar.edu>
> To: David Jones <jonesd647@yahoo.ca> 
> Cc: ncl-talk Talk <ncl-talk@ucar.edu>; wrfhelp@ucar.edu 
> Sent: Friday, November 30, 2012 5:56:38 AM
> Subject: Re: [ncl-talk] Trying to get correct map over rotated 
>lat-lon data
> 
> 
> Hi David,
> 
> 
> I think that rotated lat/lon projections in WRF is a confusing 
>subject, especially in the context of NCL, because I suspect that the 
>NCL WRF routines may not always map the attributes defined in the WRF 
>file correctly into the resource settings required for NCL (in the 
>rotated lat/lon case only). I will add the caveat that I do not fully 
>understand the WRF attributes defined in your file either. What I do 
>know is based on looking through this web page: 
>http://www.mmm.ucar.edu/wrf/users/docs/user_guide_V3/users_guide_chap3.htm#_Description_of_GEOGRID.TBL, 
>which may or may not be the most current description of the WRF 
>projections and corresponding namelist parameters.
> 
> Since I do not have your data file I cannot run your script for 
>myself, but I have attempted to approximate the map projection of 
>your data based on the WRF parameters you list in the script, the 
>documentation listed above, and the appearance of the surface 
>pressure field plot you attached. I am attaching my test script along 
>with the output. I think it is fairly close to the projection of your 
>data, but of course it would be better if I had the locations of the 
>lower left and upper right corners of your data. If you use 
>gsn_csm_contour_map instead ofwrf_map_overlays you should be able to 
>plot the correct map projection for your data. One thing the 
>documentation states is that CEN_LON and CEN_LAT do not apply for 
>rotated lat/lon projections, but I have to wonder why they are set at 
>all in that case. I am also not sure what the role of ret_lat and 
>ref_lon are. I do note that ref_lat - pole_lat is almost equal to 
>CEN_LAT and - STAND_LON is not too far
> from CEN_LON, but I don't know if that has any significance.
> 
> It would be helpful if you would upload your data file. It uses a 
>projection and a domain we have not often encountered. If you need 
>instructions on how to do that, 
>seehttp://www.ncl.ucar.edu/report_bug.shtml. Thanks,
> -dave
> 
> 
> 
> 
> 
> 
> 
> On Nov 28, 2012, at 11:50 AM, David Jones wrote:
> 
> 
>>Dear NCL talk
> help,
>>
>>
>>I would like to
> kindly ask for some advice regarding the plotting of rotated lat-lon
> output. I am currently trying to plot the output from a rotated
> lat-lon grid using NCL...
>>
>>
>>When I make the
> contour plot of my variable (psfc), all of the contoured data 
>appears
> where I expect and want it to be, the land-sea boundaries can even 
>be
> detected as the value of the variable changes. So everything seems
> fine in that respect. 
>>
>>
>>However, when I
> try to plot the national/geographical borders using a map overlay,
> the borders are all in the wrong places. Please see the attached
> image, in which it is clearly visible. 
>>
>>
>>The variable I
> am plotting (PSFC) is made up of (time, west_east, south_north), 
>with
> west_east and south_north corresponding to a point on the
> computational grid (i.e. they're defined by a lat and lon from the
> computational grid).
>>
>>
>>I also have 2D
> arrays (XLAT and XLONG) giving the geographical latitude and
> longitude for every point on the computational grid. I tried to use
> these to plot the geographical boundaries in the correct location,
> but sadly to no avail. 
>>
>>
>>I have looked
> through the NCL talk archives and on the NCL website, but
> unfortunately nothing I've tried has worked so far. 
>>
>>
>>I'll include my
> script below and also attach a copy of the output from an
> ncl_filedump on my data. You'll notice in the script below that lots
> of the map options are commented out; that's because several tries
> showed them to make no difference whatsoever to the plot I received,
> whether they were included or commented out (I did try plotting with
> them uncommented though, in several different combinations). The 
>only
> one of those options that did make a difference was setting
> mpres@mpCenterLatF   = a@CEN_LAT. But that only changed it from one
> incorrect plot to another incorrect one, albeit the new looks
> slightly better (attached). 
>>
>>
>>If you had any guidance as to where I might be going wrong I would be
> most appreciative. 
>>
>>
>>Many thanks for any advice,
>>Dave. 
>>
>>
>>----------------------SCRIPT-----------------------------
>>
>>
>>; Basic script for plotting model output
>>
>>
>>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/wrf/WRFUserARW.ncl"
>>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>
>>
>>begin
>>
>>
>>a = addfile("wrfout_d01_2009-12-01_00:00:00.nc","r")
>>
>>
>>wks = gsn_open_wks("x11","PSFC")
>>
>>
>>;---------------------------------------------------------------
>>; Set basic resources
>>res = True
>>res@MainTitle = "PSFC"			     ; Give a main title
>>pltres = True                                ; Plotting resources
>>mpres  = True                                ; Map resources
>>
>>
>>;---------------------------------------------------------------
>>; Get variables
>>;---1. WRF------------------
>>
>>
>>psfc	      = a->PSFC(0,:,:)
>>psfc_size     = dimsizes(psfc(:,:))
>>print("dimsizes(psfc) = " + psfc_size)
>>
>>
>>;---2. Lat/Lon arrays-------
>>
>>
>>lat2d  = a->XLAT(0,:,:)
>>printVarSummary(lat2d)
>>lon2d  = a->XLONG(0,:,:)
>>printVarSummary(lon2d)
>>
>>
>>;--------------------------------------------------------------
>>; Plotting options
>>opts		 = res                    ; Add basic resources
>>opts@FieldTitle = "PSFC"
>>opts@cnFillOn	 = True                           ; Shaded plot
>>opts@gsnSpreadColorEnd	= -8 ;
>>
>>
>>contour_psfc = wrf_contour(a,wks,psfc,opts) 
>>
>>
>>; Set some map options
>>
>>
>>;	mpres@tfDoNDCOverlay= True
>>
>>
>>;	pltres@sfXArray	     = lon2d
>>;	pltres@sfYArray	     = lat2d
>>;	pltres@LatLonOverlay = True
>>
>>
>>;	mpres@mpCenterLonF   = a@CEN_LON ;= -16.84654
>>mpres@mpCenterLatF   = a@CEN_LAT ;= 53.4331
>>
>>
>>;	mpres@mpProjection   = "CylindricalEquidistant"
>>;	mpres@mpLimitMode    = "Corners"
>>;	mpres@mpLeftCornerLatF = lat2d(0,0) ; bottom left
>>;	mpres@mpLeftCornerLonF = lon2d(0,0) ; bottom left
>>;	mpres@mpRightCornerLatF= lat2d(268,568) ; top right
>>;	mpres@mpRightCornerLonF= lon2d(268,568)	; top right
>>
>>
>>mpres@mpGeophysicalLineColor = "Black"
>>mpres@mpGridLineColor = "Black"
>>mpres@mpLimbLineColor = "Black"
>>mpres@mpNationalLineColor = "Black"
>>mpres@mpUSStateLineColor = "Black"
>>mpres@mpPerimLineColor = "Black"
>>mpres@mpGeophysicalLineThicknessF = 1.0
>>mpres@mpGridLineThicknessF = 1.0
>>mpres@mpLimbLineThicknessF = 1.0
>>mpres@mpNationalLineThicknessF = 1.0
>>mpres@mpUSStateLineThicknessF = 1.0
>>
>>
>>;  Namelist map projection information:
>>;  map_proj : "lat-lon" (cylindrical equidistant)
>>;  ref_lat  : 88.5 (geographical lat of grid point (ref_x,ref_y)
>>;  ref_lon  : -5.0 (geographical lon of grid point (ref_x,ref_y)
>>;  ref_x    : 285 (total x grid points = 570)
>>;  ref_y    : 265 (total y grid points = 270)
>>;  pole_lat : 35.0 (lat of geog N hem pole on comp grid)
>>;  pole_lon : 180.0(lon of geog N hem pole on comp grid)
>>;  stand_lon: 17.5 (rotation about geographic poles)
>>
>>
>>; MAKE PLOTS 
>>
>>
>>plot = wrf_map_overlays(a,wks,(/contour_psfc/),pltres,mpres)
>>
>>
>>delete(res)
>>
>>end
>>
>><ncl_filedump.txt><psfc_rotlatlon.ncl><PSFC.000001.png>_______________________________________________
>>ncl-talk mailing list
>>List instructions, subscriber options, unsubscribe:
>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>

From wrfhelp@ucar.edu Fri Dec 14 15:10:25 2012
Received: from [128.117.88.60] (account wrfhelp@mail.mmm.ucar.edu)
  by mail.ucar.edu (CommuniGate Pro WEBUSER 5.4.8)
  with HTTP id 43779956 for Christopher.Steele@uea.ac.uk; Fri, 14 Dec 2012 15:10:25 -0700
From: "WRF Help" <wrfhelp@ucar.edu>
Subject: Re: creating cross-sections in ncl
To: "Christopher Steele (ENV)" <Christopher.Steele@uea.ac.uk>
X-Mailer: CommuniGate Pro WebUser v5.4.8
Date: Fri, 14 Dec 2012 15:10:25 -0700
Message-ID: <web-43779957@mail.ucar.edu>
In-Reply-To: <7C9D2D609F8E484886E1D28BB3AFF2BA028940@ueastfexch01.UEA.AC.UK>
References: <7C9D2D609F8E484886E1D28BB3AFF2BA028940@ueastfexch01.UEA.AC.UK>
MIME-Version: 1.0
Content-Type: text/plain;charset=utf-8; format="flowed"
Content-Transfer-Encoding: 8bit

There is no degree of rotation.  Winds are still positive to the east 
and north.  However, please remember that this is map relative and not 
earth relative.

On Fri, 23 Nov 2012 13:36:44 +0000
  "Christopher Steele (ENV)" <Christopher.Steele@uea.ac.uk> wrote:
> Hi,
> 
> I'm using ncl with WRF data to calculate a vertical cross-section 
>between two points based on the cross-section script on the WRF users 
>website. The plot requires the calculation of wind speed with the 
>u-wind component being defined so that it is positive in the offshore 
>direction, where the coastline is at an angle perpendicular to the 
>cross-section. Is there any degree of rotation of the u-wind 
>component in the creation of the cross-section or is the u-wind 
>component always defined so that it is positive to towards the East 
>even after the use of the wrf_user_intrp3d function in the script? 
> 
> cheers
> 
> Chris
> 
> 
> Christopher Steele,
> School of Environmental Sciences
> University of East Anglia
> Norwich

From wrfhelp@ucar.edu Fri Dec 14 15:06:10 2012
Received: from [128.117.88.60] (account wrfhelp@mail.mmm.ucar.edu)
  by mail.ucar.edu (CommuniGate Pro WEBUSER 5.4.8)
  with HTTP id 43779877 for zhangbin06@mails.tsinghua.edu.cn; Fri, 14 Dec 2012 15:06:10 -0700
From: "WRF Help" <wrfhelp@ucar.edu>
Subject: Re: ncl wrf_vector
To: "Bin Zhang" <zhangbin06@mails.tsinghua.edu.cn>
X-Mailer: CommuniGate Pro WebUser v5.4.8
Date: Fri, 14 Dec 2012 15:06:10 -0700
Message-ID: <web-43779878@mail.ucar.edu>
In-Reply-To: <C18FE7628AA64A53BE42167674B089CA@zbPC>
References: <C18FE7628AA64A53BE42167674B089CA@zbPC>
X-Priority: 3
MIME-Version: 1.0
Content-Type: text/plain;charset=utf-8; format="flowed"
Content-Transfer-Encoding: 8bit

Could you tell me specifically what is wrong with the script.  Is 
there a plot that is wrong?  An error message?

On Thu, 13 Dec 2012 12:21:03 +0800
  "Bin Zhang" <zhangbin06@mails.tsinghua.edu.cn> wrote:
> Hi,
> 
> I am using NCL to plot wind vector. Something is wrong with my NCL 
>script. Could you please help me with this?  Thanks! 
> 
> Below 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/wrf/WRFUserARW.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
> 
> begin
>  direc="/net/9/zb/wrfout/haze_200608/ade_lin_9km/"
>  fn1=addfile(direc+"wrfout_d02_2006-08-01_00:00:00_normal.nc","r") 
>;filename
>  u10 = wrf_user_getvar(fn1,"U10",0)   ; U component of the wind at 
>10m
>  v10 = wrf_user_getvar(fn1,"V10",0)   ; V component of the wind at 
>10m
>  wks=gsn_open_wks("x11","wind_10m")
> 
>  res = True
>  res@gsnDraw  = True    ; Forces the plot to be drawn
>  res@gsnFrame = True    ; Frame advance
>  vector = wrf_vector(fn1,wks,u10,v10,res)
> 
> end
> 
> 
> Bin Zhang, Ph.D candidate
> Center for Earth System Science
> Tsinghua University, Beijing 100084, China
> Email: zhangbin06@mails.tsinghua.edu.cn

From wrfhelp@ucar.edu Tue Nov 27 14:08:43 2012
Received: from [128.117.88.60] (account wrfhelp@mail.mmm.ucar.edu)
  by mail.ucar.edu (CommuniGate Pro WEBUSER 5.4.8)
  with HTTP id 43533370 for e.rocheta@unsw.edu.au; Tue, 27 Nov 2012 14:08:43 -0700
From: "WRF Help" <wrfhelp@ucar.edu>
Subject: Re: wrf_pvo code
To: "Eytan Rocheta" <e.rocheta@unsw.edu.au>
X-Mailer: CommuniGate Pro WebUser v5.4.8
Date: Tue, 27 Nov 2012 14:08:43 -0700
Message-ID: <web-43533371@mail.ucar.edu>
In-Reply-To: <B5AFC204BEE04245AD848924DA51C89B454448A8@INFPWXM008.ad.unsw.edu.au>
References: <B5AFC204BEE04245AD848924DA51C89B454448A8@INFPWXM008.ad.unsw.edu.au>
MIME-Version: 1.0
Content-Type: text/plain;charset=utf-8; format="flowed"
Content-Transfer-Encoding: 8bit

Eytan,

This code is available in the NCL source code.  It is located in 
ncl_ncarg-*version*/ni/src/lib/nfpfort.

Abby Jaye

On Mon, 19 Nov 2012 03:58:23 +0000
  Eytan Rocheta <e.rocheta@unsw.edu.au> wrote:
> Hi,
> 
> I am applying the wrf_pvo code to wrf output and I am interested in 
>the exact formulation that is being used in this code.
> 
> Where/how can I assess the fundamental code which is used to compute 
>wrf_pvo?
> 
> Thanks
> Eytan
> 
> Eytan Rocheta
> PhD Candidate
> 
> School of Civil and Environmental Engineering
> University of New South Wales | UNSW Sydney NSW 2052
> T: 0415 935 713 | E: 
>e.rocheta@unsw.edu.au<mailto:e.rocheta@unsw.edu.au> | Web: 
>www.water.unsw.edu.au<http://www.water.unsw.edu.au>
> 
> This message is intended for the addressee named and may contain 
>confidential information.
> If you are not the intended recipient, please delete it and notify 
>the sender.
> Views expressed in this message are those of the individual sender 
>and are not necessarily the views of UNSW.
> CRICOS Provider Code 00098G
> 

From wrfhelp@ucar.edu Tue Nov 27 13:43:02 2012
Received: from [128.117.88.60] (account wrfhelp@mail.mmm.ucar.edu)
  by mail.ucar.edu (CommuniGate Pro WEBUSER 5.4.8)
  with HTTP id 43532804 for eowynb@atmos.washington.edu; Tue, 27 Nov 2012 13:43:02 -0700
From: "WRF Help" <wrfhelp@ucar.edu>
Subject: Re: wrf_map_overlays() longitude scale
To: "Eowyn Baughman" <eowynb@atmos.washington.edu>
X-Mailer: CommuniGate Pro WebUser v5.4.8
Date: Tue, 27 Nov 2012 13:43:02 -0700
Message-ID: <web-43532805@mail.ucar.edu>
In-Reply-To: <CAAq4mkL9BHiKyE4MpREvBMnP83j4yr=b2ACms+a8kw9cEra0zQ@mail.gmail.com>
References: <CAAq4mkL0v4V=3bTnDV7kE3bFhccCKT0UbmZ0OOGE+6fVp7-wMw@mail.gmail.com>
 <CAAq4mkLZNjmdTeedVSO+h2mZh=QDq30oZ8Y0S5RXREcDaX-ZDw@mail.gmail.com>
 <CAAq4mkJT-okpDb=YEFgGmCz6Yo_i5LRbhpCLZgyjeSeFXGAkVg@mail.gmail.com>
 <web-43380833@mail.ucar.edu>
 <CAAq4mkL9BHiKyE4MpREvBMnP83j4yr=b2ACms+a8kw9cEra0zQ@mail.gmail.com>
MIME-Version: 1.0
Content-Type: text/plain;charset=utf-8; format="flowed"
Content-Transfer-Encoding: 8bit

Eowyn,

The problem appears to be in the 
$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl file.  I am 
currently trying to work on a more permanent fix, but if you comment 
out a few lines of this file your file should plot without a problem. 
 Just to warn you, the line numbers I am giving you might not 
perfectly match up with your file depending on what version you have 
so I am including the text that needs to be commented out.

Comment out lines 3175-3180.  Should look like this:

;        if( isatt(in_file,"POLE_LAT") ) then
;          opts@mpCenterRotF = get_res_value_keep(opts, 
"mpCenterRotF", 90.0 - in_file@POLE_LAT)
;          delete(opts@mpCenterLonF)
;          calcen = -190.
;          opts@mpCenterLonF  = get_res_value_keep(opts, 
"mpCenterLonF", calcen )
;        end if

Comment out lines 3595-3601.  Should look like this:

;        if( isatt(opts,"pole_lat") ) then
;          delete(opts@mpCenterLonF)
;          opts@mpCenterLonF  = get_res_value_keep(opts, 
"mpCenterLonF", - 190. )
;          opts@mpCenterRotF = get_res_value_keep(opts, 
"mpCenterRotF", 90.0 - opts@pole_lat)
;        else
;          opts@mpCenterRotF = get_res_value_keep(opts, 
"mpCenterRotF", 0.0)
;        end if

Comment out lines 3997-4001.  Should look like this:

;        if( isatt(in_file,"POLE_LAT") ) then
;          map_args@mpCenterRotF = get_res_value_keep(map_args, 
"mpCenterRotF", 90.0 - in_file@POLE_LAT)
;          delete(map_args@mpCenterLonF)
;          calcen = -190.
;          map_args@mpCenterLonF  = get_res_value_keep(map_args, 
"mpCenterLonF", calcen )
;        end if

  Please let me know if you have any problems implementing this fix 
and I will let you know when I have a more permanent fix for you.

Abby Jaye


On Sat, 17 Nov 2012 09:32:32 -0800
  Eowyn Baughman <eowynb@atmos.washington.edu> wrote:
> Hi Abby,
> 
> Certainly, it is below. Thank you.
> 
> Regard,
> 
> &share
> wrf_core = 'ARW',
> max_dom = 1,
> start_date = '2008-08-16_00:00:00','2006-08-16_12:00:00',
> end_date   = '2008-09-15_18:00:00','2006-08-16_12:00:00',
> interval_seconds = 21600,
> io_form_geogrid = 2,
> /
> 
> &geogrid
> parent_id         =   1,   1,
> parent_grid_ratio =   1,   3,
> i_parent_start    =   1,  31,
> j_parent_start    =   1,  17,
> e_we              =  61, 112,
> e_sn              =  51,  97,
> geog_data_res     = '10m','2m',
> dx = 1.,
> dy = 1.,
> map_proj = 'lat-lon',
> ref_lat   =  -10.0,
> ref_lon   = 10.0,
> stand_lon = 180.0,
> geog_data_path = '/home/disk/eos10/eowynb/geog'
> /
> 
> &ungrib
> out_format = 'WPS',
> prefix = 'FILE',
> /
> 
> &metgrid
> fg_name = 'FILE'
> io_form_metgrid = 2,
> /
> 
> On Fri, Nov 16, 2012 at 10:16 AM, WRF Help <wrfhelp@ucar.edu> wrote:
> 
>> Eowyn,
>>
>> Would you be able to send me your WPS namelist so I can see how you
>> created your domain?
>>
>> Thanks!
>> Abby Jaye
>>
>>
>> On Tue, 30 Oct 2012 09:32:22 -0700
>>  Eowyn Baughman <eowynb@atmos.washington.edu> wrote:
>>
>>> Hello,
>>>
>>> I first sent the query (below) in early September. I have not yet 
>>>received
>>> a response. Please let me know if there is a better place to direct 
>>>WRF /
>>> NCL questions.
>>>
>>> Thank you.
>>>
>>> Regards,
>>> Eowyn Baughman
>>>
>>> On Wed, Sep 12, 2012 at 6:00 PM, Eowyn Baughman <
>>> eowynb@atmos.washington.edu
>>>
>>>> wrote:
>>>>
>>>
>>>  Hi again,
>>>>
>>>> I tried setting just about every relevant map resource specifying 
>>>>the
>>>> domain, but none were valid in contourplot_contour so it was
>>>> unsuccessful.
>>>> That script is also available at
>>>> http://www.atmos.washington.**edu/~eowynb/wrfnclfiles/<http://www.atmos.washington.edu/~eowynb/wrfnclfiles/>
>>>>
>>>> Thanks,
>>>> Eowyn
>>>>
>>>>
>>>> On Wed, Sep 12, 2012 at 4:33 PM, Eowyn Baughman <
>>>> eowynb@atmos.washington.edu> wrote:
>>>>
>>>>  Hi WRF-ers,
>>>>>
>>>>> I'm experienced with NCL but am using it on WRF (and regional) data 
>>>>>for
>>>>> the first time. The funny coordinates are a pain. When I use
>>>>> wrf_map_overlays(), I get the correct latitude range but the wrong
>>>>> longitude range. I've tried setting mpCenterLon but it doesn't 
>>>>>appear
>>>>> to be
>>>>> a valid resource. My domain is from -40:20 longitude, from around
>>>>> Ascension
>>>>> Island in the SE Atlantic to around Zanzibar. But wrf_map_overlay() 
>>>>>is
>>>>> doing the other circle, from Zanzibar to Ascension Island.
>>>>>
>>>>> The files are here: http://www.atmos.washington.**
>>>>> edu/~eowynb/wrfnclfiles/<http://www.atmos.washington.edu/~eowynb/wrfnclfiles/>
>>>>>
>>>>> If it was a global domain, I'd use lonflip or lonpivot, but it is 
>>>>>not.
>>>>> There's probably a simple solution; if you have one, please let me 
>>>>>know.
>>>>>
>>>>> I sent this message to ncl-talk as well, but was told this was the 
>>>>>more
>>>>> appropriate group.
>>>>>
>>>>> Thanks.
>>>>>
>>>>> Regards,
>>>>> Eowyn Baughman
>>>>>
>>>>>
>>>>
>>>>
>>

From wrfhelp@ucar.edu Fri Nov 16 11:16:15 2012
Received: from [128.117.88.60] (account wrfhelp@mail.mmm.ucar.edu)
  by mail.ucar.edu (CommuniGate Pro WEBUSER 5.4.8)
  with HTTP id 43380832 for eowynb@atmos.washington.edu; Fri, 16 Nov 2012 11:16:15 -0700
From: "WRF Help" <wrfhelp@ucar.edu>
Subject: Re: wrf_map_overlays() longitude scale
To: "Eowyn Baughman" <eowynb@atmos.washington.edu>
X-Mailer: CommuniGate Pro WebUser v5.4.8
Date: Fri, 16 Nov 2012 11:16:15 -0700
Message-ID: <web-43380833@mail.ucar.edu>
In-Reply-To: <CAAq4mkJT-okpDb=YEFgGmCz6Yo_i5LRbhpCLZgyjeSeFXGAkVg@mail.gmail.com>
References: <CAAq4mkL0v4V=3bTnDV7kE3bFhccCKT0UbmZ0OOGE+6fVp7-wMw@mail.gmail.com>
 <CAAq4mkLZNjmdTeedVSO+h2mZh=QDq30oZ8Y0S5RXREcDaX-ZDw@mail.gmail.com>
 <CAAq4mkJT-okpDb=YEFgGmCz6Yo_i5LRbhpCLZgyjeSeFXGAkVg@mail.gmail.com>
MIME-Version: 1.0
Content-Type: text/plain;charset=utf-8; format="flowed"
Content-Transfer-Encoding: 8bit

Eowyn,

Would you be able to send me your WPS namelist so I can see how you 
created your domain?

Thanks!
Abby Jaye

On Tue, 30 Oct 2012 09:32:22 -0700
  Eowyn Baughman <eowynb@atmos.washington.edu> wrote:
> Hello,
> 
> I first sent the query (below) in early September. I have not yet 
>received
> a response. Please let me know if there is a better place to direct 
>WRF /
> NCL questions.
> 
> Thank you.
> 
> Regards,
> Eowyn Baughman
> 
> On Wed, Sep 12, 2012 at 6:00 PM, Eowyn Baughman 
><eowynb@atmos.washington.edu
>> wrote:
> 
>> Hi again,
>>
>> I tried setting just about every relevant map resource specifying 
>>the
>> domain, but none were valid in contourplot_contour so it was 
>>unsuccessful.
>> That script is also available at
>> http://www.atmos.washington.edu/~eowynb/wrfnclfiles/
>>
>> Thanks,
>> Eowyn
>>
>>
>> On Wed, Sep 12, 2012 at 4:33 PM, Eowyn Baughman <
>> eowynb@atmos.washington.edu> wrote:
>>
>>> Hi WRF-ers,
>>>
>>> I'm experienced with NCL but am using it on WRF (and regional) data 
>>>for
>>> the first time. The funny coordinates are a pain. When I use
>>> wrf_map_overlays(), I get the correct latitude range but the wrong
>>> longitude range. I've tried setting mpCenterLon but it doesn't 
>>>appear to be
>>> a valid resource. My domain is from -40:20 longitude, from around 
>>>Ascension
>>> Island in the SE Atlantic to around Zanzibar. But wrf_map_overlay() 
>>>is
>>> doing the other circle, from Zanzibar to Ascension Island.
>>>
>>> The files are here: 
>>>http://www.atmos.washington.edu/~eowynb/wrfnclfiles/
>>>
>>> If it was a global domain, I'd use lonflip or lonpivot, but it is 
>>>not.
>>> There's probably a simple solution; if you have one, please let me 
>>>know.
>>>
>>> I sent this message to ncl-talk as well, but was told this was the 
>>>more
>>> appropriate group.
>>>
>>> Thanks.
>>>
>>> Regards,
>>> Eowyn Baughman
>>>
>>
>>

