National Oceanic and
Atmospheric Administration
United States Department of Commerce

working with go fland

working with go fland

Question:

How can I get `go fland' to fill the same region as the data?

Example:

yes? set data coads_climatology
yes? shade/x=-20:340 sst[l=1]
yes? go fland

working with go fland

Explanation:

The coverage issue arises because the ETOPO data sets are (presently) organizedas 20E to 20E. If you do not explicitly set the longitude region,`go fland' will default to the ETOPOdefault region, not that of your data.

This typically happens in one of two ways:

  • Using the dataset default region:
    set data my_0_to_360_dsetshade myvargo fland
  • Setting the region as a qualifier (as in the example):
    set data my_dsetshade/x=lo:hi myvargo fland

Either set of commands above results in incorrect output.

Solution:

The solution is to be explicit in the longitude region.

  • Using `go fland' with an `x=lo:hi' argument:
    set data my_0_to_360_dsetshade myvargo fland " " " " " " " " x=0:360
  • Using `set region':
    set data my_dsetset region/x=lo:hishade myvar