From 8da1f29581e6ce2ae091130ab16690b5bba9adc2 Mon Sep 17 00:00:00 2001
From: Giovanni La Mura <giovanni.lamura@inaf.it>
Date: Mon, 27 May 2024 20:25:20 +0200
Subject: [PATCH] Write CS_SPHERE, CSM_CLUSTER, CS1 and CS2 data to c_OCLU
---
src/cluster/cluster.cpp | 13 +++++++++++++
1 file changed, 13 insertions(+)
diff --git a/src/cluster/cluster.cpp b/src/cluster/cluster.cpp
index 3a73321b..d6d427f2 100644
--- a/src/cluster/cluster.cpp
+++ b/src/cluster/cluster.cpp
@@ -867,6 +867,9 @@ int cluster_jxi488_cycle(int jxi488, ScattererConfiguration *sconf, GeometryConf
output->append_line(virtual_line);
sprintf(virtual_line, " FSAS=%15.7lE%15.7lE\n\0", real(cid->c1->fsas[i]), imag(cid->c1->fsas[i]));
output->append_line(virtual_line);
+ double alamb = 2.0 * 3.141592653589793 / cid->vk;
+ sprintf(virtual_line, "INSERTION: CS_SPHERE %15.7lE%15.7lE%15.7lE%15.7lE\n\0", alamb, cid->c1->sscs[i], cid->c1->sabs[i], cid->c1->sexs[i]);
+ output->append_line(virtual_line);
csch = 2.0 * cid->vk * cid->sqsfi / cid->c1->gcsv[i];
s0 = cid->c1->fsas[i] * exri;
qschu = imag(s0) * csch;
@@ -1098,6 +1101,9 @@ int cluster_jxi488_cycle(int jxi488, ScattererConfiguration *sconf, GeometryConf
output->append_line(virtual_line);
sprintf(virtual_line, " Fk=%15.7lE\n\0", fz);
output->append_line(virtual_line);
+ double alamb = 2.0 * 3.141592653589793 / cid->vk;
+ sprintf(virtual_line, "INSERTION: CSM_CLUSTER %15.7lE%15.7lE%15.7lE%15.7lE\n\0", alamb, scasm, abssm, extsm);
+ output->append_line(virtual_line);
} // ilr210 loop
double rmbrif = (real(cid->c1ao->fsacm[0][0]) - real(cid->c1ao->fsacm[1][1])) / real(cid->c1ao->fsacm[0][0]);
double rmdchr = (imag(cid->c1ao->fsacm[0][0]) - imag(cid->c1ao->fsacm[1][1])) / imag(cid->c1ao->fsacm[0][0]);
@@ -1345,6 +1351,13 @@ int cluster_jxi488_cycle(int jxi488, ScattererConfiguration *sconf, GeometryConf
qschu, pschu, s0mag
);
output->append_line(virtual_line);
+ double alamb = 2.0 * 3.141592653589793 / cid->vk;
+ if (ilr290 == 1) {
+ sprintf(virtual_line, "INSERTION: CS1_CLUSTER %9.2lf%9.2lf%9.2lf%9.2lf%9.2lf%13.5lf%13.5lf%13.5lf\n\0", alamb, th, ph, ths, phs, scasec, abssec, extsec);
+ } else if (ilr290 == 2) {
+ sprintf(virtual_line, "INSERTION: CS2_CLUSTER %9.2lf%9.2lf%9.2lf%9.2lf%9.2lf%13.5lf%13.5lf%13.5lf\n\0", alamb, th, ph, ths, phs, scasec, abssec, extsec);
+ }
+ output->append_line(virtual_line);
bool goto190 = isam >= 0 && (jths > 1 || jphs > 1);
if (!goto190) {
cid->gapv[0] = cid->gap[0][ilr290 - 1];
--
GitLab