I want to calculate the surface exchange coefficient for heat with the WRF output variables. I think I found the file- associated with chosen surface scheme (MYJSFC, sf_sfclay_physics=2, sf_surface_physics=2) -module_sf_myjsfc.F, where the subroutine SFCDIF calculates AKHS. The only variables that I cannot identify are SIMH components like ZSL,PSIH1?PSHZL?PSHZ? Any ideas?