[Wrf-users] wrf

[Guan Song]

> Hi all
>
> I run the idealized 2D hill case. I just changed e_sn to 302 (in namelist.input) and obtain the output.
>
> I used NCL to examine output.
>
> I wrote the script as follows, it seems work, but not in the right way. I attach the script below.I also attach the figure and the error I got.
> I already got WRFHELP's response to remove the map call from my script.
> ___________________________________
> mpres=True
> map=wrf_map(wks,a,mpres)
>
> wrf_map_overlay(wks,map,(/contour/),True)
> _____________________________________
> However if I removed the following statement,it did not work. Nothing is pop up.
> So I have to keep these statement.
>
> I am not sure what is wrong with my script, and I really need your help.
>
> thanks a lot in advance.
>
> __________________________________________________________--
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> begin
>
> ;WRF ARW INPUT FILE
>
> ;a=addfile("../wrfout_d01_2003-07-15_00-00-00.nc","r")
> a=addfile("../public/WRFV2/run/wrfout_d01_0001-01-01_00:00:00.nc","r")
> type="x11"
> wks=gsn_open_wks(type,"plt_Surface")
> ;SET BASIC OPTIONS
>
> ARWres=True
> ARWres at MainTitle="REAL-TIME WRF"
>
>
> ;---------------------------------------------------------------
> ;LETS FIND OUT HOW MANY TIMES ARE IN THE DATASETS.
>
> times=wrf_user_list_times(a) ; GET TIMES IN THE FILE
> ntimes=dimsizes (times) ;NUMBER OF TIMES IN THE FILE
>
> do it=0,ntimes-1,10 ;WE ARE ONLY INTERESTED IN THE LAST TIME
> ;ARWres at TimeLabel = times (it);keep some time information
> mpres=True
> map=wrf_map(wks,a,mpres)
> u=wrf_user_getvar(a,"U",it)
> v=wrf_user_getvar(a,"V",it)
>
> hgt=wrf_user_getvar(a,"HGT",it)
>
> opts=ARWres
>
> opts at cnMonoFillPattern =True
> opts at cnFillOn=True
> opts at pmLabelBarDisplayMode = "Always" ; Turn on label bar.
> contour=wrf_contour(a,wks,u(:,:,4),opts)
> wrf_map_overlay(wks,map,(/contour/),True)
> end do
> end
> -------------------------------------------------------------------
>
> Copyright (C) 1995-2007 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 4.3.1
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
>
>
> Variable: times
> Type: string
> Total Size: 88 bytes
> 11 values
> Number of Dimensions: 1
> Dimensions and sizes: [11]
> Coordinates:
> Number Of Attributes: 2
> description : times in file
> _FillValue : missing
> (0) 0001-01-01_00:00:00
> (1) 0001-01-01_01:00:00
> (2) 0001-01-01_02:00:00
> (3) 0001-01-01_03:00:00
> (4) 0001-01-01_04:00:00
> (5) 0001-01-01_05:00:00
> (6) 0001-01-01_06:00:00
> (7) 0001-01-01_07:00:00
> (8) 0001-01-01_08:00:00
> (9) 0001-01-01_09:00:00
> (10) 0001-01-01_10:00:00
> warning:ContourPlotInitialize: scalar field is constant; ContourPlot not possible:[errno=1102]
>
> ^^^^^^^^^^^^^^^^^^^^^^^^^^
> Best regards,
> Guan
>
> ________________________________
>
> _________________________________________________________________
> Peek-a-boo FREE Tricks & Treats for You!
> http://www.reallivemoms.com?ocid=TXT_TAGHM&loc=us

_________________________________________________________________
Peek-a-boo FREE Tricks & Treats for You!
http://www.reallivemoms.com?ocid=TXT_TAGHM&loc=us