<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<?xml-model href="https://pds.nasa.gov/pds4/pds/v1/PDS4_PDS_1I00.sch" schematypens="http://purl.oclc.org/dsdl/schematron"?>

<Product_Observational xmlns="http://pds.nasa.gov/pds4/pds/v1" 
  xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" 
  xsi:schemaLocation="http://pds.nasa.gov/pds4/pds/v1 https://pds.nasa.gov/pds4/pds/v1/PDS4_PDS_1I00.xsd">
  <Identification_Area>
    <logical_identifier>urn:nasa:pds:grail_gravity_derived:data_shbdr:jggrx_1800f_shb_l100</logical_identifier>
    <version_id>1.0</version_id>
    <title>GRAIL Gravity Binary Spherical Harmonic Model: jggrx_1800f_shb_l100</title>
    <information_model_version>1.18.0.0</information_model_version>
    <product_class>Product_Observational</product_class>
    <Modification_History>
      <Modification_Detail>
        <modification_date>2026-05-16</modification_date>
        <version_id>1.0</version_id>
        <description>Initial PDS4 version</description>
      </Modification_Detail>
    </Modification_History>
  </Identification_Area>
  <Observation_Area>
    <Time_Coordinates>
      <start_date_time>2012-03-01T16:28:00.000Z</start_date_time>
      <stop_date_time>2012-12-14T20:56:00.000Z</stop_date_time>
    </Time_Coordinates>
    <Investigation_Area>
      <name>Gravity Recovery and Interior Laboratory</name>
      <type>Mission</type>
      <Internal_Reference>
        <lid_reference>urn:nasa:pds:context:investigation:mission.gravity_recovery_and_interior_laboratory</lid_reference>
        <reference_type>data_to_investigation</reference_type>
      </Internal_Reference>
    </Investigation_Area>
    <Observing_System>
      <Observing_System_Component>
        <name>Gravity Recovery and Interior Laboratory A</name>
        <type>Host</type>
        <Internal_Reference>
          <lid_reference>urn:nasa:pds:context:instrument_host:spacecraft.grail-a</lid_reference>
          <reference_type>is_instrument_host</reference_type>
        </Internal_Reference>
      </Observing_System_Component>
      <Observing_System_Component>
        <name>Gravity Recovery and Interior Laboratory B</name>
        <type>Host</type>
        <Internal_Reference>
          <lid_reference>urn:nasa:pds:context:instrument_host:spacecraft.grail-b</lid_reference>
          <reference_type>is_instrument_host</reference_type>
        </Internal_Reference>
      </Observing_System_Component>
      <Observing_System_Component>
        <name>Lunar Gravity Ranging System A for GRAIL-A</name>
        <type>Instrument</type>
        <Internal_Reference>
          <lid_reference>urn:nasa:pds:context:instrument:grail-a.lgrs-a</lid_reference>
          <reference_type>is_instrument</reference_type>
        </Internal_Reference>
      </Observing_System_Component>
      <Observing_System_Component>
        <name>Lunar Gravity Ranging System B for GRAIL-B</name>
        <type>Instrument</type>
        <Internal_Reference>
          <lid_reference>urn:nasa:pds:context:instrument:grail-b.lgrs-b</lid_reference>
          <reference_type>is_instrument</reference_type>
        </Internal_Reference>
      </Observing_System_Component>
    </Observing_System>
    <Target_Identification>
      <name>Moon</name>
      <type>Satellite</type>
      <Internal_Reference>
        <lid_reference>urn:nasa:pds:context:target:satellite.earth.moon</lid_reference>
        <reference_type>data_to_target</reference_type>
      </Internal_Reference>
    </Target_Identification>
  </Observation_Area>
  <File_Area_Observational>
    <File>
      <file_name>jggrx_1800f_shb_l100.dat</file_name>
      <creation_date_time>2024-07-18T00:00:00.000</creation_date_time>
      <comment>
        This file contains the covariance and related data for the JPL Lunar 
        gravity field GL1800F, a 1800th degree and order spherical harmonic model.
        It is a JPL gravity field that includes the entire primary
        mission (March 1, 16:30 to May 29, 16:36, 2012) and the entire extended
        mission GRAIL tracking data (August 30, 16:20 to Dec. 14, 20:56, 2012).
        This solution uses the version 4 level-1 data for the primary and
        extended missions.
        
        Some details describing this model are:
        The spherical harmonic coefficients are fully normalized.
        The reference radius = 1738.0 km
        The planetary ephemeris is de440 and defines the lunar body-fixed
        principal axes (PA) coordinate system.
        The second degree Love number solution is estimated to be k2=0.024223.
        The third degree Love number solution is estimated to be k3=0.0163.
        The second and third degree gravity coefficients of this model do not
        include the permanent tide.
        A constraint is applied starting at n=601. The a priori value is the
        topography from gravity coefficient value (i.e., not zero)
        and the a priori uncertainty equals 0.3 x the gravity from topography
        coefficient value using topography to power 23 and a variable density.
        
        The weighting of the KBRR data is 0.03 to 0.06 microns/sec for the
        primary mission and 0.05 microns/sec for the entire extended mission.
        
        This product contains the truncated n=100 covariance of the GL1800F
        gravity model or JGGRX_1800F_SHB_L100. The covariance is scaled to realistic
        uncertainties based upon differences with other solutions. This increases
        the uncertainties for degree 2 and 3 by a factor of 15 and decreases to 1 
        at degree 120. The covariance is generated from 2000 clones of GL1800F.
        
        Reference: R. S. Park, A. Berne, A. S. Konopliv, J. T. Keane,
        I. Matsuyama, F. Nimmo, M. Rovira-Navarro, M. P. Panning,
        M. Simons, D. J. Stevenson, R. C. Weber, 'Thermal asymmetry in the
        Moon's mantle inferred from monthly tidal response', Nature, 2025.
        
        This product is a set of binary tables:
        a header table, a names table, a coefficients table, and a covariance
        table. Definitions of the tables follow. This GRAIL moon gravity model
        is in the form of a Spherical Harmonics Binary Data Record (SHBDR).
      </comment>
    </File>
    <Table_Binary>
      <name>SHBDR_Header_Table</name>
      <offset unit="byte">0</offset>
      <records>1</records>
      <description>The SHBDR Header includes
        descriptive information about the spherical harmonic
        coefficients which follow in SHBDR_COEFFICIENTS_TABLE.
        The header consists of a single record of nine data
        columns requiring 56 bytes.  The Header is followed by
        a pad of binary integer zeroes.</description>
      <Record_Binary>
        <fields>9</fields>
        <groups>0</groups>
        <record_length unit="byte">56</record_length>
        <Field_Binary>
          <name>Reference_Radius</name>
          <field_number>1</field_number>
          <field_location unit="byte">1</field_location>
          <data_type>IEEE754LSBDouble</data_type>
          <field_length unit="byte">8</field_length>
          <unit>km</unit>
          <description>The assumed reference
            radius of the spherical planet.</description>
        </Field_Binary>
        <Field_Binary>
          <name>Constant</name>
          <field_number>2</field_number>
          <field_location unit="byte">9</field_location>
          <data_type>IEEE754LSBDouble</data_type>
          <field_length unit="byte">8</field_length>
          <description>For a gravity field model
            the gravitational constant GM in km cubed per second
            squared for the planet.  For a topography model, set to 1.
          </description>
        </Field_Binary>
        <Field_Binary>
          <name>Uncertainty_in_Constant</name>
          <field_number>3</field_number>
          <field_location unit="byte">17</field_location>
          <data_type>IEEE754LSBDouble</data_type>
          <field_length unit="byte">8</field_length>
          <description>For a gravity field model
            the uncertainty in the gravitational constant GM in km
            cubed per second squared for the planet.  For a topography
            model, set to 0.</description>
        </Field_Binary>
        <Field_Binary>
          <name>Degree_of_Field</name>
          <field_number>4</field_number>
          <field_location unit="byte">25</field_location>
          <data_type>SignedLSB4</data_type>
          <field_length unit="byte">4</field_length>
          <description>Degree of the model field.</description>
        </Field_Binary>
        <Field_Binary>
          <name>Order_of_Field</name>
          <field_number>5</field_number>
          <field_location unit="byte">29</field_location>
          <data_type>SignedLSB4</data_type>
          <field_length unit="byte">4</field_length>
          <description>Order of the model field.</description>
        </Field_Binary>
        <Field_Binary>
          <name>Normalization_State</name>
          <field_number>6</field_number>
          <field_location unit="byte">33</field_location>
          <data_type>SignedLSB4</data_type>
          <field_length unit="byte">4</field_length>
          <description>The normalization indicator.
            For gravity field:
            0   coefficients are unnormalized
            1   coefficients are normalized
            2   other.</description>
        </Field_Binary>
        <Field_Binary>
          <name>Number_of_Names</name>
          <field_number>7</field_number>
          <field_location unit="byte">37</field_location>
          <data_type>SignedLSB4</data_type>
          <field_length unit="byte">4</field_length>
          <description>Number of valid names in
            the SHBDR Names Table.  Also, the number of valid
            coefficients in the SHBDR Coefficients Table.</description>
        </Field_Binary>
        <Field_Binary>
          <name>Reference_Longitude</name>
          <field_number>8</field_number>
          <field_location unit="byte">41</field_location>
          <data_type>IEEE754LSBDouble</data_type>
          <field_length unit="byte">8</field_length>
          <unit>degree</unit>
          <description>The reference longitude for
            the spherical harmonic expansion; normally 0.</description>
        </Field_Binary>
        <Field_Binary>
          <name>Reference_Latitude</name>
          <field_number>9</field_number>
          <field_location unit="byte">49</field_location>
          <data_type>IEEE754LSBDouble</data_type>
          <field_length unit="byte">8</field_length>
          <unit>degree</unit>
          <description>The reference latitude for
            the spherical harmonic expansion; normally 0.</description>
        </Field_Binary>
      </Record_Binary>
    </Table_Binary>
    
    <Table_Binary>
      <name>SHBDR_Names_Table</name>
      <offset unit="byte">512</offset>
      <records>10198</records>
      <description>
        The SHBDR Names Table
        contains names for the solution parameters (including
        gravity field coefficients) which will follow in
        SHBDR_COEFFICIENTS_TABLE.  The order of the names
        in SHBDR_NAMES_TABLE corresponds identically to the
        order of the parameters in SHBDR_COEFFICIENTS_TABLE.
        Each coefficient  name is of the form  Cij  or  Sij
        where  i  is the degree of the coefficient and  j  is
        the order of the coefficient.  Both indices are three-
        digit zero-filled right-justified ASCII character strings
        (for example,  C010005  for the 10th degree 5th order C
        coefficient, or  S002001  for the 2nd degree 1st order
        S  coefficient).  The eighth byte in the table is an
        ASCII blank used to ensure that the row length
        is equal to 8 bytes.  Names of other solution
        parameters are limited to 8 ASCII characters; if less
        than 8, they will be left-justified and padded with
        ASCII blanks.  The Names Table itself will be padded
        with ASCII blanks, if necessary.</description>
      <Record_Binary>
        <fields>1</fields>
        <groups>0</groups>
        <record_length unit="byte">8</record_length>
        <Field_Binary>
          <name>Parameter_Name</name>
          <field_number>1</field_number>
          <field_location unit="byte">1</field_location>
          <data_type>ASCII_String</data_type>
          <field_length unit="byte">8</field_length>
          <description>The name of the
            coefficient or other solution parameter, left-justified and padded with ASCII blanks (if needed)
            to 8 characters.</description>
        </Field_Binary>
      </Record_Binary>
    </Table_Binary>
    
    <Table_Binary>
      <name>SHBDR_Coefficients_Table</name>
      <offset unit="byte">82432</offset>
      <records>10198</records>
      <description>The SHBDR Coefficients Table
        contains the coefficients and other solution parameters
        for the spherical harmonic model.  The order of the
        coefficients in this table corresponds exactly to the
        order of the coefficient and parameter names in
        SHBDR_NAMES_TABLE.  The SHBDR Coefficients Table will be
        padded with double precision zeroes.</description>
      <Record_Binary>
        <fields>1</fields>
        <groups>0</groups>
        <record_length unit="byte">8</record_length>
        <Field_Binary>
          <name>Coefficient_Value</name>
          <field_number>1</field_number>
          <field_location unit="byte">1</field_location>
          <data_type>IEEE754LSBDouble</data_type>
          <field_length unit="byte">8</field_length>
          <description>A coefficient Cij or
            Sij or other solution parameter as specified in the
            SHBDR Names Table.</description>
        </Field_Binary>
      </Record_Binary>
    </Table_Binary>
    
    <Table_Binary>
      <name>SHBDR_Covariance_Table</name>
      <offset unit="byte">164352</offset>
      <records>52004701</records>
      <description>The SHBDR_Covariance_Table
        contains the covariances for the spherical harmonic model
        coefficients and other solution parameters. The order of
        the covariances in this table is defined as columnwise
        vector storage of the upper triangular matrix formed by
        the product of the SHBDR Names Table with its transpose.
        For example, if the Names Table has four entries  A, B, 
        C, and D,  then the covariances are given by the column 
        vectors in the upper triangular matrix of:
        | A | [ A B C D ] = | AA AB AC AD |
        | B |               | BA BB BC BD |
        | C |               | CA CB CC CD |
        | D |               | DA DB DC DD |
        
        That is, the covariance table will list (in this order)
        AA, AB, BB, AC, BC, CC, AD, BD, CD, and DD.
        The SHBDR_Covariance_Table will be padded with double
        precision zeroes</description>
      <Record_Binary>
        <fields>1</fields>
        <groups>0</groups>
        <record_length unit="byte">8</record_length>
        <Field_Binary>
          <name>Covariance_Value</name>
          <field_number>1</field_number>
          <field_location unit="byte">1</field_location>
          <data_type>IEEE754LSBDouble</data_type>
          <field_length unit="byte">8</field_length>
          <description>The covariance value
            for the coefficients and other solution parameters
            specified by the product of SHBDR_NAMES_TABLE with
            its transpose, after omitting redundant terms.</description>
        </Field_Binary>
      </Record_Binary>
    </Table_Binary>
  </File_Area_Observational>
  
  <File_Area_Observational_Supplemental>
    <File>
      <file_name>jggrx_1800f_shb_l100.lbl</file_name>
    </File>
    <Stream_Text>
      <offset unit="byte">0</offset>
      <parsing_standard_id>PDS3</parsing_standard_id>
      <description>Original PDS3 label</description>
      <record_delimiter>Carriage-Return Line-Feed</record_delimiter>
    </Stream_Text>
  </File_Area_Observational_Supplemental>
</Product_Observational>
