[ncl-talk] Issue regridding with ESMF from unstructured

[Anil Kumar - NOAA Affiliate]

It should work actually. Are you reading a variable in correct order of
dimensions? You mention (lev, time) and your data looks (time, lev). please
check.
Maybe some reading issue of file data...

-Anil


On Fri, Jun 29, 2018 at 1:38 PM, Andrew Kren - NOAA Affiliate <
andrew.kren at noaa.gov> wrote:

> Thanks Anil. I tried a few additions from your code but still ran into
> issues. It is still giving me the same error.
>
> On Fri, Jun 29, 2018 at 1:14 PM, Anil Kumar - NOAA Affiliate <
> anil.kumar at noaa.gov> wrote:
>
>> Andrew,
>> I did similar thing to regrid unstructed data on to wrf grid (
>> geo_em_d01.nc)  here is the code and worked successfully with ncl 6.4
>> version.
>> Hope this help.
>>
>> src_file = "maxele.63.nc"
>>
>>  sfile = addfile(src_file,"r")
>>
>>
>> ;printVarSummary(sfile)
>>
>>  lat1D = sfile->y
>>
>>  lon1D = sfile->x
>>
>>  var = sfile->zeta_max(0,:)   ; zeta is function of time and variable.
>>
>>
>>  dst_file = "geo_em.d01.nc"
>>
>>  dfile    = addfile(dst_file,"r")
>>
>>
>>
>> ;---Set up regridding options
>>
>>
>>     Opt                = True
>>
>>
>> ;---"bilinear" is the default. "patch" and "conserve" are other options.
>>
>>     Opt at InterpMethod     = "neareststod"        ;;---Change (maybe)
>>
>>     Opt at WgtFileName      = "unstruct_to_rect.nc"
>>
>>
>>     Opt at SrcGridLat       = lat1D
>>
>>     Opt at SrcGridLon       = lon1D
>>
>>     Opt at SrcRegional      = False            ;;--Change (maybe)
>>
>>     Opt at SrcInputFileName = src_file          ; optional, but good idea
>>
>>
>>     Opt at SrcMask2D        = where(.not.ismissing(var),1,0) ; Necessary if
>> has
>>
>>                                                           ; missing
>> values.
>>
>>
>>     dst_lat              = dfile->XLAT_M(0,:,:)    ; Be sure to use
>> appropriate names
>>
>>     dst_lon              = dfile->XLONG_M(0,:,:)   ; here.
>>
>>     Opt at DstGridLon       = dst_lon
>>
>>     Opt at DstGridLat       = dst_lat
>>
>>     Opt at DstRegional      = True            ;;--Change (maybe)
>>
>> ;    Opt at DstMask2D = where(.not.ismissing(dst_lat).and.\
>>
>> ;                          .not.ismissing(dst_lon),1,0) ; Necessary if
>> lat/lon
>>
>>                                                       ; has missing
>> values.
>>
>>                                                       ; has missing
>> values.
>>
>>     Opt at ForceOverwrite   = True
>>
>>     Opt at Debug            = True
>>
>>     Opt at PrintTimings     = True
>>
>>
>>     var_regrid = ESMF_regrid(var,Opt)     ; Do the regridding
>>
>>
>>     printVarSummary(var_regrid)
>>
>>
>>
>>
>>
>>
>>
>>
>> On Fri, Jun 29, 2018 at 11:56 AM, Andrew Kren - NOAA Affiliate <
>> andrew.kren at noaa.gov> wrote:
>>
>>> Dear ncl-talk,
>>>
>>> I'm trying to regrid from an unstructured grid of temperature and
>>> moisture retrievals from an aircraft. The data are a function of
>>> (lev,time), with each time corresponding to a different lat/lon location. I
>>> wrote some code to regrid to a triangular mesh using ESMF, but ran into
>>> some issues. I got this error:
>>>
>>> (0) get_src_grid_info: source lat dims = (14398)
>>>
>>> (0) get_src_grid_info: source lon dims = (14398)
>>>
>>> (0) get_src_grid_info: SrcGridType and/or SrcGridLat/SrcGridLon were
>>> not set.
>>>
>>> (0)           Cannot determine the source grid type.
>>>
>>> I have attached my code and a sample file it reads. I am confused b/c I
>>> am setting the SrcGridLat and Lon explicitly.
>>>
>>> Thanks,
>>>
>>> --
>>> Andrew Kren
>>> Assistant Scientist
>>> CIMAS - NOAA/AOML
>>> 314-322-0867
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>>
>> --
>> Dr. Anil Kumar
>> NWS/NCEP/EMC, Office# 2875
>> NOAA Center for Weather and Climate Prediction (NCWCP)
>> 5830 UNIVERSITY RESEARCH CT
>> <https://maps.google.com/?q=5830+UNIVERSITY+RESEARCH+CT+College+Park,+MD+20740&entry=gmail&source=g>
>> College Park, MD 20740
>> <https://maps.google.com/?q=5830+UNIVERSITY+RESEARCH+CT+College+Park,+MD+20740&entry=gmail&source=g>
>> -3818
>> 301-683-0494
>> anil.kumar at noaa.gov
>>
>>
>
>
> --
> Andrew Kren
> Assistant Scientist
> CIMAS - NOAA/AOML
> 314-322-0867
>



-- 
Dr. Anil Kumar
NWS/NCEP/EMC, Office# 2875
NOAA Center for Weather and Climate Prediction (NCWCP)
5830 UNIVERSITY RESEARCH CT
College Park, MD 20740-3818
301-683-0494
anil.kumar at noaa.gov
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180629/12bc8fc6/attachment.html>