

    if inpr3 == 1
        FILNAME='private/gshhs_c.b';
        m_gshhs_c('save','coastl.mat');

    elseif inpr3 == 2
        FILNAME='private/gshhs_l.b';
        m_gshhs_l('save','coastl.mat');

    elseif inpr3 == 3
        FILNAME='private/gshhs_i.b';
        m_gshhs_i('save','coastl.mat');


    elseif inpr3 == 4
        FILNAME='private/gshhs_h.b';
        m_gshhs_h('save','coastl.mat');
    end
    load coastl.mat

    coastline = ncst;
    update(mainmap())
    clear  ncst coastl



