Re: ESMF gen weights problem

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Wed Nov 14 2012 - 11:55:33 MST

I checked with the ESMF folks, and they said there's no way to have this file automatically removed. They will probably create a ticket for this, but it's likely to be a low priority.

You could actually do this yourself with a system call:

system("/bin/rm PET0.RegridWeightGen.Log")

--Mary

On Nov 13, 2012, at 3:43 PM, David B. Reusch wrote:

> One last related question: is there any way to *not* have
> ESMF_regrid_gen_weights create the PET*Log file? It's not using much
> space but it is something I don't normally need and adds to file
> clutter. I didn't see anything in the available Opt attributes and I
> have Debug turned off.
>
> Thanks,
> Dave
>
> Mary Haley said the following on 11/13/12 2:59 PM:
>> On Nov 13, 2012, at 2:34 PM, David B. Reusch wrote:
>>
>>> Yes, that was it. Sigh. Thanks so much for the sharper eyes! With the
>>> power outage disruptions yesterday, I wasn't exactly at my best. Why I
>>> didn't catch it today though...who knows. Still, happier to have it be
>>> a mistake here then a nasty bug there.
>> Yeah, me too! :-)
>>
>>> There is still the question of ESMF_regrid_gen_weights saying everything
>>> is ok when it's not (i.e., n_s should *not* ever be zero, right?).
>> I would think so. I'll ask the ESMF folks about this. I've also created a ticket
>> to add a check in our ESMF_regrid_xxx code for this.
>>
>>
>>> And I also came across a super-minor type when poking around for this
>>> yesterday: at line 2063 in ESMF_regridding.ncl it says
>>> "ESMF_RegridGenWeight" when I think it should be "ESMF_RegridWeightGen".
>> Noted and corrected, thanks..
>>
>> --Mary
>>
>>> Thanks,
>>> Dave
>>>
>>> Mary Haley said the following on 11/13/12 2:23 PM:
>>>> I think the problem is with this line:
>>>>
>>>> lon2d_sub = lat2d_d(::2,::2)
>>>>
>>>> The variable on the right side should be "lon2d_d".
>>>>
>>>> Can you let me know if this fixes the issue on your end?
>>>>
>>>> --Mary
>>>>
>>>> On Nov 13, 2012, at 2:00 PM, David B. Reusch wrote:
>>>>
>>>>> Hello,
>>>>>
>>>>> I'm having a problem with ESMF regridding where ESMF_regrid_gen_weights appears to succeed but the var n_s = 0 and thus row, col and S are all empty. Not surprisingly, this leads to a segfault if I try to do a regrid_with_weights with this file. The weird part is that the target grid will work when going from one source grid but not another, as explained below.
>>>>>
>>>>> I have two potential input files/source grids:
>>>>> (1) an unmodified WRF file (technically it's been postprocessed but I haven't done any spatial subsetting yet). This grid is 357 x 381.
>>>>> (2) a spatial subset of (1) created by selecting points in the original file that lie within particular lat/lon boundaries (the "subdomain" file). This grid is 119 x 167.
>>>>>
>>>>> I have two potential target grids:
>>>>> (3) every other point from (1), 179 x 191.
>>>>> (4) every other point from (2), 60 x 84.
>>>>>
>>>>> Case 1. Input grid is lat/lon from (1) , output grid is (3). Result: OK
>>>>> Case 2. Input grid is lat/lon from (1) , output grid is (4). Result: OK
>>>>> Case 3. Input grid is lat/lon from (2) , output grid is (4). Result: n_s = 0, segfault on regrid.
>>>>>
>>>>> The clues so far point to a problem with relating source grid (2) to target grid (4), but this target grid does actually work with source grid (1).
>>>>>
>>>>> There are no messages in the PET0* log file and no obvious debug output suggesting errors.
>>>>>
>>>>> This is NCL 6.1 on OS X 10.8.2. I've put the script, input WRF files and ESMF files in file ESMF_Problem_Reusch.zip in the incoming directory at ftp.cgd.ucar.edu. For the ESMF files, prefix "Dmn" indicates case 3, prefix "Full" = case 2. Likewise, the WRF files have suffix "domain" and "full".
>>>>>
>>>>> Thanks,
>>>>> Dave
>>>>> --
>>>>> Associate Research Professor of Climatology
>>>>> Dept of Earth and Environmental Science
>>>>> MSEC 304; 801 Leroy Place
>>>>> New Mexico Tech
>>>>> Socorro, NM 87801
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>> --
>>> Associate Research Professor of Climatology
>>> Dept of Earth and Environmental Science
>>> MSEC 304; 801 Leroy Place
>>> New Mexico Tech
>>> Socorro, NM 87801
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
> --
> Associate Research Professor of Climatology
> Dept of Earth and Environmental Science
> MSEC 304; 801 Leroy Place
> New Mexico Tech
> Socorro, NM 87801
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Nov 14 11:55:43 2012

This archive was generated by hypermail 2.1.8 : Wed Nov 21 2012 - 11:16:05 MST