<div dir="ltr"><div class="gmail_default" style="font-family:"comic sans ms",sans-serif"><font color="#000000">Hello,</font></div><div class="gmail_default" style="font-family:"comic sans ms",sans-serif"><font color="#000000"><br></font></div><div class="gmail_default" style="font-family:"comic sans ms",sans-serif"><font color="#000000"> I have seen scripts regarding my issues and also tried to correct my scripts using concepts of that scripts. I am new to NCL and have tried every possible measure to plot Hovmueller diagram but still, I'm stuck into it.</font><span style="color:rgb(7,55,99)"></span></div><div class="gmail_default" style="font-family:"comic sans ms",sans-serif"><font color="#000000">It's my humble request to NCL committee that please have a look at my script and help me with my issues in plotting vector Hovmueller diagram.</font></div><div class="gmail_default" style="font-family:"comic sans ms",sans-serif"><font color="#000000">I am sending you the script upon which I have been working from last one month.</font></div><div class="gmail_default" style="font-family:"comic sans ms",sans-serif"><font color="#000000"><br></font></div><div class="gmail_default" style="font-family:"comic sans ms",sans-serif"><span style="color:rgb(0,0,255)"><br></span></div><div class="gmail_default" style="font-family:"comic sans ms",sans-serif"><span style="color:rgb(0,0,255)">begin</span><br></div><div class="gmail_default"><font face="comic sans ms, sans-serif" color="#0000ff"><div class="gmail_default"><br></div><div class="gmail_default">;---File handling</div><div class="gmail_default">   <span style="white-space:pre">     </span></div><div class="gmail_default">   <span style="white-space:pre"> </span>fn  = "<a href="http://omega.mon.mean.nc">omega.mon.mean.nc</a>" ; define filename</div><div class="gmail_default">   <span style="white-space:pre">      </span>omega  = addfile(fn,"r")<span style="white-space:pre">                  </span>   <span style="white-space:pre">        </span>; open netcdf file</div><div class="gmail_default"><span style="white-space:pre">    </span>fn  = "<a href="http://uwnd.mon.mean.nc">uwnd.mon.mean.nc</a>" ; define filename</div><div class="gmail_default">   <span style="white-space:pre">        </span>uwnd  = addfile(fn,"r")                         <span style="white-space:pre">     </span>; open netcdf file</div><div class="gmail_default"><br></div><div class="gmail_default">;---Read needed variables from file</div><div class="gmail_default">   </div><div class="gmail_default">   <span style="white-space:pre">        </span>U = uwnd->uwnd</div><div class="gmail_default">   <span style="white-space:pre">        </span>W = omega->omega</div><div class="gmail_default">   <span style="white-space:pre">      </span>P0mb = 1000.</div><div class="gmail_default"><span style="white-space:pre">  </span>hyam = uwnd->uwnd                              <span style="white-space:pre">            </span>; get a coefficiants</div><div class="gmail_default">   <span style="white-space:pre">     </span>hybm = omega->omega                              <span style="white-space:pre">  </span>; get b coefficiants</div><div class="gmail_default"><br></div><div class="gmail_default">;---Define other arguments required</div><div class="gmail_default">   interp = 2 </div><div class="gmail_default">   pnew   = (/ 1000,900,850,800,750,700,550,500,250,200,150,100/)</div><div class="gmail_default">   pnew@units = "mb"          </div><div class="gmail_default">;</div><div class="gmail_default">; Omega is significantly smaller than u, so we will scale it so that some vertical motion is visible</div><div class="gmail_default">;</div><div class="gmail_default"> <span style="white-space:pre">       </span></div><div class="gmail_default"><span style="white-space:pre">      </span>wAve   = avg(W(0,:,:,{170}))           <span style="white-space:pre">                      </span>; used for scaling</div><div class="gmail_default"> <span style="white-space:pre">  </span>uAve   = avg(U(0,:,:,{170}))</div><div class="gmail_default"> <span style="white-space:pre">      </span>scale  = fabs(uAve/wAve)</div><div class="gmail_default"> <span style="white-space:pre">   </span>wscale = W*scale                       <span style="white-space:pre">                  </span>; now scale</div><div class="gmail_default"><br></div><div class="gmail_default"> <span style="white-space:pre">        </span>copy_VarCoords(W, wscale)              <span style="white-space:pre">                       </span>; copy coordinate variables</div><div class="gmail_default"><br></div><div class="gmail_default"><br></div><div class="gmail_default">;*************************************************</div><div class="gmail_default">; calculate divergence: Use Wrap to include meta data</div><div class="gmail_default">;*************************************************</div><div class="gmail_default"><br></div><div class="gmail_default"><span style="white-space:pre">  </span>f    = addfile ("<a href="http://uwnd.mon.mean.nc">uwnd.mon.mean.nc</a>", "r")</div><div class="gmail_default">  <span style="white-space:pre"> </span>u    = f->uwnd</div><div class="gmail_default"><span style="white-space:pre">   </span>f    = addfile ("<a href="http://vwnd.mon.mean.nc">vwnd.mon.mean.nc</a>", "r")</div><div class="gmail_default">  <span style="white-space:pre"> </span>v    = f->vwnd</div><div class="gmail_default"><br></div><div class="gmail_default">  <span style="white-space:pre">       </span>div = uv2dvG_Wrap(u,v)                ; u,v ==> divergence</div><div class="gmail_default"><br></div><div class="gmail_default">;*************************************************</div><div class="gmail_default">; calculate divergent wind components </div><div class="gmail_default">;*************************************************    </div><div class="gmail_default">  <span style="white-space:pre">      </span>ud    = new ( dimsizes(u), typeof(u), u@_FillValue )</div><div class="gmail_default"><span style="white-space:pre">        </span>vd    = new ( dimsizes(v), typeof(v), v@_FillValue )</div><div class="gmail_default"><br></div><div class="gmail_default">  <span style="white-space:pre">    </span>dv2uvg(div,ud,vd) ; div  ==> divergent  wind components</div><div class="gmail_default"><br></div><div class="gmail_default">  <span style="white-space:pre">      </span>copy_VarCoords(u, ud ) </div><div class="gmail_default">  <span style="white-space:pre">   </span>copy_VarCoords(u, vd ) </div><div class="gmail_default">  <span style="white-space:pre">   </span>ud@long_name  = "Zonal Divergent Wind"</div><div class="gmail_default">  <span style="white-space:pre">  </span>ud@units      = u@units</div><div class="gmail_default">  <span style="white-space:pre"> </span>vd@long_name  = "Meridional Divergent Wind"</div><div class="gmail_default">  <span style="white-space:pre">     </span>vd@units      = v@units </div><div class="gmail_default"><br></div><div class="gmail_default">;---Create plot</div><div class="gmail_default"> </div><div class="gmail_default"><span style="white-space:pre">  </span>wks   = gsn_open_wks ("png", "Vector" )        <span style="white-space:pre">            </span>; send graphics to PNG file</div><div class="gmail_default"><br></div><div class="gmail_default"> <span style="white-space:pre">        </span>res                 = True                     <span style="white-space:pre">          </span>; plot mods desired</div><div class="gmail_default"> <span style="white-space:pre"> </span>res@tiMainString    = "Vector Hovmoller" <span style="white-space:pre">                </span>; title</div><div class="gmail_default"><br></div><div class="gmail_default"> <span style="white-space:pre">    </span>res@cnLineLabelsOn  = True                <span style="white-space:pre">          </span>; turn off line labels</div><div class="gmail_default"> <span style="white-space:pre">      </span>res@cnFillOn        = True                 <span style="white-space:pre">             </span>; turn on color fill</div><div class="gmail_default"> <span style="white-space:pre">        </span>res@cnFillPalette   = "GMT_polar"  <span style="white-space:pre">                     </span>; choose color map</div><div class="gmail_default"><br></div><div class="gmail_default"> <span style="white-space:pre"> </span>res@lbLabelStride   = 2                    <span style="white-space:pre">              </span>; every other color</div><div class="gmail_default"><br></div><div class="gmail_default"> <span style="white-space:pre">        </span>res@vcRefMagnitudeF = 3.0                <span style="white-space:pre">            </span>; define vector ref mag</div><div class="gmail_default"> <span style="white-space:pre">     </span>res@vcRefLengthF    = 0.045              <span style="white-space:pre">           </span>; define length of vec ref</div><div class="gmail_default"> <span style="white-space:pre">  </span>res@vcGlyphStyle    = "CurlyVector"      <span style="white-space:pre">             </span>; turn on curly vectors</div><div class="gmail_default"> <span style="white-space:pre">     </span>res@vcMinDistanceF  = 0.01               <span style="white-space:pre">           </span>; thin out vectors</div><div class="gmail_default"> <span style="white-space:pre">  </span>res@vcMapDirection  = True</div><div class="gmail_default"><br></div><div class="gmail_default"><span style="white-space:pre">  </span>res                 = True</div><div class="gmail_default">  <span style="white-space:pre">        </span>res@vcRefMagnitudeF = 4.                    ; make vectors larger</div><div class="gmail_default">  <span style="white-space:pre">        </span>res@vcRefLengthF    = 0.050                 ; reference vector length</div><div class="gmail_default">  <span style="white-space:pre">   </span>res@vcGlyphStyle    = "CurlyVector"         ; turn on curly vectors</div><div class="gmail_default">  <span style="white-space:pre">       </span>res@vcMinDistanceF  = 0.012                 ; thin the vectors</div><div class="gmail_default"> </div><div class="gmail_default"><br></div><div class="gmail_default"><br></div><div class="gmail_default">;---Draw plot from pole to pole at 180E</div><div class="gmail_default"> </div><div class="gmail_default"><span style="white-space:pre">    </span>plot    = gsn_csm_pres_hgt_vector(wks,U(0,:,:,{170}),W(0,:,:,{170}),wscale(0,:,:,{170}),res )  </div><div class="gmail_default"><span style="white-space:pre">   </span>plot_ov = gsn_csm_vector_map(wks,ud(0,1,:,:),vd(0,1,:,:),res)</div><div class="gmail_default"><br></div><div class="gmail_default"><span style="white-space:pre">        </span>overlay(plot,plot_ov)                       ; overlay the U-wind plot on the hovmoller plot</div><div class="gmail_default">  <span style="white-space:pre">    </span>draw(plot)                                  ; draw the temperature plot (with the U-wind plot overlaid)</div><div class="gmail_default">  <span style="white-space:pre">   </span>frame(wks)                                  ; advance the frame</div><div class="gmail_default"><br></div><div class="gmail_default">end</div><div class="gmail_default"><br></div><div class="gmail_default"><br></div><div class="gmail_default"><div class="gmail_default" style="color:rgb(34,34,34)"><font color="#000000">Thank you,</font></div><div class="gmail_default" style="color:rgb(34,34,34)"><font color="#000000">Ankita Sharma</font></div></div></font></div></div>