micapsphysics计算.docx

物理量的计算：（《天气分析预报物理量计算基础》） 干空气和水汽的比气体常数分别为Rd=287.05J/kg/K，Rv=461.5J/kg/k，其中K为开尔文温度单位； 一． 利用FNL资料计算以下物理量： （1）Richardson数 （1）修正的K指数： 地面至850hPa温度与露点平均之和减700hPa温度露点差再减500hPa温度。 mK= (Tsur+T1000+T925+T850)avg-(Td700+T500) （2）最大对流稳定度指数(?BIC ) 在边界层附近，往往有一个的极大值()，在对流层中层往往有一个的极小值()，所谓用以下表达式定义： 计算方法： 从地面开始到对流层顶，分别寻找出最大假相当位温和最小位温，求出二者之差。 （3）深对流指数 850hPa温度露点和与抬升指数LI相减之差。 诊断用深对流指数： 深对流指伸展高度具有等于或大于均质大气高度H0（与400hPa等压面高度更接近）的对流系统。利用云顶相当黑体亮温计算的深对流指数可以作为表示云顶等于或高于400hPa深对流云的指数。 预报用深对流指数DCI DCI=T850+Td850-LI （4）抬升指数LI? 平均气块从修正的低层3000英尺高度沿干绝热线上升，到达凝结高度后再沿湿绝热线上升至500hPa时所具有的温度与500hPa等压面上的环境温度(T500)的差值。 （5）风暴强度指数(P176) （6）密度加权平均风垂直切变shr(P118) (7)对流稳定度指数（P84） 假相当位温（P37） （8）全总指数 （9）强天气威胁指数（P216） （10）条件-对流稳定度指数（P84） 为饱和假相当位温，具体见P37 （11）对流有效位能（CAPE） （12）对流抑制有效位能（CIN） 水汽压 当t0时(水面)，； 当t=0时(冰面)，； （2）饱和水汽压 将上述公式中的t改为露点温度，即得饱和水汽压公式 比湿 p：气压值，Es:水汽压 饱和比湿 把上式水汽压改为饱和水汽压，即为饱和比湿。 混合比（水汽质量与干空气质量的比值） 相对湿度 实际水汽压/饱和水汽压 ６）位温  其中： ????????????? T：为温度 ????????????? P：为气压 ????????????? Rd/Cpd：一般为常数，取值0.268 ? ７）假相当位温 其中： ????????????? L：为凝结潜热，一般为常数 ????????????? Cp：为定压比热，为常数 ????????????? q：为比湿 ?(define t=TMPprs define rh=RHprs define prs=lev define es=(6.112*exp(17.67*(t-273.15)/(t-29.65))) define q=rh*(0.62197*es/(prs-es))/100. define e=prs*q/(0.62197+q)+1e-10 define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805) define theta=t*pow((1000/prs),(0.2854*(1.0-0.28*q))) define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q)) ) Or(* Script to calculate psuedoadiabatic equivalent potential temperature * (the) from cm1 GrADS-format output * Author: George H. Bryan, NCAR/MMM * Last modified: 8 August 2008 th is potential temperature;qv is water vapor mixing ratio define e=0.01*prs*qv/(0.6219718+qv)+1e-20 define t=th*pow((prs/100000.0),(287.04/1005.7)) define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805) define the=t*pow(100000.0/prs,0.2854*(1.0-0.28*qv))*exp(((3376.0/tlcl)-2.54)*qv*(1.0+0.81*qv))) ８）饱和假相当位温 将前式的温度改为露点温度即为饱和假相当位温 ? ９）CAPE 公式： 计算方法： ??? 第一步： 首先假设气块从地面开始按干绝热抬升，干绝热抬升即处处位温相等，将前面位温公式两边取对数得： ??? 逐步抬升气压（可按5m