- 博客(31)
- 收藏
- 关注
原创 grib文件转nc
使用教程参考https://www.ncl.ucar.edu/Document/Tools/ncl_convert2nc.shtml。使用ncl中的ncl_convert2nc,可以直接在命令中使用,而不用先打开ncl。
2023-12-10 21:17:26
801
原创 使用python在ECMWF的网站上批量下载数据
注意:之前装过另一个包ecmwfapi,在下载数据是报错了,可能是旧的包有问题,这个是新的包。密钥的查看地址:https://api.ecmwf.int/v1/key/1、安装包ecmwf-api-client。写入文件$HOME/.ecmwfapirc。3、下载数据,以ERA-20C为例。2、在家目录写入密钥文件。
2023-05-29 10:42:08
654
原创 CESM运行错误:Can‘t redeclare “my“ in “my“ at /~/models/drv/bld/build-namelist line 787, near “my“
【代码】CESM运行错误:Can‘t redeclare “my“ in “my“ at /~/models/drv/bld/build-namelist line 787, near “my“
2023-05-24 15:37:32
140
原创 Can‘t locate XML/LibXML.pm
一次并不能安装完perl中所有内容(pm),需要根据需要(安装的时候会报错,就知道需要什么依赖),一个一个手动下载并安装。如果机器不能联网,需要上传安装包,下载地址和安装方法,参考。该错误表示缺少perl包,或者perl里缺少部分东西。如果是管理员且机器可以联网,可以使用命令直接安装。如果不是管理员,需要安装新的perl5到家目录。这一步十分繁琐,我下载和安装内容如下。
2023-05-24 15:30:50
218
原创 CESM2中修改CAM6的初始条件
1、下载ERA5中的uwind, vwind, temperature, surface pressure和specific humidity。注意:如果直接使用线性插值可能会出现如下报错。应该重新考虑垂直方向上的插值方法,笔者采用三次样条插值后,可以正常运行。2、插值到atm_in中的ncdata文件。
2023-05-15 23:37:57
225
2
原创 CESM中PHIS与给定的边界条件中topography有不同
CESM在得到边界条件中的地形数据后,会进行“谱截断(spectral truncation)”的操作,使得大气模块在运行时有比较平滑的边界条件。具体代码位置:dynamics/eul/spetru.F90。
2023-02-06 14:11:21
108
原创 经验正交分解(EOF)的推导
Am×n表示待分析变量(anomaly),有m个格点,n个时间片,Aa1a2...anXAATXm×m表示方差-协方差方阵。Xm×mUm×mTSm×mUm×m其中S是对角方阵,Sλ是对角线上的特征值组成的向量,其中Uu1u2...um则有UAUATSUA可以写作UAua1ua2...uam。
2023-01-16 20:06:27
844
原创 ncl读取xlsx的数据
2、使用ncl函数asciiread读取:data = asciiread( file_in, (/nlat, nlon/), “float” )注意:必须先将数据格式修改为“数值”,而不是“文本”,否则在读取时会出现如下报错。1、先将xlsx的数据导出为txt格式文件。
2022-11-27 11:49:58
1378
1
原创 报错:Can‘t locate Switch.pm in @INC
注意:https://www.cpan.org/src获取最新偶数版本下载链接并替换(偶数版本为稳定版)报错内容:Can’t locate Switch.pm in @INC。2、configure,指定安装目录和使用的编译器(icc)备注:使用默认的cc编译器时可能会报错。参考和摘抄了@童心同萌的。报错原因:缺少perl。
2022-11-24 11:39:45
1666
原创 表面重力波推导
当流体作无旋运动时, 运动方程和连续方程将具有非常简单的形式. 首先, 速度旋度 为零。的条件保证在任何单连通域内存在单值速度势。
2022-10-04 14:56:09
416
原创 WaveWatch3学习笔记——安装与测试
5、设置关于netcdf的环境变量,先查看netcdf的版本(nc-config --version)备注:若出现编译错误,则运行“bin/w3_clean all”,而后重新编译。4、编译前准备,找到待使用编译器对应的comp和link文件,以Gnu为例。3、生成./model/bin/wwatch3_env,编译aux文件。
2022-10-04 14:55:52
1627
1
原创 WaveWatch3学习笔记——模拟实例(1)
wind.raw是描述风速的文本文件,首行为时间,如“20150101 000000”,然后是每一纬度的纬向风速(向东为正),从西向东写入,按照从北到南的顺序写入各个纬度上的数值,每个数值单独占一行。Pacific.bot是描述地形高度的文本文件,每一行为某一纬度的高度值,从左到右对应经度增加,从上到下各行的对应纬度值增加(默认北半球纬度为正,南半球纬度为负,东半球经度为0~180,西半球精度为180~360);Pacific.mask是描述格点类型的文本文件,文件的书写方式和Pacific.bot相同。
2022-10-04 14:55:32
2500
1
原创 SWAN学习笔记——安装与模拟实例
global_wind.dat是描述风速的文本文件,将每一纬度的纬向风速(向东为正)从西向东写入,按照从南到北的顺序写入各个纬度上的数值,每个数值单独占一行。写完纬向风速后,再写经向风速(向北为正),全部写完后再写下一个时刻的风速值。Atlantic_EW.grd是描述网格坐标的文本文件,将每一纬度的经向坐标从西向东写入,按照从南到北的顺序写入各个纬度上的数值,每个数值单独占一行。此处使用空间分辨率为0.5°×0.5°的网格,时间分辨率是1小时,空间范围是85W~40W,27N~50N。
2022-10-03 18:26:29
4203
5
原创 Linux 安装netcdf
1、安装netcdf-fortran之前需要安装netcdf-c,在安装netcdf是需要告知netcdf-c的路径。2、netcdf.inc仅在netcdf3中,netcdf4中没有。
2022-10-03 18:25:45
1649
原创 Error: ‘b‘ argument of ‘sign‘ intrinsic at (1) must be the same type and kind as ‘a‘
解决方案:在编译时将单精度数全部转为双精度数,-freal-4-real-8。原因:对于sign这个库函数,它要求两个参数必须为同样的数据类型。
2022-10-03 18:24:46
380
原创 gfortran: warning: ‘11.0‘ is not valid for ‘mmacosx-version-min
在环境变量中加入export MACOSX_DEPLOYMENT_TARGET=10.6。
2022-10-03 18:19:00
259
原创 INITIAL, FILE NOT FOUND ERROR
在读入多个hotfile时,若hotfile的命名中出现点,读取文件会失败,需要在输出文件时简单命名,如“hotfs”
2022-10-03 18:16:32
249
原创 Linux环境中非交互地安装Matlab
可以根据个人计算机下载相应系统的安装程序,打开后,点击高级选项,下载目标计算机相应系统的matlab安装包。复制license.lic到文件夹licenses,编写激活文件acitvate.ini。3、修改installer_input.txt。1、下载matlab安装程序和安装包。2、上传到目标计算机,并解压安装包。
2022-10-03 18:07:31
484
转载 MPI多线程传输大量数据
MPI多线程传输大量数据报错如下Fatal error in MPI_Waitall: Other MPI error, error stack:MPI_Waitall(392)…: MPI_Waitall(count=200, req_array=0x7ffee33bc3e0, status_array=0x7ffee33bb440) failedMPIR_Waitall_impl(221)…:MPIDI_CH3I_Progress(547)…:pkt_CTS_handler(350)…:M
2021-12-21 13:51:03
661
1
转载 Segmentation fault signal 11
Segmentation fault signal 11解决方法解决方法系统设置的堆栈大小不够,ulimit -s unlimited
2021-12-21 11:24:36
3156
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人