- 博客(19)
- 收藏
- 关注
RSS订阅原创 调式flexinvert
1、vim SETTINGS_ghg2、vim p3、4、vim prep_pathnames.f90设置GFS, fmt=’(A)’ ‘/usr/local/GFS/2014/’5、vim6、vim prep_releases.f90 mass=1. （144行）7、vim run_flexpart修改 SLURM=0 (这个=1）是用来提交大型机跑的把后面一段关于提交到大型机的段注释掉 <<!
2022-01-06 10:55:35
450
1
原创 温室气体内通量分析及排放量反演-FLEX_INVERT
tar -xvf lapack-3.1.1.tar.gzcd lapack-3.1.1cp make.inc.ifort make.incvim Makefilelib: lapacklib tmglib#lib: blaslib lapacklib tmglib改为#lib: lapacklib tmgliblib: blaslib lapacklib tmglibmake
2022-01-06 10:34:53
1009
原创 解决FLEXPART使用EC数据全球气象场的编译器问题，修改par_mod.f90，安装gfortran-5
一.先复制一个包sudo cp -r src2 src3##二.修改par_mod.f90cd src3vim par_mod.f90重新make ncf=yes但是编译通不过。这真是一个坑。后来发现发装gfortran-5sudo apt-get install gfrotran-5修改后，还要修改makefile里的参数vim makefile‘’’elseROOT_DIR = /usr/localF90 = /usr/bin/gfortran-5‘’’然后重新
2022-01-06 10:29:03
936
原创 FLEXPART安装笔记，ubantu，grib_api，eccodes，NetCDF
FLEXPART安装配置笔记 2021-1-5 操作系统版本：ubuntu 18.04 1.安装前准备说明：编译安装都安装到/usr/localpython3pip3 系统配置sudo apt-get updatesudo apt-get install g++ gfortran autoconf libtool automake flex bison cmake python3-dev python3-pip git vim curl build-essential li
2021-09-06 19:02:11
3003
8
原创 bz2,zip,txt数据解析
温室气体GC的数据经二次压缩，zip和bz2。最后为bz2文件，数据在最后解压后的txt里。一、bz2解压为zip文件#Wu dongqiao 2021.08.02import bz2import zipfileimport osimport pandas as pdimport numpy as np'''=======================================bz2解压为zip文件===========================================
2021-08-12 21:38:31
718
原创 温室气体数据订证之PICARRO（python、bz2file、pandas、matplotlib)
今年其实还是写了不少代码的，但来去还是在玩pandas，就今年重点解决的一个问题来分享。数据说明温室气体设备PICARRO打包的数据格式为bz2文件，包括了CO2,CH4,CO，水汽，切换的阀口等数据。数据处理的总设计数据的读入部分使用bz2file包来读bz2压缩文件，返回dataframe格式的数据，部分代码如下。#Author Wu Dongqiao 2020.11.18import osimport numpy as npimport pandas as pdimport ti
2020-12-31 22:40:13
790
原创 风玫瑰图的绘制，基于气象A文件（windrose）
风玫瑰图的绘制这两天研究了下python中风玫瑰图绘制，有几种方法，比较方便的是引用windrose模块。数据的获取，从气象A文件中获取，为计算年数据将每年A文件放入文件夹，首先放入效果图。代码如下：#Author:Wu Dongqiao#2020.6.6import osimport calendarimport datetimeimport numpy as npfrom matplotlib import pyplot as pltimport matplotlib as mp
2020-06-07 00:02:05
4599
4
原创 matplotlib与seaborn，箱形图绘制比较
当我开始研究seaborn时我惊呆了， seaborn之于matplotlib，相当于pandas之于numpy呀。真是强又简洁，终于不用幸苦修改参数，图画得要更漂亮了，能用一句话搞定的，python绝不说多句。下面用我绘的一个图作例子，matplotlib绘图用了30行， seaborn用了9行。原数据：2007年1月至2018年12月的平均颗粒物数据matplotlib绘图（箱型图）...
2019-08-19 22:33:47
1479
1
原创 将python代码打包成exe可执行文件(pyinstaller)
安装pyinstaller包pip install pyinstaller对编好的python代码进行打包打开cmd命令行窗口， cd到代码所在的目录，
2019-08-18 20:43:22
417
原创 文件的复制批处理（地面气象月上传数据打包程序）shutil、datetime、os、python
文件的复制批处理程序因为工作中的需要，可能需要对指定文件夹中的数据进行复制打包，其中我碰到的是在指定文件夹下打包上个月形成的数据，到分好类的13个文件夹中，有一个特点是需打包的数据为上个月的数据，文件名中都包括‘201907’这样的字符串。主要用到了三个包，os进行路径的获取与文件夹的新建,datatime用于获取时间,shutil用于复制文件。感兴趣的同志，我打包成exe测试好后可以发给你...
2019-08-13 00:32:33
363
原创 pandas对反应性气体批处理，体积浓度转质量浓度apply、lambda，小时、天统计resample
由反应性气体仪生成的臭氧、一氧化碳、二氧化硫、一氧化氮、二氧化氮、氮氧化物等数据为分钟级，从仪器上可以下载数据格式，导出的excel格式数据，一次支持最多64000行。所以我导出了几个月的数据用了4个excel表格。编写程序用统计小时平均，24小时平均值。同时处理其中的异常值，并把体积浓度转化为质量浓度。ppb、 ppm–>ug/m3。首先excel另存为csv文件，虽然pandas也...
2019-07-11 21:23:47
541
原创 年月日季平均运算，pandas中关于时间序列处理，resample
时间序列生成指定范围的日期，例：pd.date_range(‘4/1/2018’,‘4/5/2018’)指定起始或者结束日期，并指定步长(单位是“日”)，例：pd.date_range(start=‘4/1/2018’,periods=5)pd.date_range(end=‘4/5/2018’,periods=5)将时间列设为索引，时间的标准格式对应类似于20190101这样的格式，...
2019-07-10 20:23:26
16408
10
原创 pandas中关于apply、lambda的应用。
lambdalambda：输入是传入到参数列表x的值，输出是根据表达式(expression)计算得到的值。比如：lambda x, y: xy #函数输入是x和y，输出是它们的积xylambda x :x[-2:] #x是字符串时，输出字符串的后两位lambda x :func ...
2019-07-09 11:16:25
31910
2
原创 微波辐射数据绘图实现，matplotlib,pandas
pcolor彩色图绘制#Author:Wu Dongqiaofrom matplotlib import pyplot as pltimport numpy as npimport pandas as pdfrom pylab import *import tkinter.filedialog#设置微波辐射数据路径选择要读取的_lv2.CSV文件#wbfs_file_path='...
2019-07-05 11:40:24
1689
9
原创 解析微波辐射计数据，并上传到mysql数据库（pandas应用及mysql操作）
创建mysql数据库，从微波辐射数据文件_lv2中提取，廓线数据，积分数据，云底高度及雨状态，并上传数据到数据库1. 数据库创建#Author: Wu dongqiaoimport pymysql as msq# 打开数据库连接db = msq.connect("localhost", "root", "12345678", "wbfs", charset='utf8' )# 使用c...
2019-07-04 22:42:02
1373
1
原创 读取pumml文件（颗粒物数据）文件读取,numpy及pandas应用
颗粒物质量浓度数据文件名格式类似为（Z_CAWN_I_57596_20170701010000_O_AER-FLD-PMMUL.txt），为每小时上传文件，本程序可以从文件夹中读取所有pm2.5,pm10数据，去除掉异常值，并计算日平均值。#Author: Wu dongqiaoimport pandas as pdimport numpy as npimport os#定义读取文件函数...
2019-07-04 21:46:05
539
原创 提取文件夹中所有地面气象数据A文件的小时雨量数据
提取文件夹中所有A文件的小时雨量数据地面观测数据一般可以在A文件中找到小时数据，以提取雨量为例，实现从文件夹中提取所有小时雨量值和年月日时格式的时间，并形成csv文件以便后期处理 。代码如下:#从一个文件夹中读取所有A文件，提取小时雨量值和年月日时格式的时间，并形成csv文件以便后期处理#Author: Wu dongqiao July 3, 2019import pandas as ...
2019-07-04 21:04:57
1819
2
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人