[RegCNET] "LANDUSE MASK exceed the limit"

[XUNQIANG BI]

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
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   Dr. Xunqiang Bi         email:bixq at ictp.it
   Earth System Physics Group
   The Abdus Salam ICTP
   Strada Costiera, 11
   P.O. BOX 586, 34100 Trieste, ITALY
   Tel: +39-040-2240302  Fax: +39-040-2240449
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~