2 subroutine da_radzicevar_parameter_zxx(p1,para_pdx_dq,para_pdx_df,rhoa,rhox,qx,cx,dx,alphax, &
3 mm3todBZ,lambda,Kw2,pi,n0x)
4 !----------------------------------------------------------------------
5 ! purpose: calculate the constants for ice speices
6 !----------------------------------------------------------------------
10 real :: cx,dx,alphax,n0x ! for two moment, not used now
11 real :: mm3todBZ,lambda,Kw2,pi
13 real :: para_pdx_df,para_pdx_dq
16 p1=mm3todBZ*4.*gamma(7.)*lambda**4/(pi**4*Kw2)*(pi*rhox)**(-1.75)
17 para_pdx_dq=1.75*p1*n0x**(-0.75)*(rhoa*qx)**(0.75)
18 para_pdx_df=p1*n0x**(-0.75)*(rhoa*qx)**(1.75)
20 end subroutine da_radzicevar_parameter_zxx