最近工作需要用到最远点采样fps算法, 找了一些blog学习学习, 写了个原型: 测试一下, 跑的没问题, 但是速度感人 所有点的欧氏距离计算做个向量化优化: 这一下子快了1...
发简信
IP属地：浙江
-
FarthestPointSampling(FPS)算法学习&实践
-
记一个gcc7代编译器下allocate-source中有关array bounds的离谱bug
Bug情况 最近写Fortran代码的时候遇到了个怪bug, 代码如下: 然后离谱的来了, 我发现在gcc750和gcc710里面跑出来的shape_data的lbound都...
-
cmake - 编写单元测试在cmake中构建测试例就是enable_testing()+add_test()这一套。本篇文章是我在学习cmake构建单元测试的一个总结笔记。 工程概览 如下是我所定义的...
-
WPS ungrib编译报错 undefined reference to 'jpc_decode'
最近在一个新机子上用Spack配置好编译器, hdf5, netcdf-c, netcdf-fortran, jasper等一系列依赖库后编译了最新版本的WRF-4.4.2和...
-
地理坐标转换-proj4
目标：WGS-84 转为 CGCS2000坐标 使用proj4 测试数据 举例：天安门在不同类型坐标系下的经纬度值 WGS-84经纬度：116.391349,39.90737...
-
但是这不是netCDF4吗😁
利用xarray写netcdf文件
三维模式输出的数据很大，传输起来很麻烦。有时候我们只需要几个关键变量，因此把他们抽取出来写成nc文件比较方便。 例如，要单独保存megan输出的isoprene变量isop为...
-
EBAS观测数据处理
记录一下 EBAS-Thredds[https://ebas.nilu.no/thredds/]数据的处理流程, 以全站点的ozone为例 基于前一个EBAS观测数据下载[h...
-
m3diff对比分析ioapi-netcdf文件
最近要在一个新机器上跑CMAQ, 先对机器和模式做个校验, 以前都是跑出来自己写个ncl脚本比对的, 现在觉得用ioapi的m3diff也挺好, 不过还需要基于那个结果进行后...
-
EBAS观测数据下载
记录一下 下载EBAS-Thredds[https://ebas.nilu.no/thredds/]数据的流程, 以全站点的ozone为例 1. 下载ebas-thredds...
-
CFSR 6小时数据批量下载脚本-WRF模拟用
WRF官方[https://www2.mmm.ucar.edu/wrf/users/download/free_data.html]提供的CFSR数据下载地址是在UCAR[h...
-
FORTRAN调用C函数的参数不一致问题-自动额外传参特性
最近在debug的时候看到CMAQ ioapi的code, envgets.c中和ENVSTR函数和m3utilio.F中的接口参数不一致, 但是可以正常跑, 之前没注意过,...
-
CS61C-2022-su20作业记录
记录一下CS61C 20su的各个lab和proj的心得, 代码放在了这里[https://github.com/Roadelse/CS61C-2022-su20] proj...
-
ImageMagick日常使用格式记录
记录一些日常用ImageMagick处理图像的格式 1. 剪切图的一部分 magick convert src.png -gravity East -crop 200x300...
-
观测数据处理-国控环境监测站点
工作需要用到我国环境监测站点的污染物浓度数据(感谢大佬的分享), 数据很全, 不过csv格式在分析的时候尤其大量数据分析的时候并不友好, 所以一般要二次处理一下 上学的时候搞...
-
powershell数组赋值遇到的bug, 开头的逗号不是一直都能用的
今天写Powershell的时候遇到个bug, 搜了半天没搜到(不过最后解决后感觉以前似乎看到过, 可能我搜的内容不对) 简单地说, powershell弄个数组一般就直接逗...
-
numpy slice broadcast的坑, 和ncl或者xarray不同
过去用ncl作array slice习惯了, 不同的dimension之间就是直接乘起来, 例如 即可得到一个2x2的数组 结果今天在用numpy处理数据的时候发现 就只能得...
-
Sublime修改Build Python的默认指令
网上查到的配置sublime直接运行python脚本的方法大多是新建一个Build python的设置, 而直接以默认的Python build运行的话可能会报错提示找不到指...
-
好像目前是隐藏了中国站点..但是手动输入站点编号比如北京是54511, 历史数据还是可以看的
不顾20年的似乎不行
另外历史数据我之前都down了,有需要的话可以发给你下载wyoming大学的探空数据
水平有限,欢迎指正交流,共同进步! 俄怀明大学将全球探空站的数据共享使用,并通过计算了很多的变量指数.非常有用,近期有这方面的需求,暂时写了一段未经优化的代码下载.如果要大量...
-
mpi4py 快速上手
在上一篇中我们介绍了如何安装和使用 mpi4py，下面我们以几个简单的例子来展示怎么使用 mpi4py 来进行并行编程，以使读者能够快速地上手使用 mpi4py。这些例子来自...