[ncl-talk] Help Again With Plotting NARR Data

Kerwyn Texeira ktish86 at gmail.com
Tue Aug 2 15:43:49 MDT 2016


Hi ncl-talk,

I have made changes to my NARR script but it is still not plotting a zoomed
in overlay plot.  I tried several times playing around with the script by
commenting and adding stuff but I'm not seeing much of any progress.

When I commented out  all the tfDoNDCOverlay and left lat2d and lon2d
uncommented, I got the entire USA plotted with both temp, geopotential
height and winds overlayed.

When I commented out lat2d and lon2d and left tfDoNDCOverlay uncommented, I
got mainly a plain figure with one or two wind barbs. (temp_narr) figure
attached.

Kindly help, please.  I do appreciate the help that I'm getting but I'm
just not quite reach my ultimate goal which is a zoomed in overlay plot of
temps, winds, and geopotential height.

Thanks

Some outputs:

Variable: temp
Type: float
Total Size: 386692 bytes
            96673 values
Number of Dimensions: 2
Dimensions and sizes:   [gridx_221 | 277] x [gridy_221 | 349]
Coordinates:
Number Of Attributes: 15
  lv_ISBL3 :    600
  sub_center :  The North American Regional Reanalysis (NARR) Project
  center :      US National Weather Service - NCEP (WMC)
  long_name :   Temperature
  units :       K
  _FillValue :  1e+20
  coordinates : gridlat_221 gridlon_221
  level_indicator :     100
  grid_number : 221
  parameter_table_version :     131
  parameter_number :    11
  model :       North American Regional Reanalysis (NARR)
  forecast_time :       0
  forecast_time_units : hours
  initial_time :        01/28/2014 (00:00)

Variable: u_wind
Type: float
Total Size: 386692 bytes
            96673 values
Number of Dimensions: 2
Dimensions and sizes:   [gridx_221 | 277] x [gridy_221 | 349]
Coordinates:
Number Of Attributes: 15
  lv_ISBL3 :    600
  sub_center :  The North American Regional Reanalysis (NARR) Project
  center :      US National Weather Service - NCEP (WMC)
  long_name :   u-component of wind
  units :       m/s
  _FillValue :  1e+20
  coordinates : gridlat_221 gridlon_221
  level_indicator :     100
  grid_number : 221
  parameter_table_version :     131
  parameter_number :    33
  model :       North American Regional Reanalysis (NARR)
  forecast_time :       0
  forecast_time_units : hours
  initial_time :        01/28/2014 (00:00)

Variable: v_wind
Type: float
Total Size: 386692 bytes
            96673 values
Number of Dimensions: 2
Dimensions and sizes:   [gridx_221 | 277] x [gridy_221 | 349]
Coordinates:
Number Of Attributes: 15
  lv_ISBL3 :    600
  sub_center :  The North American Regional Reanalysis (NARR) Project
  center :      US National Weather Service - NCEP (WMC)
  long_name :   v-component of wind
  units :       m/s
  _FillValue :  1e+20
  coordinates : gridlat_221 gridlon_221
  level_indicator :     100
  grid_number : 221
  parameter_table_version :     131
  parameter_number :    34
  model :       North American Regional Reanalysis (NARR)
  forecast_time :       0
  forecast_time_units : hours
  initial_time :        01/28/2014 (00:00)

Variable: hgt
Type: float
Total Size: 386692 bytes
            96673 values
Number of Dimensions: 2
Dimensions and sizes:   [gridx_221 | 277] x [gridy_221 | 349]
Coordinates:
Number Of Attributes: 15
  lv_ISBL3 :    600
  sub_center :  The North American Regional Reanalysis (NARR) Project
  center :      US National Weather Service - NCEP (WMC)
  long_name :   Geopotential height
  units :       gpm
  _FillValue :  1e+20
  coordinates : gridlat_221 gridlon_221
  level_indicator :     100
  grid_number : 221
  parameter_table_version :     131
  parameter_number :    7
  model :       North American Regional Reanalysis (NARR)
  forecast_time :       0
  forecast_time_units : hours
  initial_time :        01/28/2014 (00:00)

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/csm/contributed.ncl"

begin
a = addfile("./narr-a_221_20140128_0000_000.grb", "r")

wks_type = "png"
wks_type at wkWidth = 2500
wks_type at wkHeight = 2500

wks= gsn_open_wks(wks_type, "temp_narr")
gsn_define_colormap(wks,"matlab_jet")

lat2d = a->gridlat_221
lon2d = a->gridlon_221
hgt = a->HGT_221_ISBL({600},:,:)
temp = a->TMP_221_ISBL({600},:,:)
pres = a->PRES_221_HYBL
u_wind = a->U_GRD_221_ISBL({600},:,:)
v_wind = a->V_GRD_221_ISBL({600},:,:)
lev    = a->lv_ISBL3            ; 29 levels

printVarSummary(temp)
printVarSummary(u_wind)
printVarSummary(v_wind)
printVarSummary(hgt)

spd = (u_wind*u_wind + v_wind*v_wind)^(0.5) ; m/s

pres = pres*0.01  ; converts to hpa

u_wind = u_wind*1.94384449
v_wind = v_wind*1.94384449

temp1 = temp - 273

temp1 at lat2d = lat2d
temp1 at lon2d = lon2d

hgt at lat2d = lat2d
hgt at lon2d = lon2d

;pres at lon2d = lon2d
;pres at lat2d = lat2d

u_wind at lon2d = lon2d
u_wind at lat2d = lat2d

spd = spd*1.94384449

spd at lon2d = lon2d
spd at lat2d = lat2d


res = True
res at gsnDraw = False
res at gsnFrame = False
res at cnLineLabelsOn = False
res at cnFillOn = True
res at tiMainString = "Temperature (C) at 600hpa on Jan 28 at 00:00UTC"
res at tiMainFont = "helvetica"
res at tiMainFontHeightF = 0.02
res at pmTickMarkDisplayMode = "Always"
res at mpProjection = "CylindricalEquidistant"
res at mpDataBaseVersion = "MediumRes"
res at mpOutlineOn = "True"
res at lbOrientation = "Vertical"
res at tiMainOffsetYF = -0.03
res at mpFillOn  = False
res at mpOutlineOn = True
;res at tfDoNDCOverlay = True
res at mpMinLatF = 37.85
res at mpMaxLatF = 38.50
res at mpMinLonF = -120.0
res at mpMaxLonF = -119.0
res at gsnLeftString = "Geopotential Height (m) and Winds (kts)"
res at gsnStringFontHeightF = 0.020
res at gsnRightString = ""
res at gsnMaximize = True
res at mpShapeMode = "FreeAspect"
res at lbTitleString = "Temperature (C)"
res at lbTitleDirection = "Across"
res at lbTitleAngleF = 90
res at lbTitlePosition = "Right"
res at lbTitleFontHeightF = 0.020
res at gsnAddCyclic = False
;res at cnLevelSelectionMode = "ManualLevels"
;res at cnMinLevelValF =
;res at cnMaxLevelValF =
;res at cnLevelSpacingF =

res1 = True
res1 at gsnDraw = False
res1 at gsnFrame = False
res1 at vcWindBarbLineThicknessF = 7.0
res1 at vcRefLengthF =0.018
res1 at vcRefMagnitudeF = 10
;res1 at tfDoNDCOverlay = True
res1 at vcMinDistanceF = 0.017
res1 at vcGlyphStyle = "WindBarb"
res1 at gsnRightString =""
res1 at gsnLeftString = ""
res1 at vcRefAnnoOn = False
res1 at gsnAddCyclic = False

res2 = True
res2 at gsnDraw = False
res2 at gsnFrame = False
res2 at cnLineColor = "Black"
res2 at gsnContourLineThicknessesScale = 5.0
res2 at gsnAddCyclic = False
;res2 at tfDoNDCOverlay = True
res2 at gsnLeftString = ""
res2 at gsnRightString = ""
res2 at cnInfoLabelOn = False
;res2 at cnLevelSpacingF = 2
res2 at cnLineLabelBackgroundColor = "white"
res2 at cnLineLabelDensityF = 1

geo = gsn_csm_contour(wks, hgt, res2)

t600 = gsn_csm_contour_map(wks, temp1, res)

vector = gsn_csm_vector(wks, u_wind, v_wind, res1)

overlay(t600, geo)
overlay(t600, vector)

draw(t600)
frame(wks)

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

end
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160802/c5a9e509/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: temp_narr.png
Type: image/png
Size: 199471 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160802/c5a9e509/attachment.png 


More information about the ncl-talk mailing list