Re: EASE grid and esmf warning messages

Hi Mary,

Yes, that was exactly it: this grid has missing values in the first/last
rows/cols. Technically speaking, I did "see" this Tuesday when I was
looking at my new ncfile in ncview, or at least I saw that there was a
FillValue along the edge. Didn't think much of it other than briefly
wondering why. Certainly didn't connect that with the ESMF warning
messages.

Rather than going unstructured (I really, really want to keep it
curvilinear if I can) I first changed my ease_to_nc.ncl script to first
look for totally empty missing cols/rows. Then after I confirmed the
above situation, I changed it to drop the first/last row/col before
proceeding. Not very flexible but it does the job for this file.

I'm happy to share with the ESMF regridding page. My plan is to regrid
CMIP5 GCMs to this EASE grid so any variable from any of those files
should work as an example. I just created a 1981-2000 average of daily
January tas from CCSM4 in /glade/scratch/dbr/tas_avg_1981-2000_01.nc.
Haven't tried regridding it but have no reason to think it won't work.

I really do try to follow Dennis' maxim of "know your data" but somehow
didn't quite see this one. Thanks for doing the hunting down I should
have done myself.

Cheers,
Dave

Mary Haley said the following on 5/22/2014 2:47 PM:
> Dave,
>
> I think the issue is that your lat/lon arrays contain missing values all along two rows and two columns. Even though you are setting the SrcMask array, it still has trouble because of the nature of the missing values.
>
> One possibility is to treat this as an unstructured grid, by throwing away all the missing points, and saving the valid points in 1D arrays for a call to unstructured_to_ESMF. See the attached script.
>
> This is a unique example that could be useful for our ESMF regridding page. If you could provide me with the variable you want to regrid, and what you are regridding to (or from), then I can try to come up with a complete example for this page.
>
> Thanks,
>
> --Mary
>
>
>
> On May 22, 2014, at 9:21 AM, David B. Reusch <dreusch@ees.nmt.edu> wrote:
>
>> Hi,
>>
>> I am getting a huge number of warning messages when doing curvilinear_to_SCRIP for an EASE grid with lots of FillValues. I recognize they are *warning* mesages so possibly nothing to worry about, but I wanted to be sure.
>>
>> I'm using the southern hemisphere EASE grid from NSIDC so the square input array has a circle of valid data in the middle surrounded by FillValues. I seem to be getting two of the message below for each FillValue point:
>> warning:avg: Entire input array contains missing values, can't compute average
>>
>> Assuming it would work, and as long as I'm not hiding a real problem, I'm willing to use the "how to suppress warnings" code (from the FAQ).
>> err = NhlGetErrorObjectId()
>> setvalues err
>> "errLevel" : "Fatal" ; only report Fatal errors
>> end setvalues
>> I've put esmf_def_curvi.tgz in the incoming ftp. It has the ncl script (define_grid_curvilinear.ncl), input (EASE_125KM_lat_lon.nc) & output (EASE_125KM_esmf_grid_reg.nc) files and output messages (def_EASE125.out) from the script.
>>
>> 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 Pl
New Mexico Tech
Socorro, NM 87801

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu May 22 15:19:07 2014