GRIB-659: Provide key to get angular precision

definitions/grib1/boot.def

@@ -29,9 +29,6 @@ meta section0Pointer section_pointer(offsetSection0,section0Length,0);
 # of the message length and of the section4 lenth, so instead of
 # section_length[3] totalLength                     ;
 # we get:
-
-
-
 g1_message_length[3] totalLength(section4Length) ;
 position startOfHeaders;
 unsigned[1] editionNumber = 1 : edition_specific,dump;
@@ -48,6 +45,7 @@ meta GDSPresent gds_is_present(gridDescriptionSectionPresent,gridDefinition,bitm
 flagbit bitmapPresent(section1Flags,6) :dump;
 alias bitmapSectionPresent=bitmapPresent;
 alias geography.bitmapPresent=bitmapPresent;
+transient angularPrecision=1000;  # milli degrees
 
 if(gridDescriptionSectionPresent){
    template section_2 "grib1/section.2.def"     ;

definitions/grib2/boot.def

@@ -12,6 +12,7 @@ constant million = 1000000 : hidden;
 constant grib2divider   = 1000000;
 alias extraDimensionPresent=zero;
 alias is_tigge = zero;
+transient angularPrecision=grib2divider;  # micro degrees
 
 meta gts_header gts_header() : no_copy,hidden,read_only;
 meta gts_TTAAii gts_header(20,6) : no_copy,hidden,read_only;

definitions/grib2/template.3.grid.def

@@ -37,8 +37,10 @@ unsigned[4] subdivisionsOfBasicAngle = missing() : can_be_missing;
 transient angleDivisor = 1000000;
 when (missing(subdivisionsOfBasicAngle) || subdivisionsOfBasicAngle == 0) {
   set angleDivisor  = 1000000;
+  set angularPrecision = 1000000;
 } else {
   set angleDivisor  = subdivisionsOfBasicAngle;
+  set angularPrecision = subdivisionsOfBasicAngle;
 }
 
 #  La1 - latitude of first grid point