Commit 8da1f295 authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Write CS_SPHERE, CSM_CLUSTER, CS1 and CS2 data to c_OCLU

parent f987530e

src/cluster/cluster.cpp

+13 −0
Original line number Diff line number Diff line
@@ -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];