parameter (IUNIT_OUT = 83)

      character*1 BC_W, BC_L, BC_P, BC_0, BC_I(0:9), bci
      parameter (BC_W='W',BC_L='L',BC_P='P',BC_0='0')

      common /masks/BC_I

      parameter(MAX_NX = 1000,  MAX_NY = 1000)

      parameter (PREC = 1.e-6)
      parameter (PI_MATH=3.14159265)
      parameter (R_EARTH = 6378000)

c-----------------------------------------------------------------
      common /MOD_CONST/ CNST_EPS, CNST_NORM, GLUBINA, CNST_EPT
c----------------------------------------------------------

      character*1 mask(1)
      integer list(1), ilst(1)
      pointer (pmask,mask), (plist,list), (pilst, ilst)
      common /MOD_INDICES/ pmask, plist, pilst

      common /MOD_GRID/ NX, NY, NXY, NPACK , if_per, iper
      common /MOD_GRID2/ X_MAX, Y_MAX, X_MIN, Y_MIN

      real deph(1),gradHx(1),gradHy(1)
      real fcor(1), bemx(1), bemy(1),  bemxx(1), bemyy(1)
      real relx_m(1), relx_p(1), rely_m(1), rely_p(1), bcosli(1)
      pointer (pdeph, deph), (pfcor,fcor)
      pointer (pgradHx, gradHx),(pgradHy, gradHy)
      pointer (pbemx, bemx), (pbemy, bemy)
      pointer (pbemxx, bemxx), (pbemyy, bemyy)
      pointer (pbcosli, bcosli)
      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, pbemxx, pbemyy, pbcosli
      common /MOD_DT/  ptaux, ptauy, pgradHx, pgradHy

      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 rhsi(1), rhs(1), rhs1(1)
      real bound_rhs(1),psim(1)
      real rhs_bc(1), rhs_bci(1), sol(1)
      pointer (paa, aa)
      pointer (prhsi, rhsi), (prhs, rhs), (prhs1,rhs1)
      pointer (pbound_rhs, bound_rhs), (ppsim, psim)
      pointer (prhs_bc, rhs_bc), (psol, sol), (prhs_bci, rhs_bci)

      common /ISLAND_1/ prhs_bci, psol, prhsi, prhs, ib, prhs1
      common /ISLAND_2/ b_island(0:9), a_island(0:9,0:9)
     *                 ,i_min1,i_max1,j_min1,j_max1
     *                 ,i_mins(10),i_maxs(10),j_mins(10),j_maxs(10)
      common /MOD_MATRX/ paa, prhs_bc, pbound_rhs, ppsim

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, nonlin_baro

      logical use_per_island, use_stan_island
      common /baro_island/
     *             alons_min(10),alons_max(10),alats_min(10),alats_max(10)
     *            ,alon1_min,alon1_max,alat1_min,alat1_max
     *            ,per_lat,n_sunk,use_per_island,use_stan_island
     *            ,itrans,ftrans

      common /new_island/n_islands,in_rhs(0:9),ityp_island(0:9),
     *               alon_island(9),alat_island(9),,btrans(0:9)
     *               n_poly(9),alon_poly(100,9),alat_poly(100,9),
     *               i_poly(100,9),j_poly(100,9)
c--------------------------------------------------------------


        common /sink_em/ nsunk(9),isunk(1000,9)