Add automatic gestion of the case where you want to estimate spatial grid with a certain radius but coeffs radius do not correspond so internal change_reference() is applied

Author: hulecom

doc/tutorials/Advanced_lnharmo_validation.ipynb

@@ -22,7 +22,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": 2,
    "id": "08455b0e-32fb-49fd-b247-79a8067cb24c",
    "metadata": {},
    "outputs": [],
@@ -75,7 +75,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 3,
    "id": "07df9789-d606-489b-9a71-1a424f806c91",
    "metadata": {},
    "outputs": [],
@@ -84,18 +84,15 @@
     "\n",
     "ds_ref = xr.open_dataset(file_model, engine='lenapyGfc', no_date=True)\n",
     "\n",
-    "# Convert radius and GM of reference to be coherent with ICGEM computation\n",
-    "ds_ref = ds_ref.lnharmo.change_reference(new_radius=6378137.0, new_earth_gravity_constant=3.986005e14)\n",
-    "\n",
-    "# Default ds.attrs radius and earth_gravity_constant are used for the computation (set with change_reference())\n",
-    "grid_lenapy = ds_ref.lnharmo.to_grid(unit='potential', ellipsoidal_earth=False)\n",
+    "# Inside convertion of SH coefficients by giving a_earth and GM of the new reference to be coherent with ICGEM computation\n",
+    "grid_lenapy = ds_ref.lnharmo.to_grid(unit='potential', ellipsoidal_earth=False, a_earth=6378137.0, earth_gravity_constant=3.986005e14)\n",
     "\n",
-    "grid_lenapy_ellps = ds_ref.lnharmo.to_grid(unit='potential', ellipsoidal_earth=True, f_earth=1/298.257222101)"
+    "grid_lenapy_ellps = ds_ref.lnharmo.to_grid(unit='potential', ellipsoidal_earth=True, f_earth=1/298.257222101, a_earth=6378137.0, earth_gravity_constant=3.986005e14)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 4,
    "id": "28ee5716-c77e-4a0e-a872-b97184a7b0b5",
    "metadata": {},
    "outputs": [
@@ -150,7 +147,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 5,
    "id": "9f129e2a-c5a6-4869-abcd-b4e8628dae41",
    "metadata": {},
    "outputs": [
@@ -214,7 +211,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 6,
    "id": "8821900e-a794-4ba1-a1ce-edc9b5ff9ace",
    "metadata": {},
    "outputs": [],
@@ -229,7 +226,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 7,
    "id": "7dd1e41a-d705-44b8-a17b-ae21e3e9bd51",
    "metadata": {},
    "outputs": [
@@ -282,7 +279,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": 8,
    "id": "ef5824e1-54e2-4007-af08-065457822d3c",
    "metadata": {},
    "outputs": [