Hi,
There is the same/similar problem on the Matlab side (as there is no
writing and reading of xml-files including the time group).
So a general hint for all. If you just want to set sensor_time to some
value, you can do this by these method calls
timeNow
ArrayOfTimeSetConstant( sensor_time, 1, time )
Bye,
Patrick
On 2021-02-16 14:53, eric.sauvag...@iap.unibe.ch wrote:
Dear ARTS community,
In order to try a recent fix suggested and implemented (commit n.
5951a72fbef6d7ac9a3de8d2c503b58ef5af7d17, fromTyphon mailing list) on
the development version of ARTS I came into another problem trying to
implement my ground-based retrievals through PyARTS. In fact, I can't
succeed anymore in setting the sensor_time variable through PyARTS,
which is then required to perform the retrievals.
I'm not sure if something got mixed maybe between the 2 ARTS version
installed on my laptop and the 2 attached PyARTS versions (each in a
separate conda env) or if this is really a bug from the development
version.
As a very short example, please find below a minimalist python script
that reproduces the problem. It seems to work fine if I use ARTS2.4 and
produces the following error with ARTS Dev version:
###############################################################################
import os
import numpy as np
import xarray as xr
import matplotlib.pyplot as plt
from pyarts.workspace import Workspace, arts_agenda
from dotenv import load_dotenv
#load_dotenv('/home/eric/Documents/PhD/ARTS/arts-examples/.env.t490-arts2.5')
load_dotenv('/home/eric/Documents/PhD/ARTS/arts-examples/.env.t490-arts2.4')
ARTS_DATA_PATH = os.environ['ARTS_DATA_PATH']
ARTS_BUILD_PATH = os.environ['ARTS_BUILD_PATH']
ARTS_INCLUDE_PATH = os.environ['ARTS_INCLUDE_PATH']
def test_sensor():
# Initializing Workspace object
ws = Workspace(verbosity=0, agenda_verbosity=0)
ws.execute_controlfile("general/general.arts")
ws.execute_controlfile("general/agendas.arts")
ws.execute_controlfile("general/continua.arts")
ws.execute_controlfile("general/planet_earth.arts")
ws.sensor_los = np.array([50])
ws.sensor_pos = np.array([1e3])
ws.sensor_time = np.array([1])
print('sensor_time is :', ws.sensor_time.value)
if __name__=="__main__":
test_sensor()
###############################################################################
_______________________________________________
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
_______________________________________________
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
