请问如何关闭wrf的感热通量？

月贝凡

请问

做热力敏感试验如何关闭wrf初始场的感热通量？wrfinput里面没有这个量，改变相关的量对其他的计算会有影响的，那如何实现呢？

白雪红梅

请问楼主最后解决了没？

ztftom

你是说sensible heat flux么？ 这个量是计算出来的，不是输入的，所以无论wrfinput还是wrfbdy你都不会找到这个量的。你可以考虑isfflx=0， 但这样的话humidity flux也会为0

月贝凡

ztftom 发表于 2012-4-3 13:20

                               
登录/注册后可看大图

你是说sensible heat flux么？ 这个量是计算出来的，不是输入的，所以无论wrfinput还是wrfbdy你都不会找到这 ...

谢谢你的回复，isfflx=0会使水汽通量为0，但只想关闭感热通量，也就是sensible heat flux，琢磨好久就是没找到起作用的那部分，包括phys下面的一些程序，请问版主有没有好的建议给我？

wangmol

只关掉感热通量貌似比较困难吧，通量是从近地层方案里面计算出来的，如果所有通量都不考虑，您可以把所有与边界层相关的选项都关掉试试。sf_sfclay_physics=0，sf_surface_physics=0，bl_pbl_physics=0，不考虑近地层、路面过程和边界层闭合。
如果真想只让感热为0，或许动sf_sfclay_physics对应的近地层程序可以,把HFX赋值0（比如module_sf_sfclay.f sf_sfclay_physics=1)，重新编译一下。
我说的不一定对，您可以试试。

月贝凡

wangmol 发表于 2012-4-4 10:08

                               
登录/注册后可看大图

只关掉感热通量貌似比较困难吧，通量是从近地层方案里面计算出来的，如果所有通量都不考虑，您可以把所有与 ...

很感谢您的回复，按道理的确应该是修改sf_sfclay_physics对应的近地层程序，我几乎修改了m0dule_sf_*等一系列的程序，重新编译后仍然对降水没有任何影响，wrfout里面的HFX也基本没有变化，头大了

wangmol

我只仔细看过module_sf_sfclay.f 感热通量是由Ｈ ＝－ρCｐU＊θ＊最后计算出来的，这里面应该有迭代，但是貌似不在这个子程序里面，你可以看到有这么一段：
!-----COMPUTE SURFACE HEAT FLUX:                                                
!                                                                                
  390 CONTINUE                                                                 
      DO 400 I=its,ite
        IF(XLAND(I)-1.5.GT.0.)THEN                                          
          HFX(I)=FLHC(I)*(THGB(I)-THX(I))                                
          IF ( PRESENT(ISFTCFLX) ) THEN
             IF ( ISFTCFLX.NE.0 ) THEN
! AHW: add dissipative heating term
                HFX(I)=HFX(I)+RHOX(I)*USTM(I)*USTM(I)*WSPDI(I)
             ENDIF
          ENDIF
        ELSEIF(XLAND(I)-1.5.LT.0.)THEN                                       
          HFX(I)=FLHC(I)*(THGB(I)-THX(I))                                
          HFX(I)=AMAX1(HFX(I),-250.)
        ENDIF
  400 CONTINUE
这里面也讨论了是陆地还是水面。FLHC在前面已经算好了
FLHC(I)=CPM(I)*RHOX(I)*UST(I)*MOL(I)/(THX(I)-THGB(I))
CPM应该是在当前湿度条件下的定压比容，RHOX空气密度，UST摩擦速度，MOL温度尺度，THX应该为参考高度的位温，THGB地面位温。XLAND=1陆面XLAND=2水面。
所以你要改感热HFX的话，基本就在这里改就行了，至于能不能对降水产生影响，我也不清楚了。

月贝凡

wangmol 发表于 2012-4-5 09:49

                               
登录/注册后可看大图

我只仔细看过module_sf_sfclay.f 感热通量是由Ｈ ＝－ρCｐU＊θ＊最后计算出来的，这里面应该有迭代，但是 ...

非常感谢您的耐心回复，前面也是修改的这部分，对降水没有产生影响，现在在尝试修改其他计算程序中与感热相关的部分，基本上现在就是挨个试了，呵呵。

wangmol

不用谢，很高兴大家互相交流学习。

lgg00828

月贝凡 发表于 2012-4-5 11:17

                               
登录/注册后可看大图

非常感谢您的耐心回复，前面也是修改的这部分，对降水没有产生影响，现在在尝试修改其他计算程序中与感热 ...

请问大神做出来了吗？

小小reed

呃，好复杂的赶脚，看来楼主对wrf运用到一定程度了。。。我想请教楼主，初始场是FNL文件包，我想把初始场里的风切变去掉~该如何是好呢？？