Commit c4efca0b authored by jonas.sinjan's avatar jonas.sinjan
Browse files

defines.h change cmilos, hrt pipe pmilos change

parent d8814c86
......@@ -65,7 +65,7 @@
#define CUANTIC_JUPI 0
#define NOISE_SIGMA 0.004
#define NOISE_SIGMA 0.001
#define CLASSICAL_ESTIMATES_SAMPLE_REF 4 //Muestra referencia para cambio de cuadrante de azimuth. Depende del numero de muestras y posicion Continuo
......
......@@ -779,7 +779,7 @@ def pmilos(data_f, wve_axis_arr, data_shape, cpos_arr, data, rte, field_stop, st
printc(err.args[1],color=bcolors.FAIL)
return
wavelength = 6173.3356
wavelength = 6173.3354
for scan in range(int(data_shape[-1])):
......@@ -832,7 +832,7 @@ def pmilos(data_f, wve_axis_arr, data_shape, cpos_arr, data, rte, field_stop, st
hdul.writeto(f'./P-MILOS/run/data/input_tmp.fits', overwrite=True)
if rte == 'RTE':
cmd = "mpiexec -n 16 ../pmilos.x pmilos.minit" #../milos.x pmilos.mtrol" ##
cmd = "mpiexec -n 64 ../pmilos.x pmilos.minit" #../milos.x pmilos.mtrol" ##
if rte == 'CE':
cmd = "mpiexec -np 16 ../pmilos.x pmilos_ce.minit"
......@@ -841,6 +841,9 @@ def pmilos(data_f, wve_axis_arr, data_shape, cpos_arr, data, rte, field_stop, st
print("CE+RTE not possible on PMILOS, performing RTE instead")
cmd = "mpiexec -np 16 ../pmilos.x pmilos.minit"
if rte == 'RTE_seq':
cmd = '../milos.x pmilos.mtrol'
del sdata
#need to change settings for CE or CE+RTE in the pmilos.minit file here
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment