如何在WRF中增加物理方案中的输出

不吃鱼的小虾米

想在wrfout中增加mp_physics中诊断量的输出，比如Thompson方案中QVAPOR的source和sink项
请问这个具体步骤是什么呢？如果有相关贴也可以~在论坛中没有找到

youxia1873

修改register和thompson方案的代码，要是搞定了，教教我可以不

xinshou6

state   real    qv     ikjftb   moist     1       -      i0rh01usdf=(bdy_interp:dt) "QVAPOR"         "Water vapor mixing ratio"      "kg kg-1"
microphysics_driver ,SELECT CASE(mp_physics),CASE (THOMPSON),CALL mp_gt_driver,搜thompson，module_mp_thompson ,SUBROUTINE mp_gt_driver(qv, qc, qr, ，This is a wrapper routine designed to transfer values from 3D to 1D，REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT):: qv, qc, 水汽混合比是湿空气内水汽与干空气的质量比。Source/sink terms.  First 2 chars: "pr" represents source/sink of mass while "pn" represents source/sink of number.，"v" for water vapor，"de" for sublimation/deposition, "ev" for evaporation, "fz" for freezing, "ml" for melting, "au" for autoconversion, "nu" for ice nucleation, "hm" for Hallet/Mossop secondary ice production, and "c" for collection，就看到prv_rev(k)，DOUBLE PRECISION, DIMENSION(kts:kte)，registry没有所以得加state,用real行不，公式是啥 啥是源汇项 Calculate tendencies of all species but constrain the number of ice to reasonable values，Water vapor tendency是这个吗 qvten(k) = qvten(k) + (-pri_inu(k) - pri_iha(k) - pri_ide(k)  - prs_ide(k) - prs_sde(k) - prg_gde(k))  * orho,pri_inu是啥意思 cloud ice nucleation产生的cloud ice 质量吗，后面还有个qv(k) = MAX(1.E-10, qv1d(k) + DT*qvten(k))，要是输出是不前面那些也得*dt啊