Skip to content

Commit

Permalink
Merge pull request #136 from GESS-research-group/joss-review-7550
Browse files Browse the repository at this point in the history 
Joss review 7550
  • Loading branch information
@lsmvivek
lsmvivek authored Jan 27, 2025
2 parents 219c8f6 + 29376b3 commit a7cdfb0
Show file tree
Hide file tree
Showing 8 changed files with 216 additions and 5,739 deletions.
13 changes: 6 additions & 7 deletions docs/theory.md
View file
Original file line number Diff line number Diff line change
Expand Up @@ -157,28 +157,27 @@ $k_n$ refers to the load Love number of degree $n$.
The `GRACE` data products are being developed, processed and archieved in a shared Science Data System between the `Jet Propulsion Laboratory(JPL)`, the `University of Texas Center for Space Research (UT-CSR)` and `GeoForschungsZentrum Potsdam (GFZ)`.

+ Level 0:<br>
The level-0 data are the result of the data reception, collection and decommutation by the Raw Data Center (RDC) of the Mission Operation System (MOS) located in Neustrelitz, Germany. The MOS receives twice per day using its Weilheim and Neustrelitz tracking antennae the science instrument and housekeeping data from each GRACE satellite which will be stored in two appropriate files in the level-0 rolling archive at DFD/Neustrelitz. The SDS retrieves these files and extracts and reformats the orresponding instrument and ancillary housekeeping data like GPS navigation solutions,space segment temperatures or thruster firing events. Level-0 products are available 24-hours after data reception.
The level-0 data are the result of the data reception, collection and decommutation by the Raw Data Center (RDC) of the Mission Operation System (MOS) located in Neustrelitz, Germany. The MOS receives twice per day using its Weilheim and Neustrelitz tracking antennae the science instrument and housekeeping data from each GRACE satellite which will be stored in two appropriate files in the level-0 rolling archive at DFD/Neustrelitz. The SDS retrieves these files and extracts and reformats the corresponding instrument and ancillary housekeeping data like GPS navigation solutions, space segment temperatures or thruster firing events. Level-0 products are available 24-hours after data reception.

+ Level 1:<br>
The level-1 data are the preprocessed, time-tagged and normal-pointed instrument data. These are the K-band ranging, accelerometer, star camera and GPS data of both satellites. Additionally the preliminary orbits of both GRACE satellites will be generated. Level-1 data processing software is developed by JPL with support from GFZ (e.g. accelerometer data preprocessing). Processing of level-1 products is done primarily at JPL. An identical processing system (hardware/software) is installed at GFZ to serve as a backup system in case of hardware or network problems. This double implementation is necessary to guarantee the envisaged level-1 product delay of 5 days. All level-1 products are archived at JPL’s Physical Oceanography Distributed Active Data Center(PODAAC) and at GFZ’s Integrated System Data Center (ISDC) . Both archives are harmonized on a sub-daily timeframe.
The level-1 data are the preprocessed, time-tagged and normal-pointed instrument data. These are the K-band ranging, accelerometer, star camera and GPS data of both satellites and laser ranging data for GRACE-FO. Additionally the preliminary orbits of both GRACE satellites will be generated. Level-1 data processing software is developed by JPL with support from GFZ (e.g. accelerometer data preprocessing). Processing of level-1 products is done primarily at JPL. An identical processing system (hardware/software) is installed at GFZ to serve as a backup system in case of hardware or network problems. This double implementation is necessary to guarantee the envisaged level-1 product delay of 5 days. All level-1 products are archived at JPL’s Physical Oceanography Distributed Active Data Center(PODAAC) and at GFZ’s Integrated System Data Center (ISDC) . Both archives are harmonized on a sub-daily timeframe.

+ Level 2: `Spherical Harmonic Coefficients` for the geopotential

Level-2 data include the short term (30 days) and mean gravity field derived from calibrated and validated GRACE level-1 data products. This level also includes ancillary data sets (temperature and pressure fields, ocean bottom pressure, and hydrological data) which are necessary to eliminate time variabilities in gravity field solutions. Additionally the precise orbits of both GRACE satellites are generated. All level-2 products are archived at JPL’s PODAAC and at GFZs ISDC and are available 60 days after data taking. The level-2 processing software were developed independently by all three processing centres using already existing but completely independent software packages which were upgraded for GRACE specific tasks. Common data file interfaces guarantees a strong product validation. Routine processing is done at UTCSR and GFZ, while JPL only generate level-2 products at times for verification purposes.
Level-2 data include the short term (30 days) and mean gravity field derived from calibrated and validated GRACE level-1 data products. This level also includes ancillary data sets (temperature and pressure fields, ocean bottom pressure, and hydrological data) which are necessary to eliminate time variabilities in gravity field solutions. Additionally the precise orbits of both GRACE satellites are generated. All level-2 products are archived at JPL’s PODAAC and at GFZs ISDC and are available 60 days after data taking. The level-2 processing software were developed independently by all three processing centres using already existing but completely independent software packages which were upgraded for GRACE specific tasks. Common data file interfaces guarantees a strong product validation. Routine processing is done at JPL, UTCSR and GFZ.

+ Level 3: `Mascons`<br>
+ Level 3: `Spatial fields`<br>
consists of mass anomalies or other standardized products such as Monthly Ocean/Land Water Equivalent Thickness, Surface-Mass Anomaly. Similarly mass concentration blocks or `mascons` are also available.

+ Level 4: `Time Series`<br>
Time-series of catchment level hydrological estimates of TWSA

`PySHBundle` provides the capability to obtain grided Total Water Storage Anomaly(TWSA) from Level 2 data.
`PySHBundle` provides the capability to obtain grided anoamlies (like TWSA, geopotential height etc) from Level 2 data.

## Spherical Harmonics Data Formats

Spherical harmonic functions or coefficients, Legendre functions and their derivatives can be arranged in different ways. There are multiple functions in SHBundle for reordering from one format to another. Some of them have been translated to Python in PySHBundle. Couple of new ones have also been added.

## Spherical Harmonics Data Formats

### clm-format

This is a standard format to store spherical harmonic coefficients in the indexed column-vector-format (abbreviatedL clm-format)
Expand Down
73 changes: 34 additions & 39 deletions examples/01_intro_spherical_harmonics_with_pyshbundle.ipynb
View file

Large diffs are not rendered by default.

82 changes: 24 additions & 58 deletions examples/02_load_data.ipynb
View file
Original file line number Diff line number Diff line change
Expand Up @@ -98,7 +98,7 @@
},
{
"cell_type": "code",
"execution_count": 14,
"execution_count": 3,
"id": "e2a0677d-3b5e-47b7-ae6b-58bda3731b4f",
"metadata": {},
"outputs": [],
Expand Down Expand Up @@ -142,7 +142,7 @@
},
{
"cell_type": "code",
"execution_count": 6,
"execution_count": 4,
"id": "a92d2ff6",
"metadata": {},
"outputs": [],
Expand All @@ -152,7 +152,7 @@
},
{
"cell_type": "code",
"execution_count": 7,
"execution_count": 5,
"id": "9e1ad2ad",
"metadata": {},
"outputs": [
Expand Down Expand Up @@ -193,7 +193,7 @@
},
{
"cell_type": "code",
"execution_count": 8,
"execution_count": 6,
"id": "245a71a8",
"metadata": {},
"outputs": [
Expand All @@ -206,7 +206,7 @@
" 'Slm_sdev': np.float64(0.0)}"
]
},
"execution_count": 8,
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
Expand Down Expand Up @@ -242,7 +242,7 @@
},
{
"cell_type": "code",
"execution_count": 9,
"execution_count": 7,
"id": "d373cd59",
"metadata": {},
"outputs": [
Expand All @@ -262,7 +262,7 @@
},
{
"cell_type": "code",
"execution_count": 10,
"execution_count": 8,
"id": "69b2ce55",
"metadata": {},
"outputs": [],
Expand Down Expand Up @@ -297,7 +297,7 @@
},
{
"cell_type": "code",
"execution_count": 11,
"execution_count": 9,
"id": "d117b268",
"metadata": {},
"outputs": [
Expand Down Expand Up @@ -357,7 +357,7 @@
},
{
"cell_type": "code",
"execution_count": 15,
"execution_count": 10,
"id": "137425a9",
"metadata": {},
"outputs": [],
Expand All @@ -371,7 +371,7 @@
},
{
"cell_type": "code",
"execution_count": 16,
"execution_count": 11,
"id": "9f83a5ec",
"metadata": {},
"outputs": [],
Expand Down Expand Up @@ -416,7 +416,7 @@
},
{
"cell_type": "code",
"execution_count": 19,
"execution_count": 12,
"id": "22288c18-1c81-4c26-b9c7-762126898822",
"metadata": {},
"outputs": [],
Expand All @@ -429,7 +429,7 @@
},
{
"cell_type": "code",
"execution_count": 20,
"execution_count": 13,
"id": "13f777d8",
"metadata": {},
"outputs": [],
Expand All @@ -439,7 +439,7 @@
},
{
"cell_type": "code",
"execution_count": 21,
"execution_count": 14,
"id": "ebe9f813",
"metadata": {},
"outputs": [
Expand Down Expand Up @@ -469,7 +469,7 @@
},
{
"cell_type": "code",
"execution_count": 22,
"execution_count": 15,
"id": "a4398d89",
"metadata": {},
"outputs": [],
Expand All @@ -479,7 +479,7 @@
},
{
"cell_type": "code",
"execution_count": 23,
"execution_count": 16,
"id": "cdbbc9b4",
"metadata": {},
"outputs": [
Expand Down Expand Up @@ -522,43 +522,14 @@
},
{
"cell_type": "code",
"execution_count": null,
"execution_count": 17,
"id": "ca759125",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "2a767edf",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "d571b9aa",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "62fdbe39",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "bf8d72ab",
"metadata": {},
"outputs": [],
"source": []
"source": [
"##### This section has been replaced with a simple for loop with explicit description of the replacement procedure.\n",
"##### Please see examples/04_TWSA_time_series.ipynb"
]
},
{
"cell_type": "markdown",
Expand All @@ -578,7 +549,7 @@
},
{
"cell_type": "code",
"execution_count": 24,
"execution_count": 18,
"id": "8faf394f",
"metadata": {},
"outputs": [],
Expand All @@ -603,7 +574,7 @@
},
{
"cell_type": "code",
"execution_count": 25,
"execution_count": 19,
"id": "23f99565",
"metadata": {},
"outputs": [],
Expand All @@ -619,7 +590,7 @@
},
{
"cell_type": "code",
"execution_count": 26,
"execution_count": 20,
"id": "f4028dcb",
"metadata": {},
"outputs": [
Expand Down Expand Up @@ -688,7 +659,7 @@
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"display_name": "pysh",
"language": "python",
"name": "python3"
},
Expand All @@ -703,11 +674,6 @@
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.9.19"
},
"vscode": {
"interpreter": {
"hash": "c942a33ffa182542ea7eb705673c4cde756b073012d32ad559a923b584d30dd7"
}
}
},
"nbformat": 4,
Expand Down
Loading

0 comments on commit a7cdfb0

Please sign in to comment.