Here is a quora answer about the slab ocean
Slab ocean models consider this whole upper layer as one single body of water, or a 'slab' with one single velocity vector (u,v), and density. The Navier-Stokes equation is modified accordingly to compute the resultant velocity of that slab based on wind forcing and coriolis force. The mixed-layer depth in such models is usually kept fixed.
Here is a FAQ about the slab ocean model. Two things to note:
SOM is a kind of data model, can be switched in a namelist variable F compset, env_run.xml.
Some standard SOM forcing is available in the inputdata repository.
Now we may start to test the slab ocean. compset choose: E_1850_CAM5
ERROR comes:
svn: URL 'https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/ocn/docn7/SOM/UNSET' non-existent in that revision
There should be something wrong with the set of the ocean forcing files. So come to this dir, no problem.
We need to set the specific foring files in the env_run.xml file. (check FAQ)
We choose pop_frc.b.c40.B1850CN.f19_g16.100105.nc for test.
In env_run.xml, only give the name of the file is okay to continue.
After downloading the file, you can build and submit the case. Yet there is no problem even after altering the physics into cam5 -chem none.
Perfect!
突然想到能否有什么方法知道发出去的某封邮件是否被别人查阅，查了下，有人提到图片挂参的方法可以办到，回想之前IP验证和php自动发邮件，突然想到可以做这么一套东西，每天固定时间统计一下今天发出去的邮件多少被查阅了，然后做个报告发到邮箱。将来撒海网套磁的时候，绝对是利器啊。LOL
说搞就搞：
首先，测试一下php返回图片
参照这个帖子，bingo，原来是改一下http的header就可以。
然后引入访问IP的检测，直接抄13年检测服务器访问沿海自动站IP的程序，调用IP138结果记录。
居然就可以了～给php文件挂个参数标志是哪一个邮件，嗯嗯
最后设置下crontab，每天固定时间把log发到邮箱，bingo～再写一个可以随时访问的查看页面，done～
Some points to be aware:
2D fields, not same by binary output and netcdf (high efficiency) output, one point value missing, might be the polar point. Therefore, *.bin control and *.nc sensitive is not under the same structure.
2D fields work fine under v1.9 of grads. However, 3D fields might meet some trouble. Try to use the “standard” attribution, dimention to tackle this problem. An example from NCEP-DOE dataset:
float level(level) ; level:units = “millibar” ; level:long_name = “Level” ; level:positive = “down” ; level:GRIB_id = 100s ; level:GRIB_name = “hPa” ; level:actual_range = 1000.f, 10.f ; level:axis = “Z” ;