module const
    integer, parameter::    r8=selected_real_kind(13)
    integer, parameter::    i2t=320
    integer, parameter::    j2t=288
    integer, parameter::    k0=56

    integer, parameter::    i0t=i2t+2
    integer, parameter::    j0t=j2t+2
    integer, parameter::    i1t=i0t-1
    integer, parameter::    j1t=j0t-1 

    integer, parameter::    ngx0=32
    integer, parameter::    ngy0=9
    integer, parameter::    n0=4
    integer, parameter::    ng0=ngx0*ngy0 
    integer, parameter::    ngx1=ngx0-1
    integer, parameter::    ngy1=ngy0-1
    integer, parameter::    nb0=j2t/ngy0/n0
    integer, parameter::    nb1=nb0-1 
    integer, parameter::    i0=i2t/ngx0+2
    integer, parameter::    i1=i0-1
    integer, parameter::    i2=i0-2 
    integer, parameter::    j0=j2t/ngy0+2
    integer, parameter::    j1=j0-1
    integer, parameter::    j2=j0-2 
    integer, parameter::    k1=k0-1
    integer, parameter::    k2=k0-2

    integer, parameter::    ncopt=20
    integer, parameter::    ini_clim=1

    integer, parameter::   nbg0=nb0*ngy0
    integer, parameter::   nbg1=nbg0-1
    integer, parameter::   i01=i0+1
    integer, parameter::   ibir=i2t/ngy0
    integer, parameter::   jg0=j2*ngy0+2
    integer, parameter::   j01=j0+1
    integer, parameter::   jg1=jg0-1
    integer, parameter::   jg2=jg0-2
    integer, parameter::   ig0=i2*ngx0+2
    integer, parameter::   ig1=ig0-1
    integer, parameter::   ig2=ig0-2

    real*8, parameter :: o12=1./12.
    real*8, parameter :: o24=1./24.
    real*8, parameter :: cn0=0.
    real*8, parameter :: cn1=1.
    real*8, parameter :: cn2=2.
    real*8, parameter :: cn10=10.
    real*8, parameter :: cn1000=1000.  
    real*8, parameter :: grav=980.6  !gravit. accel. (cm/s^2)
    real*8, parameter :: eps=1.0e-10
    real*8, parameter :: eps2=1.0e-20
    real*8, parameter :: pi=4.d0*atan(1.d0)

    real*8, parameter :: ppt_to_salt=1.d-3 ! from psu to msu. namely, ppt to g/g.
    real*8, parameter :: salt_to_ppt=1000.d0 ! from msu to psu
    real*8, parameter :: sea_ice_salinity=4.0  !salinity of sea ice formed (psu)
    real*8, parameter :: ocn_ref_salinity=34.7 !ocean reference salinity (psu)
    real*8, parameter :: rho_sw=4.1d0/3.996d0 !density of seawater in (g/cm^3)
    real*8, parameter :: cp_sw=3.996e7    ! specific heat salt water
    real*8, parameter :: hflux_factor=1000.d0/(rho_sw*cp_sw) !(W/m^2)to(Ccm/s)
   
    

end module