!--------------------------------------------------------------------- ! Copyright (C) GFD Dennou Club, 2005. All rights reserved. !--------------------------------------------------------------------- !=begin != Subroutine DissipativeHeating ! ! * Developer: KITAMORI Taichi ! * Version: $Id: dissipativeheating.f90,v 1.1.2.2 2005/10/11 04:57:06 kitamo Exp $ ! * Tag Name: $Name: arare3m-20051114 $ ! * Change History: ! !== Overview ! !»¶°ï²ÃÇ®¤ò·×»»¤¹¤ë. ! !== Error Handling ! !== Known Bugs ! !== Note ! ! ! !== Future Plans ! !=end subroutine DissipativeHeating(ss_Km_bl, ss_QDis) !=begin !==Dependency use dc_trace, only: BeginSub, EndSub use gridset, only: DimXMin, DimXMax, DimZMin, DimZMax, RegZMin, RegZMax, DelXZ use basicset, only: ss_CpBasicZ, ss_PotTempBasicZ, ss_ExnerBasicZ use arareset, only: Cm !=end !==°ÅÌۤη¿Àë¸À¤ò¶Ø»ß implicit none !=begin !==Input real(8), intent(in) :: ss_Km_bl(DimXMin:DimXMax, DimZMin:DimZMax) ! ±¿Æ°Î̤ËÂФ¹¤ëÍðή³È»¶·¸¿ô !==Output real(8), intent(out) :: ss_QDis(DimXMin:DimXMax, DimZMin:DimZMax) ! »¶°ï²ÃÇ®¥¨¥Í¥ë¥®¡¼ !=end call BeginSub("DissipativeHeating", & & fmt="%c", & & c1="Calculate dissipative heating (ss_QDis)." ) ss_QDis = ss_Km_bl ** 3.0d0 / ( ss_ExnerBasicZ * ss_CpBasicZ * Cm**2.0d0 * DelXZ**4.0d0 ) call EndSub("DissipativeHeating") end subroutine DissipativeHeating