How to calculate PET, SPI, SPEI and Palmer drought indices in Python? | Jalal | Research Society

[@DrJalal90]

Data download link: www.drought.gov/data-maps-tools/climate-and-drought-indices-python-spi-spei-pet Alternative link: github.com/monocongo/climate_indices

[@travelwonder]

Hi Dr. Jalal, The link provided for data have been obsolete. Can you provide a new link where the data for tutorial can be found? Thanks

[@supercomedy86]

thank you and much appreciated Dr.Lal

[@DrJalal90]

Thanks

[@khairullahkhairullah205]

Thank you very much Jalal for the video. It's useful for me.

[@agirl3517]

Hey did your command run?

[@SANTHOSHKUMAR-rx3eo]

Sir the zip file is not showing in that link, it was removed

[@MrMandarpriya]

Sir the data link is not available, so how to download the data.

[@agirl3517]

Hello sir, can we open this output in ArcMap and if yes, how will it look in spatial format, since we can open netcdf in ArcMap

[@lixiangAn-t8q]

Hello Dr. Jalal.Thank you for this great tutorial.Can you make a video how to claculate sri in netcdf.

[@DrJalal90]

Thanks. I will try.

[@lixiangAn-t8q]

Thank you! Looking forward to your update, and thank you very much for your video@@DrJalal90

[@Hydroclim]

Thank you Dr. Jalal for the great video. I was just wondering about the calibration start year and calibration end year. Could you please provide a bit more insight on them? As an example, I want to calculate spei for 1970-2010 dataset. So my calibration start and end year would be 1970 and 2010. I wonder if it would be needed to use some of this data for calibration, such as 1970-1980.

[@mayanktyagi9022]

Respected sir, Sir I want to calculate the SPEI of a basin . Should it be necessary to calculate the grided SPEI of that basin, or any another method of calculation of the SPEI.

[@DrJalal90]

You may use station data. Here is a software so called climpact2: climpact-sci.org/assets/climpact2-user-guide.pdf

[@林语丛]

hi can I ask why I get the error? UserWarning: resource_tracker: There appear to be 8 leaked semaphore objects to clean up at shutdown warnings.warn('resource_tracker: There appear to be %d '

[@DrJalal90]

Hello dear, I have no idea.

[@olgasolvas7751]

Hello! Thank you for sharing this video. I found it really informative and I appreciate the effort you put into creating it. However, I wanted to let you know that the links you have provided in the video description do not seem to be working properly. When I clicked on the links, they did not download the same file as the one you mentioned in the video, and they did not work as expected. I hope you share the source file, as I am really interested in accessing the resources you mentioned. Thank you again for sharing your expertise with us.

[@DrJalal90]

climate-indices.readthedocs.io/en/latest/

[@hemraj977]

So, when we calculate SPEI and just give one input (TEMP), but not precipitation, how would it calculate drought just from temperature...or do we need to give bunch of variables?

[@Mr.Goetia]

The (Data Download Link) is not working! The page is not available.

[@DrJalal90]

github.com/monocongo/climate_indices

[@Mr.Goetia]

@@DrJalal90 Thank you very much!

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@muhammadasyrofialkindy8066]

Thank you sir

[@christophergama1588]

The dataset containing precipitation values should it be only in txt file format.? or other format can also work such as CSV file etc

[@travelwonder]

Hi Mr Jalal, can you again check if this code is working fine. I am using the procedure for CHIRPS netcdf data and code is reading data and starting calculation but then giving error messages. I have checked the data carefully and netcdf data is all fine. I will be thankful if you could check again and guide me. Thanks

[@Chandramani106]

Please make of video calculating the same parameters using Jupyter notebok platform

[@DrJalal90]

I will try.

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@mahsabozorgi]

Thank you very much for the useful video and channel. I have a question. In command line I put the full address of input and output data. I faced this error "--output_file_base: command not found" Do you know the solution?Thanks

[@ibrahimerkol4511]

Thanks for the video, it's super helpful. Yet, I couldn't reach out to the zip file. The link you provided for data download in the description says page or link not found. Could anyone please share with the zip file? I couldn't find it anywhere at all. Thanks in advance

[@DrJalal90]

github.com/monocongo/climate_indices

[@ibrahimerkol4511]

Thank you so much for your return, yet the files in the link do not match with yours in the video :(((

[@DrJalal90]

climate-indices.readthedocs.io/en/latest/

[@musafirgabut8053]

@@DrJalal90 ​ @Dr. Jalal Sorry mr Jalal, i couldn't reach out the zip file from your current link. Please share again. Thanks a lot.

[@samirahasani9508]

hi jalal, when i use the given example data the code work, but when i use my own precipitation data (cmip6) i face with error. moreover cdo doesn't know the structure of given precip data in example, would you pls explain?

[@DrJalal90]

Please send me the screenshot of CDO.

[@tadessemosissa7226]

I am very interested ur video. could send me the codes used for the process climate indices calculation in your video.

[@DrJalal90]

Please check the video description.

[@agirl3517]

In your " A comprehensive statistical assessment of drought indices to monitor drought status in Bangladesh" u have used spei using R . Why didn't you use it using this python software when it was spei was available? Also I am unable to download this python software for spi and spei by NIDIS now, has it been removed or is the link change? Please help me with this

[@DrJalal90]

www.drought.gov/data-maps-tools/climate-and-drought-indices-python-spi-spei-pet

[@DrJalal90]

github.com/monocongo/climate_indices

[@agirl3517]

Sir thankyou for the package and video. I am getting this error, process_climate_indices is not recognised, please sir help me with this problem I have tried changing path. But it is still showing same error

[@yahayaibrahim6790]

Hi, Mr Jalal can u please show us how to calculate PET,SPEI using penman monteith method. Tnx

[@ibrahimerkol4511]

Thanks again for the video, I tried to run the code but I get "FileNotFoundError: [Errno 2] No such file or directory" type of an error. Could you please help me solve the problem dear professor. I already changed directory for the file in which the NC file is located. Thanks in advance

[@DrJalal90]

You need to check the working directory carefully.

[@riezkymaisyar4038]

hallo sir, i have question, how i change the dimention my file nc using nco on windows? because my data have parameter like precip(time, lat, lon) instead like precip(lat, lon, time), how i can change the order? thank you

[@cherryleheu3985]

Hello! Do you know of any Python packages that takes a .csv file as an input and calculates SPI using a calibration period?

[@DrJalal90]

The following video might help you. How to calculate drought index SPI | Jalal | Research Society kzbin.info/www/bejne/jYXPoqqnZrBlb8U

[@yahayaibrahim6790]

Is very helpful and educative, The link you put to download the indices_python_beta_v1.2.zip is not responding? Any other alternative. Tnx

[@DrJalal90]

Yes, please check the highlighted comment.

[@lankz2050]

Thank you sir...👌

[@DrJalal90]

Welcome

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@valeriabedoyapineda7574]

Hello, thank you very much for the explanation. I have a problem I get the error "ValueError: Unable to calculate Pearson Type III parameters due to invalid L-moments", do you know how to solve it?

[@DrJalal90]

Sorry, I don’t know.

[@songtanglei7539]

great video，but when my nc file not the same as that in the example，error occurs，anyone can tell how can i change the attributes in my pre data。

[@DrJalal90]

No ideas.

[@Hydroclim]

I have faced the same issue. Please let me know if you found a solution.

[@Hydroclim]

I managed to process by changing attributes of my netcdf file (chirps monthly precipitation) using jupyter notebook. #open and slice your netcdf file pr_month = xr.open_dataset('/Precipitation/chirps-v2.0.monthly.nc') pr_month = pr_month['precip'] pr_month = pr_month.sel(longitude=slice(71,73),latitude=slice(31,33),time=slice('1981-01-01','2010-01-01')) #changing attributes pr_month = pr_month.rename({'latitude':'lat','longitude':'lon','time':'time'}) pr_month.attrs["units"]='millimeter' preferred_dims = ('lat', 'lon', 'time') pr_month = pr_month.transpose(*preferred_dims) #save update netcdf file pr_month.to_netcdf('Precipitation/processed_pr_monthly.nc')

[@jeromeloukou8059]

allo I have one problem with my precipitation data can you help me

[@DrJalal90]

What is your problem?

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@agumasekindie7672]

Dear i can not get indices_python-v1.2-beta.zip. how i can get it ?

[@DrJalal90]

github.com/monocongo/climate_indices

[@agumasekindie7672]

@@DrJalal90 thank you so much, i can't get exactly 'indices_python--beta_v1.2zip ' . climate indices found in 'github.com/monocongo/climate_indices' is not similar to what you have shown in your video.i kindly ask you to share 'indices_python--beta_v1.2zip,

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@chandrakantkadam1953]

if SPI calculator in Excel format is available? Pls, share just to understand the calculations.

[@DrJalal90]

I will inform you if I find it.

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@hilmaniipaya258]

How do I convert my timeseries from a csv file to a netcdf file?

[@DrJalal90]

stackoverflow.com/questions/56293164/convert-time-series-data-from-csv-to-netcdf-python

[@DrJalal90]

You can do it in NCL. www.ncl.ucar.edu/Training/Workshops/KState/Scripts/CSV_Files/

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@local9692]

I keep getting this error. ImportError: DLL load failed while importing _arpack: The specified procedure could not be found.

[@DrJalal90]

What was your command?

[@local9692]

@@DrJalal90 $ process_climate_indices --index palmers --periodicity monthly --netcdf_precip ./nclimgrid_lowres_prcp.nc --var_name_precip prcp --netcdf_pet ./nclimgrid_lowres_pet.nc --var_name_pet pet --netcdf_awc ./nclimgrid_lowres_soil.nc --var_name_awc awc --output_file_base ../results/nclimgrid_lowres --calibration_start_year 1951 --calibration_end_year 2010 Traceback (most recent call last): File "c:\users\bikra\anaconda3\lib unpy.py", line 194, in _run_module_as_main return _run_code(code, main_globals, None, File "c:\users\bikra\anaconda3\lib unpy.py", line 87, in _run_code exec(code, run_globals) File "C:\Users\bikra\anaconda3\Scripts\process_climate_indices.exe\__main__.py", line 4, in File "c:\users\bikra\anaconda3\lib\site-packages\climate_indices\__main__.py", line 14, in from climate_indices import compute, indices, utils File "c:\users\bikra\anaconda3\lib\site-packages\climate_indices\compute.py", line 9, in import scipy.stats File "c:\users\bikra\anaconda3\lib\site-packages\scipy\stats\__init__.py", line 388, in from .stats import * File "c:\users\bikra\anaconda3\lib\site-packages\scipy\stats\stats.py", line 180, in from . import distributions File "c:\users\bikra\anaconda3\lib\site-packages\scipy\stats\distributions.py", line 8, in from ._distn_infrastructure import (entropy, rv_discrete, rv_continuous, File "c:\users\bikra\anaconda3\lib\site-packages\scipy\stats\_distn_infrastructure.py", line 23, in from scipy import optimize File "c:\users\bikra\anaconda3\lib\site-packages\scipy\optimize\__init__.py", line 387, in from .optimize import * File "c:\users\bikra\anaconda3\lib\site-packages\scipy\optimize\optimize.py", line 36, in from ._numdiff import approx_derivative File "c:\users\bikra\anaconda3\lib\site-packages\scipy\optimize\_numdiff.py", line 6, in from scipy.sparse.linalg import LinearOperator File "c:\users\bikra\anaconda3\lib\site-packages\scipy\sparse\linalg\__init__.py", line 114, in from .eigen import * File "c:\users\bikra\anaconda3\lib\site-packages\scipy\sparse\linalg\eigen\__init__.py", line 9, in from .arpack import * File "c:\users\bikra\anaconda3\lib\site-packages\scipy\sparse\linalg\eigen\arpack\__init__.py", line 20, in from .arpack import * File "c:\users\bikra\anaconda3\lib\site-packages\scipy\sparse\linalg\eigen\arpack\arpack.py", line 43, in from . import _arpack ImportError: DLL load failed while importing _arpack: The specified procedure could not be found.

[@Chandramani106]

What does the 'scale 30 90' or 'scales 6 12' signify? Does it signify the spatial grid or the time scale?

[@DrJalal90]

30-day and 90-day timescales climate-indices.readthedocs.io/en/latest/#spi-monthly

[@Chandramani106]

@@DrJalal90 thank you bhai

[@jorgearmandohernandezlopez8080]

Thank Jalal's, I have information monthly in file Raster but i dont know convertion in nc. Can you help me?. I work in arcgis but i cant transformation file. Thank

[@DrJalal90]

In ArcGIS, there is an option to save raster to netcdf file.

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@tugbaozturk1664]

Thank you very much for the video. What about calculating PET by using Hargreaves equation?

[@DrJalal90]

Sorry, I am not familiar with Hargreaves equation.

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@shammunulislam4235]

If you use R, you can use CDT package to calculate PET using Hargreaves equation

[@arkchannel7830]

please continue this video and done SPEI calculation for us

[@DrJalal90]

Everything is there. Please follow the instructions.

[@arkchannel7830]

@@DrJalal90 thank you very much

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@rohitprasadshaw1032]

Sir can I calculate ARIDITY INDEX using python ? help me plz🙏🙏 from INDIA

[@DrJalal90]

I have no idea.

[@mohd.sayeedulhasan8902]

great

[@DrJalal90]

Thank you.

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@pawansoyam6702]

Please make video on extraction and export into csv of MODIS level 2 AOD data in Python

[@sin168gh]

sir can you please explain what are times scales 1,3,6,12 months..please solve my problem

[@DrJalal90]

doi.org/10.1175/2009JCLI2909.1

[@DrJalal90]

Please read the paper carefully.

[@DrJalal90]

Do you understand Bengali?

[@DrJalal90]

মডেল ভিত্তিক Drought গবেষণা এবং দক্ষিন এশিয়ার উপর এর ভবিষ্যত ক্ষতিকর প্রভাব kzbin.info/www/bejne/qnnWaGqFadhsrKc

[@sin168gh]

@@DrJalal90no sir I don't understand Bengali

[@Mr.2-3]

Mr, what you have shown here?? you just run some codes and showing result as a file. where is the plot, how to visualize, what are the values in results, how it distributed, how to show these results in a map???? there are many unanswered questions in your video. without all of these, how it become helpful for a researcher or student to learn about climate indices?????

[@kipishism]

This is really dope! Unfortunately the zip file has been removed.......is there a way you can share the file you downloaded?

[@DrJalal90]

github.com/monocongo/climate_indices

[@kipishism]

@@DrJalal90 thanks dude!

[@kipishism]

Unfortunately I am only able to find data on the notebooks directory, this does not conform to the example you shared.....I am at a loss ;(

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@ibrahimerkol4511]

@@DrJalal90 could you please share the zip file with us?

[@rosscherin7051]

Hi Jalal! Thank you for this great tutorial. However, I am getting an error while calculating the SPI at the end. "File "c:\users\dell\anaconda3\lib\site-packages\xarray\core\variable.py", line 1068, in chunk data = da.from_array(data, chunks, name=name, lock=lock, **kwargs) TypeError: from_array() got an unexpected keyword argument 'meta'" Can you please suggest to me what I should do to fix it?

[@DrJalal90]

Please check the following link. www.gitmemory.com/issue/monocongo/climate_indices/435/866770257

[@mohammedkareem6476]

thanks

[@DrJalal90]

Welcome

[@DrJalal90]

Python free online training for international participants The training will start on 20 September 2021. Deadline for registration: 15 September 2021 Registration link: researchsociety20.org/registration/ Download link for training outline: www.dropbox.com/s/aybmuw0rrxu7zjb/Outline_python_International.png?dl=0 Event link: fb.me/e/XMyJYOwn Please share the event with your friends or colleagues.

[@mohammedkareem6476]

@@DrJalal90 thanks

[@agumasekindie7672]

how i can solve it, i am using Get bash in window 10 Traceback (most recent call last): File "process_grid.py", line 10, in from climate_indices import indices, utils File "C:\Users\TOSHIBA\anaconda3\lib\site-packages\climate_indices\indices.py", line 8, in from climate_indices import compute, eto, palmer, utils File "C:\Users\TOSHIBA\anaconda3\lib\site-packages\climate_indices\compute.py", line 8, in import scipy.special File "C:\Users\TOSHIBA\anaconda3\lib\site-packages\scipy\special\__init__.py", line 633, in from . import _ufuncs ImportError: DLL load failed while importing _ufuncs: The specified module could not be found.

[@DrJalal90]

Please check the following link. github.com/marcelotduarte/cx_Freeze/issues/355

[@ElisaFrankBuss]

Hello, I tried this easy way to calculate spi, but y received this error: "File "/home/usuario/anaconda3/envs/chirps/lib/python3.11/site-packages/xarray/core/utils.py", line 281, in either_dict_or_kwargs raise ValueError(f"the first argument to .{func_name} must be a dictionary") ValueError: the first argument to .drop_sel must be a dictionary". Someone can tell me what is wrong? Thank you!

[@changliu6799]

Hi! Thanks for sharing this video. Can you tell me what wrong is with my code? Thanks! 2022-11-02 19:28:59 INFO Start time: 2022-11-02 19:28:59.991679 Error in calling operator ncpdq with: >>> D:\anaconda3\envs\indices_env\Library\bin cpdq -a lat,lon,time -O --output=C:\Users\Chang\AppData\Local\Temp cpdq_pr_Amon_FGOALS-f3-L_ssp370_r1i1p1f1_gr_201901-205012_china.ncveh08oxy.tmp pr_Amon_FGOALS-f3-L_ssp370_r1i1p1f1_gr_201901-205012_china.nc

[@DrJalal90]

I have no idea.

[@thiandoumadama7438]

@@DrJalal90 process_climate_indices --index spei --periodicity monthly --netcdf_precip /media/espace-dev/CrucialX6/CMIP6-Merge/princeton-choix/prcp_daily_princeton_2C_WA_1981-2016.mon.nc --var_name_precip prcp --netcdf_pet /media/espace-dev/CrucialX6/CMIP6-Merge/princeton-choix/pet_data_princeton_1981-2016.mon.nc --var_name_pet pet --calibration_start_year 1981 --calibration_end_year 2016 --output_file_base /media/espace-dev/CrucialX6/CMIP6-Merge/princeton-choix/spei_output --scales 1 3 6 12 2024-10-23 15:31:09 INFO Start time: 2024-10-23 15:31:09.302110 2024-10-23 15:31:09 INFO Computing 1-month SPEI/Pearson 2024-10-23 15:31:09 ERROR Invalid dimensions for variable 'prcp': ('time', 'lat', 'lon') 2024-10-23 15:31:09 ERROR Failed to complete Traceback (most recent call last): File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 1717, in process_climate_indices _compute_write_index(kwrgs) File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 803, in _compute_write_index output_shape = _drop_data_into_shared_arrays_grid( File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 581, in _drop_data_into_shared_arrays_grid raise ValueError(message) ValueError: Invalid dimensions for variable 'prcp': ('time', 'lat', 'lon') 2024-10-23 15:31:09 ERROR Failed to complete Traceback (most recent call last): File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 1589, in main process_climate_indices(arguments=arguments) File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 1717, in process_climate_indices _compute_write_index(kwrgs) File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 803, in _compute_write_index output_shape = _drop_data_into_shared_arrays_grid( File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 581, in _drop_data_into_shared_arrays_grid raise ValueError(message) ValueError: Invalid dimensions for variable 'prcp': ('time', 'lat', 'lon') Traceback (most recent call last): File "/home/espace-dev/anaconda3/envs/mon_env/bin/process_climate_indices", line 8, in sys.exit(main()) File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 1589, in main process_climate_indices(arguments=arguments) File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 1717, in process_climate_indices _compute_write_index(kwrgs) File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 803, in _compute_write_index output_shape = _drop_data_into_shared_arrays_grid( File "/home/espace-dev/anaconda3/envs/mon_env/lib/python3.8/site-packages/climate_indices/__main__.py", line 581, in _drop_data_into_shared_arrays_grid raise ValueError(message) ValueError: Invalid dimensions for variable 'prcp': ('time', 'lat', 'lon')

[@cmontes]

does anyone have this error trying to run for SPEI? File "/var/folders/d6/m2bgqkw971v4qb3h1cq2tmmm0000gn/T/ipykernel_8047/4073413229.py", line 1 process_climate_indices --index spei --periodicity monthly --netcdf_precip ./data/nclimgrid_lowres_prcp.nc --var_name_precip prcp --netcdf_pet ./data/nclimgrid_lowres_pet.nc --var_name_pet pet --output_file_base ./output/nclimgrid_lowres --scales 9 18 --calibration_start_year 1951 --calibration_end_year 2010 ^ SyntaxError: invalid syntax

[@DrJalal90]

No