Re: wrf_user_latlon_to_ij usage

I need to also warn that wrf_user_latlon_to_ij uses a different algorithm than wrf_ll_to_ij, and that it may be incorrect for some cases. I would be careful about using this old version.

I've CC'ed wrfhelp in case they have some tips on how to use wrf_user_latlon_to_ij to get correct results.

--Mary

On Sep 2, 2011, at 1:43 PM, Paula Doubrawa Moreira wrote:

> Chris-
>
> I think it's a bad idea to replace anything that is on your $NCARG_ROOT directory.
>
> Download it, and save it as a different name wherever you want. For example, call it "WRFUserARW_v6.ncl".
>
> Then, when loading the functions in your script, do not load the one located in $NCARG_ROOT/blablabla... but only the new one, with the full path.
>
> On Fri, Sep 2, 2011 at 11:33 AM, Chris Klich <caklich@gmail.com> wrote:
> I have tried this and it still seems to not work properly, and I'm not sure why. I still get the statement that the function is undefined or cannot be found. Did you just replace the old WRFUserARW with a new one or did you make a duplicate?
>
>
> On Fri, Sep 2, 2011 at 3:29 PM, Paula Doubrawa Moreira <pmoreira@alaska.edu> wrote:
> I also often need to use a function that is not available on the version of ncl that I run on the university cluster.
>
> If the function is available in an ncl script (like this one is) I just download the latest version of the script and save it in my own root directory. Then I just load it as
>
> load "~/WRFUserARW_v6.ncl" and the function is there.
>
>
> On Fri, Sep 2, 2011 at 11:16 AM, Chris Klich <caklich@gmail.com> wrote:
> Hi, I am currently using an older version of NCL (version 5.0.0) due to library issues on my network cluster. Therefore, I am not able to use wrf_ij_to_ll yet, which would most likely do what I need to do. However, I am trying to find the i,j point closest to my ref_lat and ref_lon for my inner (second) domain. Because I am using WRF/Chem, I need to know the center point in domain 2 to build my chemical emissions. As of right now, I am just guessing lat/lon points and running the 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/WRF_contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
>
> begin
>
> a = addfile("/sully/r0/cklich/WRF/wrfout_d02_2008-04-18_06:00:00.nc","r") ;Open file for reading
>
> loc = wrf_user_latlon_to_ij(a, 40, 145)
> print("Y/X location is: " +loc)
>
>
> end
>
> until I receive the i,j that is exactly in the middle of the domain. For example, if the grid size is 250x250, I guess points until I receive an i,j that is 125, 125. However, I have also read that this function's array is base 0, not base 1. Does this mean when I use this function I have to guess points until it is 124,124? Also, is there any easier way this can be done until I am able to upgrade NCL?
>
> Thank you.
>
> --
> Christopher Klich
> Graduate Student
> Florida State University
> (908) 208-9743
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>
> --
> Paula Doubrawa Moreira
> Graduate Student - Research Assistant
> International Arctic Research Center
> University of Alaska Fairbanks
>
>
>
> --
> Christopher Klich
> Graduate Student
> Florida State University
> (908) 208-9743
>
>
>
>
> --
> Paula Doubrawa Moreira
> Graduate Student - Research Assistant
> International Arctic Research Center
> University of Alaska Fairbanks
> _______________________________________________
> 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
Received on Tue Sep 6 08:48:54 2011