Geopointset is the format used by Metview to combine a set of Geopoints variables into a single entity for ease of processing. Thus, a set of observations can be grouped in the same way that fields are grouped into a fieldset variable.
Creating a geopointset
A geopointset can be created with the create_geo_set() function, which takes any number of geopoints variables as arguments, or none. Both geopoints and geopointset variables can be concatenated to a geopointset.
set1 = create_geo_set() # creates an empty set set2 = create_geo_set(g1, g2, g3) # assuming that g1,g2,g3 are geopoints variables set3 = set1 & g1 & g2 # set3 has 2 geopoints set4 = set2 & set3 # set4 has 5 geopoints
The Geopointset file format
The format for a geopointset file is very simply a header followed by a contatenation of geopoints files - see Geopoints for details of the format. The overall header is this:
#GEOPOINTSET
The subsequent geopoints structures should all share the same format as each other. Here's an example with 3 geopoints files inside the set:
#GEOPOINTSET #GEO # lat lon height date time value # Missing values represented by 3e+38 (not user-changeable) #DATA 69.6523 18.9057 0 20130512 0 100869.8625 63.4882 10.8795 0 20130512 0 100282.3392 63.5657 10.694 0 20130512 0 100241.1666 61.2928 5.0443 0 20130512 0 99852.18932 #GEO # lat lon height date time value # Missing values represented by 3e+38 (not user-changeable) #DATA 60.82 23.5 0 20130512 600 101045.8 #GEO # lat lon height date time value # Missing values represented by 3e+38 (not user-changeable) #DATA 55.01 8.41 0 20130513 0 100949.1809 54.33 8.6 0 20130513 0 101027.9101 53.71 7.15 0 20130513 0 100846.619