[ncl-talk] lcmask_4.ncl_error

[Beata Szabo]

Hi Mary,

Thank you very  much for your suggestion. It works properly.

Beata
> 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
>>
>>
>


Dr. Beata Szabo-Takacs
Global Change Research Institute CAS
Bělidla 986/4a
60300 Brno
Czech Republic