[ncl-talk] Problem with mask_9.ncl example

[Lyndz]

Dear Sir Dennis,

Thank you so much for help (always).

I have a question about the output of this sample script.

I am not sure why the masked data only have a few gridpoints compared to
the original data overlayed with the shapefile?

The output of the script is attached to this email.

I'll appreciate any comment on this.

Lyndz

On Mon, Apr 6, 2020 at 9:48 AM Dennis Shea <shea at ucar.edu> wrote:

> Change:
>    u          = fin->precip(:,:,:)
> To
>    u          = fin->precip(0,:,:)
>
> On Sun, Apr 5, 2020 at 11:37 AM Lyndz via ncl-talk <ncl-talk at ucar.edu>
> wrote:
>
>> Dear NCL-experts,
>>
>> I am trying to replicate the mask_9.ncl example but using a new shapefile.
>> When I run the script, I got the following errors:
>>
>>  Copyright (C) 1995-2019 - All Rights Reserved
>>
>>  University Corporation for Atmospheric Research
>>
>>  NCAR Command Language Version 6.6.2
>>
>>  The use of this software is governed by a License Agreement.
>>
>>  See http://www.ncl.ucar.edu/ for more details.
>>
>> (0) shapefile_mask_data: Error: not a valid rectilinear, curvilinear, or
>> unstructured grid
>>
>>
>>
>> Attached are the shapefile and script.
>> The ncdump of the input netcdf file is below:
>>
>> *dimensions:*
>>
>> * time = UNLIMITED ; // (1 currently)*
>>
>> * bnds = 2 ;*
>>
>> * longitude = 360 ;*
>>
>> * latitude = 280 ;*
>>
>> *variables:*
>>
>> * double time(time) ;*
>>
>> * time:standard_name = "time" ;*
>>
>> * time:bounds = "time_bnds" ;*
>>
>> * time:units = "days since 1951-1-1 0" ;*
>>
>> * time:calendar = "standard" ;*
>>
>> * time:axis = "T" ;*
>>
>> * double time_bnds(time, bnds) ;*
>>
>> * double longitude(longitude) ;*
>>
>> * longitude:standard_name = "longitude" ;*
>>
>> * longitude:long_name = "Longitude" ;*
>>
>> * longitude:units = "degrees_east" ;*
>>
>> * longitude:axis = "X" ;*
>>
>> * double latitude(latitude) ;*
>>
>> * latitude:standard_name = "latitude" ;*
>>
>> * latitude:long_name = "Latitude" ;*
>>
>> * latitude:units = "degrees_north" ;*
>>
>> * latitude:axis = "Y" ;*
>>
>> * float precip(time, latitude, longitude) ;*
>>
>> * precip:long_name = "daily precipitation analysis interpolated onto
>> 0.25deg grids [mm/day]" ;*
>>
>> * precip:units = "mm/day" ;*
>>
>> * precip:_FillValue = -99.9f ;*
>>
>> * precip:missing_value = -99.9f ;*
>>
>> * precip:level_description = "Earth surface" ;*
>>
>> Any advise on how to do this correctly in NCL?
>>
>> I'll appreciate any help.
>>
>> Sincerely,
>> Lyndz
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200406/79d50754/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2020-04-06 at 9.34.01 PM.png
Type: image/png
Size: 32887 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200406/79d50754/attachment.png>