Fixed logic in editmask

[arango]

Coded a generic logic in matlab/landmask/editmask.m to process the coastline data used to determine land/sea masking in realistic applications. The internal switch got_coast was not initialized when the coastline file was provided as a second argument.

For the last few years, we have been adding the coastline data to the GRID NetCDF file into variables lon_coast and lat_coast:

dimensions:

        xi_rho = 258 ;
        xi_u = 257 ;
        xi_v = 258 ;
        xi_psi = 257 ;
        eta_rho = 130 ;
        eta_u = 130 ;
        eta_v = 129 ;
        eta_psi = 129 ;
        bath = UNLIMITED ; // (1 currently)
        coast = 11908 ;

variables:

        double lon_coast(coast) ;
                lon_coast:long_name = "Coastline longitude" ;
                lon_coast:units = "degree_east" ;
        double lat_coast(coast) ;
                lat_coast:long_name = "Coastline latitude" ;
                lat_coast:units = "degree_north" ;

This is a very useful strategy: the coastline data used to compute the land/sea masking for an application is in its GRID NetCDF file. This data can be also used for plotting.

Many thanks to Evridiki Chrisagi for ​reporting this.