Seawater Sound Speed Calculation

import asyncio
import json
from mcp.client.streamable_http import streamablehttp_client
from mcp import ClientSession

class OceanClient:
    """OceanGSW-Tool MCP Client"""

    def __init__(self, server_url: str, api_key: str):
        self.server_url = server_url
        self.api_key = api_key
        self.session = None

    async def connect(self):
        """Establish connection and initialize session"""
        try:
            self.transport = streamablehttp_client(
                url=self.server_url,
                headers={"SCP-HUB-API-KEY": self.api_key}
            )
            self.read, self.write, self.get_session_id = await self.transport.__aenter__()
            self.session_ctx = ClientSession(self.read, self.write)
            self.session = await self.session_ctx.__aenter__()
            await self.session.initialize()
            return True
        except Exception as e:
            print(f"✗ connect failure: {e}")
            return False

    async def disconnect(self):
        """Disconnect from server"""
        try:
            if self.session:
                await self.session_ctx.__aexit__(None, None, None)
            if hasattr(self, 'transport'):
                await self.transport.__aexit__(None, None, None)
        except Exception as e:
            print(f"✗ disconnect error: {e}")

    def parse_result(self, result):
        """Parse MCP tool call result"""
        try:
            if hasattr(result, 'content') and result.content:
                content = result.content[0]
                if hasattr(content, 'text'):
                    return json.loads(content.text)
            return str(result)
        except Exception as e:
            return {"error": f"parse error: {e}", "raw": str(result)}

## Initialize client
client = OceanClient(
    "https://scp.intern-ai.org.cn/api/v1/mcp/34/OceanGSW-Tool",
    "<your-api-key>"
)

if not await client.connect():
    print("connection failed")
    exit()

## Input: Seawater properties
input_params = {
    'SP': [35.0, 5.0],           # Practical salinity (PSU)
    't': [15.0, 10.0],           # In-situ temperature (°C)
    'p': [1000.0, 1000.0],       # Pressure (dbar)
    'lon': [120.0, 165.0],       # Longitude (degrees East)
    'lat': [30.0, 45.0]          # Latitude (degrees North)
}

## Step 1: Calculate absolute salinity (SA)
result = await client.session.call_tool(
    "gsw_example_absolute_salinity",
    arguments={
        "SP": input_params['SP'],
        'p': input_params['p'],
        'lon': input_params['lon'],
        'lat': input_params['lat']
    }
)

result_data = client.parse_result(result)
SA_result = result_data["st"]
print("Absolute Salinity:")
for i, sa in enumerate(SA_result):
    print(f"  SP={input_params['SP'][i]} → SA={sa:.4f} g/kg")

## Step 2: Calculate conservative temperature (CT)
result = await client.session.call_tool(
    "gsw_example_conservative_temperature",
    arguments={
        "SA": SA_result,
        't': input_params['t'],
        'p': input_params['p']
    }
)

result_data = client.parse_result(result)
CT_result = result_data["st"]
print("\nConservative Temperature:")
for i, ct in enumerate(CT_result):
    print(f"  t={input_params['t'][i]}°C → CT={ct:.4f}°C")

## Step 3: Calculate sound speed
result = await client.session.call_tool(
    "gsw_example_sound_speed",
    arguments={
        "SA": SA_result,
        'CT': CT_result,
        'p': input_params['p']
    }
)

result_data = client.parse_result(result)
sound_speed_result = result_data["st"]["sound_speed"]

print("\nSound Speed Results:")
for i, speed in enumerate(sound_speed_result):
    print(f"{i+1}. SA={SA_result[i]:.2f} g/kg, CT={CT_result[i]:.2f}°C, p={input_params['p'][i]} dbar")
    print(f"   Sound speed: {speed:.2f} m/s\n")

await client.disconnect()