[RegCNET] Terrain setup-Exceed the limit

[Rauscher Sara A]

Hi Satyaban,
Bi answered this question a while back - I've copied his answer below.
cheers
sara
On Fri, 8 Jun 2007, Rauscher Sara A wrote:

> What can I do about a "LANDUSE MASK exceed the limit" error in terrain?
>

There are some reports about the above problem.

"LANDUSE MASK            0 exceed the limit"

This usually happens at coastal points. Due to the missing value
(-129) in the global datasets (GLCC????_BATS.CDF files).

Before we provide the updated PreProc and updated global LANDUSE
datasets, the following tricky way can solve this problem temporarily.

1. edit terrain.f

search the places have this message

                 write(*,*) 'LANDUSE MASK', nint(lndout(i,j)),
     &                            ' exceed the limit'
                 STOP

and replace the above lines with

                ch(i,j) = 'X'
c                  write(*,*) 'LANDUSE MASK', nint(lndout(i,j)),
c      &                            ' exceed the limit'
c                  STOP

There are two places like this which should be changed.

2. Then run (you should have  FUDGE_LND=.false. in domain.param
    ./terrain.x

3. once you have LANDUSE file, edit it and replace 'X' with the type of
LANDUSE of neighbouring points, set FUDGE_LND=.true. in domain.param, rerun
    ./terrain.x


Best regards,
Xunqiang Bi

Satyaban Bishoyi Ratna wrote:
>
>  
> Dear RegCm Users,
> I am making a terrain setup over Indian region with 30 km
> resolution having 249X288 grid points. This is working fine.
> But when I increased the number of grid points to 329X350 with
> same 30 km resolution, it is not working and showing following
> error.
> *******************************************************
> ***** Terrain resolution (min):    30.00000
> after calling SETUP
> after calling MAP PROJECTION
> xminlat,xmaxlat,xminlon,xmaxlon,ntypec=    -21.67    55.09    30.5 4 
>   124.69        30
> after calling MXMNLL
>       29140 terrain heights read from land use volume
> after calling RDLDTR
> after calling XYOBSLL
>   rin,ds(m) =  0.150E+01  0.300E+05
>   rin,ds(m) =  0.150E+01  0.300E+05
> after calling ANAL2
> after calling SURF
> LANDUSE MASK            0 exceed the limit
> FORTRAN STOP
> *****************************************************
> Plese suggest me what went wrong when incresed the number of
> grid. I am using 'BATS' as LSMTYP and '10' as ntypec.
>
> Thanking you.
> Best regards
> Satyaban Bishoyi Ratna
>
> ****************************************************************************************************
>
> Satyaban Bishoyi Ratna
> Research Scholar
> Department of Meteorology and Oceanography
> Andhra University
> Visakhapatnam - 530003
> Andhra Pradesh, India
> Cell: 94403 30367
> E-mail: satyaban_2001 at rediffmail.com
>        satyaban14 at yahoo.co.in
>        
> ****************************************************************************************************
>
> Samsung
> <http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-home.htm/1050715198@Middle5/2047694_2040388/2047704/1?PARTNER=3&OAS_QUERY=null>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
>   

-- 

Dr. Sara A. Rauscher
Earth System Physics Section 
The Abdus Salam International Centre for Theoretical Physics
Strada Costiera 11 
34014 Trieste ITALY
ph: +39 040 2240 225 fax: +39 040 2240 449
email: srausche at ictp.it