tsutterley / gravity-toolkit Goto Github PK

where to change or reshape the input

OSError Traceback (most recent call last)
Cell In [12], line 22
20 end_mon = np.max(months)
21 missing = sorted(set(np.arange(start_mon,end_mon+1)) - set(months))
---> 22 Ylms = gravtk.grace_input_months(widgets.base_directory, PROC, DREL, DSET,
23 LMAX, start_mon, end_mon, missing, SLR_C20, DEG1, MMAX=MMAX,
24 SLR_21=SLR_21, SLR_22=SLR_22, SLR_C30=SLR_C30, SLR_C40=SLR_C40,
25 SLR_C50=SLR_C50, POLE_TIDE=POLE_TIDE, ATM=ATM)
26 # create harmonics object and remove mean
27 GRACE_Ylms = gravtk.harmonics().from_dict(Ylms)

File ~/Seasonal/read-GRACE-harmonics-main/notebooks/gravity_toolkit/grace_input_months.py:781, in grace_input_months(base_dir, PROC, DREL, DSET, LMAX, start_mon, end_mon, missing, SLR_C20, DEG1, **kwargs)
779 if months_test:
780 gm = ','.join(f'{gm:03d}' for gm in months_test)
--> 781 raise IOError(f'No Matching Geocenter Months ({gm})')
782 # for each considered date
783 for i,grace_month in enumerate(months):

OSError: No Matching Geocenter Months (248,249)

Hello, when I read your file named gen_averaging_kernel.py, I need to input (gclm, gslm, eclm, eslm, sigma, hw) 6 parameters. I followed the method in Swenson et al. (2002) and made a little modification to the file named gen_harmonics.py to perform spherical harmonic expansion on the basin $\vartheta_{lm}^c,\vartheta_{lm} ^s$ . Here are some of my questions:
(1) When I use gen_harmonics.py to perform spherical harmonic expansion of the basin, the parameter is selected as METHOD='integration', and coeff = 1.0/(4.0*np.pi) is defined in your file, but it seems that this step is not required in the introduction of Swenson et al. (2002), so I removed it, and based on this, I obtained $(\vartheta_{lm}^c ,\vartheta_{lm}^s)$, I think this should be reasonable. Is this available as an input parameter (gclm, gslm) in the gen_averaging_kernel.py file?
(2) You don't seem to have made too many introductions and notes on the gen_harmonics function, and I haven't found it used in other scripts, so I don't know how to calculate your three input parameters (eclm, eslm, sigma)?
(3) By the way, both gen_disc_load.py and gen_spherical_cap.py seem to be closely related to Jacob et al. (2012). I briefly read this paper, but I can’t seem to find a complete derivation formula in the text to support my understanding of the code in your two files and the method used by your code. Can you provide me with some help, or briefly explain the functions of these two files? Where can I find a more complete description of these methods? These two files do not seem to affect my main direction of learning, but I still want to understand the principle.
(4) On the basis of (3), there is a question about mass coefficients and geoid coefficients. I seem to confuse these two. I always think that there is the following relationship between them. It is also possible that I put these two Both definitions are wrong, and the following formula refers to the notation in Wahr et al. (1998).

But the following code is used in your code to convert mass coefficients to geoid coefficients

# calculate SH degree dependent factors to convert from coefficients
# of mass into normalized geoid coefficients
dfactor = 4.0*np.pi*factors.cmwe/(1.0 + 2.0*factors.l)

This code does not seem to correspond to the above formula. Is it where I confused the concept, and other more forms of spherical harmonic coefficients? Can you help me answer it?

Hi there,

I am using VS Code in Windows 10, and when I run the setup.py file, I end up with the following issue regarding what looks like decoding problems? I am somewhat inexperienced with Python, so any help would be great.

../read-GRACE-harmonics/setup.py install
Traceback (most recent call last):
  File "..\read-GRACE-harmonics\setup.py", line 10, in <module>
    long_description = fh.read()
  File "C:\_Miniconda3\envs\gracefo\lib\encodings\cp1252.py", line 23, in decode
    return codecs.charmap_decode(input,self.errors,decoding_table)[0]
UnicodeDecodeError: 'charmap' codec can't decode byte 0x90 in position 3118: character maps to <undefined>

cnes_grace_sync.py worked with the prior download service access point (which is now non-operational as of 2023/05/04). Need to update to work with the new data access catalogue.

TODO:

• authentication
• parsing or querying the catalogue
• download

Hello, when I read your file named gen_averaging_kernel.py, I need to input (gclm, gslm, eclm, eslm, sigma, hw) 6 parameters. I followed the method in Swenson et al. (2002) and made a little modification to the file named gen_harmonics.py to perform spherical harmonic expansion on the basin $\vartheta_{lm}^c,\vartheta_{lm} ^s$ . Here are some of my questions:
(1) When I use gen_harmonics.py to perform spherical harmonic expansion of the basin, the parameter is selected as METHOD='integration', and coeff = 1.0/(4.0*np.pi) is defined in your file, but it seems that this step is not required in the introduction of Swenson et al. (2002), so I removed it, and based on this, I obtained $(\vartheta_{lm}^c ,\vartheta_{lm}^s)$, I think this should be reasonable. Is this available as an input parameter (gclm, gslm) in the gen_averaging_kernel.py file?
(2) You don't seem to have made too many introductions and notes on the gen_harmonics function, and I haven't found it used in other scripts, so I don't know how to calculate your three input parameters (eclm, eslm, sigma)?
(3) By the way, both gen_disc_load.py and gen_spherical_cap.py seem to be closely related to Jacob et al. (2012). I briefly read this paper, but I can’t seem to find a complete derivation formula in the text to support my understanding of the code in your two files and the method used by your code. Can you provide me with some help, or briefly explain the functions of these two files? Where can I find a more complete description of these methods? These two files do not seem to affect my main direction of learning, but I still want to understand the principle.
(4) On the basis of (3), there is a question about mass coefficients and geoid coefficients. I seem to confuse these two. I always think that there is the following relationship between them. It is also possible that I put these two Both definitions are wrong, and the following formula refers to the notation in Wahr et al. (1998).

But the following code is used in your code to convert mass coefficients to geoid coefficients

# calculate SH degree dependent factors to convert from coefficients
# of mass into normalized geoid coefficients
dfactor = 4.0*np.pi*factors.cmwe/(1.0 + 2.0*factors.l)

This code does not seem to correspond to the above formula. Is it where I confused the concept, and other more forms of spherical harmonic coefficients? Can you help me answer it?

I would like to calculate the mass change of the mascon, but I am not very familiar with how to generate the necessary mascon file. I am also not entirely clear about the data input required in the mascon.py program. Would you mind helping me understand it better and provide an explanation?

Hello, when I use the GRACE-Harmonic-Plots.ipynb file under the gravity_toolkit/notebooks folder, when I go to Section 1.7 (Set Parmeter to Visualize Harmonics), I need to set four parameters in this section for importing and formatting some files (the comment is "removed from the GRACE/ GRACE-FO data", I don't understand what was removed), these parameters include widgets.remove_file, widgets.remove_format, widgets.redistribute_removed and widgets.mask, what these files are and what they do, this seems to be related to the recovery signal leakage in section 1.8, but I can't fully understand the role of these files, especially the lines 37~70 in section 1.8, can you help me understand what these files are and how to set these parameters? If so, can you provide some supporting literature related to it, thank you!

https://podaac.jpl.nasa.gov/cloud-datasets/migration

Hello, in the grace_input_months function, why do CS21/CS22/C30/C40/C50 set the grace_month greater than 179 (in the source code, grace_month > 176) to replace the GRACE observations with SLR observations? Where is it? It is mentioned in the literature, can you answer or provide a link to the literature? Thanks!

Hello, why is there a crash when running import gravity_toolkit as gravtk, Process finished with exit code -1073741819 (0xC0000005)

ModuleNotFoundError: No module named 'read_GRACE_geocenter
everytime i run the code .... i get the above error.
when I searched for "read_GRACE_geocenter" it was not their in your gravity toolkit folder... i think that is why it is giving error.
please have look over your code and pls instruct me how to run it. my mtech disseratation depends on it and I am close to the deadline

Hello! I tried to stock up on demo examples, but didn't get far: !podaac_cumulus.py --catalog=$widgets.the base directory
Do I understand correctly that this script is not working at the moment?

The link at the end:Tracking (last last call):
File "/usr/local/bin/podaac_cumulus.py ", line 444, in
main()
File "/usr/local/bin/podaac_cumulus.py ", line 436, in the main
podaac_cumulus(client, args.directory, PROC=args.center,
file "/usr/local/bin/podaac_cumulus.py ", line 135, in podaac_cumulus
url, = [url for the url in urls if R1.search(url)]
ValueError error: not enough values to decompress (1 expected, 0 received)
INFO: root:Grace Date Program: CSR RL06 GAC
Backtracking (last call):
File "/usr/local/bin/run_grace_date.py ", line 165, in
main()
File "/usr/local/bin/run_grace_date.py ", line 160, in the main
file run_grace_date(args.directory, args.center, args.release,
file "/usr/local/bin/run_grace_date.py ", line 112, in the run_grace_date file
gravtk.grace_date(base index, PROC=p, DREL=r, DSET=d,
File "/usr/local/lib/python3.10/dist-packages/gravity_toolkit/grace_date.py ", line 165, in grace_date
causes the FileNotFoundError error(f'{str(index_file)} not found")
FileNotFoundError error: /content/CSR/RL06/GAC/index.txt not found

temp.month = np.int64(calendar_to_grace(temp.time))