分享10个时间-高度剖面图及gs（温度平流、涡度平流、水汽通量散度、垂直速度......）

飞飞

这个帖子太实用了

筱筱筱筱筱雅痴

请问一下楼主，为什么我水汽通量输出的图风向杆跟时间刻度不对应，错开来了。附上我的程序和出图
'reinit'
'sdfopen J:\gz\2019ktsb\grads\data\uwnd.2015.nc'
'sdfopen J:\gz\2019ktsb\grads\data\vwnd.2015.nc'
'sdfopen J:\gz\2019ktsb\grads\data\shum.2015.nc'
'sdfopen J:\gz\2019ktsb\grads\data\air.2015.nc'
'sdfopen J:\gz\2019ktsb\grads\data\rhum.2015.nc'

'set t 179 191'
'set lev 1000 500'
'define p=lev'
'define e=exp(17.67*(air.4-273.15)/(air.4-29.65))*(rhum.5/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=uwnd.1*q/9.8'
'define q2=vwnd.2*q/9.8'

'set grads off'
'set grid off'
'set csmooth on'
'set parea 1.5 10 1.5 8'
'set xlopts 1 4 0.25'
'set ylopts 1 4 0.25'

'set lev 1000 500'
'set lon 120.2'
'set lat 33.43'
'set zlog on'
'set ylevs 1000 975 950 925 900 850 800 750 700 650 600 550 500'
'set xlabs 20|02|08|14|20|02|08|14|20|02|08|14|20|02'
'set gxout contour'
'set clopts 0'
'set gxout shaded'
'set cint 0.6'
'd mag(q1,q2)'
'set gxout barb'
'set ccolor 2'
'd uwnd.1*2.5;vwnd.2*2.5'
'run cbarn.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 0.6'
'set cthick 1'
'set clskip 1'
'd mag(q1,q2)'

'writehz 0.7 0.8 2月14日 1 2 1 1.3 0.6 0 7'
'writehz 1.8 0.8 2月15日 1 2 1 1.3 0.6 0 7'
'writehz 4.3 0.8 2月16日 1 2 1 1.3 0.6 0 7'
'writehz 7.0 0.8 2月17日 1 2 1 1.3 0.6 0 7'
'printim J:\gz\2019ktsb\grads\pic\z-sqtl.png white'
C:\Users\Administrator\Desktop

qiuba

你检查一下时间设置起点“set t 179“是对应20时吗？

linyu临宇

楼主你好，使用你的程序绘画涡度平流提示错误是什么问题呢？

Adiligirl

感谢楼主分享，mark一下

qiuba

linyu临宇 发表于 2020-8-14 22:21
楼主你好，使用你的程序绘画涡度平流提示错误是什么问题呢？

检查一下你的原始资料顺序是否和我的有变化，应该完全按我的资料顺序列举，后续的算法中涉及到了各原始资料的位置（顺序）

linyu临宇

qiuba 发表于 2020-8-16 07:05
检查一下你的原始资料顺序是否和我的有变化，应该完全按我的资料顺序列举，后续的算法中涉及到了各原始资 ...

我用fnl的资料，变量名都改过了，nc资料太大了老是下载不下来

qiuba

linyu临宇 发表于 2020-8-18 15:42
我用fnl的资料，变量名都改过了，nc资料太大了老是下载不下来

没关系的，只要变量名和资料里的变量一致，且计算式中引用位置正确就不会出错，反之则很容易出错。建议你在家园里可下载与fnI资料对应的计算公式计算

linyu临宇

qiuba 发表于 2020-8-18 22:51
没关系的，只要变量名和资料里的变量一致，且计算式中引用位置正确就不会出错，反之则很容易出错。建议你 ...

是的，我就是用fnl资料，计算式中引用的位置也正确，画涡度平流和假相当位温温度平流，提示图片上显示问题，没有看懂是什么意识

qiuba

linyu临宇 发表于 2020-8-19 11:54
是的，我就是用fnl资料，计算式中引用的位置也正确，画涡度平流和假相当位温温度平流，提示图片上显示问 ...

我用的nc资料是2.5x2.5的，如果用1x1的fnI资料，变量的表达式就变了，你检一下对应的表达式书写是否正确，引用位置是否正确