allow incidence angle to be in a column

Author: scottprahl

docs/manual.pdf

@@ -11,13 +11,15 @@ typedef struct {
 typedef struct Cache Cache;
 
 int cache_init(Cache * c);
+
 void cache_free(Cache * c);
 
 int cache_get(const Cache * c, double a, double b, double g, double *ur1, double *ut1, double *uru, double *utu);
 
 void cache_put(Cache * c, double a, double b, double g, double ur1, double ut1, double uru, double utu);
 
 size_t cache_size(const Cache * c);
+
 int cache_get_by_index(const Cache * c, size_t i,
     double *a, double *b, double *g, double *ur1, double *ut1, double *uru, double *utu);
 

src/cache.h

@@ -206,6 +206,17 @@ int Read_Data_Line_Per_Labels(FILE *fp, struct measure_type *m, struct invert_ty
         case 'g':
             r->default_g = x;
             break;
+        case 'i':
+            m->slab_cos_angle = cos(x * M_PI / 180.0);
+            if (x < 0 || x > 90) {
+                fprintf(stderr, "Error\n");
+                fprintf(stderr, "    incident angle '%8.4f'\n", x);
+                fprintf(stderr, "    must be between 0 and 90 degrees\n");
+                exit(EXIT_FAILURE);
+            }
+            if (r->method.quad_pts % 12 != 0)
+                r->method.quad_pts = ((r->method.quad_pts + 11) / 12) * 12;
+            break;
         case 'L':
             m->lambda = x;
             break;

src/iad_io.c

@@ -247,6 +247,17 @@ int Read_Data_Line_Per_Labels(FILE *fp, struct measure_type *m, struct invert_ty
             case 'g':
                 r->default_g = x;
                 break;
+            case 'i':
+                m->slab_cos_angle = cos(x*M_PI/180.0);
+                if (x<0 || x>90) {
+                    fprintf(stderr, "Error\n");
+                    fprintf(stderr, "    incident angle '%8.4f'\n", x);
+                    fprintf(stderr, "    must be between 0 and 90 degrees\n");
+                    exit(EXIT_FAILURE);
+                }
+                if (r->method.quad_pts % 12 != 0)
+                    r->method.quad_pts = ((r->method.quad_pts + 11) / 12) * 12;
+                break;
             case 'L':
                 m->lambda = x;
                 break;

src/iad_io.w

@@ -27,7 +27,7 @@
 #define BIG_A_CALC_VALUE 999999.0
 #define SMALL_A_CALC_VALUE 0.00001
 #define MAX_ABS_G 0.999999
-#define GRID_SIZE 41 \
+#define GRID_SIZE 101 \
 
 #define IAD_NO_ERROR 0 \
 

src/iad_type.h

@@ -36,7 +36,7 @@ the header file \.{iad\_type.h} and nothing else.
 @d BIG_A_CALC_VALUE              999999.0
 @d SMALL_A_CALC_VALUE            0.00001
 @d MAX_ABS_G                     0.999999
-@d GRID_SIZE                     41
+@d GRID_SIZE                     101
 
 @ Need error codes for this silly program
 

src/iad_type.w

@@ -1,5 +1,5 @@
 /* Autogenerated v3-16-3 from https://github.com/scottprahl/iad */
 
 The version number obtained from <.. / Makefile >
-    The date is the last change to < version.h > */char *Version = "3-16-3 (19 Jun 2025)";
+    The date is the last change to < version.h > */char *Version = "3-16-3 (18 Jul 2025)";
 char *VersionShort = "3-16-3";