fix angles bug

magpowell · magpowell · commit 56ff4dd2bbc2 · 2025-04-08T10:23:13.000-07:00
diff --git a/src_test/test_rte_rrtmgp_rt.cu b/src_test/test_rte_rrtmgp_rt.cu
@@ -449,11 +449,11 @@ void solve_radiation(int argc, char** argv)
     mu0 = input_nc.get_variable<Float>("mu0", {n_col_y, n_col_x});
     azi = input_nc.get_variable<Float>("azi", {n_col_y, n_col_x});
 
-    tica_sza = acos(mu0.v()[0]) * 180.0 / M_PI;
-    tica_azi = azi.v()[0];
     
     if (switch_tica)
     {
+        tica_sza = acos(mu0.v()[0]);
+        tica_azi = azi.v()[0];
         for (int icol=1; icol<=n_col; ++icol)
         {
             mu0({icol}) = 1.0;

Original file line number	Diff line number	Diff line change
`@@ -449,11 +449,11 @@ void solve_radiation(int argc, char** argv)`
`449`	`449`	`mu0 = input_nc.get_variable<Float>("mu0", {n_col_y, n_col_x});`
`450`	`450`	`azi = input_nc.get_variable<Float>("azi", {n_col_y, n_col_x});`
`451`	`451`
`452`		`- tica_sza = acos(mu0.v()[0]) * 180.0 / M_PI;`
`453`		`- tica_azi = azi.v()[0];`
`454`	`452`
`455`	`453`	`if (switch_tica)`
`456`	`454`	`{`
	`455`	`+ tica_sza = acos(mu0.v()[0]);`
	`456`	`+ tica_azi = azi.v()[0];`
`457`	`457`	`for (int icol=1; icol<=n_col; ++icol)`
`458`	`458`	`{`
`459`	`459`	`mu0({icol}) = 1.0;`