<div dir="ltr"><br>[0] <br>You can advect <b>any </b>quantity. I thought it is best to advect relative vorticity and then add the coriolis force to get the 'new' absolute vorticity. However, you could advect absolute vorticity for you needs. Good question for your advisor.  :-)<br><br><div>[1]<br></div><div>You are starting with a WRF grid. The 'advect_variable' function requires a rectilinear grid. Hence, you need to generate a weight file that takes data on a WRF grid to a rectilinear grid. Any name can be used for the weight file. This is fine:  <br>     "WRF_to_Rect.WgtFile_"+InterpMethod+"."+nlat+"_x_"+nlon+".nc"<br></div><div><br>[2]<br>Now you want to take the advected vorticity on a rectilinear grid and interpolate it to the original WRF grid. The following name is fine.<br><br>    Opt@WgtFileName     = "Rect_to_WRF.WgtFile_"+InterpMethod+"."+nlat+"_x_"+nlon+".nc"<br><br>---<br>You wrote: "I'm still obtaining errors." However, you did not state what the errors are.<br></div><div> I think the following  <br><br>    adv_wrf                         =  ESMF_regrid(adv_vrt,<span style="color:rgb(255,0,0)"><b> Opt)<br><br></b></span></div><div><span style="color:rgb(255,0,0)"><b>should be<br><br></b></span></div><div>    adv_wrf                         =  ESMF_regrid(adv_vrt,<span style="color:rgb(255,0,0)"><b> Opt@WgtFileName</b></span>)<br><br>---<br></div><div>You did not "ftp a WRF file"<br><br></div><div>Good Luck<br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 23, 2018 at 4:30 AM, Lara Quitián Hernández <span dir="ltr"><<a href="mailto:laraquitianhernandez@gmail.com" target="_blank">laraquitianhernandez@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thank you so much. I was absolutely wrong and I did not realize my mistake.<div><br></div><div>I attach here my script because I'm still obtaining errors and I'm afraid I do not fully understand the procedure of the different calculations. For example, I do not know if I have to generate a weight file for regridding the relative vorticity (wr) or which name give to that weight file (Opt@WgtFileName). Moreover, I do not know why do I have to use the relative vorticity and not the absolute vorticity to study the vorticity advection at the end. </div><div><br></div><div>On the other hand, I will also ftp a WRF file and my script.</div><span class=""><div><br></div><div>Thanks in advance.</div><div><br></div><div>Regards, </div><div>Lara</div><div><br></div></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">2018-05-22 17:40 GMT+02:00 Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Lara,<div>You are passing in your original data into advect_variable instead of the regridded variables. Change this:</div><span><div><span style="color:rgb(39,78,19);font-family:arial,sans-serif;font-size:12.8px;font-style:italic;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"> adv_vrt = advect_variable(u, v, wr, 1, "vorticity advection", "(m/s)(s-1)",0) </span><br></div></span><div><span style="color:rgb(39,78,19);font-family:arial,sans-serif;font-size:12.8px;font-style:italic;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">to this:</span></div><div><span style="color:rgb(39,78,19);font-family:arial,sans-serif;font-size:12.8px;font-style:italic;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="color:rgb(39,78,19);font-family:arial,sans-serif;font-size:12.8px;font-style:italic;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"> adv_vrt = advect_variable(u_reclin, v_reclin, wr_reclin, 1, "vorticity advection", "(m/s)(s-1)",0) </span><br></span></div><div><span style="color:rgb(39,78,19);font-family:arial,sans-serif;font-size:12.8px;font-style:italic;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span style="color:rgb(39,78,19);font-family:arial,sans-serif;font-size:12.8px;font-style:italic;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Gaussian grids are indeed rectlinear grids. See the definition here:</span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><a href="https://www.ncl.ucar.edu/Document/glossary.shtml#RectilinearGrid" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/glossary.shtml#Rectilinear<wbr>Grid</a> </span><br></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">If you have any further questions please respond to ncl-talk.</span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Adam</span></div></div><div class="gmail_extra"><div><div class="m_7015166853841739983h5"><br><div class="gmail_quote">On Tue, May 22, 2018 at 7:24 AM, Lara Quitián Hernández <span dir="ltr"><<a href="mailto:laraquitianhernandez@gmail.com" target="_blank">laraquitianhernandez@gmail.co<wbr>m</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">First of all, I apologize for asking you so many questions and for this long email.<div><br></div><div>I'm trying to interpolate the WRF u,v, relative vorticity to a rectilinear grid and use the <b>advect_variable</b> on the rectilinear grid. However, I have obtained this warning and I do not know how to fixed it: <b><font color="#660000">advect_var: Warning: No named dimensions. Can not check grid order.</font></b></div><div><br></div><div>After defining the diverse variables, I have performed the regrid in this way:</div><div><br></div><div><div><i><font color="#274e13">;=============================<wbr>==============================<wbr>=====================</font></i></div><div><i><font color="#274e13">; Part A: WRF (curvilinear) to Rectilinear regrid variable and generate  weights </font></i></div><div><i><font color="#274e13">;=============================<wbr>==============================<wbr>=====================</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">;Interpolate the WRF u, v and relative vorticity to a rectilinear grid with as fine [or finer] resoultion as the WRF grid</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">;---Set up options for regridding</font></i></div><div><i><font color="#274e13">    </font></i></div><div><i><font color="#274e13">    InterpMethod= "bilinear"                    ; define interpolation method</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">;---Title and filename options</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">    Opt                = True</font></i></div><div><i><font color="#274e13">    Opt@SrcTitle       = "WRF grid"   ; optional</font></i></div><div><i><font color="#274e13">    Opt@WgtFileName    = "WRF_to_Rect.WgtFile_"+InterpM<wbr>ethod+"."+nlat+"_x_"+nlon+".nc<wbr>"</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">;---Source grid options</font></i></div><div><i><font color="#274e13">    Opt@SrcFileName    = "WRF.SCRIP_grid_description."+<wbr>nlat+"_x_"+nlon+".nc"  ; Name of source and</font></i></div><div><i><font color="#274e13">    Opt@SrcRegional    = True</font></i></div><div><i><font color="#274e13">    Opt@SrcGridLat     = lat2d</font></i></div><div><i><font color="#274e13">    Opt@SrcGridLon     = lon2d</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">;---Destination grid options</font></i></div><div><i><font color="#274e13">    Opt@DstFileName    = "Rectilinear.SCRIP_grid_descri<wbr>ption."+nlat+"_x_"+nlon+".nc" ; destination files</font></i></div><div><i><font color="#274e13">    Opt@DstGridType    = "rectilinear"</font></i></div><div><i><font color="#274e13">    Opt@DstGridLat     = lat</font></i></div><div><i><font color="#274e13">    Opt@DstGridLon     = lon</font></i></div><div><i><font color="#274e13">    Opt@DstRegional    = True</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">;---Specify other options  </font></i></div><div><i><font color="#274e13">    Opt@ForceOverwrite = True</font></i></div><div><i><font color="#274e13">    Opt@InterpMethod   = InterpMethod</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">;---Perform the regrid: WRF ==> rectilinear (_reclin)</font></i></div><div><i><font color="#274e13">    u_reclin = ESMF_regrid(u, Opt)     ; Do the regridding for u   </font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">;---Use the just generated weight file</font></i></div><div><i><font color="#274e13">    v_reclin = ESMF_regrid_with_weights(v,Opt<wbr>@WgtFileName,False)</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">;---Perform the regrid for wr: WRF ==> rectilinear (_reclin)</font></i></div><div><i><font color="#274e13">    wr_reclin = ESMF_regrid_with_weights(wr, Opt@WgtFileName,False)  </font></i></div></div><div><i><font color="#274e13"><br></font></i></div><div><font color="#000000">Afterwards, I use the <b>advect_variable:</b></font></div><div><font color="#000000"><b><br></b></font></div><div><div style="color:rgb(39,78,19);font-style:italic">;=============================<wbr>==============================<wbr>===========</div><div style="color:rgb(39,78,19);font-style:italic">;Use 'advect_variable' on the rectilinear grid</div><div style="color:rgb(39,78,19);font-style:italic">;=============================<wbr>==============================<wbr>===========</div><div style="color:rgb(39,78,19);font-style:italic"><br></div><div style="color:rgb(39,78,19);font-style:italic">   adv_vrt = advect_variable(u, v, wr, 1, "vorticity advection", "(m/s)(s-1)",0)   ;I use gridType=1. I suppose rectilinear and fixed grid are quite similar, or even, the same.</div><div style="color:rgb(39,78,19);font-style:italic"><br></div><div style="color:rgb(39,78,19);font-style:italic">  printVarSummary(adv_vrt) </div><div style="color:rgb(39,78,19);font-style:italic">  printMinMax(adv_vrt,0)</div><div style="color:rgb(39,78,19);font-style:italic">  print("----")</div><div style="color:rgb(39,78,19);font-style:italic"><br></div><div><font color="#000000">Finally, I obtain this results and I do not know if they are correctly calculated:</font></div></div><div><font color="#000000"><br></font></div><div><div><i><font color="#274e13">Variable: <b>u_reclin</b></font></i></div><div><i><font color="#274e13">Type: float</font></i></div><div><i><font color="#274e13">Total Size: 5972400 bytes</font></i></div><div><i><font color="#274e13">            1493100 values</font></i></div><div><i><font color="#274e13">Number of Dimensions: 3</font></i></div><div><i><font color="#274e13">Dimensions and sizes:<span style="white-space:pre-wrap">        </span>[bottom_top | 35] x [lat | 180] x [lon | 237]</font></i></div><div><i><font color="#274e13">Coordinates: </font></i></div><div><i><font color="#274e13">            lat: [25.832275390625..40.830703735<wbr>35156]</font></i></div><div><i><font color="#274e13">            lon: [-34.179443359375..-9.30551147<wbr>4609375]</font></i></div><div><i><font color="#274e13">Number Of Attributes: 8</font></i></div><div><i><font color="#274e13">  missing_value :<span style="white-space:pre-wrap">     </span>9.96921e+36</font></i></div><div><i><font color="#274e13">  remap :<span style="white-space:pre-wrap">      </span>remapped via ESMF_regrid_with_weights: Bilinear</font></i></div><div><i><font color="#274e13">  _FillValue :<span style="white-space:pre-wrap">     </span>9.96921e+36</font></i></div><div><i><font color="#274e13">  stagger :<span style="white-space:pre-wrap">    </span> </font></i></div><div><i><font color="#274e13">  units :<span style="white-space:pre-wrap">       </span>m s-1</font></i></div><div><i><font color="#274e13">  description :<span style="white-space:pre-wrap">      </span>x-wind component</font></i></div><div><i><font color="#274e13">  MemoryOrder :<span style="white-space:pre-wrap">   </span>XYZ</font></i></div><div><i><font color="#274e13">  FieldType :<span style="white-space:pre-wrap">  </span>104</font></i></div><div><i><font color="#274e13">(0)<span style="white-space:pre-wrap">     </span>u_reclin: min=-21.454   max=49.9299</font></i></div><div><i><font color="#274e13">(0)<span style="white-space:pre-wrap">   </span>----</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">Variable: <b>v_reclin</b></font></i></div><div><i><font color="#274e13">Type: float</font></i></div><div><i><font color="#274e13">Total Size: 5972400 bytes</font></i></div><div><i><font color="#274e13">            1493100 values</font></i></div><div><i><font color="#274e13">Number of Dimensions: 3</font></i></div><div><i><font color="#274e13">Dimensions and sizes:<span style="white-space:pre-wrap">  </span>[bottom_top | 35] x [lat | 180] x [lon | 237]</font></i></div><div><i><font color="#274e13">Coordinates: </font></i></div><div><i><font color="#274e13">            lat: [25.832275390625..40.830703735<wbr>35156]</font></i></div><div><i><font color="#274e13">            lon: [-34.179443359375..-9.30551147<wbr>4609375]</font></i></div><div><i><font color="#274e13">Number Of Attributes: 8</font></i></div><div><i><font color="#274e13">  _FillValue :<span style="white-space:pre-wrap">        </span>9.96921e+36</font></i></div><div><i><font color="#274e13">  stagger :<span style="white-space:pre-wrap">    </span> </font></i></div><div><i><font color="#274e13">  units :<span style="white-space:pre-wrap">       </span>m s-1</font></i></div><div><i><font color="#274e13">  description :<span style="white-space:pre-wrap">      </span>y-wind component</font></i></div><div><i><font color="#274e13">  MemoryOrder :<span style="white-space:pre-wrap">   </span>XYZ</font></i></div><div><i><font color="#274e13">  FieldType :<span style="white-space:pre-wrap">  </span>104</font></i></div><div><i><font color="#274e13">  remap :<span style="white-space:pre-wrap">      </span>remapped via ESMF_regrid_with_weights: Bilinear</font></i></div><div><i><font color="#274e13">  missing_value :<span style="white-space:pre-wrap">  </span>9.96921e+36</font></i></div><div><i><font color="#274e13">(0)<span style="white-space:pre-wrap">     </span>v_reclin: min=-51.6305   max=46.8953</font></i></div><div><i><font color="#274e13">(0)<span style="white-space:pre-wrap">  </span>----</font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">Variable: <b>wr_reclin</b></font></i></div><div><i><font color="#274e13">Type: float</font></i></div><div><i><font color="#274e13">Total Size: 5972400 bytes</font></i></div><div><i><font color="#274e13">            1493100 values</font></i></div><div><i><font color="#274e13">Number of Dimensions: 3</font></i></div><div><i><font color="#274e13">Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[35] x [lat | 180] x [lon | 237]</font></i></div><div><i><font color="#274e13">Coordinates: </font></i></div><div><i><font color="#274e13">            lat: [25.832275390625..40.830703735<wbr>35156]</font></i></div><div><i><font color="#274e13">            lon: [-34.179443359375..-9.30551147<wbr>4609375]</font></i></div><div><i><font color="#274e13">Number Of Attributes: 3</font></i></div><div><i><font color="#274e13">  _FillValue :<span style="white-space:pre-wrap">     </span>9.96921e+36</font></i></div><div><i><font color="#274e13">  remap :<span style="white-space:pre-wrap">      </span>remapped via ESMF_regrid_with_weights: Bilinear</font></i></div><div><i><font color="#274e13">  missing_value :<span style="white-space:pre-wrap">  </span>9.96921e+36</font></i></div><div><i><font color="#274e13">(0)<span style="white-space:pre-wrap">     </span>wr_reclin: min=-34.6481   max=63.9581</font></i></div><div><i><font color="#274e13">(0)<span style="white-space:pre-wrap"> </span>----</font></i></div><div><i><font color="#274e13">(0)<span style="white-space:pre-wrap">    </span>RectilinearGrid: nmsg=180075</font></i></div><div><i><font color="#274e13">(0)<span style="white-space:pre-wrap">    </span>----</font></i></div><div><i><font color="#274e13">(0)</font><span style="color:rgb(39,78,19);white-space:pre-wrap">   </span><font color="#660000"><b>advect_var: Warning: No named dimensions. Can not check grid order.</b></font></i></div><div><i><font color="#274e13"><br></font></i></div><div><i><font color="#274e13">Variable: <b>adv_vrt</b></font></i></div><div><i><font color="#274e13">Type: float</font></i></div><div><i><font color="#274e13">Total Size: 5972400 bytes</font></i></div><div><i><font color="#274e13">            1493100 values</font></i></div><div><i><font color="#274e13">Number of Dimensions: 3</font></i></div><div><i><font color="#274e13">Dimensions and sizes:<span style="white-space:pre-wrap">   </span>[35] x [180] x [237]</font></i></div><div><i><font color="#274e13">Coordinates: </font></i></div><div><i><font color="#274e13">Number Of Attributes: 2</font></i></div><div><i><font color="#274e13">  long_name :<span style="white-space:pre-wrap">  </span>vorticity advection</font></i></div><div><i><font color="#274e13">  units :<span style="white-space:pre-wrap">      </span>(m/s)(s-1)</font></i></div><div><i><font color="#274e13">(0)<span style="white-space:pre-wrap">      </span>vorticity advection ((m/s)(s-1)) : min=-0.0107484   max=0.0109576</font></i></div><div style="color:rgb(0,0,0)"><br></div></div><div><font color="#000000">Unfortunately, I also obtain this errors:</font></div><div><font color="#660000"><br></font></div><div><div><font color="#660000">(0)<span style="white-space:pre-wrap">    </span>is_SCRIP: '<a href="http://Rectilinear.SCRIP_grid_description.nc" target="_blank">Rectilinear.SCRIP_grid_descri<wbr>ption.nc</a>' doesn't exist.</font></div><div><font color="#660000">(0)<span style="white-space:pre-wrap"> </span>ESMF_regrid_gen_weights: invalid or missing SCRIP source grid file.</font></div></div><div><br></div><div>And I obtain this three files:</div><div><b><br></b></div><div><b><a href="http://Rectilinear.SCRIP_grid_description.180_x_237.nc" target="_blank">Rectilinear.SCRIP_grid_descrip<wbr>tion.180_x_237.nc</a><br></b></div><div><b><a href="http://WRF.SCRIP_grid_description.180_x_237.nc" target="_blank">WRF.SCRIP_grid_description.180<wbr>_x_237.nc</a><br></b></div><div><b><a href="http://WRF_to_Rect.WgtFile_bilinear.180_x_237.nc" target="_blank">WRF_to_Rect.WgtFile_bilinear.1<wbr>80_x_237.nc</a></b><br></div><div><b><br></b></div><div>Which file do I have to use to display the vorticity advection? None, I suppose. They are all description files, except one.</div><div> </div><div>If you need my script, just tell me. </div><div><br></div><div>Thanks in advance.</div><div><br></div><div>Regards, </div><div>Lara</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-05-22 11:08 GMT+02:00 Lara Quitián Hernández <span dir="ltr"><<a href="mailto:laraquitianhernandez@gmail.com" target="_blank">laraquitianhernandez@gmail.co<wbr>m</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Good morning, <div><br></div><div>Thank you very much. There is only one doubt left unanswered: the <b>advect_function </b>requires a <u>Gaussian, regular o fixed grid</u>, but not <u>rectilinear</u>. Is there any relation between them?</div><div><br></div><div>Regards, </div><span class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492HOEnZb"><font color="#888888"><div>Lara</div></font></span></div><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492HOEnZb"><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492h5"><div class="gmail_extra"><br><div class="gmail_quote">2018-05-21 19:21 GMT+02:00 Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>ESMF 37 was <br>(a)  interpolating u and v onto a rectilinear grid; <br>(b)  calculating divergence using uv2dv_cfd which requires u and v to be on a rectilinear grid; <br>(c)  interpolating the derived divergence on a rectilinear grid back onto the original WRF grid<br><br>---<br></div>In your case:<br><div><br>[1] Interpolate the WRF u, v and relative vorticity to a rectilinear grid with as fine [or finer] resoultion as the WRF grid<br></div><div>[2] Use 'advect_variable' on the rectilinear grid<br>[3] interpolate the advected quantity on the rectilinear grid back onto the original WRF grid<br><br>---<br></div><div>Caveat: Technically, on a sphere,  u and v must be interpolated together using the 'rotation' angle.<br></div><div>Why? because the new location has a 'slightly' different orientation relative to the pole(s) <br></div><div><br>  unew        = u*cos(rot) - v*sin(rot)  <br>  vnew        = u*sin(rot) + v*cos(rot)<br><br></div><div>The above approach, independently interpolating u and v ignores the rotation. It is assumed to be very small.<br><br> <br></div><div><br></div></div><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251HOEnZb"><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, May 21, 2018 at 3:45 AM, Lara Quitián Hernández <span dir="ltr"><<a href="mailto:laraquitianhernandez@gmail.com" target="_blank">laraquitianhernandez@gmail.co<wbr>m</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Good morning,</div><div><br></div>Thank you so much. Unfortunately, I have some doubts. <div><br></div><div>1.- I have read the information for Example 37 but I do not understand the need of computing the divergence ("uv2dv_cfd"). I just need my variables on a Gaussian, regular o fixed grid in order to apply the <b>advect_variable</b> function</div><div>2.- I need to convert my curvilinear WRF grid to a Gaussian, regular o fixed grid. However, in both Examples (29a and 37) a regrid from a WRF grid to a rectilinear grid is carried out. Are those grids (Gaussian, regular o fixed grid) also rectilinear grids?</div><div>3.-Do I need to regrid back to the WRF curvilinear grid as it is done in both examples?</div><span><div><br></div><div>Thanks in advance. </div><div><br></div><div>Regards, </div><div>Lara</div><div><br></div><div><br></div></span></div><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767HOEnZb"><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767h5"><div class="gmail_extra"><br><div class="gmail_quote">2018-05-20 1:10 GMT+02:00 Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div>NCL has no builtin function to advect a field (eg: vorticity) on a curvilinear grid (eg: your WRF Grid)..<br></div>Maybe ask wrfhelp?<br><br>===<br></div>That said<br><br><a href="https://www.ncl.ucar.edu/Applications/ESMF.shtml" target="_blank">https://www.ncl.ucar.edu/Appli<wbr>cations/ESMF.shtml</a><br></div>Example 29a<br></div><div>or, maybe better,<br></div><div>Example 37<br></div><div><br></div>outlines a method that might be used for your needs.<br><br></div>===<br></div>Good luck<br><div><br></div></div><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275HOEnZb"><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 16, 2018 at 7:40 AM, Lara Quitián Hernández <span dir="ltr"><<a href="mailto:laraquitianhernandez@gmail.com" target="_blank">laraquitianhernandez@gmail.co<wbr>m</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello, <div><br></div><div>Thank you. In the case of calculating the horizontal advection for <b>vorticity</b>, it would be better to use the direct WRF variable ( vo = wrf_user_getvar(a,"<b class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275m_-4196797236910994329m_-4307035958762947733gmail-">avo</b>",-1)) or define it via different functions: <a href="https://www.ncl.ucar.edu/Applications/vort.shtml" target="_blank">https://www.ncl.uca<wbr>r.edu/Applications/vort.shtml</a></div><div><br></div><div>The only example in the link defines the Vorticity like this:</div><div><pre style="color:rgb(0,0,0);font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial;word-wrap:break-word;white-space:pre-wrap">;*****************************<wbr>*******************
; calculate vorticity on a Gaussian Grid
; scale for visual convenience
;*****************************<wbr>*******************
  scale = 1.e05                                  
  vrt   = u                                      ; retain coordinates
  vrt   = uv2vrG_Wrap(u,v) * scale
  vrt@long_name = "vorticity"
  vrt@units     = "scaled"</pre></div><div><br></div><div>However, my simulations are generated on a Lambert Conformal grid projection, not Gaussian one. Does this change the procedure to calculate it?</div><div><br></div><div>Regards, </div><span class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275m_-4196797236910994329HOEnZb"><font color="#888888"><div>Lara</div></font></span></div><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275m_-4196797236910994329HOEnZb"><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275m_-4196797236910994329h5"><div class="gmail_extra"><br><div class="gmail_quote">2018-05-16 14:28 GMT+02:00 Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Maybe:<br><br><a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/advect_variable.shtml" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/Functions/Contributed/adve<wbr>ct_variable.shtml</a><br><br></div>'x' would be vorticity<br><br></div>Good Luck<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275m_-4196797236910994329m_-4307035958762947733h5">On Wed, May 16, 2018 at 4:01 AM, Lara Quitián Hernández <span dir="ltr"><<a href="mailto:laraquitianhernandez@gmail.com" target="_blank">laraquitianhernandez@gmail.co<wbr>m</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275m_-4196797236910994329m_-4307035958762947733h5"><div dir="ltr">Good morning, <div><br></div><div>I would like to know if there is a possibility of calculating the Vorticity Advection with NCL.</div><div><br></div><div>Thanks in advance.</div><div><br></div><div>Regards,  </div><span class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275m_-4196797236910994329m_-4307035958762947733m_2745945185124896099HOEnZb"><font color="#888888"><div>Lara</div><div><br clear="all"><div><br></div>-- <br><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275m_-4196797236910994329m_-4307035958762947733m_2745945185124896099m_-1448946543173881759gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><img src="https://docs.google.com/uc?export=download&id=0B5glieJcogGocjdGWHA3NWIwdmc&revid=0B5glieJcogGocU45NUpENzV4UWhCWk9tb1ZGdXB5MEtGYjJjPQ"><br><div><font color="#000000"><br></font></div><div><b><font color="#000000"><span style="font-size:12.8px">Lara Quitián Hernández, </span></font><span style="color:rgb(0,0,0);font-size:12.8px">PhD Student</span></b><font color="#000000"><b><br style="font-size:12.8px"></b><span style="font-size:12.8px">Departamento de Física de la Tierra, Astronomía y Astrofísica</span><br style="font-size:12.8px"><span style="font-size:12.8px">Facultad de Ciencias Físicas</span><br style="font-size:12.8px"><span style="font-size:13px"><a href="https://maps.google.com/?q=Ciudad+Universitaria,+Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">Ciudad Universitaria, Plaza Ciencias, 1</a></span></font></div><div><font color="#000000"><span style="font-size:12.8px"><a href="https://maps.google.com/?q=Ciudad+Universitaria,+Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">28040 Madrid. Spain</a></span><br style="font-size:12.8px"><span style="font-size:12.8px">Email: </span><a href="mailto:mlunar@aemet.es" style="font-size:12.8px" target="_blank">lquitian@ucm.es</a><br style="font-size:12.8px"></font><br></div></div></div></div></div>
</div></font></span></div>
<br></div></div>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275m_-4196797236910994329m_-4307035958762947733gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><img src="https://docs.google.com/uc?export=download&id=0B5glieJcogGocjdGWHA3NWIwdmc&revid=0B5glieJcogGocU45NUpENzV4UWhCWk9tb1ZGdXB5MEtGYjJjPQ"><br><div><font color="#000000"><br></font></div><div><b><font color="#000000"><span style="font-size:12.8px">Lara Quitián Hernández, </span></font><span style="color:rgb(0,0,0);font-size:12.8px">PhD Student</span></b><font color="#000000"><b><br style="font-size:12.8px"></b><span style="font-size:12.8px">Departamento de Física de la Tierra, Astronomía y Astrofísica</span><br style="font-size:12.8px"><span style="font-size:12.8px">Facultad de Ciencias Físicas</span><br style="font-size:12.8px"><span style="font-size:13px"><a href="https://maps.google.com/?q=Ciudad+Universitaria,+Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">Ciudad Universitaria, Plaza Ciencias, 1</a></span></font></div><div><font color="#000000"><span style="font-size:12.8px"><a href="https://maps.google.com/?q=Ciudad+Universitaria,+Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">28040 Madrid. Spain</a></span><br style="font-size:12.8px"><span style="font-size:12.8px">Email: </span><a href="mailto:mlunar@aemet.es" style="font-size:12.8px" target="_blank">lquitian@ucm.es</a><br style="font-size:12.8px"></font><br></div></div></div></div></div>
</div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251m_3603104038903624767m_-1059488269469768275gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><img src="https://docs.google.com/uc?export=download&id=0B5glieJcogGocjdGWHA3NWIwdmc&revid=0B5glieJcogGocU45NUpENzV4UWhCWk9tb1ZGdXB5MEtGYjJjPQ"><br><div><font color="#000000"><br></font></div><div><b><font color="#000000"><span style="font-size:12.8px">Lara Quitián Hernández, </span></font><span style="color:rgb(0,0,0);font-size:12.8px">PhD Student</span></b><font color="#000000"><b><br style="font-size:12.8px"></b><span style="font-size:12.8px">Departamento de Física de la Tierra, Astronomía y Astrofísica</span><br style="font-size:12.8px"><span style="font-size:12.8px">Facultad de Ciencias Físicas</span><br style="font-size:12.8px"><span style="font-size:13px">Ciudad Universitaria, <a href="https://maps.google.com/?q=Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">Plaza Ciencias, 1</a></span></font></div><div><font color="#000000"><span style="font-size:12.8px"><a href="https://maps.google.com/?q=Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">28040 Madrid. Spain</a></span><br style="font-size:12.8px"><span style="font-size:12.8px">Email: </span><a href="mailto:mlunar@aemet.es" style="font-size:12.8px" target="_blank">lquitian@ucm.es</a><br style="font-size:12.8px"></font><br></div></div></div></div></div>
</div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><span class="m_7015166853841739983m_-3511754566990517289HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492m_2749377207748208251gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><img src="https://docs.google.com/uc?export=download&id=0B5glieJcogGocjdGWHA3NWIwdmc&revid=0B5glieJcogGocU45NUpENzV4UWhCWk9tb1ZGdXB5MEtGYjJjPQ"><br><div><font color="#000000"><br></font></div><div><b><font color="#000000"><span style="font-size:12.8px">Lara Quitián Hernández, </span></font><span style="color:rgb(0,0,0);font-size:12.8px">PhD Student</span></b><font color="#000000"><b><br style="font-size:12.8px"></b><span style="font-size:12.8px">Departamento de Física de la Tierra, Astronomía y Astrofísica</span><br style="font-size:12.8px"><span style="font-size:12.8px">Facultad de Ciencias Físicas</span><br style="font-size:12.8px"><span style="font-size:13px">Ciudad Universitaria, <a href="https://maps.google.com/?q=Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">Plaza Ciencias, 1</a></span></font></div><div><font color="#000000"><span style="font-size:12.8px"><a href="https://maps.google.com/?q=Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">28040 Madrid. Spain</a></span><br style="font-size:12.8px"><span style="font-size:12.8px">Email: </span><a href="mailto:mlunar@aemet.es" style="font-size:12.8px" target="_blank">lquitian@ucm.es</a><br style="font-size:12.8px"></font><br></div></div></div></div></div>
</font></span></div><span class="m_7015166853841739983m_-3511754566990517289HOEnZb"><font color="#888888">
</font></span></div></div></blockquote></div><span class="m_7015166853841739983m_-3511754566990517289HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div class="m_7015166853841739983m_-3511754566990517289m_5759546487367864492gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><img src="https://docs.google.com/uc?export=download&id=0B5glieJcogGocjdGWHA3NWIwdmc&revid=0B5glieJcogGocU45NUpENzV4UWhCWk9tb1ZGdXB5MEtGYjJjPQ"><br><div><font color="#000000"><br></font></div><div><b><font color="#000000"><span style="font-size:12.8px">Lara Quitián Hernández, </span></font><span style="color:rgb(0,0,0);font-size:12.8px">PhD Student</span></b><font color="#000000"><b><br style="font-size:12.8px"></b><span style="font-size:12.8px">Departamento de Física de la Tierra, Astronomía y Astrofísica</span><br style="font-size:12.8px"><span style="font-size:12.8px">Facultad de Ciencias Físicas</span><br style="font-size:12.8px"><span style="font-size:13px">Ciudad Universitaria, <a href="https://maps.google.com/?q=Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">Plaza Ciencias, 1</a></span></font></div><div><font color="#000000"><span style="font-size:12.8px"><a href="https://maps.google.com/?q=Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g" target="_blank">28040 Madrid. Spain</a></span><br style="font-size:12.8px"><span style="font-size:12.8px">Email: </span><a href="mailto:mlunar@aemet.es" style="font-size:12.8px" target="_blank">lquitian@ucm.es</a><br style="font-size:12.8px"></font><br></div></div></div></div></div>
</font></span></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br></div></div><span class="m_7015166853841739983HOEnZb"><font color="#888888"><div class="m_7015166853841739983m_-3511754566990517289gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist,  </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a>   </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_7015166853841739983gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><img src="https://docs.google.com/uc?export=download&id=0B5glieJcogGocjdGWHA3NWIwdmc&revid=0B5glieJcogGocU45NUpENzV4UWhCWk9tb1ZGdXB5MEtGYjJjPQ"><br><div><font color="#000000"><br></font></div><div><b><font color="#000000"><span style="font-size:12.8px">Lara Quitián Hernández, </span></font><span style="color:rgb(0,0,0);font-size:12.8px">PhD Student</span></b><font color="#000000"><b><br style="font-size:12.8px"></b><span style="font-size:12.8px">Departamento de Física de la Tierra, Astronomía y Astrofísica</span><br style="font-size:12.8px"><span style="font-size:12.8px">Facultad de Ciencias Físicas</span><br style="font-size:12.8px"><span style="font-size:13px">Ciudad Universitaria, <a href="https://maps.google.com/?q=Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g">Plaza Ciencias, 1</a></span></font></div><div><font color="#000000"><span style="font-size:12.8px"><a href="https://maps.google.com/?q=Plaza+Ciencias,+1+28040+Madrid.+Spain&entry=gmail&source=g">28040 Madrid. Spain</a></span><br style="font-size:12.8px"><span style="font-size:12.8px">Email: </span><a href="mailto:mlunar@aemet.es" style="font-size:12.8px" target="_blank">lquitian@ucm.es</a><br style="font-size:12.8px"></font><br></div></div></div></div></div>
</div>
</div></div></blockquote></div><br></div>