[ncl-talk] lcmask_4.ncl_error

[Mary Haley]

Hi Beata,

This script was written for a larger lat/lon area.

If you simply change these two lines:

lat_values = ispan(toint(minlat),toint(maxlat),10) * 1.
lon_values = ispan(toint(minlon+10),toint(maxlon-10),10) * 1.

to:

lat_values = ispan(toint(minlat),toint(maxlat),1) * 1.
lon_values = ispan(toint(minlon+10),toint(maxlon-10),1) * 1.

it should work.

--Mary




On Thu, Sep 22, 2016 at 5:57 AM, Beata Szabo <szabo.b at czechglobe.cz> wrote:

> Dear All, I applied lcmask_4.ncl script from
> https://www.ncl.ucar.edu/Applications/Scripts/lcmask_4.ncl. It is able to
> create lambert conic projected map if I use large European area
> geographical values :
> minlat = 30
> maxlat = 75
> minlon = -20
> maxlon = 40
>
> I tried to run this script to create lambert projected map about Iberian
> Peninsula with:
> minlat = 36
> maxla = 44
> minlon = -10
> maxlon = 3
>
> but I got error messages
> fatal:divide: Division by 0, Can't continue
> fatal:Div: operator failed, can't continue
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 43 in
> file lambert_IP.ncl
>
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 198 in
> file lambert_IP.ncl
>
> where line 43 is:
>   slope_lft = (lat2_ndc-lat1_ndc)/(lon2_ndc-lon1_ndc)
> and line 198 is:
>  add_lc_labels(wks,plot(0),minlat,maxlat,minlon,maxlon)
>
> I attached the script. Can someone suggest me how I can fix this error?
> Thank you for your help in advance!
> Beata
>
>
> Dr. Beata Szabo-Takacs
> Global Change Research Institute CAS
> Bělidla 986/4a
> 60300 Brno
> Czech Republic
> _______________________________________________
> 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/20160926/a38ba979/attachment.html