Re: about wrf_real.ncl

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Tue, 6 Feb 2007 10:04:17 -0700 (MST)

On Mon, 5 Feb 2007, Sir.Kadarsah wrote:

>
> Dear Dr.Mary Haley,
>
> When I execute the command
> ncl < wrf_real.ncl
> I get error:
>
> Copyright (C) 1995-2006 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 4.2.0.a034
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> warning:WorkstationClassInitialize:Unable to access rgb color database - named colors unsupported:[errno=2]:No such file or directory
> fatal:NhlCvtStringGenArrayToFloatGenArray:Unable to Convert "white"
> warning:Error retrieving resource wkBackgroundColor from DB - Using default value
> fatal:NhlCvtStringGenArrayToFloatGenArray:Unable to Convert "black"
> warning:Error retrieving resource wkForegroundColor from DB - Using default value
>
>
> Variable: times
> Type: string
> Total Size: 72 bytes
> 9 values
> Number of Dimensions: 1
> Dimensions and sizes: [9]
> Coordinates:
> Number Of Attributes: 2
> description : times in file
> _FillValue : missing
> (0) 2005-08-28_00:00:00
> (1) 2005-08-28_03:00:00
> (2) 2005-08-28_06:00:00
> (3) 2005-08-28_09:00:00
> (4) 2005-08-28_12:00:00
> (5) 2005-08-28_15:00:00
> (6) 2005-08-28_18:00:00
> (7) 2005-08-28_21:00:00
> (8) 2005-08-29_00:00:00
> fatal:MapV40DHClassPartInit: cannot open map data file: NhlMapData
> fatal:MapV40DHClassPartInit: error initializing map outline records
> fatal:Unable to initialize layer-Can't Create
> fatal:Unable to access object with id:-4
> fatal:PID #-4 can't be found in NhlSetValues
> fatal:Invalid plot ID=-4 passed to NhlGetBB
> fatal:NhlGetValues:PID #-4 is invalid
> fatal:Execute: Error occurred at or near line 1616 in file gsn_code.ncl
>
> fatal:Execute: Error occurred at or near line 2256 in file gsn_code.ncl
>
> fatal:Execute: Error occurred at or near line 7518 in file gsn_code.ncl
>
> fatal:Execute: Error occurred at or near line 818 in file WRFPlot.ncl
>
> fatal:Execute: Error occurred at or near line 66
>

Dear Kadarsah,

I think the problem is that the environment variable NCARG_ROOT is
either not set, or it is set incorrectly. It needs to be set to the
root directory of where NCL was installed on your system. It is
necessary to set this so NCL can find things like the map databases
and the fonts, which reside in the directory $NCARG_ROOT/lib/ncarg.

See:

   http://www.ncl.ucar.edu/Download/install.shtml#SetNCARG_ROOT

for more information about setting this environment variable.

Good luck,

--Mary

>
> and
> this is wrf_real.ncl
>
>
> ; Script to produce standard plots for a WRF real-data run,
> ; with the ARW coordinate dynamics option.
>
> load "WRFOptions.ncl" ; set basic plot options here
> ;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "gsn_code.ncl"
> load "WRFPlot.ncl"
> load "WRFUserARW.ncl"
> load "SkewTFunc.ncl"
>
> begin
> ;
> ; The WRF ARW input file.
> ; This needs to have a ".nc" appended, so just do it.
>
> a = addfile("wrfout.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_real_kadarsah")
>
> ; Debug information.
>
> debug = False
> ; debug = True
> icount = 0
>
> ; Basic Plot Information
>
> res_at_MainTitle = "REAL-TIME WRF"
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> ; What times and how many time steps are in the data set?
>
> FirstTime = True
> times = wrf_user_list_times(a) ; get times in the file
> ntimes = dimsizes(times) ; number of times in the file
>
> ; The specific pressure levels that we want the data
> ; interpolated to.
>
> pressure_levels = (/ 850., 700., 500., 300./) ; pressure levels to plot
> nlevels = dimsizes(pressure_levels) ; number of pressure levels
>
> ; This is the big loop over all of the time periods to process.
>
> do it = 0,ntimes-1,2
> ; do it = 0,2,2
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> ;Save some Time Information and create the Map Background
>
> time = it
> if (FirstTime) then
> times_sav = times(it)
> end if
> res_at_TimeLabel = times(it)
>
> map = wrf_map(wks,a,res)
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> ; First get the variables we will need
>
> slvl = wrf_user_getvar(a,"slvl",time) ; psl
> wrf_user_filter2d(slvl, 3) ; smooth psl
> qv = wrf_user_getvar(a,"QVAPOR",time) ; Qv
> tc = wrf_user_getvar(a,"tc",time) ; T in C
> td = wrf_user_getvar(a,"td",time) ; Td in C
> if ( it .eq. 0 ) then
> tc2 = tc(0,:,:) ; Use lowest T at time zero
> else
> tc2 = wrf_user_getvar(a,"T2",time) ; T2 in Kelvin
> tc2 = tc2-273.16 ; T2 in C
> end if
> if (it .eq. 0) then
> td2 = td(0,:,:) ; Use lowest Td at time zero
> else
> td2 = wrf_user_getvar(a,"td2",time) ; Td2 in C
> end if
> u = wrf_user_getvar(a,"ua",time) ; u averaged to mass points
> v = wrf_user_getvar(a,"va",time) ; v averaged to mass points
> uvm = wrf_user_getvar(a,"umeta",time) ; umet,vmet averaged to mass points
> if (it .eq. 0) then
> u10 = u(0,:,:) ; Use lowest level at time 0
> v10 = v(0,:,:)
> else
> u10 = wrf_user_getvar(a,"U10",time) ; u at 10 m, mass point
> v10 = wrf_user_getvar(a,"V10",time) ; v at 10 m, mass point
> end if
> w = wrf_user_getvar(a,"wa",time) ; vertical velocity, averaged to half levels
> w = 100.*w ; vert vel in cm/s
> p = wrf_user_getvar(a, "p",time) ; pressure is our vertical coordinate
> z = wrf_user_getvar(a, "Z",time) ; grid point height
> rh = wrf_user_getvar(a,"rh",time) ; relative humidity
> ter = wrf_user_getvar(a,"HGT",time) ; terrain height
>
> ; Get non-convective, convective and total precipitation
> ; Calculate tendency values
> rain_exp = wrf_user_getvar(a,"RAINNC",time)
> rain_con = wrf_user_getvar(a,"RAINC",time)
> rain_tot = rain_exp + rain_con
> if( FirstTime ) then
> rain_exp_save = rain_exp
> rain_con_save = rain_con
> rain_tot_save = rain_tot
> end if
> rain_exp_tend = rain_exp - rain_exp_save
> rain_con_tend = rain_con - rain_con_save
> rain_tot_tend = rain_tot - rain_tot_save
> ; Bookkeeping, just to allow the tendency at the next time step
> rain_exp_save = rain_exp
> rain_con_save = rain_con
> rain_tot_save = rain_tot
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> ; ************************************************************
>
> ; 1) Surface Temp + SLP + wind barbs
> ; Surface Td + wind barbs
>
>
> ; Plotting options for T
>
> tc_plane = 1.8*tc2+32. ; Turn temperature into Fahrenheit - Typical of a surface map
> opts_tc = res
> opts_tc_at_FieldTitle = "Surface Temperature"
> opts_tc_at_UnitLabel = "F"
> opts_tc_at_ContourParameters = (/ -20., 90., 5./)
> opts_tc_at_cnFillOn = True
> opts_tc_at_gsnSpreadColorEnd = -3 ; End third from the last color in color map
> contour_tc = wrf_contour(a,wks,tc_plane,opts_tc)
> print_opts("opts_tc", opts_tc, debug)
>
>
> ; Plotting options for Td
>
> td_plane = 1.8*td2+32. ; Turn temperature into Fahrenheit - Typical of a surface map
> opts_td = res
> opts_td_at_FieldTitle = "Surface Dew Point Temp"
> opts_td_at_UnitLabel = "F"
> opts_td_at_ContourParameters = (/ -20., 90., 5./)
> opts_td_at_cnFillOn = True
> opts_td_at_cnLinesOn = True
> opts_td_at_cnLineLabelsOn = True
> opts_td_at_cnLineLabelBackgroundColor = -1
> opts_td_at_gsnSpreadColorEnd = -3 ; End third from the last color in color map
> contour_td = wrf_contour(a,wks,td_plane,opts_td)
> print_opts("opts_td", opts_td, debug)
>
> ; Plotting options for SLP
>
> opts_psl = res
> ;opts_psl_at_FieldTitle = "Sea Level Pressure"
> ;opts_psl_at_UnitLabel = "mb"
> opts_psl_at_ContourParameters = (/ 900., 1100., 4. /)
> opts_psl_at_cnLineColor = "NavyBlue"
> opts_psl_at_cnHighLabelsOn = True
> opts_psl_at_cnLowLabelsOn = True
> opts_psl_at_cnLineLabelBackgroundColor = -1
> opts_psl_at_gsnContourLineThicknessesScale = 2.0
> contour_psl = wrf_contour(a,wks,slvl,opts_psl)
> print_opts("opts_psl", opts_psl, debug)
>
> ; Plotting options for Wind Vectors
>
> u_plane = u10*1.94386 ; Turn wind into knots - Typical of a surface map
> v_plane = v10*1.94386
> opts_vct = res
> opts_vct_at_FieldTitle = "Winds"
> opts_vct_at_UnitLabel = "kts"
> opts_vct_at_NumVectors = 47
> vector = wrf_vector(a,wks,u_plane, v_plane,opts_vct)
> print_opts("opts_vct", opts_vct, debug)
>
> ; MAKE PLOTS
>
> wrf_map_overlay(wks,map,(/contour_tc,contour_psl,vector/),False)
> print_header(icount,debug)
>
> wrf_map_overlay(wks,map,(/contour_td,vector/),False)
> print_header(icount,debug)
>
> ; Delete options and fields, so we don't have carry over
> delete(opts_tc)
> delete(opts_td)
> delete(opts_vct)
> delete(opts_psl)
> delete(tc_plane)
> delete(td_plane)
> delete(u_plane)
> delete(v_plane)
>
> ; ************************************************************
>
> ; 2) 2d precip plots
>
> ; Plotting options for Sea Level Pressure
>
> opts_psl = res
> opts_psl_at_ContourParameters = (/ 900., 1100., 2. /)
> opts_psl_at_cnLineColor = "Blue"
> opts_psl_at_cnInfoLabelOn = False
> opts_psl_at_cnHighLabelsOn = True
> opts_psl_at_cnLowLabelsOn = True
> opts_psl_at_gsnContourLineThicknessesScale = 1.0
> contour_psl = wrf_contour(a,wks,slvl,opts_psl)
> print_opts("opts_psl",opts_psl, debug)
>
>
> ; Plotting options for Precipitation
>
> opts_r = res
> opts_r_at_UnitLabel = "mm"
> opts_r_at_cnLevelSelectionMode = "ExplicitLevels"
> opts_r_at_cnLevels = (/ .1, .2, .4, .8, 1.6, 3.2, 6.4, \
> 12.8, 25.6, 51.2, 102.4/)
> opts_r_at_cnFillColors = (/"White","White","DarkOliveGreen1", \
> "DarkOliveGreen3","Chartreuse", \
> "Chartreuse3","Green","ForestGreen", \
> "Yellow","Orange","Red","Violet"/)
> opts_r_at_cnInfoLabelOn = False
> opts_r_at_cnConstFLabelOn = False
> print_opts("opts_r", opts_r, debug)
>
>
> ; MAKE PLOTS
> opts_r_at_cnFillOn = True
>
> ; Total Precipitation (color fill)
> opts_r_at_FieldTitle = "Total Precipitation"
> contour_tot = wrf_contour(a,wks, rain_tot, opts_r)
>
> ; Precipitation Tendency (color fill)
> opts_r_at_FieldTitle = "Precipitation Tendency"
> opts_r_at_SubFieldTitle = "from " + times_sav + " to " + times(it)
> contour_tend = wrf_contour(a,wks, rain_tot_tend,opts_r)
>
> ; Non-Convective Precipitation Tendency (color fill)
> opts_r_at_FieldTitle = "Explicit Precipitation Tendency"
> opts_r_at_SubFieldTitle = "from " + times_sav + " to " + times(it)
> contour_res = wrf_contour(a,wks,rain_exp_tend,opts_r)
>
> ; Convective Precipitation Tendency (contour lines)
> opts_r_at_FieldTitle = "Param Precipitation Tendency"
> opts_r_at_SubFieldTitle = "from " + times_sav + " to " + times(it)
> opts_r_at_cnFillOn = False
> opts_r_at_cnLineColor = "Red4"
> contour_prm = wrf_contour(a,wks,rain_con_tend,opts_r)
>
>
> ; Total Precipitation Tendency + SLP
>
> wrf_map_overlay(wks,map,(/contour_tend,contour_psl/),False)
> print_header(icount,debug)
>
> ; Total Precipitation
>
> wrf_map_overlay(wks,map,contour_tot,False)
> print_header(icount,debug)
>
> ; Non-Convective and Convective Precipiation Tendencies
>
> wrf_map_overlay(wks,map,(/contour_res,contour_prm/),False)
> print_header(icount,debug)
>
>
> ; Delete options, so we don't have carry over
> delete(opts_r)
> delete(opts_psl)
>
> ; ************************************************************
>
> ; 3) 3d plots: height, rh, vert vel, temp
>
> ; Loop over the selected pressure levels to plot
>
> do level = 0,nlevels-1
>
> pressure = pressure_levels(level)
>
> tc_plane = wrf_user_intrp3d(tc,p,ter,"h",pressure,0.)
> z_plane = wrf_user_intrp3d( z,p,ter,"h",pressure,0.)
> rh_plane = wrf_user_intrp3d(rh,p,ter,"h",pressure,0.)
> u_plane = wrf_user_intrp3d( u,p,ter,"h",pressure,0.)
> v_plane = wrf_user_intrp3d( v,p,ter,"h",pressure,0.)
> spd = (u_plane*u_plane + v_plane*v_plane)^(0.5) ; m/sec
> u_plane = u_plane*1.94386 ; kts
> v_plane = v_plane*1.94386 ; kts
>
>
> ; Plotting options for T
>
> opts_tc = res
> opts_tc_at_cnLineColor = "Red"
> opts_tc_at_ContourParameters = (/ 5.0 /)
> opts_tc_at_cnInfoLabelOrthogonalPosF = 0.06
> opts_tc_at_gsnContourLineThicknessesScale = 2.0
>
> ; Plotting options for Geopotential Heigh
>
> opts_z = res
> opts_z_at_cnLineColor = "NavyBlue"
> opts_z_at_gsnContourLineThicknessesScale = 3.0
>
> ; Plotting options for RH
>
> opts_rh = res
> opts_rh_at_ContourParameters = (/ 10., 90., 10./)
> opts_rh_at_cnFillOn = True
> opts_rh_at_cnFillColors = (/"White","White","White", \
> "White","Chartreuse","Green",\
> "Green3","Green4", \
> "ForestGreen","PaleGreen4"/)
>
>
>
> ; Plotting options for Wind Speed
>
> opts_spd = res
> opts_spd_at_FieldTitle = "Wind Speed"
> opts_spd_at_UnitLabel = "m/s"
> opts_spd_at_cnLineColor = "MediumSeaGreen"
> opts_spd_at_ContourParameters = (/ 10. /)
> opts_spd_at_cnInfoLabelOrthogonalPosF = 0.06
> opts_spd_at_gsnContourLineThicknessesScale = 3.0
>
> ; Plotting options for Wind Vectors
>
> opts_vct = res
> opts_vct_at_FieldTitle = "Winds"
> opts_vct_at_UnitLabel = "kts"
> opts_vct_at_NumVectors = 47
>
> ; MAKE PLOTS
>
> if ( pressure .eq. 850 ) then
> ; temp, rh, height, barbs
>
> opts_z_at_ContourParameters = (/ 20.0 /)
>
> contour_rh = wrf_contour(a,wks,rh_plane,opts_rh)
> print_opts("opts_rh", opts_rh, debug)
> contour_tc = wrf_contour(a,wks,tc_plane,opts_tc)
> print_opts("opts_tc", opts_tc, debug)
> contour_height = wrf_contour(a,wks, z_plane,opts_z)
> print_opts("opts_z", opts_z, debug)
> vector = wrf_vector(a,wks,u_plane, v_plane,opts_vct)
> print_opts("opts_vct", opts_vct, debug)
>
> wrf_map_overlay(wks,map,(/contour_rh,contour_tc,contour_height, \
> vector/),False)
> print_header(icount,debug)
>
> end if
>
> if ( pressure .eq. 700 ) then
> ; temp, height, barbs
>
> opts_z_at_ContourParameters = (/ 30.0 /)
>
> contour_tc = wrf_contour(a,wks,tc_plane,opts_tc)
> print_opts("opts_tc", opts_tc, debug)
> contour_height = wrf_contour(a,wks, z_plane,opts_z)
> print_opts("opts_z", opts_z, debug)
> vector = wrf_vector(a,wks,u_plane, v_plane,opts_vct)
> print_opts("opts_vct", opts_vct, debug)
>
> wrf_map_overlay(wks,map,(/contour_tc,contour_height, \
> vector/),False)
> print_header(icount,debug)
>
> end if
>
>
> if ( pressure .eq. 500 ) then
> ; temp, height, barbs
>
> opts_z_at_ContourParameters = (/ 60.0 /)
>
> contour_tc = wrf_contour(a,wks,tc_plane,opts_tc)
> print_opts("opts_tc", opts_tc, debug)
> contour_height = wrf_contour(a,wks, z_plane,opts_z)
> print_opts("opts_z", opts_z, debug)
> vector = wrf_vector(a,wks,u_plane, v_plane,opts_vct)
> print_opts("opts_vct", opts_vct, debug)
>
> wrf_map_overlay(wks,map,(/contour_tc,contour_height, \
> vector/),False)
> print_header(icount,debug)
>
> end if
>
> if ( pressure .eq. 300 ) then
> ; windspeed, height, barbs
>
> opts_z_at_ContourParameters = (/ 60.0 /)
>
> contour_spd = wrf_contour(a,wks,spd,opts_spd)
> print_opts("opts_spd", opts_spd, debug)
> contour_height = wrf_contour(a,wks, z_plane,opts_z)
> print_opts("opts_z", opts_z, debug)
> vector = wrf_vector(a,wks,u_plane, v_plane,opts_vct)
> print_opts("opts_vct", opts_vct, debug)
>
> wrf_map_overlay(wks,map,(/contour_spd,contour_height, \
> vector/),False)
> print_header(icount,debug)
>
> end if
>
> end do ; END OF 3D LOOP
>
> ; Delete options and fields, so we don't have carry over
> delete(opts_tc)
> delete(opts_z)
> delete(opts_rh)
> delete(opts_spd)
> delete(opts_vct)
> delete(z_plane)
> delete(tc_plane)
> delete(rh_plane)
> delete(u_plane)
> delete(v_plane)
> delete(spd)
>
> ; ************************************************************
>
> ; 4) xsection plots
>
> do ip = 1, 2
> ; |
> ; |
> ; Two different xsections, angle=90 is |, angle=0 is ------
> ; |
> ; |
> if(ip .eq. 1) then
> angle = 90.
> Orientation = "Cross-Section Orientation : West-East"
> else
> angle = 0.
> Orientation = "Cross-Section Orientation : South-North"
> end if
>
> ; Options for XY Plots
> opts_xy = res
> opts_xy_at_tiXAxisString = "Number of Grid Points"
> opts_xy_at_tiYAxisString = "Height (km)"
> opts_xy_at_PlotType = "XY"
> opts_xy_at_AspectRatio = 0.75
> opts_xy_at_PlotOrientation = Orientation
> opts_xy_at_cnMissingValPerimOn = True
> opts_xy_at_cnMissingValFillColor = 0
> opts_xy_at_cnMissingValFillPattern = 11
> opts_xy_at_tmYLMode = "Explicit"
> opts_xy_at_tmYLValues = ispan(0,100,10)
> opts_xy_at_tmYLLabels = ispan(0,20,2)
> opts_xy_at_tiXAxisFontHeightF = 0.020
> opts_xy_at_tiYAxisFontHeightF = 0.020
> opts_xy_at_tmXBMajorLengthF = 0.02
> opts_xy_at_tmYLMajorLengthF = 0.02
>
>
> ; Build planes of data for temp and rh
>
> dimsrh = dimsizes(rh)
> plane = new(2,float)
> plane = (/ dimsrh(2)/2, dimsrh(1)/2 /)
>
> rh_plane = wrf_user_intrp3d(rh,z,ter,"v",plane,angle)
> tc_plane = wrf_user_intrp3d(tc,z,ter,"v",plane,angle)
>
> ; Plotting options for RH
>
> opts_rh = opts_xy
> opts_rh_at_ContourParameters = (/ 10., 90., 10. /)
> opts_rh_at_cnFillOn = True
> opts_rh_at_cnFillColors = (/"White","White","White", \
> "White","Chartreuse","Green", \
> "Green3","Green4", \
> "ForestGreen","PaleGreen4"/)
>
> ; Plotting options for Temperature
>
> opts_tc = opts_xy
> opts_tc_at_ContourParameters = (/ 5. /)
>
> ; Get the contour info for the rh and temp
>
> contour_tc = wrf_contour(a,wks,tc_plane,opts_tc)
> print_opts("opts_tc", opts_tc, debug)
> contour_rh = wrf_contour(a,wks,rh_plane,opts_rh)
> print_opts("opts_rh", opts_rh, debug)
>
> ; Overlay the fields
>
> wrf_overlay(wks,(/contour_rh,contour_tc/),False)
> print_header(icount,debug)
>
> ; Delete options and fields, so we don't have carry over
> delete(opts_tc)
> delete(opts_rh)
> delete(tc_plane)
> delete(rh_plane)
> end do ; end of xsection plot rotation
>
> ; ************************************************************
>
> ; 5) skew-T plots
>
> u = uvm(0,:,:,:)*1.94386
> v = uvm(1,:,:,:)*1.94386
>
> do ip = 12, 13
>
> ; Define a few skew-T plotting options
>
> skewtOpts = True
> dataOpts = True
> dataOpts_at_Parcel = 0
> dataOpts_at_WspdWdir = False ; wind speed and dir [else: u,v]
> dataOpts_at_PlotWindH = False ; plot wind barbs at h lvls [pibal; special]
> dataOpts_at_HspdHdir = True ; wind speed and dir [else: u,v]
>
> ; Choose the lat/lon of the skew-T plots
>
> ip_lats = (/ 39.77, 37.77, 41.12, 41.30, 39.07, 35.22, \
> 37.23, 35.22, 31.95, 32.50, 43.56, 38.94, \
> 35.27 /)
> ip_lons = (/-104.87, -99.97, -100.67, -95.90, -95.62, -97.45, \
> -93.38, -101.72, -102.20, -97.30, -116.24, -77.44, \
> -75.55 /)
> ip_locs = (/"Denver (DEN)", "Dodge City (DDC)", "North Platte (LBF)", \
> "Omaha (OMA)", "Topeka (TOP)", "Norman (OUN)", \
> "Springfield (SGF)", "Amarillo (AMA)", "Midland (MAF)", \
> "Fort Worth (FWD)", "Boise (BOI)", "Washington DC (IAD)", \
> "Cape Hatteras (HAT)"/)
>
> ip_locs = ip_locs + ", Valid at " + times(it)
> locr = wrf_user_find_ij_lat_long(a, ip_lats(ip-1), ip_lons(ip-1))
> loc = floattointeger(locr)
>
> ; Get the skew-T background
>
> skewtOpts_at_tiMainString = ip_locs(ip-1)
> skewt_bkgd = skewT_BackGround (wks, skewtOpts)
>
> ; Draw the skew-T plot
>
> draw (skewt_bkgd)
> skewT_data = skewT_PlotData(wks, skewt_bkgd, p(:,loc(0), loc(1)), \
> tc(:,loc(0), loc(1)), \
> td(:,loc(0), loc(1)), \
> z(:,loc(0), loc(1)), \
> -u(:,loc(0), loc(1)), \
> -v(:,loc(0), loc(1)), \
> dataOpts)
>
> ; Close the frame
>
> frame(wks)
> print_header(icount,debug)
>
> ; Delete the temporary arrays
>
> delete(skewtOpts)
> delete(dataOpts)
> delete(skewT_data)
> delete(skewt_bkgd)
>
> end do
>
> ; ************************************************************
>
> times_sav = times(it)
> FirstTime = False
> end do ; end of the time loop
>
> end
>
>
> I am attaching the script wrf_real.ncl for your reveiw.
> I am still in learning mode and want to understand how to resolve this
> error.
>
> Thank you for your help
> Kadarsah
>
>
> *Cogito Ergo Sum*
> ----Kadarsah---
> ---------------------------------------------
> Meteorologist
> Kanda Laboratory
> Tokyo Institute of Technology
>
> **************
> Homepage:http:/www.kadarsah.50megs.com
>
>
>
>
>
>
>
> ---------------------------------
> It's here! Your new message!
> Get new email alerts with the free Yahoo! Toolbar.
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 06 2007 - 10:04:17 MST

This archive was generated by hypermail 2.2.0 : Tue Feb 06 2007 - 15:01:55 MST