<table style="width: 99.8%; "><tbody><tr><td id="QQMAILSTATIONERY" style="background:url(https://res.mail.qq.com/zh_CN/htmledition/images/xinzhi/bg/a_04.jpg) repeat-x #cdede2; min-height:550px; padding: 100px 55px 200px "><div><span style="font-family: 'lucida Grande', Verdana, 'Microsoft YaHei'; font-size: 14px; background-color: rgb(255, 255, 255);"><font size="4" style="line-height: 30px;">Hi,all:</font></span><div style="font-family: 'lucida Grande', Verdana, 'Microsoft YaHei'; font-size: 14px; line-height: 23px; background-color: rgb(255, 255, 255);"><span style="line-height: 24px;"><font size="4" style="line-height: 30px;">&nbsp; &nbsp; I am calculating a index and try to plot it.</font></span></div><div style="background-color: rgb(255, 255, 255);"><font face="lucida Grande, Verdana, Microsoft YaHei" size="4"><span style="line-height: 30px;">&nbsp; &nbsp; &nbsp;lindex=abs(sreh)*(qgsum+qisum+qssum)<span id="_editor_bookmark_start_1" style="display: none; line-height: 0px;">‍</span></span></font></div><div style="font-family: 'lucida Grande', Verdana, 'Microsoft YaHei'; font-size: 14px; line-height: 23px; background-color: rgb(255, 255, 255);"><span style="line-height: 24px;"><font size="4" style="line-height: 30px;">&nbsp; &nbsp; but after I write the script according the function above and try to run it,it just calculated the sreh and has no errors warning and wanted pictures appeared.</font></span></div><div style="font-family: 'lucida Grande', Verdana, 'Microsoft YaHei'; font-size: 14px; line-height: 23px; background-color: rgb(255, 255, 255);"><span style="line-height: 24px;"><font size="4" style="line-height: 30px;">The bash window and script:</font></span></div><div style="font-family: 'lucida Grande', Verdana, 'Microsoft YaHei'; font-size: 14px; line-height: 23px; background-color: rgb(255, 255, 255);"><br></div><div style="font-family: 'lucida Grande', Verdana, 'Microsoft YaHei'; font-size: 14px; line-height: 23px; background-color: rgb(255, 255, 255);"><img src="cid:D21A0106@AB695143.71342154.jpg" filesize="75277"></div><div style="background-color: rgb(255, 255, 255);"><font face="lucida Grande, Verdana, Microsoft YaHei" size="4"><span style="line-height: 30px;"><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"&nbsp;</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"&nbsp;</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div><br></div><div>begin</div><div>;</div><div>; The WRF ARW input file. &nbsp;</div><div>; This needs to have a ".nc" appended, so just do it.</div><div>&nbsp; a = addfile("/home/Huanglei/data/d032"+".nc","r") &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;</div><div><br></div><div><br></div><div>; We generate plots, but what kind do we prefer?</div><div>&nbsp; type = "pdf"</div><div>; type = "pdf"</div><div>; type = "ps"</div><div>; type = "ncgm"</div><div>&nbsp; wks = gsn_open_wks(type,"plt_lightningindex")</div><div><br></div><div>&nbsp; gsn_define_colormap(wks,"precip_11lev")</div><div>; Set some basic resources</div><div>&nbsp; res = True</div><div>&nbsp; res@MainTitle = "REAL-TIME WRF"</div><div>&nbsp;; res@gsnDraw &nbsp; &nbsp; &nbsp;= &nbsp;False &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;</div><div>&nbsp; ;res@gsnFrame &nbsp; &nbsp; = &nbsp;False</div><div><br></div><div>&nbsp; mpres &nbsp;= True &nbsp;; Map resources</div><div>&nbsp; mpres@mpOutlineOn = False &nbsp;; Turn off map outlines</div><div>&nbsp; mpres@mpFillOn &nbsp; &nbsp;= False &nbsp;; Turn off map fill</div><div>&nbsp; mpres@mpGridAndLimbOn = True</div><div>&nbsp;;res@mpProjection &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= "Lambert"</div><div>&nbsp; pltres = True ; Plot resources</div><div>&nbsp; pltres@PanelPlot &nbsp;= True &nbsp; ; Tells wrf_map_overlays not to remove overlays</div><div><br></div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>; What times and how many time steps are in the data set?</div><div>&nbsp; times = wrf_user_getvar(a,"times",-1) &nbsp;; get all times in the file</div><div>&nbsp; ntimes = dimsizes(times) &nbsp; &nbsp; &nbsp; &nbsp; ; number of times in the file</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>&nbsp; do it =1, ntimes-1,1 &nbsp; &nbsp; &nbsp; &nbsp;; TIME LOOP</div><div><br></div><div>&nbsp; &nbsp; print("Working on time: " + times(it) )</div><div>&nbsp; &nbsp; res@TimeLabel = times(it) &nbsp; ; Set Valid time to use on plots</div><div>&nbsp; &nbsp; sreh = wrf_user_getvar(a,(/"helicity","3000"/),it) &nbsp; &nbsp;<span class="Apple-tab-span" style="white-space:pre">        </span>; here 3km is specifically set - same as above</div><div><br></div><div>&nbsp; &nbsp; cnres &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= res</div><div>&nbsp; &nbsp; cnres@cnFillOn &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = True</div><div>&nbsp; &nbsp; cnres@cnSmoothingOn &nbsp; &nbsp; &nbsp; &nbsp;= True</div><div>&nbsp; &nbsp; cnres@cnSmoothingDistanceF = .005</div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div>; First get the variables we will need &nbsp; &nbsp; &nbsp; &nbsp;</div><div>&nbsp; &nbsp; if(isfilevar(a,"QGRAUP"))</div><div><span class="Apple-tab-span" style="white-space:pre">        </span> &nbsp;qgraup = wrf_user_getvar(a,"QGRAUP",it )</div><div>&nbsp; &nbsp; &nbsp; qgraup = qgraup*1000.</div><div>&nbsp; &nbsp; &nbsp; qgraup@units = "g/kg" &nbsp;&nbsp;</div><div>&nbsp; &nbsp; end if</div><div><br></div><div>&nbsp;if(isfilevar(a,"QICE"))</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre">        </span> &nbsp;qice = wrf_user_getvar(a,"QICE",it )</div><div>&nbsp; &nbsp; &nbsp; qice = qice*1000.</div><div>&nbsp; &nbsp; &nbsp; qice@units = "g/kg" &nbsp;&nbsp;</div><div>&nbsp; &nbsp; end if</div><div><br></div><div>&nbsp;if(isfilevar(a,"QSNOW"))</div><div><span class="Apple-tab-span" style="white-space:pre">        </span> &nbsp;qsnow = wrf_user_getvar(a,"QSNOW",it )</div><div>&nbsp; &nbsp; &nbsp; qsnow = qsnow*1000.</div><div>&nbsp; &nbsp; &nbsp; qsnow@units = "g/kg" &nbsp;&nbsp;</div><div>&nbsp; &nbsp; end if</div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>&nbsp; &nbsp;</div><div>&nbsp; &nbsp; &nbsp; opts = res</div><div>&nbsp; &nbsp; &nbsp; opts@cnFillOn &nbsp; &nbsp; &nbsp; &nbsp; = True</div><div>&nbsp; &nbsp; &nbsp; opts@gsnSpreadColors &nbsp;= False</div><div>&nbsp; &nbsp; &nbsp; opts@ContourParameters &nbsp; &nbsp; &nbsp; = (/ 1, 19, 2 /)</div><div><span class="Apple-tab-span" style="white-space:pre">        </span> &nbsp;opts@gsnDraw &nbsp; &nbsp; &nbsp;= &nbsp;False &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;</div><div>&nbsp; &nbsp; &nbsp; opts@gsnFrame &nbsp; &nbsp; = &nbsp;False</div><div><br></div><div><br></div><div>&nbsp; &nbsp; &nbsp; if (isvar("qgraup"))</div><div>&nbsp; &nbsp; &nbsp; &nbsp; qgs &nbsp;= qgraup(13 + (/0,1,2,3,4,5,6/),:,:)</div><div>&nbsp; &nbsp; &nbsp; &nbsp; qgsum = dim_sum_n_Wrap(qgs, 0)</div><div>&nbsp; &nbsp; &nbsp; end if</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre">        </span> &nbsp; if (isvar("qice"))</div><div>&nbsp; &nbsp; &nbsp; &nbsp; qis &nbsp;= qice(13 + (/0,1,2,3,4,5,6/),:,:)</div><div>&nbsp; &nbsp; &nbsp; &nbsp; qisum = dim_sum_n_Wrap(qis, 0)</div><div>&nbsp; &nbsp; &nbsp; end if</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre">        </span> &nbsp; if (isvar("qsnow"))</div><div>&nbsp; &nbsp; &nbsp; &nbsp; qss &nbsp;= qsnow(13 + (/0,1,2,3,4,5,6/),:,:)</div><div>&nbsp; &nbsp; &nbsp; &nbsp; qssum = dim_sum_n_Wrap(qss, 0)</div><div>&nbsp; &nbsp; &nbsp; end if</div><div><br></div><div>&nbsp; &nbsp;lindex=abs(sreh)*(qgsum+qisum+qssum)</div><div>&nbsp; &nbsp;contour = wrf_contour(a,wks,lindex,opts)</div><div>&nbsp; &nbsp;plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)</div><div>;&gt;============================================================&lt;</div><div>; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;add China map</div><div>;&gt;------------------------------------------------------------&lt;</div><div><br></div><div>&nbsp; &nbsp; &nbsp;</div><div>&nbsp; shp_name1 &nbsp; &nbsp;= "/home/Huanglei/map/China/diquJie_polyline.shp"</div><div><br></div><div>&nbsp; lnres &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= True</div><div>&nbsp; lnres@gsLineColor &nbsp; &nbsp; &nbsp;= "gray25"</div><div>&nbsp; lnres@gsLineThicknessF = 0.5 &nbsp;&nbsp;</div><div><br></div><div>&nbsp;id = gsn_add_shapefile_polylines(wks,plot,shp_name1,lnres)</div><div>&nbsp; shp_name2 &nbsp; &nbsp;= "/home/Huanglei/map/China/cnmap/cnhimap.shp"</div><div><br></div><div>&nbsp; prres=True</div><div>&nbsp; prres@gsLineThicknessF = 2.0 &nbsp; &nbsp; &nbsp;&nbsp;</div><div>&nbsp; prres@gsLineColor = "black"</div><div>&nbsp; plotcn3 = gsn_add_shapefile_polylines(wks,plot,shp_name2,prres)</div><div><br></div><div>&nbsp; &nbsp;txres2 &nbsp;= True</div><div>&nbsp; &nbsp;txres2@txFont &nbsp;= 10</div><div>&nbsp; &nbsp;txres2@txFontHeightF =0.01</div><div>&nbsp; &nbsp;txres2@txFontColor = "Blue"</div><div>&nbsp; &nbsp;txdum1 =gsn_add_text(wks, plot, "Chengdu", 104.06,30.67, txres2)</div><div><br></div><div>&nbsp; draw(plot) &nbsp; &nbsp; &nbsp; ; This will draw the map and the shapefile outlines.</div><div>&nbsp; frame(wks)</div><div>&nbsp; &nbsp;delete(opts)</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>&nbsp; end do &nbsp; &nbsp; &nbsp; &nbsp;; END OF TIME LOOP &nbsp; &nbsp;&nbsp;</div><div><br></div><div>end<span id="_editor_bookmark_start_2" style="display: none; line-height: 0px;">‍</span></div></span></font></div><div style="font-family: 'lucida Grande', Verdana, 'Microsoft YaHei'; font-size: 14px; line-height: 23px; background-color: rgb(255, 255, 255);"><span style="line-height: 24px;"><font size="4" style="line-height: 30px;"><br></font></span></div><div style="font-family: 'lucida Grande', Verdana, 'Microsoft YaHei'; font-size: 14px; line-height: 23px; background-color: rgb(255, 255, 255);"><span style="line-height: 24px;"><font size="4" style="line-height: 30px;">&nbsp; &nbsp; &nbsp;Do you guys have some advice?</font></span></div><div style="font-family: 'lucida Grande', Verdana, 'Microsoft YaHei'; font-size: 14px; line-height: 23px; background-color: rgb(255, 255, 255);"><span style="line-height: 24px;"><font size="4" style="line-height: 30px;">&nbsp; &nbsp; &nbsp;any information will be appreciated</font></span></div></div><div><span style="line-height: 24px;"><font size="4" style="line-height: 30px;"><span id="_editor_bookmark_start_0" style="display: none; line-height: 0px;">‍</span><br></font></span></div></td></tr></tbody></table>