grisinv issueshttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues2021-06-24T08:02:44Zhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/12Fix systematic error in time coordinate2021-06-24T08:02:44ZVigeesh GangadharanFix systematic error in time coordinate- [ ] Deal with `CSYER` for the time coordinate.
- [ ] Also, what's with `LONPOLE`, `LATPOLE`, `MJDREF` added by `astropy.WCS`.
- [ ] Also, change the units of WCS keywords to be consistent with L1.- [ ] Deal with `CSYER` for the time coordinate.
- [ ] Also, what's with `LONPOLE`, `LATPOLE`, `MJDREF` added by `astropy.WCS`.
- [ ] Also, change the units of WCS keywords to be consistent with L1.Vigeesh GangadharanVigeesh Gangadharanhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/17Option for adding new spectral lines by the user2021-03-24T10:03:10ZVigeesh GangadharanOption for adding new spectral lines by the userIt would be nice to have an option to add new spectral lines on the fly. <br>
Need to think of how to deal with it in the CLI.It would be nice to have an option to add new spectral lines on the fly. <br>
Need to think of how to deal with it in the CLI.Vigeesh GangadharanVigeesh Gangadharanhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/18Sunpy Map doesn't understand Wavelength info2021-03-24T13:59:17ZCarl SchafferSunpy Map doesn't understand Wavelength info```
In [45]: m[0]
Out[45]:
<sunpy.map.mapbase.GenericMap object at 0x7f55d741f0d0>
SunPy Map
---------
Observatory: Teide Observatory
Instrument: GRIS
Detector:
Measurement: 0.0 Angstrom
Wavelength: 0.0 Angstrom
Observation ...```
In [45]: m[0]
Out[45]:
<sunpy.map.mapbase.GenericMap object at 0x7f55d741f0d0>
SunPy Map
---------
Observatory: Teide Observatory
Instrument: GRIS
Detector:
Measurement: 0.0 Angstrom
Wavelength: 0.0 Angstrom
Observation Date: 2020-09-01 08:57:01
Exposure Time: 0.000000 s
Dimension: [10. 10.] pix
Coordinate System: helioprojective
Scale: [-3.75e-05 -3.75e-05] deg / pix
Reference Pixel: [-1. -1.] pix
Reference Coord: [ 0.02302444 -0.03129083] deg
array([[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
[0., 0., 0., 0., 0., 0., 0., 0., 0., 0.]])
```
Printing a sunpy Map instance on the terminal seems like it has problem understanding some of the information from the header. the `Measurement`, `Wavelength` and `Exposure Time` fields are not filled in the instance, I suspect they haven't been understood.
Reproduce this issue by using the file from importer_test_data as a map inputhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/21Fix image flip2022-02-24T07:52:36ZVigeesh GangadharanFix image flipVigeesh GangadharanVigeesh Gangadharanhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/22ds9 error messages2021-03-25T16:17:06ZVigeesh Gangadharands9 error messagesdifferent versions of ds9 throw different errors for the same file.
**Verison 7.3.2**
```
!! AST: Error at line 3509 in file fitsimage.C.
! astRead(FitsChan): A longitude axis ('HPLN-TAN') was found without a corresponding latitude ax...different versions of ds9 throw different errors for the same file.
**Verison 7.3.2**
```
!! AST: Error at line 3509 in file fitsimage.C.
! astRead(FitsChan): A longitude axis ('HPLN-TAN') was found without a corresponding latitude axis.
```
**Versionn 8.1**
```
!! AST: Error in routine wcsSkyFrame at line 37 in file frame/wcsast.C.
! astSet(FrameSet): Cannot maintain FrameSet integrity.
! Unable to accommodate the attribute setting "System=FK5, Equinox=J2000".
```Vigeesh GangadharanVigeesh Gangadharanhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/32Issue with static linking of mpi2021-04-22T09:27:48ZVigeesh GangadharanIssue with static linking of mpiby default, I see the following links on the complied code.
```
libopenblas.so.0 => /home/vigeesh/conda/envs/gris_env/lib/libopenblas.so.0 (0x00007fbbe7119000)
libmpi.so.12 => /home/vigeesh/conda/envs/gris_env/lib/libmpi.so.12 (0x00007fb...by default, I see the following links on the complied code.
```
libopenblas.so.0 => /home/vigeesh/conda/envs/gris_env/lib/libopenblas.so.0 (0x00007fbbe7119000)
libmpi.so.12 => /home/vigeesh/conda/envs/gris_env/lib/libmpi.so.12 (0x00007fbbe65ad000)
libmpifort.so.12 => /home/vigeesh/conda/envs/gris_env/lib/libmpifort.so.12 (0x00007fbbe9240000)
libgfortran.so.5 => /home/vigeesh/conda/envs/gris_env/lib/libgfortran.so.5 (0x00007fbbe640d000)
```
but if I have mpi loaded, even before I compile.
```
module load mpi/openmpi-x86_64
```
then, `mpi` seems to not use the conda path.
```
libopenblas.so.0 => /home/vigeesh/conda/envs/gris_env/lib/libopenblas.so.0 (0x00007effb110c000)
libmpi.so.12 => /usr/lib64/openmpi/lib/libmpi.so.12 (0x00007effb0e28000)
libmpifort.so.12 => /home/vigeesh/conda/envs/gris_env/lib/libmpifort.so.12 (0x00007effb3233000)
libgfortran.so.5 => /home/vigeesh/conda/envs/gris_env/lib/libgfortran.so.5 (0x00007effb0c88000)
```
I guess, this a problem with `mpif90` provided by `mpi4py`, although it seems to use proper `rpath`
the default, shows,
```
$ /usr/lib64/openmpi/bin/mpif90 -showme
gfortran -I/usr/include/openmpi-x86_64 -pthread -m64 -I/usr/lib64/openmpi/lib -Wl,-rpath -Wl,/usr/lib64/openmpi/lib -Wl,--enable-new-dtags -L/usr/lib64/openmpi/lib -lmpi_usempi -lmpi_mpifh -lmpi
```
the conda one, shows,
```
$ /home/vigeesh/conda/envs/gris_env/bin/mpif90 -show
x86_64-conda-linux-gnu-gfortran -I/home/vigeesh/conda/envs/gris_env/include -L/home/vigeesh/conda/envs/gris_env/lib -Wl,-rpath,/home/vigeesh/conda/envs/gris_env/lib -I/home/vigeesh/conda/envs/gris_env/include -I/home/vigeesh/conda/envs/gris_env/include -L/home/vigeesh/conda/envs/gris_env/lib -lmpifort -Wl,-rpath -Wl,/home/vigeesh/conda/envs/gris_env/lib -Wl,--enable-new-dtags -lmpi
```Vigeesh GangadharanVigeesh Gangadharanhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/44[IFU-integration] Header does not change for timeseries data2022-02-02T16:48:34ZVigeesh Gangadharan[IFU-integration] Header does not change for timeseries dataFor timeseries data, the header is temporarily fixed to the first map. The axes of the preview (animated) gif are fixed and does not change with the FOV. This needs to be corrected.For timeseries data, the header is temporarily fixed to the first map. The axes of the preview (animated) gif are fixed and does not change with the FOV. This needs to be corrected.Vigeesh GangadharanVigeesh Gangadharanhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/46add proper BTYPE label for uncertainties2022-01-11T10:31:49ZVigeesh Gangadharanadd proper BTYPE label for uncertaintiesVigeesh GangadharanVigeesh Gangadharanhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/50Add an option to apply user-supplied weights2022-03-10T11:38:55ZVigeesh GangadharanAdd an option to apply user-supplied weightsRelated to `inv_utils.f90`:
```fortran
42 WEIGHTS(:)=NOISE
43 !-----------------------------------------------------------
44 WEIGHTS(1)=(2.0*(1.0+4.0*RATUMQS-ICVAL/ICMEAN)/(1.0+RATUMQS*ICVAL/ICMEAN))
45...Related to `inv_utils.f90`:
```fortran
42 WEIGHTS(:)=NOISE
43 !-----------------------------------------------------------
44 WEIGHTS(1)=(2.0*(1.0+4.0*RATUMQS-ICVAL/ICMEAN)/(1.0+RATUMQS*ICVAL/ICMEAN))
45 IF (WEIGHTS(1) > 2.5) WEIGHTS(1)=2.5
46 IF (WEIGHTS(1) < 0.5) WEIGHTS(1)=0.5
47 WEIGHTS(1)=WEIGHTS(1)*NOISE
```
Try to get this out of Fortran into PythonVigeesh GangadharanVigeesh Gangadharanhttps://gitlab.leibniz-kis.de/sdc/gris/grisinv/-/issues/51Array size mismatch2022-07-27T09:45:54ZVigeesh GangadharanArray size mismatch```
SDC:GRIS Inversion Pipeline: v0.6.0+4.gf83d1e4
Reading fits files 0/128: | | 0% (0.0s)
Traceback (most recent call last):
File "/home/vigeesh/conda/envs/gris_env_test/bin/vfisv", line 8, in <module>
...```
SDC:GRIS Inversion Pipeline: v0.6.0+4.gf83d1e4
Reading fits files 0/128: | | 0% (0.0s)
Traceback (most recent call last):
File "/home/vigeesh/conda/envs/gris_env_test/bin/vfisv", line 8, in <module>
sys.exit(main())
File "/home/vigeesh/conda/envs/gris_env_test/lib/python3.9/site-packages/click/core.py", line 829, in __call__
return self.main(*args, **kwargs)
File "/home/vigeesh/conda/envs/gris_env_test/lib/python3.9/site-packages/click/core.py", line 782, in main
rv = self.invoke(ctx)
File "/home/vigeesh/conda/envs/gris_env_test/lib/python3.9/site-packages/click/core.py", line 1066, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/vigeesh/conda/envs/gris_env_test/lib/python3.9/site-packages/click/core.py", line 610, in invoke
return callback(*args, **kwargs)
File "/home/vigeesh/conda/envs/gris_env_test/lib/python3.9/site-packages/grisinv/invert.py", line 2198, in main
data, stokes, header = vfisv(path, id, line, width,numproc,weights=weights)
File "/home/vigeesh/conda/envs/gris_env_test/lib/python3.9/site-packages/grisinv/invert.py", line 214, in vfisv
vfisv_data = VFISVpackage(path, id, line, width)
File "/home/vigeesh/conda/envs/gris_env_test/lib/python3.9/site-packages/grisinv/invert.py", line 1669, in __init__
self.get_data_slit()
File "/home/vigeesh/conda/envs/gris_env_test/lib/python3.9/site-packages/grisinv/invert.py", line 1921, in get_data_slit
self.SI[iff, ...],
ValueError: could not broadcast input array from shape (424,132) into shape (424,201)
```Vigeesh GangadharanVigeesh Gangadharan