天天编程je - 博客园

wrfout 变量提取

摘要： load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/ 阅读全文

posted @ 2024-01-01 22:40 天天编程je 阅读(139) 评论(0) 推荐(0)

WRF------------单层水汽通量

摘要： `; WRFOUT D01 a = addfile(".nc","r") u = wrf_user_getvar(a,"ua",-1) v = wrf_user_getvar(a,"va",-1) p = wrf_user_getvar(a,"pressure",-1) ; z = wrf_user 阅读全文

posted @ 2023-12-20 21:07 天天编程je 阅读(159) 评论(0) 推荐(0)

如何求wrfout日降水量？

摘要： wrfout日降水量如何处理？rainnc+rainc是累计降水利用ncl如何循环读取每日降水量 ; 计算24h累计降水，因为逐1h输出一次，所以3-7月一共3672个文件 do i=0,ntimes-1 rainc(i,:,:) = wrf_user_getvar(f, "RAINC", i + 阅读全文

posted @ 2023-06-14 17:53 天天编程je 阅读(337) 评论(2) 推荐(0)

win10出现0x80071ac3无法完成操作因为卷有问题怎么解决?

摘要： 1.win+r 2.输入 cmd 3.输入 chkdsk 你的符盘:/f 阅读全文

posted @ 2023-04-23 21:32 天天编程je 阅读(182) 评论(0) 推荐(0)

WRFdomain选取

摘要：可以使用该网站WRF Domain Wizard (jiririchter.github.io) 直接选定模拟区域，相关参数会出来，点击open可以保存该方案至本地阅读全文

posted @ 2023-04-23 21:20 天天编程je 阅读(181) 评论(0) 推荐(0)

ncl 常用绘图属性之常规设置

摘要： wks_type = "png" wks_type@wkWidth = 2500 ；设置分辨率 wks_type@wkHeight = 2500 res = True res@tmXTOn = False res@tmYROn = False ;;;;;;;;;;;;;;;;;;;;;;;;;;;; 阅读全文

posted @ 2022-11-09 21:51 天天编程je 阅读(1077) 评论(0) 推荐(0)

ncl中使用mask 但最后mask全是缺测解决方法

摘要：；当我在使用mask提取我所需要的区域时，发现提取出来的全是缺测值；主要问题出现在：；变量名var，维度属性为 time lon lat ；当我调整维度为： time lat lon 便解决问题 printVarsummaty( var ) ；通过该函数查看变量维度属性 var1= var({t 阅读全文

posted @ 2022-11-09 10:09 天天编程je 阅读(501) 评论(0) 推荐(0)

ncl批量读取多个nc文件

摘要： diri = "/mnt/g/data/3/" files = systemfunc("ls " + diri + "*nc4") ;批量读取文件路径 f1 = addfiles(files,"r") ListSetType (f1, "cat") lat = f1[0]->lat({20:30}) 阅读全文

posted @ 2022-11-07 21:08 天天编程je 阅读(2123) 评论(0) 推荐(0)

ncl中截取所需要的时间段数据

摘要： ;;;;;;;;;;;装换时间 t = f->time YYYY =cd_calendar(t, -1) yrstrt =201901 ；起始年、月 yrlast =201912 ；结束年、月 ist =ind(yrstrt.eq.YYYY) iet =ind(yrlast.eq.YYYY) pre 阅读全文

posted @ 2022-10-21 21:34 天天编程je 阅读(663) 评论(0) 推荐(0)

在某路径下建文件夹

摘要： os.chdir('/mnt/g/DA/1') if not os.path.exists('day'): os.mkdir('day') if not os.path.exists('month'): os.mkdir('month') if not os.path.exists('year'): 阅读全文

posted @ 2022-10-13 21:44 天天编程je 阅读(32) 评论(0) 推荐(0)

PREjelly

公告