Seawater Freezing Temperature Calculation

Usage

1. MCP Server Definition

import asyncio
import json
from contextlib import AsyncExitStack
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):
        try:
            self.transport = streamablehttp_client(
                url=self.server_url,
                headers={"SCP-HUB-API-KEY": self.api_key}
            )
            self._stack = AsyncExitStack()
            await self._stack.__aenter__()
            self.read, self.write, self.get_session_id = await self._stack.enter_async_context(self.transport)
            self.session_ctx = ClientSession(self.read, self.write)
            self.session = await self._stack.enter_async_context(self.session_ctx)
            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 hasattr(self, '_stack'):
                await self._stack.aclose()
            print("✓ already disconnect")
        except Exception as e:
            print(f"✗ disconnect error: {e}")
    def parse_result(self, 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)}

2. Freezing Temperature Workflow

Calculate the freezing point of seawater based on salinity and pressure.

Workflow Steps:

1. Calculate Absolute Salinity - Convert practical salinity
2. Calculate Freezing Temperature - Compute freezing point

Implementation:

## 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 parameters
input_params = {
    'SP': [35.0, 5.0],
    'p': [1000.0, 1000.0],
    'lon': [120.0, 165.0],
    'lat': [30.0, 45.0],
    'saturation_fraction': [0.0, 0.0]  # 0=air-free, 1=air-saturated
}

## Step 1: Calculate absolute salinity
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']
    }
)

SA_result = client.parse_result(result)["st"]

## Step 2: Calculate freezing temperature
result = await client.session.call_tool(
    "gsw_example_freezing_temp",
    arguments={
        "SA": SA_result,
        "p": input_params['p'],
        "saturation_fraction": input_params['saturation_fraction']
    }
)

t_freeze_result = client.parse_result(result)["st"]["t_freeze"]

print("Freezing Temperature Results:")
for i, t_freeze in enumerate(t_freeze_result):
    print(f"{i+1}. SA={SA_result[i]:.2f} g/kg, p={input_params['p'][i]} dbar")
    print(f"   Freezing temp: {t_freeze:.3f}°C\n")

await client.disconnect()

Tool Descriptions

OceanGSW-Tool Server:

• gsw_example_absolute_salinity: Convert practical to absolute salinity
• gsw_example_freezing_temp: Calculate freezing temperature
  • Args:
    • SA (list): Absolute salinity (g/kg)
    • p (list): Pressure (dbar)
    • saturation_fraction (list): Air saturation (0-1)
  • Returns: Freezing temperature (°C)

Input/Output

Input:

• SA: Absolute salinity (0-42 g/kg)
• p: Sea pressure (0-11000 dbar)
• saturation_fraction: 0 (air-free) to 1 (air-saturated)

Output:

• Freezing temperature in °C (typically -2 to 0°C)

Use Cases

• Sea ice formation prediction
• Polar oceanography
• Marine engineering in cold regions
• Climate modeling

Performance Notes

• Standards: TEOS-10
• Accuracy: ±0.001°C