本项目提供了一套基于ISQG模型的MATLAB代码，旨在通过二维海表面的数据（如高度和温度），推算出三维海洋的状态。为研究海洋动力学提供了有力工具。
在MATLAB环境中使用QG（拟地转）方法的用户可以利用invert.py脚本作为黑盒工具来计算从一组数据得出的isQG解。以下是所需输入参数：
- ssd(Ny,Nx)：海面密度，单位为kg/m^3。
- ssh(Ny,Nx)：海表面高度，单位为米。
- z(Nz)：垂直坐标数组，表示从海底到表层的位置（例如z=[-3000...-10] 米）。请注意不要将表层设在z=0处。因为无论怎样，在z=0位置的N^2都没有物理意义。
- n2(Nz)：分层点上的n2值，即重力加速度g除以海水密度ρ_0乘以dρ/dz（单位为/s）。
- lat(Ny): 纬度信息
- lon(Nx): 经度信息，使用度数表示
此外还有一个可选参数：
- useanomaly：布尔值，默认设置为True。当设为True时，在反演过程中会应用ssd和ssh的异常场（即原始2D场与最佳拟合二次曲面之间的差异）。
确保上述变量已准备好并保存在名为datain.mat的MATLAB文件中，且以小写形式命名相同名称的变量。假设数据已经存储在名为datain.mat的文件里，则可以通过以下命令行运行invert.py脚本：
```
python invert.py datain.mat dataout.mat
```
或者直接调用 `invert.py` 并提供相应的输入和输出文件名：
```
invert.py datain.mat datao
```
5星
