qq_38607066的博客

原创 文章记录笔记

Ocean surface energy balance allows a constraint on the sensitivity of precipitation to global warming背景：气候模型预测，未来气候变暖，降水增加。但降水增加是否是对气温增加的响应仍然在争论。量化降水的历史变化是为了更好地预测未来降水的变化。但是气候模型对历史降水气温敏感性的符号和量级均存在差异。温度上升，大气水汽含量上升。水汽一方面吸收地面长波辐射，另一方面反射太阳短波辐射。大气水汽含量上升可以解

原创 MATLAB 绘图

文章目录1. 图形对象的理解1.1 对象间的层次结构1.2 对象的查询和修改1. plotSyntaxInput ArgumentsOutput ArgumentsExample2. 线属性 Line PropertiesLineMarkersGeographic Coordinate Data 地理坐标数据LegendAnnotation — Control for including or excluding object from legend其它3. 标签、标题、文本、注释、图例xlabelylab

原创 CMD命令

文章目录文件目录参考：WINDOW常用命令https://www.cnblogs.com/kekec/p/3662125.html文件目录cdcd …

原创 VPD-Forcing10年循环时间步长匹配问题

10. 测试forcing循环时的时间步长问题问题：forcing是1901-1910年循环。运行1911-1920年时，发现1920-12-31日的最后一个步长（1920-12-31 23:30:00-1921——01-01 00:00:00）之后应当是1921-01-01日的第一个步长（1921-01-01 00:00:00——1921-01-01 00:30:00），发现时间戳对应不上。由于是10年循环，1921-01-01日的第一个步长（1921-01-01 00:30:00）的时间戳是1901-

原创 LINUX VIM 笔记

跳转到指定的行命令行模式下输入（n为指定的行号）：（1）ngg / nG（2）:n（3）vim +n filename（注意这里要输入 + 号）三者的区别在于，ngg/nG输入以后不需要按下回车键，而输入:n之后还要按下回车键才会跳转，vim +n filename只能在打开文件时跳转到相应的行。...

原创 CESM2笔记

文章目录1. 自定义Namelist和Output fields2. 复制案例create_clone3. RESTART重启和BRANCH 分支运行1. 自定义Namelist和Output fields参考链接：CESM 1.2 用户指南：CESM 1.2 CAM Namelist Variablesavgflag_pertape （平均值、瞬时值、最大值、最小值）Sets the averaging flag for all variables on a particular his

原创 Numpy基础

文章目录ndarray1.1 创建numpy数组 ndarray1.2 ndarray的数据类型1.3 数组与标量之间的运算1.4 索引和切片 0-based1.5 切片索引1.6 布尔型索引ndarrayndarray：numpy的N维数组对象，其中的所有元素必须是相同的数据类型属性和方法的区别：属性没有参数，不需要括号；方法有参数，有括号2. 属性ndarray.shapendarray.dtype方法ndarray.astype()1.1 创建numpy数组 nda

原创 Markdown Latex

文章目录Markdown 首行缩进Markdown 首行缩进Markdown原生功能就不支持缩进找了很多方法，最终选择了 因为一个 的宽度刚好是一个中文字符的宽度两个 即可替代传统的缩进两个字符，效果如下  成功缩进两个字符...

原创 MATLAB常见函数

文章目录分组1.tabulate分组1.tabulateFrequency tableSyntaxtbl = tabulate(x)tabulate(x)Descriptiontbl = tabulate(x) creates a frequency table of data in vector x. Information in tbl is arranged as follows:1st column — The unique values of x2nd column — Th

原创 多变量微积分

112

原创 LINUX 常见操作

1. Linux统计文件夹、文件数量的命令查看当前目录下的文件数量（不包含子目录中的文件）ls -l|grep "^-"| wc -l查看当前目录下的文件数量（包含子目录中的文件） 注意：R，代表子目录ls -lR|grep "^-"| wc -l查看当前目录下的文件夹目录个数（不包含子目录中的目录），同上述理，如果需要查看子目录的，加上Rls -l|grep "^d"| wc -l查询当前路径下的指定前缀名的目录下的所有文件数量例如：统计所有以“20161124”开头的

原创 TIFF NC HDF数据的正反问题

%% TIFF NC HDF数据的正反问题clcclear alllat1=zeros(360,1);for i=1:360 lat1(i,1)=90+90.0/360.0-(180.0/360.0)*i;end% 数组为正，从上到下为90：-90，从左至右为-180:180data = repmat(lat1,1,720);data(:,1) = -180;data(:,end) = 180;%% TIFF% write 保证数组为正，lat*lonR=geo

原创 ERA5批量下载

教程链接：https://cds.climate.copernicus.eu/api-how-to#install-the-cds-api-keyWINDOWS前提条件：有anaconda，注册CDS或者ADS，有账号获取API KEY新建.cdsapirc文件，内容为你的API KEY。一般位于C:\Users\Username folderpip install cdsapi准备脚本下载我的.cdsapirc文件，url: https://cds.climate.copernicu

原创 NCO笔记

调用nco包module load nco/4.6.0-icc-14.0.2ncdumpncdump -h **.ncncdump -v a variable name **.nc | lessncks从nc文件中抽取部分变量，保存为新的nc文件。ncks -v lat,lon,time,landmask,area,QBOT FHIST_BGC_hist10.clm2.h0.1901-01-01-03600.nc 1901.ncncrcat将多个nc文件按公共字段时间顺序连接ncr

原创 NetCDF读写

文章目录简介1.1 文件格式1.2 限制文件的结构与性能2.1 经典格式的部分数据类型3.1 外部数据类型3.2 用户自定义类型组件4.1 数据模型4.2 Enhanced Data Model in NetCDF-4/HDF5 Filesnetcdf的使用简介netcdf：数组形式，数据访问、存储、检索、自描述、可移植、通用1.1 文件格式CDF-1：经典格式，默认。适用范围最广，单个文件不超过2GBCDF-2: 64-bit offset formatHDF5 format / netCD

原创 陆地生态系统生态学原理

1. 生态系统概念1.1 引言生态系统生态学将生物有机体及其环境作为一个整体，强调它们之间的相互作用。1.2 生态系统综述经由有机体和物理环境的能量流和物质流为我们理解地球物理环境和生物过程的形式和功能的多样性提供框架生态系统生态学试图了解调节生态系统中物质和能量库（总量）和通量（流）的因子。这些物质包括碳、水、氮、磷、杀虫剂、放射性核素等。库分为生物库和非生物库：生物库：植物、动物、土壤微生物非生物库：土壤、岩石、水、大气生态系统过程是指能量和物质从一个库向另一个库的转移。生态系统过

原创 气孔导度

原创 论文笔记：气候变化、水循环和植被动态

植被与水文，水循环植被影响植物蒸腾，影响蒸散发，影响降水，影响径流和水储量，进而影响水循环，水资源植被动态和水循环的交互：目标：未来气候条件下，植被如何调节水通量？在气候引起的水胁迫下，植被如何变化？水分传输：水力过程-水的潜在梯度驱动干旱对水分传输的影响：干旱导致土壤水分下降，栓塞效应，导管输送水分的能力下降，目前，土壤湿度和大气湿度VPD是高度相关的。但在未来气候情境下，土壤湿度有升有降，但大气干旱普遍上升（气温上升）。如何优化参数？植被影响气候植被通过

原创 FHIST_BGC方案错误

1. read_variable_2d ERROR:sum of PCT_NAT_PFT not 1.0 at nl =12原因：这个可能是模式读取的地表资料存在问题，从 lnd_in 里找出来哪个文件里包括 PCT_NAT_PFT 变量，看看是否满足要求利用ncdump -h ****.nc查找lnd_in中所有的输入文件，发现landuse.timeseries文件中存在PCT_NAT_PFT计算PCT_NAT_PFT的和，发现在1893年并不等于1。即该输入文件下载错误。解决方法：重新下载替

原创 我的Fortran之路

文章目录环境入门例子规范数据类型环境ubantu服务器自带的gfortran本地电脑的VS code入门例子编程：新建一个addnumbers.f90文件，键入一下内容program addnumbers ! This is sample program adds two numbers implicit none ! Type declaration statements real ::a,b,result ! executable stat

原创 CESM2：基本框架

文章目录1. 基本框架2. CESM components3. grid4. compset4.1 命名规范：4.2 分类4.3 options4.4 Examples of simulations using different compsets:1. 基本框架CouplerAtmosphere ModelsCAM,CAM-CHEM,WACCM Community Atmosphere ModelDATM Climatological Data ModelLand Mo

原创 天河2的基本操作

编译程序编译在登陆节点上进行，程序运行在计算节点上进行。一个case的编译进程不超过4。不要同时编译多个进程。在计算节点编译目前普通计算节点不支持编译的，如您需要在计算节点上编译的话，可使用rhenv分区的节点， 可使用命令“ yhrun -n 1 -N 1 -p rhenv ./case.build --skip-provenance-check &> log.case.build ”提交到rhenv分区下编译（rhenv分区节点是更完整的操作系统），提交作业时可使用“-p rh

原创 CESM2 实验笔记

文章目录1. 重新启动1.1 关于重启文件2. RUN_TYPE 初始化运行方式startupbranchhybrid3. branch 或者 hybrid运行的示例4. 其他声明5. 关于RESUBMIT1. 重新启动参考：https://bb.cgd.ucar.edu/cesm/threads/how-to-restart-cesm2-1-3-fhist-from-a-backup.5245/#post-36275适用情况：模型运行因时间不足或其它原因中断，需要继续运行方法1：将CONTIN

原创 FHIST_BGC方案

任务：月值输出：径流量、土壤湿度、gpp、npp、nep、lai、总辐射、净辐射、植物蒸腾、土壤蒸发、植被冠层截留蒸发、感热半小时输出：湿度、温度case: FHIST-BGC1890-1900：预热10年，输出年值1901-1920：历史输出前20年，月值和半小时1921-2000：历史输出后80年，月值变量名是否有标准输出径流量yesQOVER总径流， QH2OSFC地表径流土壤湿度yesH2OSOIgppyesGPPnppyes

原创 天河安装CESM2.1.3

[sysu_syli_3@lon26%tianhe2-H FHIST_BGCtest]$ ./preview_namelists - Prestaging REFCASE (/BIGDATA1/sysu_syli_3/CESM2/inputdata/cesm2_init/b.e20.BHIST.f09_g17.20thC.297_01_v3/1979-01-01) to /BIGDATA1/sysu_syli_3/CESM2/scratch/FHIST_BGCtest/runCopy rpoi...

原创 CESM2.1.3移植 XML文件的填写

config_machines.xml <machine MACH="tianhe2"> <DESC>NSCC-gz Tianhe2, os is RedHat/Centos, 24 pes/node, batch system is Slurm</DESC> <NODENAME_REGEX>tianhe2</NODENAME_REGEX> # 此处的应该填写hostname.*， # 例如 运行hostname，为

原创 验证是否正确迁出CESM2

劭怡老师的CESM2迁出似乎有问题验证方法1参考：https://escomp.github.io/CESM/versions/cesm2.1/html/downloading_cesm.htmlTo confirm a successful download of all components, you can run checkout_externals with the status flag to show the status of the externals:./manage_exter

原创 FWHIST 实验

大气化学模块[sysu_syli_3@ln42%tianhe2-H FWHIST-BGC]$ ./xmlquery CAM_CONFIG_OPTS CAM_CONFIG_OPTS: -phys cam6 -age_of_air_trcs -chem waccm_tsmlt_mam4./xmlchange --append "CAM_CONFIG_OPTS=-phys cam6 -chem none"[sysu_syli_3@ln42%tianhe2-H FWHIST-BGC]$ ./xmlque

原创 BHIST实验

1. create_newcasevim B1850.shsource BHIST.sh./create_newcase --compset BHIST --res f09_g17 --case /BIGDATA1/sysu_syli_3/CESM2/scratch/BHISTtest &> BHISTtest.create_newcase.log查询PE配置：[sysu_syli_3@ln42%tianhe2-H BHISTtest]$ ./preview_runCASE

原创 F 实验运行过程

三个case：FWHIST_BGC : HIST_CAM60%WCTS_CLM50%BGC-CROP_CICE%PRES_DOCN%DOM_MOSART_CISM2%NOEVOLVE_SWAVFWmaHIST_BGC : HIST_CAM60%WCCM_CLM50%BGC-CROP_CICE%PRES_DOCN%DOM_MOSART_CISM2%NOEVOLVE_SWAVFW1850 : 1850_CAM60%WCTS_CLM50%SP_

原创 Linux wget 下载

wget 下载# "nohup+&"后台不挂断的运行nohup wget xxxx &nohup wget https://zenodo.org/record/400950/files/agx716_undetermined_reads.fastq.gz & vim temp # 全选网页文件,复制后粘贴至vim编辑器grep '.fastq.gz' temp > temp1cat temp1 | perl -ne 'print "wget https://z

原创 运行时的错误

错误1：Building cesm with output to /BIGDATA1/sysu_syli_3/CESM2/scratch/mycase/bld/cesm.bldlog.201008-213936 ERROR: Error gathering provenance information from manage_externals.manage_externals error message: ERROR:root:SVN returned invalid XML messag

原创 Linux 笔记：文本编辑器

文章目录vi or vim另外的两款文字编辑软件Using 'pico'Using 'joe'vi or vim两种模式：文本模式 （按i进入）——更改文本内容命令模式 （按Esc进入）——输入命令创建一个新的文本文件vi example.txt退出与保存:w # write 保存:q # quit:wq # write and quit:q! # 强制退出另外的两款文字编辑软件Using ‘pico’Using ‘joe’...

原创 英语笔记：好词好句

尽己所能做某事 do everything in their power to do sth.For this reason they have done everything in their power to ‘warn’ the public that Linux is something that you probably don’t want and don’t need.

原创 CESM2环境搭建

文章目录学院服务器自带学院服务器自带(base) [19fush@login sources]$ perl --versionThis is perl, v5.10.1 (*) built for x86_64-linux-thread-multiCopyright 1987-2009, Larry WallPerl may be copied only under the terms of either the Artistic License or theGNU General Pub

原创 IDL编程基础

文章目录快捷键快捷键按键功能$换行;注释ctrl+;批量注释ALT+/补全函数ctrl+shift+F自动对齐（搜狗输入法快捷键冲突）alt+?联想查词、工作空间工程过程函数帮助：帮助-帮助内容； 选中，双击，F1； 控制台，输入?；demo...

原创 张扬建：遥感与生态学的结合

张扬建,范春捆,黄珂,刘瑶杰,俎佳星,朱军涛.遥感在生态系统生态学上应用的机遇与挑战[J].生态学杂志,2017,36(03):809-823.遥感作为生态学研究从单点尺度向区域尺度扩展的关键技术手段，在当今的生态学研究中，不仅可以为地面监测提供补充资料，甚至具备替代实地测量的潜力。特定尺度：遥感技术的特点决定了其在生态学特定尺度能发挥更佳研究效用。在个体生态学，群落生态学和种群生态学尺度上...

原创 PM2.5与城市化的关系

一、 吴浪,周廷刚,温莉,刘晓璐,朱晓波.基于遥感数据的PM2.5与城市化的时空关系研究——以成渝城市群为例[J].长江流域资源与环境,2018,27(09):2142-2152.研究数据：1998-2013年夜间灯光、NASA PM2.5栅格数据集、统计年鉴（）研究方法：GWR、 地理探测器研究区：成渝城市群研究内容：基于GWR模型，探讨了PM2.5浓度与夜间灯光（一定程度上表示城市...

原创 ggplot2 颜色

查看ggplot2 支持的色带library(RColorBrewer)display.brewer.all(type = "seq")display.brewer.all(type = "div")display.brewer.all(type = "qual")翻转色带# define a color rmap palettecolr = colorRampPalett...

原创 R 地图

参考：https://stackoverflow.com/questions/33227182/how-to-set-use-ggplot2-to-map-a-rasterlevelplot方法：https://www.rdocumentation.org/packages/rasterVis/versions/0.46/topics/levelplot-methods利用rasterVis:...