/ Recent content in Home on Shu Hugo -- gohugo.io zh-cn Tue, 16 Nov 2021 00:00:00 +0000 /zh/eosdata/ Tue, 16 Nov 2021 00:00:00 +0000 /zh/eosdata/ 首先你要确保有已经拥有earth Data的账户。如果没有，请先申请，申请链接在此https://disc.gsfc.nasa.gov/#top 方法1：wget在Linux/Mac上的使用 本节基本内容来自NASA的数据说明 https://disc.gsfc.nasa.gov/data-access 1.1 安装wget 在mac上，需要事先安装homebrew，然后再安装wget。 安装wget方法只需一行命令： brew install wget 1.2 下载 1.2.1 明文密码下载 wget --user=<uid> --password=<password> <url> <uid>，<password>和<url>分别为用户名，密码和欲下载的文件/文件夹链接。 明文下载是最简单粗暴的方式，但是用户名和密码都是明文的，如果写成脚本存储在本地也可以，但就是对自己的账号不太安全。但如果你对自己的密码没那么在乎，你可能不需要看下面一节的内容。直接跳过下一节，查看更复杂的wget参数一节。 1.2.2 用Cookie方式下载 1.2.2.1 创建.netrc文件 步骤为： cd ~ touch .netrc chmod 0600 .netrc echo "machine urs.earthdata.nasa.gov login <uid> password <password>" >> .netrc 此处<uid>是你的用户名，<password>是你的密码。至此，你已经准备好了.netrc文件，保存好了你的用户名和密码。 此处我们距离使用了远程主机名称为urs.earthdata.nasa.gov，请根据你下载数据的具体服务器名称修改。 1.2.2.2 创建.urs_cookies文件 cd ~ touch .urs_cookies 如果你有过错误的登录，此处可能需要重新创建.urs_cookies文件。 1.2.2.3 用cookie下载文件 注意：<url>就是你要下载的目标文件或者文件夹。 单个文件/文件夹 wget --load-cookies ~/.urs_cookies --save-cookies ~/.urs_cookies --auth-no-challenge=on --keep-session-cookies <url> 如果你的<url>是文件名字，则以上命令下载该文件，并存储到当前命令行所在的目录。 /zh/20200328_naturepapers/ Sat, 28 Mar 2020 00:00:00 +0000 /zh/20200328_naturepapers/ 反复检测才是战胜新冠病毒(2019 nCoVID)的关键 /zh/20201128_netcdf/ Sat, 28 Mar 2020 00:00:00 +0000 /zh/20201128_netcdf/ The following packages are required to build netCDF-C using CMake. netCDF-C Source Code CMake version 2.8.12 or greater. Optional Requirements: HDF5 Libraries for netCDF4/HDF5 support. libcurl for DAP support. PnetCDF libraries for parallel I/O support to classic netCDF files If you have libraries installed in a custom directory, you may need to specify the CMAKE_PREFIX_PATH variable to tell cmake where the libraries are installed. /zh/20200205_rscopus/ Wed, 05 Feb 2020 22:05:38 +0000 /zh/20200205_rscopus/ 以获取武汉病毒所所有博导的发文章数据为例。 一 获取名单。 列出你需要获得数据的人名和研究所。 人名和单位名字都必须是英文。 从此链接进入，武汉病毒所导师介绍。 可以看到31位博士生导师的姓名，通过手动或者软件获得这些导师的人名拼音——即英文名。将这些英文名写入一个txt文件。下面切图我隐去了他们的姓名。 如： 从scopus获取名单上教授的发文章数据。 直接使用一下代码，其中需要路径，使得程序可以找到第一步中存储的名字列表。 affil_id = ‘60027289’ 是武汉病毒所的单位代码。如果需要其他研究所，请自行查找单位代码。 aff = ‘virology’ 是病毒所的英文名关键字。 api_key =‘xxx’ 是你需要获取scopus一个帐号，以此来获取数据。 rm(list = ls()) library('rvest') library('dplyr') library('scholar') # install.packages('scholar') library(rscopus) # source('Rcode/function/function.GS.R') dir.data = 'wuhan' fn.namelist = 'wuhan/BoDao.txt' outdir = 'wuhan/out' # Read names from fn.namelist # download the profile html pages from GS, save them to local drive. dir.create(outdir, showWarnings = F, recursive = T) dir.create(dir.data, showWarnings = F, recursive = T) tmp = readLines(fn. /zh/20200203_wushui/ Mon, 03 Feb 2020 22:05:38 +0000 /zh/20200203_wushui/ 现在几个批评文章意思差不多，基本主旨是这样： 学术能力不够，却当所长。证据：只有24篇文章，11篇通讯文章，说不够分量。 年龄太小，资历不够。 证据：1981年生，2005年大学毕业，2010年博士毕业。 王延轶和导师舒红兵是夫妻关系，上位过程都是靠她老公。 可能是大学一毕业就确立了夫妻关系。 给了一个“双黄连之母”的标签。 我逐条批驳。 第一条，文章数量 真实数据：47篇文章，18篇是顶刊(PNAS及以上期刊)。19篇为一作或通讯作者，19篇中，6篇为顶刊。 1.1 她的文章本身并不差。 请看SCOPUS上对她文章的统计：[Author details (Wang, Yanyi)](https://www.scopus.com/authid/detail.uri?authorId=7601493171） SCOPUS上的数据，标志了重要时间点 王延轶个人网页自述有28篇文章（一作、通讯、N作）引用1700；实际数据应该是47篇，引用3410次。有三个文章数：24， 28， 47，分别来自病毒所介绍、批评文章和SCOPUS，SCOPUS数据更新最快，应该以此为准。朋友圈的批评文章有意还是无意采用的是百度的页面，24篇，不知道是那一年的数据了。 顶刊有：6篇PNAS, 4篇Cell Host and Microbe, 3篇Cell Research，3篇Immunity，图上缺一篇NC，一篇Molecular Cell。一共18篇大于9影响因子的文章，PNAS刚好就是9.5，所以我选9作为阈值。\ 批评文章里面说她作为通讯作者只有11篇，她作为导师尤其是中科院的导师，因为中科院的学生数量少，所以11篇通讯文章已经不少了。实际上她一作和通讯一共有19篇。 拿科学院研究的文章数量和任何一个大学同级别教授比较，科学院都要输——因为科学院学生数量少。 scopus下载的详细数据， 到2020为止。后面影响因子是我查询的数据 我们一边批评学术评价体系只看发文章的八股式科研制度， 同时又以文章数量少来批评或者指责一个人学术不合格。即便是做吃瓜群众也不可以这么双标吧——何况人家的文章不少。 /en/20200201_nature/ Sat, 01 Feb 2020 00:00:00 +0000 /en/20200201_nature/ Stop using the offensive virus names Inappropriate naming of the new human disease inevitably causes prejudice and discrimination upon particular people and undue fear of the diseases on global social media, although unintentionally. After the outbreak of Novel Coronavirus (2019-nCoV) causing a severe respiratory illness, some of the research papers and news on nature (such as Nature 577, 450 (2020), Nature 577, 605-607 (2020), doi: 10.1038/d41586-020-00146-w, doi: 10.1038/d41586-020-00180-8, doi: 10.1038/d41586-020-00253-8 and doi: 10. /zh/20200201_nature/ Sat, 01 Feb 2020 00:00:00 +0000 /zh/20200201_nature/ Nature在报道新冠病毒时候，使用的名字都很刺耳，于是我写了一篇短评论给Nature。 Stop using the offensive virus names Inappropriate naming of the new human disease inevitably causes prejudice and discrimination upon particular people and undue fear of the diseases on global social media, although unintentionally. After the outbreak of Novel Coronavirus (2019-nCoV) causing a severe respiratory illness, some of the research papers and news on nature (such as Nature 577, 450 (2020), Nature 577, 605-607 (2020), doi: 10.1038/d41586-020-00146-w, doi: 10.1038/d41586-020-00180-8, doi: 10. /en/20200110_calender/ Fri, 10 Jan 2020 12:05:38 +0000 /en/20200110_calender/ Is the current calendar is the best or only way to represent days? The days we experienced, undergoing, and expecting are recorded a piece of paper; we call it a calendar, where marks days of years, months, weeks, and holidays. The current calendar appearance that organized with month and weeks are very efficient to display the days of week and month. But, the calendar drops the original meaning of “Month”, which was related to the moon and moon phase, while the lunar calendar and lunisolar calendar (such as Chinese Calendar) still conserve the sense in the month. /zh/20191202_xueshuwubi/ Mon, 02 Dec 2019 22:05:38 +0000 /zh/20191202_xueshuwubi/ 我经常在想，金庸笔下的挥刀自宫究竟是什么意思？ 金庸笔下的江湖。并不是一个空穴来风的虚构世界，而是将现实社会抽象化之后的江湖。金庸笔下说的江湖说的每一个门派、每一个大侠、每一个普通人，在我们现实中都可以找到类似的例子。 在金庸笔下，江湖中有若干绝世武功，都可以通过艰苦卓绝的努力而成为一代大侠。比如，岳不群说的紫霞神功，比如风清扬的独孤九剑，比如黄药师的九阴白骨爪，比如武当山的九阳神功。这些都是绝世武功，但想要获得这样的绝世武功和成为一代大侠，一需要有名师传授，二需要刻苦练习的毅力，第三还需要有天赋，有悟性。所以这些条件限制了一个人获得成功、获得绝世武功的能力。 青城山的余沧海，塞北明驼木高峰，君子剑岳不群，嵩山左冷禅，他们都没有这个啊，机会、毅力、条件通过这条路径获得他们想要的名和利。于是他们想获得一个可以获得绝世武功的秘籍，这个秘籍并不需要名师指导，似乎也不需要刻苦练习，只需要获得这个秘籍就可以获得成功，他们都走上了追求武功秘籍的道路，放弃了勤勤恳恳而获得成功的途径。 那么，什么是挥刀自宫呢？ 江湖中每一个人都有武功门派，也都有师傅指导，只要通过自己艰苦努力，或多或少总会进步，也有机会获得更大的成功，也可能成为一代大侠，比如天资并不聪明，江南七怪的徒弟郭靖。但是这样的过程，需要艰苦的努力，需要很长的时间，甚至还要赌上一些运气。并不是每一个年轻人都愿意走这一趟充满不确定性的艰苦道路。 对成功的渴望，与现实的焦虑压迫着每一个江湖中的人。于是有的年轻人——甚至是门派掌门，就开始思考用其他更简便易行、离经叛道的方式来获得他们的名利。他们首先相信一点要获得成功，首先要获得《葵花宝典》或者《辟邪剑谱》这样的绝世武功典籍。他们相信，获得《葵花宝典》和《辟邪剑谱》是传说中的人物成功的基本条件。于是他们千方百计要首先获得剑谱。 于是江湖中有很多人都拿到了《葵花宝典》和《辟邪剑谱》。比如百年前华山派、魔教、任我行都拿到了《葵花宝典》，但是魔教、华山派、任我行都没有真正修炼这部剑谱，为什么？那是因为《葵花宝典》的开篇第一句是“欲练神功，必先自宫”。 为什么会拒绝挥刀自宫呢？因为挥刀自宫就意味着要放弃自己，要放弃作为人所最基本的那些东西；放弃了自己对自己的认同，也放弃了自己底线的底线。在金庸笔下，一个人若能够若敢于挥刀自宫，那么人世间一切的亲情、恩情、正义都不存在都不重要。只要你敢挥刀自宫那你就放弃了，除了名利之外这个世界上一切的东西。 那么现实中，谁是那个挥刀自宫的人呢？ 在学术界那些剽窃、抄袭、篡改数据、篡改结果的人，不就是挥刀自宫的人吗？选择舞弊的人，首先放弃了通过艰苦卓绝努力而获得学术进步和影响力，而选择那种最轻易，却不道德的方式，来获得他所想要的名利。当一个做学术的人，愿意抄袭的时候，也就放弃了作为学术工作者的底线。 在学术界那些剽窃、抄袭、篡改数据、篡改结果的人，不就是挥刀自宫的人吗？放弃了一个学术人为学术人的最基本底线，当然，这样的挥刀自宫，可以获得成功可以发表很多文章，可以发表重要文章，一段时间成为一代学术上的武功高手。 但是和那些真的拿起刀子切鸡鸡的人一样。这样在学术上挥刀自宫的人，同样不会再有后代——他如何去教他的学生，他的门徒重复和他一样的实验，发一样的论文？ 当有一天，那些挥刀自宫的学术工作者发现竟然切了鸡鸡，那他就丢尽了脸面，人人都知道他是没有底线的败类。那未来他在学术的江湖中，人人都会躲着他，唯恐与他有任何瓜葛，他在学术界将无立锥之地，他将会被人唾骂，遗臭万年。 现实中想成为学术大侠和江湖中的武功大侠一样，需要艰苦卓绝的努力才可获得。但时常艰苦的努力，并不能让你心想事成，但是挥刀自宫却可以让你轻易获得那份忐忐忑忑的荣耀。 《葵花宝典》放在你面前，你愿意挥刀自宫吗？ 我知道这世界上有人愿意，比如他们： /en/20191119_sundials/ Tue, 19 Nov 2019 22:05:38 +0000 /en/20191119_sundials/ SUNDIALS: SUite of Nonlinear and DIfferential/ALgebraic Equation Solvers SUNDIALS (https://computing.llnl.gov/projects/sundials) is a very powerful mathematical library that is helpful in solving engineering and scientific problems efficiently. CVODE is one of the tools, aiming to solve the Ordinary Differential Equation, in C language. In the Simulator for Hydrologic Unstructured Domains (SHUD) model, the system of ODEs describing the hydrological processes is fully coupled and solved simultaneously at each time step ($ \Delta t=t_{n}-t_{n-1}$) using CVODE, a stiff solver based on Newton-Krylov iteration. /zh/20191119_sundials/ Tue, 19 Nov 2019 22:05:38 +0000 /zh/20191119_sundials/ SUNDIALS: SUite of Nonlinear and DIfferential/ALgebraic Equation Solvers SUNDIALS (https://computing.llnl.gov/projects/sundials) 是一个很强大的数学库，用于求解工程和科学上的数学问题。CVODE是其中重要工具之一，用于求解常微分方程(Ordinary Differential Equation, ODE)，由C语言写成. SHUD模型(Simulator of Hydrologic Unstructured Domains)使用有限体积法求解流域控制方程，即理查德方程(Richard Equation)和圣维南方程(Saint Venant Equation)。 当前最新的SUNDIALS是其v5.x版本，2019年11月发布. 本文介绍如何在Mac/Linux上安装SUNDIALS/CVODE。 本安装过程以及在 Mac 10.13 and **Ubuntu 18.04.3 LTS (GNU/Linux 4.15.0-66-generic x86_64)**上进行测试。 下载 SUNDIALS，地址：https://computing.llnl.gov/projects/sundials/sundials-software. 或者也可以从github上获取最新源代码： git clone https://github.com/LLNL/sundials.git 此刻，你应该已经在你的目录中有sundials的源码了，请进入该目录。假设你已经进入了 sundials所在的目录. 创建安装目录 安装目录名称为：installSundials 执行以下命令： mkdir installSundials cd installSundials rm -rf builddir instdir srcdir mkdir builddir mkdir instdir mkdir srcdir cd builddir/ 利用CCMAKE配置SUNDIALS安装参数. 有两种方式配置其安装参数，一种是图像界面(GUI)，另一种是直接命令行参数赋值。 方法1：, 键入命令： ccmake . /zh/20191121_quadraticequation/ Tue, 19 Nov 2019 22:05:38 +0000 /zh/20191121_quadraticequation/ 我忘了写了~~~ 一定补上，一定补上 /en/20191003_r_animation/ Thu, 03 Oct 2019 00:15:07 +0000 /en/20191003_r_animation/ This script is to demonstrate the capability and scheme to realise the animation. You have to prepare your own data for this animation script. Script for Global NDVI rm(list=ls()) library(RoundAndRound) # devtools::install_github('shulele/RoundAndRound') col.fun=function(x){ colorspace::terrain_hcl(x, l=c(65,95), power=c(1/3, 1.5)) } plotStar <- function(x, y, radius, angle=180, xyratio=1, ...){ i <- 1:11 alpha <- 2*pi / 10 r <- radius * (i %% 2 + 1)/2 omega <- alpha * i + angle * pi /180 invisible(mapply(function(a, b, . /zh/20191003_r_animation/ Thu, 03 Oct 2019 00:15:07 +0000 /zh/20191003_r_animation/ This script is to demonstrate the capability and scheme to realise the animation. You have to prepare your own data for this animation script. Script for Global NDVI rm(list=ls()) library(RoundAndRound) # devtools::install_github('shulele/RoundAndRound') col.fun=function(x){ colorspace::terrain_hcl(x, l=c(65,95), power=c(1/3, 1.5)) } plotStar <- function(x, y, radius, angle=180, xyratio=1, ...){ i <- 1:11 alpha <- 2*pi / 10 r <- radius * (i %% 2 + 1)/2 omega <- alpha * i + angle * pi /180 invisible(mapply(function(a, b, . /zh/20191001_nature_greening/ Tue, 01 Oct 2019 00:00:00 +0000 /zh/20191001_nature_greening/ 先上结论：这篇文章是一个非专业记者的不专业报道。 这是篇新闻( Nature News)，并不是篇科学论文，所以见到《Nature》就膝盖软的先平身。 为什么说这篇文章非常不靠谱？ 分四节来说 这个文章写了什么，怎么写的 究竟植树有没有好处？ 几个扩展问题。 这个文章怎么写才靠谱？ 第一节： 这个文章写了什么，怎么写的 文章标题： China’s tree-planting drive could falter in a warming world。 全球变暖或将削弱中国的植树动机。 副标题：Researchers warn that the country’s push to hold back its deserts could strain water resources. 研究者警告，中国的治沙行动可能使水资源不堪重负。 文章标题单看没什么问题。但副标题说治沙行动导致水资源紧张，在文章中缺乏足够的证据。 文章的主标题和副标题分别用了“种树”和“治沙”两个相关却不同的概念，文章内部指责的基本上是“种树”。 文章配图的图1 ​ 配图文字写：中国正通过种树来阻止沙漠扩展 上面这个配图写的是中国种树来阻止沙漠扩张，图片上沙漠和“树“泾渭分明，广袤的沙漠和”树“一线之隔。但实际上，地图上并不是树，甚至很难说是植被，它叫”草方格“，是中国沙漠地区防止沙漠扩展非常有效的手段。 作者先入为主告诉你中国大面积种的就是树、树、树、树，一切都以种树为靶子开始，但是配图是草方格。本文后面有专门对治沙的草方格有介绍。 第10段The growth of forests is significant and necessary progress in the fight against desertification, says Jianping Huang, a climate researcher at Lanzhou University (兰州大学 黄建平，种树可以抵御沙漠化). /zh/20190916_iahs_hydrology/ Mon, 16 Sep 2019 00:00:00 +0000 /zh/20190916_iahs_hydrology/ The 23 unsolved problems in hydrology identified by the community process in 2018. Link of the paper:link Time variability and change Is the hydrological cycle regionally accelerating/decelerating under climate and environmental change, and are there tipping points (irreversible changes)?
 How will cold region runoff and groundwater change in a warmer climate (e.g. with glacier melt and permafrost thaw)?
 What are the mechanisms by which climate change and water use alter ephemeral rivers and groundwater in (semi-) arid regions? /zh/20190627_ndvi/ Thu, 27 Jun 2019 00:00:00 +0000 /zh/20190627_ndvi/ The normalized difference vegetation index (NDVI) is a simple graphical indicator that reflect the greenness and the heath condition of vegetations. The raw data is from NASA MOD13A. The global NDVI variation from 2002 to 2017. The NDVI varies from 2002 to 2017 in US (Continental) How to download the rawdata The code of the monthly NDVI from MODIS dataset is MOD13C2, so type the code below in your terminal and the wget will download all availalble data from server to you local matchine. /en/20190627_ndvi/ Thu, 27 Jun 2019 00:00:00 +0000 /en/20190627_ndvi/ The normalized difference vegetation index (NDVI) is a simple graphical indicator that reflect the greenness and the heath condition of vegetations. The raw data is from NASA MOD13A. The global NDVI variation from 2002 to 2017. The NDVI varies from 2002 to 2017 in US (Continental) How to download the rawdata The code of the monthly NDVI from MODIS dataset is MOD13C2, so type the code below in your terminal and the wget will download all availalble data from server to you local matchine. /zh/20190611_wrfhydro/ Tue, 11 Jun 2019 02:05:38 +0000 /zh/20190611_wrfhydro/ WRF-Hydro做得非常好。 能够完成大尺度(全美国900万平方公里)实时的洪水预报，这个就已经超出所有分布式水文模型的能力了。而且他们在推广上面作了非常好的工作，数据的预处理后处理都有完整的解决办法，还免费开班授课推广WRF-Hydro。 下图用WRF-Hydro做出来的全美国河流流量图。这个图可以实时预测每条河流的流量，非常强大。 它最重要的突破应该是在耦合上面，给WRF的模拟一个更精细的边界层模拟，WRF和NOAH LSM的一层耦合，NOAH与Hydro又一层耦合，这个在时间和空间分辨率上的复杂耦合关系它都处理得非常得当。 也许它最初的名字应该叫WRF-NOAH-Hydro模型，而不是WRF-Hydro。WRF-Hydro模型耦合上做得非常非常好的案例。 我们组上周（2019.06）的在北卡罗来纳Ashevill的项目总结会也是讨论模型耦合，耦合流域、湖泊、农业和社会经济，叫Coupled-Human-Nature System，是NSF的资助，几个大学合作做了3年，有很多进步但还算不上“耦合”，只能到数据单向连接的程度，用于科学研究尚可，不能做应用预报。 美国DARPA还有一个叫MINT的项目，基本思路也是耦合，也还没有到达“耦合”的程度，将来做预报和决策。项目里面做技术的是南加州大学做人工智能方向的人，但是他们不懂也不愿意懂模型，搞不清楚各种模型的数据意义和时空错位问题。所以我挺担心将来他们暴力链接数据，弄出搞笑结果来。 当然了，这种把自然和人类社会耦合模拟，的确比多种自然过程的耦合模拟更难。未来前景未知。 不过WRF-Hydro在水文层面上的处理相当简单了，我觉得像是从大气-陆面-水的设计思路，也很像VIC-Routing的解决方式。 蒸发、下渗、产流、地下水补给都是由NOAH在较粗的分辨率上完成，Hydro部分计算时空分辨率更高，要实现Overland Routing, River Routing，Groundwater baseflow和湖泊。Routing都是一维、单向流动，地下水的处理是不具有物理意义的Bucket Model。湖泊也是个Bucket Model。 Noah和Hydro的交互是通过前时间步长的水量做权重分配。 下面两张图是WRF-Hydro的水文过程的思路和计算办法(左)和VIC Routing(右)的办法。 WRF-Hydro的思路能很好的模拟大气、水、能量在时空上的变化，但是在水文学家来看，水的模拟有很多改进空间（当然是根据具体的需求）。 WRF-Hydro更重视大尺度的洪水问题，主要就是关心河道洪水，所以它的处理办法速度快、计算稳定、模型预热快、参数也容易调。通过计算出来河道的流量、对比测量的河道Rating Curve，可以推算洪水发生的位置和可能性，做出预报。 WRF-Hydro的限制就是除了河道、积雪、土壤水量之外的水文变量——地下水、Baseflow、Percolation等等都不够可靠，更进一步也限制了它与其他水资源、环境、水质、地貌、生态、BioBGC的耦合。所以它关心的终极目标可以说是河道中水的流量； 它可以预报洪水流量，但是对于淹没区等等的信息就不够可靠了。 如果使用isotope和Water Age等数据来挑战它的计算，很可能出现大的偏差。 当然了，能解决上述这些限制的水文模型，必然是求解速度很慢的模型了。 再讲讲水文模型。 对水文模型而言，**物理性(1)、准确性(2)、求解速度(3)**应该是个不可能三角。统计模型和数据挖掘在2/3上做得不错，1就很差了。PIHM，PAWS，tRIBs, ModFlow, ParFlow,GSflow等等在1/2都不错，但是3非常差。1/2/3都比较平衡的也有，简单的如Topmodel，新安江模型，复杂的如PRMS和SWAT。 以上速度慢的都使用了有限体积(FV), 有限元(FE)、有限差分(FD)等等数值求解办法，FV/FE/FD方法利用非常广泛，但是那些最经典应用都得益于其问题设计，求解域内存量和流量在时间、空间上的变化具有稳定性或者惯性，于是理论求解矩阵比较规矩（稀疏矩阵、对称矩阵、条形矩阵等等柔性矩阵），例如工程上的热传递和单纯的河道水动力学问题。 相比之下，水文中将地表、地下和河道的计算同时包含的矩阵就非常刚性(Stiff Matrix)，求解非常困难。这一困难主要是由水文问题的空间特征决定的，空间的不同计算单元的形状，参数的差异性太大了。 即便是Modflow, M3D这种专门针对地下水的简单模型(问题求解使用同一公式)，应对空间异质性都有很高的挑战。 水文使用FE/FD另外一个挑战是变量时间上的变化速度，计算都要满足CFL Condition才能保证收敛性；CFL与变化速度相关，但是河道、地面、地下水等等的流速都有数量级上的巨大差异，于是在求解的时间步长上的要求非常高。当地下水流动可以在小时甚至天间隔上收敛的同时，河道需要在分钟甚至秒的步长上保证收敛，于是同时求解的步长就需要在分钟和秒的级别计算。所以现阶段的最常见提速办法都是分别计算域，然后同步——也就是相当于依次计算多个域，然后不同域在特定时间步长上交换变量，成为各自的边界条件。PAWS和WRF-Hydro的是这样处理时间步长的。 用FE/FD方法水文建模的人，与其说是水文学家，不如说是水文-计算-软件学家。 回到WRF-Hydro， 如果关心实时洪水和大气耦合，WRF-Hydro是非常棒的。 给WRF-Hydro做个广告，想上课的朋友搜“WRF-Hydro Training”，春秋两季都开课，研究生和博士后都可以申请到500美元的Travel Funding。 /en/20190129roundandround/ Tue, 29 Jan 2019 00:00:00 +0000 /en/20190129roundandround/ I planed to write an R package and submit on CRAN, that can plot the planets orbiting in Solar System or satellite/moon orbiting around Earth. That package must plot the moving planets, the moon, moon phase and satellite in 2D or 3D canvas. In the future, it may calculate the orbit based on orbit mechanics. My son loves the song “The Wheels on the bus go round and round“, and he has an obsession with planets. /en/20190129kids-ucd/ Tue, 29 Jan 2019 00:00:00 +0000 /en/20190129kids-ucd/ external_link of the project Motivation of the project This is fantastic program to teach the local kids and lead them into the discovering sciences. Panoview of the Math&Science classroom. Week 03: Week 03: Jessica is writing students observation and discussion of plants in two different soils (Sperpentine and sandy loam). Week 05: Week 05: The growth condition of their beans. /zh/20190129roundandround/ Tue, 29 Jan 2019 00:00:00 +0000 /zh/20190129roundandround/ I planed to write an R package and submit on CRAN, that can plot the planets orbiting in Solar System, or satellite/moon orbiting around earth. That package must plot the moving planets, the moon, moon phase and satellite in 2D or 3D canvas. In future, it may calculate the orbit based on orbit mechanics. My son loves the song “The wheels on the bus go round and round“, and he has an obsession on planets. /zh/20190129kids-ucd/ Tue, 29 Jan 2019 00:00:00 +0000 /zh/20190129kids-ucd/ 项目外部链接 Motivation of the project This is fantastic program to teach the local kids and lead them into the discovering sciences. Panoview of the Math&Science classroom. Week 03: Week 03: Jessica is writing students observation and discussion of plants in two different soils (Sperpentine and sandy loam). Week 05: Week 05: The growth condition of their beans. /about/ Mon, 01 Jan 0001 00:00:00 +0000 /about/ 这里是我的个人博客，内容都繁杂，都是些鸡毛蒜皮的事情。 有关工作的事情或者联系方式，请访问我的工作网址：https://shud.xyz/