水汽稳定度修正函数\Psi_q对潜热通量影响--模式验证工作

我之前提出了一个水汽通量廓线关系，这项工作偏理论，如果对下面说的背景不了解的话可以看下

https://agupubs.onlinelibrary.wiley.com/share/YNSG74MV8B8BAAUMCHN3?target=10.1029/2022JD036708

那会没把提出的水汽稳定度修正函数加到CAS-ESM,当时对CAS-ESM模式还没这么熟悉,也想着师兄能帮我，但是师兄太忙了。现在我对模式已经比较熟悉啦，之前又把师兄加COARE的工作重复了一下，万事俱备只欠东风了。

试验设计：

为了探究温度稳定度修正函数psit和比湿稳定度修正函数对湍流潜热通量的模拟差异，我们采用CAS-ESM2.0全耦合模式，设置了两组实验，除了湍流通量方案中稳定度修正函数以外，两组实验的其他设置都是相同的。我们将采用Psit的COARE算法的CAS-ESM耦合模式作为控制实验（CTRL），我们将采用Psiq的COARE算法的CAS-ESM耦合模式定义为实验一（EXP1）。两组实验都是初始积分（startup）开始,为了防止非平衡初值或扰动的条件下模式调整到平衡态的过程对实验结果影响，我们的两组实验均积分10年，从1995年1月1日到2004年12月31日，并选取2000年到2004年的数据进行分析。

技术细节写在这里：

我已经编写了模式的运行脚本，我也将COARE算法放入到了CAS-ESM中，现在只要在module_coare_model.f90中增加我的psiq的计算子程序再调用就可以了。

 

 real function psiq_26(zet)
!% computes specific humidity  structure function (ma,2022)real(r8) :: zetreal(r8) :: dzet, psi, x, psik, psic, fdzet = min(50.,0.35*zet) !% stablepsiq_26 = -0.6*zetif(zet < 0.) then !% unstablex = (1.-15.*zet)**0.5psik = 2.*log((1.+x)/2.)x = (1.-34.15*zet)**0.3333psic = 1.5*log((1.+x+x*x)/3.) - sqrt(3.)*atan((1.+2.*x)/sqrt(3.)) + &4.*atan(1.)/sqrt(3.)f = zet*zet/(1.+zet*zet)psiq_26 = (1.-f)*psik + f*psicendifend function psiq_26

在稳定的时候调用我们的比湿稳定度修正函数，计算潜热通量。下面是我的调试用的脚本，每次只要./tiaoshi.sh就可以自动进行模式运行试验了。我把这次试验叫HIST_coare_psiq_ri。

#case试验名称
CASE_NAME='HIST_coare_psiq_ri' 
#试验开始年月日YYYYMMDD
start_ymd='20000101'
stop_option='nyears'
stop_n='1'#判断和删除残留文件夹if [ -d "/data/chengxl/CAS-ESM2.0-test1/scripts/${CASE_NAME}" ]; thenrm -r /data/chengxl/CAS-ESM2.0-test1/scripts/${CASE_NAME}
fiif [ -d "/data/chengxl/CAS-ESM2.0-test1/run/${CASE_NAME}" ]; thenrm -r /data/chengxl/CAS-ESM2.0-test1/run/${CASE_NAME}
fi#创建case试验和编译
cd /data/chengxl/CAS-ESM2.0-test1/scripts./create_newcase -case ${CASE_NAME} -compset HIST_C6_C -res fd14_licom -mach huan_defaultcd /data/chengxl/CAS-ESM2.0-test1/scripts/${CASE_NAME}./configure -case./${CASE_NAME}.huan_default.build#运行试验前的设置cd /data/chengxl/CAS-ESM2.0-test1/run/${CASE_NAME}/run #cp /data/chengxl/CAS-ESM2.0-test1/run_demo/atm_in . 
#cp /data/chengxl/CAS-ESM2.0-test1/run_demo/drv_in . #输出变量管理
sed -i "12 a\\nhtfrq=0,0" atm_in
sed -i "13 a\\fincl1 = 'lhf_tbf','tau_tbf' ,'shf_tbf','n2','pbl_h','h_wave','tau_coare','hsb_coare','hlb_coare', 'LHFLX', 'SHFLX','qsss'" atm_in
sed -i "14 a\\fincl2 = 'lhf_tbf','tau_tbf','shf_tbf','tau_coare','hsb_coare','hlb_coare','h_wave' ,'h_wave_ln','pbl_h','ustar','obklen','n2','LHFLX','SHFLX','ri'"  atm_in#运行时间控制
sed -i "1,$ s/continue/startup/g" drv_in
sed -i "1,$ s/start_ymd      =  00010101/start_ymd      =${start_ymd}/g" drv_in 
sed -i "1,$ s/ndays/${stop_option}/g" drv_in
sed -i "1,$ s/stop_n         = 5/stop_n         =${stop_n}/g" drv_in
sed -i "1,$ s/restart_n      = 5/restart_n      =${stop_n}/g" drv_incp /data/chengxl/CAS-ESM2.0-test1/run_demo/run.slurm .
sbatch run.slurm
sleep 2
squeue
sleep 2

经过二十四小时左右的运行，模式结果全部都输出好了，没有错误。我们对模式的输出进行分析。

我将我的边界层模块的变量都放在了h1文件中。我们首先查看12月的都在。

(base) [chengxl@login01 run]$ ls B20TR_C35_y100.iap.h1.2000-*

B20TR_C35_y100.iap.h1.2000-01.nc B20TR_C35_y100.iap.h1.2000-05.nc B20TR_C35_y100.iap.h1.2000-09.nc

B20TR_C35_y100.iap.h1.2000-02.nc B20TR_C35_y100.iap.h1.2000-06.nc B20TR_C35_y100.iap.h1.2000-10.nc

B20TR_C35_y100.iap.h1.2000-03.nc B20TR_C35_y100.iap.h1.2000-07.nc B20TR_C35_y100.iap.h1.2000-11.nc

B20TR_C35_y100.iap.h1.2000-04.nc B20TR_C35_y100.iap.h1.2000-08.nc B20TR_C35_y100.iap.h1.2000-12.nc

 我们首先将其转化为一年平均的，因为是地球系统模式，所以分析年平均是有意义的，时间太短反而科学性有待考量。我们用一句语句来解决：

ls B20TR_C35_y100.iap.h1.2000-* | xargs -I{} cdo yearmean {} CASESM-COARE-psiq.2000ym.nc
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo).nc
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.27s 50MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.31s 49MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.39s 49MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.31s 50MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.30s 49MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.44s 49MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.37s 49MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.19s 49MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.32s 49MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.49s 49MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.31s 50MB].
 
cdo: /public/software/apps/anaconda3/5.3.0/lib/libuuid.so.1: no version information available (required by cdo)
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable date_written!
 
Warning (cdfCheckVars): Unsupported data type (char/string), skipped variable time_written!
 
cdo yearmean: Processed 2720139 values from 29 variables over 1 timestep [0.23s 49MB].
 

 

 成功得到了我们想要的数据CASESM-COARE-psiq.2000ym.nc
 

因为我比较习惯于使用本地主机上的Python还有panoply来分析我的结果，所以我先用ftp把我的资料从超算上下载到本地。
 

 

我们用panoply来简单看一下潜热的计算结果。
 

 

先保存一下，hlb_coare_in_CASESM-COARE-psiq.2000ym.png
 

然后我们再看下如果是原先用psit计算出来的潜热通量
 

 

我们看到计算出来的结果是合理的。比原先的潜热模拟值要大，这也是我们预料之内的。因为稳定度参数随着稳定度的增加减少量是变小的，这样的话分母变小，系数是相对变大的。
 

 

哦对了，因为极地区域没有经过处理，所以也是偏大的。这一点还存在一点问题。下次修改代码的时候把极地的通量采用原先的计算结果。就是在海冰的时候不该使用COARE算法。
 

除了这个问题，其他没有问题了。我们开始用Python来进行一些简单的数据处理。就是和OAflux潜热，以及Psiq,Psit,之间的对比
 


ctrl_data = xr.open_dataset(r'H:\fluxdeepl\CAS_ESM_psi_test\CASESM-COARE-psit.2000ym.nc', engine='netcdf4')exp1_data = xr.open_dataset(r'H:\fluxdeepl\CAS_ESM_psi_test\CASESM-COARE-psiq.2000ym.nc', engine='netcdf4')oaflux_data = xr.open_dataset(r'H:\fluxdeepl\CAS_ESM2_flux_test\latent256X128.nc', engine='netcdf4')lhf_ctrl = ctrl_data.hlb_coare.data
lhf_ctrl = np.squeeze(lhf_ctrl, axis = 0)lhf_exp1 = exp1_data.hlb_coare.data
lhf_exp1 = np.squeeze(lhf_exp1, axis = 0)lhf_oa   = oaflux_data.lhtfl.data[192:204,:,:]
lhf_oa = np.nanmean(lhf_oa,axis = 0)exp1_oa = lhf_exp1 - lhf_oa
ctrl_oa = lhf_ctrl - lhf_oa 
exp1_ctrl = lhf_exp1 - lhf_ctrl
write_to_nc_1(lhf_oa,'oa.nc')
write_to_nc_1(lhf_exp1,'psiq.nc')
write_to_nc_1(lhf_ctrl,'psit.nc')
write_to_nc_1(exp1_oa,'psiq_oa.nc')
write_to_nc_1(ctrl_oa,'psit_oa.nc')
write_to_nc_1(exp1_ctrl,'psiq_psit.nc')

 

 然后我们就用panpoly来画图了，懒得用NCL或者python画了。
 

 

我们可以看见： 

 

 

 
 
1.潜热通量对水汽稳定度修正函数响应较大。
 
2.热带海气边界层湍流潜热通量的模拟有所改进。
 
 












本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.hqwc.cn/news/344135.html


如若内容造成侵权/违法违规/事实不符，请联系编程知识网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！












相关文章










Unity游戏开发面试知识点全解读








Unity游戏开发面试知识点全解读




Unity游戏开发面试知识点全解读 
在数字化世界中&#xff0c;Unity游戏开发不仅是一种艺术形式和商业活动&#xff0c;而且已成为推动整个数字娱乐产业创新和进步的重要力量。Unity以其易用性、高效性和灵活性&#xff0c;赢得了全球开发者的青睐&#xff0c;从独立开发者到大型…




阅读更多...

















伐木工 - 华为OD统一考试








伐木工 - 华为OD统一考试




OD统一考试 题解： Java / Python / C++ 题目描述 
一根X米长的树木，伐木工切割成不同长度的木材后进行交易，交易价格为每根木头长度的乘积。规定切割后的每根木头长度都为正整数,也可以不切割，直接拿整根树木进行交易。请问伐木工如何尽量少的切割，才能使收益最大化? 
输…




阅读更多...

















mysql8.0 1055报错问题 [安装包]安装修改








mysql8.0 1055报错问题 [安装包]安装修改




mysql1055报错问题 this is incompatible with sql_modeonly_full_group_by 刚从5.7 升级到8.0 再运行之前的代码就给我搞事情 
例如 select a , b from table group by b 之前没问题 8.0默认有个only_full_group_by 这个是强效验证的  网上看了看资料 发现我的sql_Mode里面 没…




阅读更多...

















EOCR电机保护器485通讯协议概念








EOCR电机保护器485通讯协议概念




Modbus是由Modicon&#xff08;现为施耐德电气公司的一个品牌&#xff09;在1979年发明的&#xff0c;是全球第一个真正用于工业现场的总线协议。为更好地普及和推动Modbus在基于以太网上的分布式应用&#xff0c;目前施耐德公司已将Modbus协议的所有权移交给IDA&#xff08;In…




阅读更多...

















无法解析的外部符号 “public: virtual void * __cdecl MyTcpsocket::qt_metaca








无法解析的外部符号 “public: virtual void * __cdecl MyTcpsocket::qt_metaca




问题&#xff1a;严重性    代码    说明    项目    文件    行    禁止显示状态 错误    LNK2001    无法解析的外部符号 "public: virtual void * __cdecl MyTcpsocket::qt_metacast(char const *)" (?qt_metacastMyTcpsocketUEAAPEAXPEBDZ)    SmartTool    D:\…




阅读更多...

















使用 CMake 和 Ninja 构建 C/C++ 项目的教程








使用 CMake 和 Ninja 构建 C/C++ 项目的教程




使用 CMake 和 Ninja 构建 C/C 项目的教程 
CMake 是一个跨平台的开源构建工具&#xff0c;它简化了项目的构建过程。而 Ninja 是一个快速、轻量级的构建系统&#xff0c;与 CMake 配合使用可以提高项目的构建效率。本教程将向你介绍如何使用 CMake 和 Ninja 来构建你的 C/C 项…




阅读更多...

















matlab/simulink风电光伏储能虚拟同步机VSG下垂控制虚拟惯量控制，光伏变压减载控制一次调频二次调频研究，储能下垂控制SOC








matlab/simulink风电光伏储能虚拟同步机VSG下垂控制虚拟惯量控制，光伏变压减载控制一次调频二次调频研究，储能下垂控制SOC




非无穷大系统&#xff01; 
第一个介绍风光储VSG算例研究。 
其中风机采用虚拟惯性控制&#xff0c;下垂控制&#xff0c;网侧采用VSG虚拟同步机控制。储能下垂控制。 上述结果为整个系统的仿真波形&#xff0c;包括频率&#xff0c;pmsg永磁风机出力&#xff0c;储能出力&…




阅读更多...

















无重复字符的最长字串








无重复字符的最长字串




题目 给定一个字符串 s &#xff0c;请你找出其中不含有重复字符的 最长子串 的长度。 示例 
示例 1: 
输入: s  "abcabcbb"
输出: 3 
解释: 因为无重复字符的最长子串是 "abc"&#xff0c;所以其长度为 3。示例 2: 
输入: s  "bbbbb"
输出: 1
解…




阅读更多...

















揭露欧拉骗局4.“Σ1/n²=π²/6”里的猫腻








揭露欧拉骗局4.“Σ1/n²=π²/6”里的猫腻




自然数平方倒数求和Σ1/n是一个并不复杂的问题&#xff0c;但它困扰了欧洲大陆整整90年&#xff0c;在欧系数学里它被称为“巴塞尔级数”。 
解决巴塞尔级数让欧拉一战成名&#xff0c;然而欧拉采用的方法对数学这门学问是严重的侮辱。数学是工具学科&#xff0c;数学的宗旨是化…




阅读更多...

















聚焦老年生活与健康，“老有所依·情暖夕阳”元岗街社区微型养老博览会顺利开展








聚焦老年生活与健康，“老有所依·情暖夕阳”元岗街社区微型养老博览会顺利开展




尊老敬老是中华民族的传统美德&#xff0c; 
爱老助老是全社会的共同责任。 
家有一老&#xff0c;如有一宝&#xff0c; 
长者的生活情况是一个家庭的头等大事&#xff0c; 
做好长者服务是街道和社区的重要工作。 2024年1月6日&#xff0c;由元岗街道党工委、元岗街道办事处、…




阅读更多...

















挑选知识付费平台，只需掌握这几招








挑选知识付费平台，只需掌握这几招




明理信息科技知识付费saas租户平台  
随着知识经济的兴起&#xff0c;越来越多的人开始重视知识付费&#xff0c;并希望通过打造自己的知识付费平台来实现自己的知识变现。本文将介绍如何打造自己的知识付费平台&#xff0c;并从定位、内容制作、渠道推广、运营维护四个方面进行…




阅读更多...

















【PyQt小知识 - 8】：QFileDialog — 文件选择对话框（选择文件夹/文件）








【PyQt小知识 - 8】：QFileDialog — 文件选择对话框（选择文件夹/文件）




文章目录 QFileDialog的介绍和使用 QFileDialog的介绍和使用 
QFileDialog是Qt框架中提供的一个文件选择对话框。它可以让用户通过图形界面的方式选择文件或目录&#xff0c;并提供了许多与文件操作相关的实用功能&#xff0c;例如文件过滤、文件类型限制、文件名默认值、默认目…




阅读更多...


















推荐文章












最新文章


















Copyright @ 2022~2023