Re: NCL script for Idealized cases

This is from old web page and the script is from old NCL tar file
that we provide. You may download
this tar from the site where you download the WRF code (http://
www.mmm.ucar.edu/wrf/users/download/get_source.html),
and look for 'Older WRF NCL scripts'. These idealized scripts have
not been ported to NCL 4.3.1 at this time.

wrfhelp

On Oct 16, 2007, at 5:23 PM, Mary Haley wrote:

> Hi Guan,
>
> You ask a good question, but I'm afraid I don't know the answer to
> this one. I looked around on the mmm website for the scripts, but
> couldn't find a link to them.
>
> I've cc'ed wrfhelp_at_ucar.edu as they should be able to point you to the
> scripts.
>
> --Mary
>
>
> On Tue, 16 Oct 2007, Guan Song wrote:
>
>>
>> Hi,all
>>
>> I am running newest WRF with idealized cases (such as 2dhill).
>> Does anyone have the NCL scripts for wrf output?
>>
>> I looked at a website, which supposes to have scripts for all
>> idealized cases.
>>
>> http://www.mmm.ucar.edu/mm5/WRF_post/WRF_NCL/b_wave.html
>>
>> however,it does not have any link for this script.
>>
>> I found the old version of NCL scripts for WRF output, however, it
>> did not work due to Z variable is available anymore from WRF output.
>>
>>
>> The old script (show as follows) to plot standard 2dhill show
>> undefined variable Z (grid point height) so it cannot plot
>> anything. Does anyone know what
>> variables we probably need to replace Z and let the script work
>> for WRF output? or do you have any other scripts which should
>> work in this case.
>>
>>
>> ; Script to produce standard plots for a WRF hill2d run
>>
>> load "WRFOptions.ncl"
>> ;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "gsn_code.ncl"
>> load "WRFPlot.ncl"
>> load "WRFUserARW.ncl"
>>
>> begin
>> ;
>> ; The WRF ARW input file.
>> ; This needs to have a ".nc" appended, so just do it.
>>
>> a = addfile("../WRFV2/run/wrfout_d01_0001-01-01_00:00:00.nc","r")
>>
>> ; We generate plots, but what kind do we prefer?
>>
>> type = "x11"
>> ; type = "pdf"
>> ; type = "ps"
>> ; type = "ncgm"
>>
>> wks = gsn_open_wks(type,"wrf_hill2d")
>>
>> ; Debug information.
>>
>> debug = False
>> ; debug = True
>> icount = 0
>>
>> ; Basic Plot options, including header information
>>
>> res_at_MainTitle = "WRF HILL2x"
>>
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>
>> ; What times and how many time steps are in the data set?
>>
>> times = wrf_user_list_times(a) ; get times in the file
>> ntimes = dimsizes(times) ; number of times in the file
>>
>> ; The specific plane we want to plot data on
>>
>> plane = (/ 20., 0./) ; (x,y) point for vertical plane
>> angle = 90.0
>> pii = 3.14159
>> aspect_ratio = .7
>>
>> ; This is the big loop over all of the time periods to process.
>>
>> do it = 0,ntimes-1,30
>>
>> time = it
>> res_at_TimeLabel = times(it)
>> res_at_AspectRatio = aspect_ratio
>> res_at_PlotOrientation = "(" + plane(0)+","+plane(1) +
>> ")" + " angle " + angle
>>
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>> ; First get the variables we will need
>>
>> u_in = wrf_user_getvar(a,"ua",time) ; ua is u averaged to
>> mass points
>> v_in = wrf_user_getvar(a,"va",time) ; va is v averaged to
>> mass points
>> w_in = wrf_user_getvar(a,"wa",time) ; vertical velocity
>> th_in = wrf_user_getvar(a,"th",time) ; get temperature (C)
>> z_in = wrf_user_getvar(a, "Z",time) ; grid point height
>> ter = wrf_user_getvar(a,"HGT",time) ; need terrain height
>> sometimes
>>
>> u = u_in(0:34,:,:)
>> v = v_in(0:34,:,:)
>> w = w_in(0:34,:,:)
>> th = th_in(0:34,:,:)
>> z = z_in(0:34,:,:)
>>
>> u_plane = wrf_user_intrp3d( u,z,ter,"v",plane,angle)
>> v_plane = wrf_user_intrp3d( v,z,ter,"v",plane,angle)
>> w_plane = wrf_user_intrp3d( w,z,ter,"v",plane,angle)
>> th_plane = wrf_user_intrp3d(th,z,ter,"v",plane,angle)
>>
>> vel_normal = u_plane*cos(2.*pii*angle/360.) - v_plane*sin
>> (2.*pii*angle/360.)
>> vel_tangent = u_plane*sin(2.*pii*angle/360.) + v_plane*cos
>> (2.*pii*angle/360.)
>> vel_tangent = vel_tangent - 10.
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>
>> ; Theta
>> opts_th = res
>> opts_th_at_FieldTitle = th_in_at_description
>> opts_th_at_cnFillOn = True
>> opts_th_at_gsnSpreadColorEnd = -10
>> contour_th = wrf_contour(a,wks,th_plane,opts_th)
>> print_opts("opts_th", opts_th, debug)
>>
>> ; Vertical Velocity
>> opts_w = res
>> opts_w_at_FieldTitle = w_in_at_description
>> contour_w = wrf_contour(a,wks, w_plane,opts_w)
>> print_opts("opts_w", opts_w, debug)
>>
>> ; Vel Tangent
>> opts_vt = res
>> opts_vt_at_FieldTitle = "Perturbation u"
>> opts_vt_at_UnitLabel = "m/s"
>> contour_vt = wrf_contour(a,wks,vel_tangent,opts_vt)
>> print_opts("opts_vt", opts_vt, debug)
>>
>>
>> wrf_overlay(wks,(/contour_th, contour_vt/),False)
>> print_header(icount,debug)
>> wrf_overlay(wks,(/contour_th, contour_w/),False)
>> print_header(icount,debug)
>>
>> ; ************************************************************
>>
>> end do ; end of the time loop
>>
>> end
>>
>> ;*****************************************************************8
>> Thanks in advance.
>>
>> Guan Song
>> Environmental sciences department
>> University of virginia
>> Charlottesville,VA,
>>
>>
>>> From: ncl-talk-request_at_ucar.edu
>>> Subject: ncl-talk Digest, Vol 47, Issue 19
>>> To: ncl-talk_at_ucar.edu
>>> Date: Tue, 16 Oct 2007 12:00:06 -0600
>>>
>>> Send ncl-talk mailing list submissions to
>>> ncl-talk_at_ucar.edu
>>>
>>> To subscribe or unsubscribe via the World Wide Web, visit
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>> or, via email, send a message with subject or body 'help' to
>>> ncl-talk-request_at_ucar.edu
>>>
>>> You can reach the person managing the list at
>>> ncl-talk-owner_at_ucar.edu
>>>
>>> When replying, please edit your Subject line so it is more specific
>>> than "Re: Contents of ncl-talk digest..."
>>>
>>>
>>> Today's Topics:
>>>
>>> 1. newbie documentation and julian date question (Mary Jo
>>> Brodzik)
>>> 2. Re: newbie documentation and julian date question (Dennis
>>> Shea)
>>> 3. Re: newbie documentation and julian date question
>>> (Mary Jo Brodzik)
>>> 4. Floattoint truncating issue. (Wang, Yi)
>>> 5. Re: ncl error (Bernie Baer)
>>> 6. Re: ncl error (bbaer_at_internode.on.net)
>>> 7. Re: Problems with Stereographic projection plots (Mary Haley)
>>>
>>>
>>> --------------------------------------------------------------------
>>> --
>>>
>>> Message: 1
>>> Date: Mon, 15 Oct 2007 12:35:15 -0600 (MDT)
>>> From: Mary Jo Brodzik <brodzik_at_wuzzles.colorado.edu>
>>> Subject: newbie documentation and julian date question
>>> To: ncl-talk_at_ucar.edu
>>> Message-ID: <Pine.LNX.4.64.0710151210300.24524_at_wuzzles.colorado.edu>
>>> Content-Type: TEXT/PLAIN; charset=US-ASCII
>>>
>>>
>>> Hello, everyone,
>>>
>>> I'm new to NCL, and with Dennis Shea's help (thanks, Dennis!)
>>> I've gotten
>>> a little script working that converts a couple of flat binary
>>> files into
>>> netCDF. I have three questions:
>>>
>>> 1) I'll bet that there are julian date format convertors out there,
>>> somewhere. Can someone let me know where they are?
>>>
>>> 2) As a general logistical question, is there some way to query
>>> for help
>>> on a general topic at the command line? Something like "help
>>> date"
>>> If not, then how do NCL users find out about routines they
>>> might be
>>> looking for, and/or how to call those routines once they know
>>> they're
>>> there?
>>>
>>> 3) I'm an emacs user, and have found and installed the ncl.el
>>> file for
>>> formatting, but I'm wondering if there's anything like an
>>> ncl-shell-mode, where I can do line-by-line debugging with the
>>> script
>>> in one window, tied to an ncl execution shell in the other?
>>>
>>> Thanks,
>>> Mary Jo Brodzik
>>>
>>> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
>>> ^^^
>>> Mary Jo Brodzik, Lead Scientific Programmer, 303-492-8263
>>> NSIDC/CIRES, Univ. of Colo. at Boulder, 449 UCB, Boulder, CO
>>> 80309-0449
>>> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
>>> ^^^
>>>
>>>
>>>
>>> ------------------------------
>>>
>>> Message: 2
>>> Date: Mon, 15 Oct 2007 14:08:41 -0600
>>> From: Dennis Shea <shea_at_ucar.edu>
>>> Subject: Re: newbie documentation and julian date
>>> question
>>> To: brodzik_at_nsidc.org
>>> Cc: ncl-talk_at_ucar.edu
>>> Message-ID: <4713C8C9.7050104_at_ucar.edu>
>>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>>>
>>> Mary Jo Brodzik wrote:
>>>> Hello, everyone,
>>>>
>>>> I'm new to NCL, and with Dennis Shea's help (thanks, Dennis!)
>>>> I've gotten
>>>> a little script working that converts a couple of flat binary
>>>> files into
>>>> netCDF. I have three questions:
>>>>
>>>> 1) I'll bet that there are julian date format convertors out there,
>>>> somewhere. Can someone let me know where they are?
>>>>
>>> [a] http://www.ncl.ucar.edu/
>>> [b1] Click "Documentation"; Click "Functions";
>>> Under "Earth Science" Click "Date"
>>>
>>> greg2jul,
>>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/greg2jul.shtml>
>>> jul2greg, day_of_year, ut_calendar, ut_inv_calendar
>>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/jul2greg.shtml>
>>>
>>> [b2] In the upper left of *any* NCL page type what
>>> you want
>>> in the Search box. eg: julian
>>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/jul2greg.shtml>
>>>> 2) As a general logistical question, is there some way to query
>>>> for help
>>>> on a general topic at the command line? Something like "help
>>>> date"
>>>>
>>> Nope .... bummer
>>>
>>>> If not, then how do NCL users find out about routines they
>>>> might be
>>>> looking for, and/or how to call those routines once they know
>>>> they're
>>>> there?
>>>>
>>> use the "search" box in the upper left. Usually, there are
>>> examples of
>>> usage at the bottom of the page.
>>>
>>>
>>>
>>>> 3) I'm an emacs user, and have found and installed the ncl.el
>>>> file for
>>>> formatting, but I'm wondering if there's anything like an
>>>> ncl-shell-mode, where I can do line-by-line debugging with
>>>> the script
>>>> in one window, tied to an ncl execution shell in the other?
>>>>
>>>
>>> I don't use EMACS ... somebody else must answer.
>>>
>>> Regards
>>> D
>>>
>>>> Thanks,
>>>> Mary Jo Brodzik
>>>>
>>>> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
>>>> ^^^^
>>>> Mary Jo Brodzik, Lead Scientific Programmer, 303-492-8263
>>>> NSIDC/CIRES, Univ. of Colo. at Boulder, 449 UCB, Boulder, CO
>>>> 80309-0449
>>>> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
>>>> ^^^^
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk_at_ucar.edu
>>>> 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 |
>>> ======================================================
>>>
>>>
>>>
>>> ------------------------------
>>>
>>> Message: 3
>>> Date: Mon, 15 Oct 2007 14:28:25 -0600 (MDT)
>>> From: Mary Jo Brodzik <brodzik_at_wuzzles.colorado.edu>
>>> Subject: Re: newbie documentation and julian date
>>> question
>>> To: Dennis Shea <shea_at_ucar.edu>
>>> Cc: ncl-talk_at_ucar.edu
>>> Message-ID: <Pine.LNX.4.64.0710151423570.24524_at_wuzzles.colorado.edu>
>>> Content-Type: TEXT/PLAIN; charset=US-ASCII
>>>
>>>
>>> Thanks, Dennis, the search facility and documentation pages are
>>> exactly
>>> what I'm looking for.
>>>
>>> (Now *that's* embarrassing: I've been using the red nav bar on
>>> the left
>>> side of those pages, and never saw the yellow nav across the
>>> top. I don't
>>> think there's such a thing as "red-yellow" color-blind that I
>>> could even
>>> use as an excuse! Thank you for being so polite about it,
>>> too. ;-) )
>>>
>>> --Mary Jo
>>>
>>> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
>>> ^^^
>>> Mary Jo Brodzik, Lead Scientific Programmer, 303-492-8263
>>> NSIDC/CIRES, Univ. of Colo. at Boulder, 449 UCB, Boulder, CO
>>> 80309-0449
>>> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
>>> ^^^
>>>
>>>
>>>
>>> ------------------------------
>>>
>>> Message: 4
>>> Date: Mon, 15 Oct 2007 14:01:27 -0700
>>> From: "Wang, Yi" <yi.wang_at_pnl.gov>
>>> Subject: Floattoint truncating issue.
>>> To: "ncl-talk_at_ucar.edu" <ncl-talk_at_ucar.edu>
>>> Message-ID: <C3392337.5AC8%yi.wang_at_pnl.gov>
>>> Content-Type: text/plain; charset="US-ASCII"
>>>
>>> Hi Folks:
>>>
>>> For example, I have a time series from 19961010 to 19961020.
>>>
>>> When reading in as float, I get this (print-out):
>>>
>>> (0) 1.996101e+07
>>> (1) 1.996101e+07
>>> (2) 1.996101e+07
>>> (3) 1.996101e+07
>>> (4) 1.996101e+07
>>> (5) 1.996102e+07
>>> (6) 1.996102e+07
>>> (7) 1.996102e+07
>>> (8) 1.996102e+07
>>> (9) 1.996102e+07
>>> (10) 1.996102e+07
>>>
>>>
>>> To do a "floattoint" on above series, I get this (print-out):
>>>
>>> (0) 19961010
>>> (1) 19961012 "" (should be 19961011)""
>>> (2) 19961012
>>> (3) 19961012 ""(should be 19961013)""
>>> (4) 19961014
>>> (5) 19961016
>>> (6) 19961016
>>> (7) 19961016
>>> (8) 19961018
>>> (9) 19961020
>>> (10) 19961020
>>>
>>> The actual date should go from 19961010 to 19961020. How do I fix
>>> this
>>> floattoint trancating issue?
>>>
>>> Thanks,
>>>
>>> Yi
>>>
>>>
>>> ------------------------------
>>>
>>> Message: 5
>>> Date: Tue, 16 Oct 2007 07:27:41 +1000
>>> From: "Bernie Baer" <bbaer_at_internode.on.net>
>>> Subject: Re: ncl error
>>> To: "Mary Haley" <haley_at_ucar.edu>
>>> Cc: ncl-talk_at_ucar.edu
>>> Message-ID: <01f001c80f72$386e24d0$6400a8c0_at_amd>
>>> Content-Type: text/plain; format=flowed; charset="iso-8859-1";
>>> reply-type=response
>>>
>>> Hi Mary,
>>> thanks for your response. The script is part of a suite I got
>>> from here:
>>> http://www.drjack.info/RASP/DOWNLOAD/
>>> and the original version is WRF/NCL/skewt_func.ncl from the tarball
>>> rasp_ncl.tar.gz (I have to double check that is the
>>> tarball when I get to work) and the failing line was at 770 in
>>> the original.
>>> I'm not sure of the NCL version right now.
>>> I'm not an NCL programmer, just a Unix hacker! :-)
>>> More info: I am trying to plot three meteorological soundings,
>>> the first
>>> one, Camden, plots OK but the next two do not.
>>> Hope this is sufficient info if not please let me know.
>>> Kind regards, Bernie in Australia.
>>>
>>>
>>>
>>>
>>> ------------------------------
>>>
>>> Message: 6
>>> Date: Tue, 16 Oct 2007 14:00:00 +0930
>>> From: "bbaer_at_internode.on.net" <bbaer_at_internode.on.net>
>>> Subject: Re: ncl error
>>> To: ncl-talk_at_ucar.edu
>>> Message-ID: <47143e48.53.11fe.10128_at_internode.on.net>
>>> Content-Type: text/plain; charset="iso-8859-1"
>>>
>>> Mary et al,
>>> my problem has been fixed by DrJack :
>>>
>>> http://www.drjack.info/twiki/bin/view/RASPop/AlternateVersions
>>>
>>> Thank you for your help!
>>> Regards, Bernie.
>>>
>>>
>>> ------------------------------
>>>
>>> Message: 7
>>> Date: Tue, 16 Oct 2007 09:17:50 -0600 (MDT)
>>> From: Mary Haley <haley_at_ucar.edu>
>>> Subject: Re: Problems with Stereographic projection plots
>>> To: ncl-talk_at_ucar.edu
>>> Message-ID: <Pine.SOC.4.61.0710160906290.20660_at_windom.ucar.edu>
>>> Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed
>>>
>>>
>>> If you are using the wrf_xxxx scripts to do the plotting, then a fix
>>> was issued in which the wrong value was being used for the center
>>> latitude for the stereographic projection case. It was changed from
>>> TRUELAT1 to CEN_LAT (both attributes of the variable on the WRF
>>> file).
>>>
>>> This fix should make it into the next release of NCL, hopefully
>>> the end of this month.
>>>
>>> If you have V4.3.1, you can easily make the change yourself. You
>>> may
>>> need "root" permissions to do this:
>>>
>>> 1. Edit the file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/
>>> WRFUserARW.ncl
>>>
>>> 2. Change the line:
>>> opts_at_mpCenterLatF = get_res_value_keep(opts,
>>> "mpCenterLatF", in_file_at_TRUELAT1)
>>>
>>> to:
>>> opts_at_mpCenterLatF = get_res_value_keep(opts,
>>> "mpCenterLatF", in_file_at_CEN_LAT)
>>>
>>> --Mary
>>>
>>>
>>> On Fri, 12 Oct 2007, Andrea Hahmann wrote:
>>>
>>>> Hi!
>>>>
>>>> I have been having problems trying to plot a land mask (or land
>>>> use)
>>>> field on a stereographic projection. The data comes from a WRF
>>>> "geo"
>>>> file. When I use res_at_tfDoNDCOverlay = True, then the coastline
>>>> of the
>>>> data does not match the drawn geographic coastline. When I turn
>>>> it off,
>>>> the coastlines line up but I get a strange white band on the
>>>> right side
>>>> of the plot (also a small one on the bottom). I have
>>>> experimented with
>>>> various parameters to set the projection (CEN_LAT, CEN_LON, and
>>>> STAND_LON), but nothing seems to help.
>>>>
>>>> Any other ideas of what might be going wrong?
>>>>
>>>> The data is on an ftp site
>>>> (ftp://ftp.rap.ucar.edu/pub/hahmann/geo_em.d01.nc) if anybody is
>>>> interested in trying the script. The script and plot are attached.
>>>>
>>>> Thank you so much!
>>>> Andrea
>>>>
>>>> --
>>>> ----------------------------------------------------------------
>>>> Andrea N. Hahmann, Ph.D.
>>>> Research Applications Laboratory
>>>> Natl. Center for Atmospheric Research Phone: 1-303-497-8383
>>>> PO BOX 3000 Fax: 1-303-497-8401
>>>> Boulder, CO 80301 hahmann_at_ucar.edu
>>>> ----------------------------------------------------------------
>>>>
>>>>
>>>>
>>>
>>>
>>> ------------------------------
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk_at_ucar.edu
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>> End of ncl-talk Digest, Vol 47, Issue 19
>>> ****************************************
>>
>> _________________________________________________________________
>> Boo! Scare away worms, viruses and so much more! Try Windows Live
>> OneCare!
>> http://onecare.live.com/standard/en-us/purchase/trial.aspx?
>> s_cid=wl_hotmailnews
Received on Tue Oct 16 2007 - 17:56:09 MDT