Bug with Bottles matchup

[gbolzon]

La profilazione genera profili per il 2019, ma NON per la data 2019-03-18

Poi il NutrientsReader.Selector(None, 2019,Med) trova profili in quella data ==> crash

Il profilatore originale è
/g100_scratch/userexternal/camadio0/RUN_FIUMI_ATM/SIMU_8_v1/wrkdir/POSTPROC/ogstm_postproc/Rivers/profiler_RA_N_bottle.py
che fa riferimento a
/g100_work/OGS_devC/camadio/RUN_FIUMI_ATM/VALIDATION_DATASETS_VDB/BOTTLE_ONLY/MedBGCins_nut_v2.nc

BITSEA=
/g100_work/OGS_test2528/camadio/RUN_FIUMI_ATM/archieve/CLIM/wrkdir/POSTPROC/bit.sea_BOTTLE

Fino a 15-20 giorni fa funzionava, non dava errori.
L'errore si genera a riga 62 di
/g100_scratch/userexternal/camadio0/RUN_FIUMI_ATM/SIMU_8_v1/wrkdir/POSTPROC/ogstm_postproc/Rivers/river_valid_BOTTLE.sh

L'errore è generato da

bit.sea/src/bitsea/instruments/matchup_manager.py

Lines 136 to 152 in 98a64bf

	def modeltime(self,profile):
	'''
	Argument:
	* profile * a profile object.
	E.g. an element of list provided by FloatSelector(), or MooringSelector()
	Returns:
	* datetime object * , the time of model corresponding to the selected profile
	'''
	assert isinstance(profile, Profile)
	for Model_time,INTERESTED_Indices in self.Coupled_List:
	if profile in [self.PROFILE_LIST[k] for k in INTERESTED_Indices]:
	break
	else:
	raise ValueError("No model profile corresponds to %s. Please check time interval in your profiler." % ( profile.ID(), ) )
	return Model_time

[gbolzon]

Sospetti: python diverso tra generazione profili e validazione

[gbolzon]

Lavori in corso in /g100_scratch/userexternal/gbolzon0/debug

La sequenza col vecchio python funziona

source /g100_work/OGS23_PRACE_IT/COPERNICUS/sequence3.sh
export PYTHONPATH=/g100_work/OGS_test2528/camadio/RUN_FIUMI_ATM/archieve/CLIM/wrkdir/POSTPROC/bit.sea_BOTTLE/src
python profiler_RA_N_bottle.py

ma con il nuovo no ( non dipende dal PYTHONPATH)

source /g100_work/OGS23_PRACE_IT/COPERNICUS/sequence.sh
PYTHONPATH=/g100_work/OGS_test2528/camadio/RUN_FIUMI_ATM/archieve/CLIM/wrkdir/POSTPROC/bit.sea_BOTTLE/src
python profiler_RA_N_bottle.py

Traceback (most recent call last):
  File "/g100_scratch/userexternal/gbolzon0/debug/profiler_RA_N_bottle.py", line 28, in <module>
    ALL_PROFILES = N.Selector(None,T_INT, Rectangle(-6,36,30,46))
  File "/g100_work/OGS_test2528/camadio/RUN_FIUMI_ATM/archieve/CLIM/wrkdir/POSTPROC/bit.sea_BOTTLE/src/bitsea/static/Nutrients_reader.py", line 87, in Selector
    sublist=self.DataExtractor.selector(myvar, T_int, region)
  File "/g100_work/OGS_test2528/camadio/RUN_FIUMI_ATM/archieve/CLIM/wrkdir/POSTPROC/bit.sea_BOTTLE/src/bitsea/static/DatasetExtractor.py", line 126, in selector
    return self.profileGenerator(Time, Lon, Lat, values, depth, dataset)
  File "/g100_work/OGS_test2528/camadio/RUN_FIUMI_ATM/archieve/CLIM/wrkdir/POSTPROC/bit.sea_BOTTLE/src/bitsea/static/DatasetExtractor.py", line 49, in profileGenerator
    Values=values[inthisprofile]
  File "/g100_work/OGS23_PRACE_IT/COPERNICUS/py_env_3.9.18_new/lib/python3.9/site-packages/numpy/ma/core.py", line 3254, in __getitem__
    dout = self.data[indx]
IndexError: too many indices for array: array is 1-dimensional, but 2 were indexed

[gbolzon]

Di conseguenza rigiro tutto con il vecchio python - e funziona - e ottengo
/g100_scratch/userexternal/gbolzon0/debug/output/DENSITY_PLOTS/N3n/densplot_N3n_lgm.png

@CarolinaAmadio , vedo che il tuo
/g100_scratch/userexternal/camadio0/RUN_FIUMI_ATM/SIMU_8_v1/wrkdir/POSTPROC/ogstm_postproc/Rivers/river_valid_BOTTLE.sh
fa riferimento al nuovo python, e questo è un errore.

Tuttavia, come fai col nuovo python a eseguire
ALL_PROFILES = N.Selector(None,T_INT, Rectangle(-6,36,30,46)) ?
Perché non ti si pianta?