<div dir="ltr">Dennis,<div>I am still not able to mask out the WRF selected region from a large covered region on unstructured data. </div><div>any advice or suggestion to try. Thanks Anil</div><div><br></div><div>*********************************</div><div> src_file = "../<a href="http://fort.63.nc">fort.63.nc</a>"<br> sfile = addfile(src_file,"r")<br><br>;printVarSummary(sfile)<br> lat1D = sfile->y<br> lon1D = sfile->x<br> var = sfile->zeta(89,:)<br></div><div><div>dst_file = "../<a href="http://geo_em.d01.nc/" target="_blank">geo_em.d01.nc</a>"  //WRF grid <span class="gmail-im" style="color:rgb(80,0,80)"><br> dfile    = addfile(dst_file,"r")<br><br><br>    dst_lat              = dfile->XLAT_M(0,:,:)    ; Be sure to use appropriate names<br>    dst_lon              = dfile->XLONG_M(0,:,:)   ; here.<br></span></div><div><br></div><div><span class="gmail-im" style="color:rgb(80,0,80)">   Opt@SrcGridLat       = lat1D<br>    Opt@SrcGridLon       = lon1D<br></span>    Opt@SrcRegional      = True            ;;--Change (maybe)<span class="gmail-im" style="color:rgb(80,0,80)"><br>    Opt@SrcInputFileName = src_file          ; optional, but good idea<br><br></span>    Opt@SrcGridCornerLat = (/ min(dst_lat),  max(dst_lat)/)   ;sfile->lat_vertices    ; corners are necessary<br>    Opt@SrcGridCornerLon =  (/ min(dst_lon),  max(dst_lon)/)  ;sfile->lon_vertices    ; for "conserve" method</div></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Apr 1, 2022 at 11:04 AM Anil Kumar - NOAA Affiliate <<a href="mailto:anil.kumar@noaa.gov">anil.kumar@noaa.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Dennis,<div>Thanks for getting back to me. </div><div>Opt@SrcGridCornerLat and Opt@SrcGridCornerLon points provided through WRF max and min coordinates. <br></div><div>I don't know Is this way we can do it.</div><div>Src grid (Unstructured) is covering the Gulf and Atlantic and I want to mask out the region where my wrf grid is located. For example WRF grid is set up over the Houston Coastal region. </div><div>Is this correct?</div><div>Thanks</div><div><br></div><div> dst_file = "../<a href="http://geo_em.d01.nc" target="_blank">geo_em.d01.nc</a>"  //WRF grid <br> dfile    = addfile(dst_file,"r")<br><br><br>    dst_lat              = dfile->XLAT_M(0,:,:)    ; Be sure to use appropriate names<br>    dst_lon              = dfile->XLONG_M(0,:,:)   ; here.<br></div><div><br></div><div>   Opt@SrcGridLat       = lat1D<br>    Opt@SrcGridLon       = lon1D<br>    Opt@SrcRegional      = True            ;;--Change (maybe)<br>    Opt@SrcInputFileName = src_file          ; optional, but good idea<br><br>    Opt@SrcGridCornerLat = (/ min(dst_lat),  max(dst_lat)/)   ;sfile->lat_vertices    ; corners are necessary<br>    Opt@SrcGridCornerLon =  (/ min(dst_lon),  max(dst_lon)/)  ;sfile->lon_vertices    ; for "conserve" method<br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 28, 2022 at 9:22 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Given no responses, I speculate no user knows how to generate the edge vertices.</div><div><br></div><div><a href="https://www.ncl.ucar.edu/Applications/ESMF.shtml" target="_blank"><b>https://www.ncl.ucar.edu/Applications/ESMF.shtml</b></a></div><div><br></div><div>Example 6 reads the vertices from a file.</div><div><br></div><div>ESMF_regrid_6.ncl:    Opt@SrcGridCornerLat = sfile->lat_vertices    ; corners are necessary</div><div>ESMF_regrid_6.ncl:    Opt@SrcGridCornerLon = sfile->lon_vertices    ; for "conserve" method</div><div><br></div><div>===</div><div>The text above the examples contains:</div><div><br></div><div>"

You can (and are highly recommended to) provide your own grid corners
via the special Src/Dst/GridCornerLat/Lon resources. They must be
provided as an <i>N</i> x 4 array, where <i>N</i> represents the
dimensionality of the center lat/lon grid. For example, if the center
lat/lon grid is dimensioned 256 x 220, then the corner arrays must be
256 x 220 x 4.  Sometimes the corner grids are provided on the file
along with the center grids, with names like
"lat_vertices"/"lon_vertices" or "lat_bounds"/"lon_bounds"."</div><div><br></div><div>Good Luck<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 23, 2022 at 7:17 AM Anil Kumar - NOAA Affiliate via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi All,<div>I am using ncl script for esmf regrid and converting unstructured onto WRF grid. </div><div>There is nothing wrong with the script, it was working fine with the previous unstructured grid nodes ( node = 1813443). </div><div>Now I am facing issues running these scripts for new unstructured data nodes ( node = 6057307). It is taking a long time and CPU time terminates the script and even job_script takes more than 8 hours and time limitation issues with job_scripts. </div><div>Any suggestions to resolve this issue? Below is my script<br></div><div><br></div><div>However, when I run script, it says </div><div> This can be VERY slow if you have a large array.<br>(0)                           To speed up code, provide your own edge vertices via the<br>(0)                           'ElementVertices' option.<br></div><div>How to use ElementVertices?</div><div><br></div><div>*******************************************************</div><div> src_file = "../<a href="http://fort.63.nc" target="_blank">fort.63.nc</a>"<br> sfile = addfile(src_file,"r")<br> lat1D = sfile->y<br> lon1D = sfile->x<br> var = sfile->zeta(89,:)<br><br> dst_file = "../<a href="http://geo_em.d01.nc" target="_blank">geo_em.d01.nc</a>"<br> dfile    = addfile(dst_file,"r")<br>  Opt                = True<br><br>;---"bilinear" is the default. "patch" and "conserve" are other options.<br>    Opt@InterpMethod     = "neareststod"        ;;---Change (maybe)<br>    Opt@WgtFileName      = "unstruct_to_rect_2906_D1.nc"<br>    Opt@SrcGridLat       = lat1D<br>   Opt@SrcGridLon       = lon1D<br>    Opt@SrcRegional      = False            ;;--Change (maybe)<br>    Opt@SrcInputFileName = src_file          ; optional, but good idea<br><br>    Opt@SrcMask2D        = where(.not.ismissing(var),1,0) ; Necessary if has<br>                                                          ; missing values.<br></div><div><br>    dst_lat              = dfile->XLAT_M(0,:,:)    ; Be sure to use appropriate names<br>    dst_lon              = dfile->XLONG_M(0,:,:)   ; here.<br>    Opt@DstGridLon       = dst_lon<br>    Opt@DstGridLat       = dst_lat<br>    Opt@DstRegional      = True            ;;--Change (maybe)<br>    Opt@ForceOverwrite   = True<br>    Opt@Debug            = True<br>    Opt@PrintTimings     = True<br>    var_regrid = ESMF_regrid(var,Opt)     ; Do the regridding<br>    printVarSummary(var_regrid)<br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dr. Anil Kumar<div>Environmental Modeling Center/NCEP/NWS/NOAA<br><div cols="72">College Park, MD, USA  </div></div><div><br></div><div><br></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dr. Anil Kumar<div>Environmental Modeling Center/NCEP/NWS/NOAA<br><div cols="72">College Park, MD, USA  </div></div><div><br></div><div><br></div></div></div></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dr. Anil Kumar<div>Environmental Modeling Center/NCEP/NWS/NOAA<br><div cols="72">College Park, MD, USA  </div></div><div><br></div><div><br></div></div></div></div></div></div></div></div>