c  ISLANDS 
c--------------------------------

#ifdef ANTARCTICA
#define ISLAND         ! is there a periodic island? (antarctica)
#endif

#ifdef OTHER_ISLAND
#define ISLAND1        ! is there a standard kind of island? (e.g., New Zealand)
#endif
c#define ISLAND_SINK    ! are there islands to be sunk?

c--------------------------------
      parameter (IUNIT_OUT = 18)

      character*1 BC_W, BC_L, BC_P, BC_0, BC_1, BC_S
      parameter (BC_W='W',BC_L='L',BC_P='P',BC_0='0',BC_1='1',BC_S='S')

      integer    SCH_5, SCH_9
      parameter (SCH_5=5, SCH_9=9)

      parameter(MAX_NX = 1000,  MAX_NY = 1000)

      parameter (PREC = 1.e-6)
      real*8    MATH_PI,MATH_E
      parameter (MATH_PI=3.14159265358979323846)
      parameter (MATH_E= 2.7182818284590452354)

#define to_RAD(v) ((v)*MATH_PI/180.)
#define to_DEG(v) ((v)*180./MATH_PI)
c-----------------------------------------------------------------
      common /MOD_STATE/ mod_cosys
      common /MOD_CONST/ CNST_EPS, CNST_2OMEGA, R_EARTH
     *           , CNST_BETA, CNST_F0, CNST_NORM, GLUBINA
     *           , mod_island, CNST_EPT
c----------------------------------------------------------
      character*1 mask(1)
      real*8  hx, hy, R_EARTH
      integer list(1)
      pointer (pmask,mask), (plist,list)

      real*4 xx(MAX_NX), yy(MAX_NY)
      real*4 hx_cos(MAX_NY), r_tan(MAX_NY)
      common /MOD_GRID/ NX_MODEL, NY_MODEL, NX, NY, NXY, NPACK
     *                 ,pmask, plist, if_per, iper
      common /MOD_GRID1/ hx, hy, hx_s, hy_s
      common /MOD_GRID2/ X_MAX, Y_MAX, X_MIN, Y_MIN, xx, yy
      common /MOD_GRID3/ hx_cos, r_tan

      real deph(1)
      real fcor(1), bemx(1), bemy(1), bemxy(1), bemxx(1), bemyy(1)
      real relx_m(1), relx_p(1), rely_m(1), rely_p(1)
      pointer (pdeph, deph), (pfcor,fcor)
      pointer (pbemx, bemx), (pbemy, bemy), (pbemxy, bemxy) 
      pointer (pbemxx, bemxx), (pbemyy, bemyy)
      pointer (prelx_m,relx_m),(prelx_p,relx_p)
      pointer (prely_m,rely_m),(prely_p,rely_p)
      real taux(1), tauy(1)
      pointer (ptaux, taux), (ptauy, tauy)
      common /MOD_DATA/ pdeph, prelx_p, prelx_m, prely_p, prely_m, pfcor
     *                  , pbemx, pbemy, pbemxy, pbemxx, pbemyy
      common /MOD_DT/  ptaux, ptauy

      integer iro(1), ico(1), sn(1), ha(1)
      pointer (piro, iro), (pico, ico), (pha, ha), (psn, sn)
      real a1(1), pivot(1)
      pointer (pa1, a1), (ppivot, pivot)
      common /MOD_GRAPH/ NONZ, NGRAPH
      common /MOD_Y1/  NN12, pha, psn, ifail
      common /MOD_Y2/ piro, pico, pa1, ppivot
c-----------------------------------------------------------
      real aa(1)
      real rhs0(1), rhs1(1), rhs(1)
      real bound_rhs(1)
      real rhs_bc(1), rhs_bc0(1), sol(1), rhs_bc1(1)
      pointer (paa, aa)
      pointer (prhs0, rhs0), (prhs1, rhs1), (prhs, rhs)
      pointer (pbound_rhs, bound_rhs)
      pointer (prhs_bc, rhs_bc), (psol, sol), (prhs_bc0, rhs_bc0)
      pointer (prhs_bc1, rhs_bc1)
      real    b_island(0:2), a_island(0:2,0:2)

      common /ISLAND_1/ prhs_bc0, psol, prhs0, prhs_bc1, ib, prhs1, prhs
      common /ISLAND_2/ b_island, a_island
     *                 ,i_min1,i_max1,j_min1,j_max1
      common /MOD_MATRX/ paa, prhs_bc, pbound_rhs
c--------------------------------------------------------------
      real    rtmp(1)
      integer ilst(1)
      pointer (prtmp, rtmp), (pilst, ilst)
      common /MOD_SPACE/ prtmp, pilst
c--------------------------------------------------------------


c--------------------------------------------------------------
c   new common blocks to be added also to model_input subroutine
c--------------------------------------------------------------
      character*80       f_bar
      common /baro_files/ n_bar, f_bar
      integer iflag(12), RM12_NN 
      real    aflag(12)
      common /y12m_input/ RM12_NN, aflag, iflag
      common /baro_input/ n_def_cor, mod_scheme, mod_solver, BAR_DELTA,
     *                    BAR_DSINK, ibar_key, nbaro, rayl
c--------------------------------------------------------------