Hello,
I read on the forum about spin-up time that WRF is quite rapid to stabilize : for mesoscale, 12h is commonly advised. My case runs nested domains (7km and 2.333 km) at time-step for d01 at 40 seconds (on climate modeling time-scale).
I've been advised to look at TKE in the first (few) month(s) of my simulation to determine ideal spin-up time. I wonder however how to proceed using U10 and V10 (and other Us and Vs at lower pressure levels for example).
My research weren't successful on that matter. I looked at Skamarock et al., 2004 (mentionned by Kelly somewhere on the forum) to plot log-log TKE vs. wavelength but it seems to focus on free troposphere and above ; hence I'm not quite sure that I can use it for surface winds.
I've been told to compute TKE as (u'² + v'²), with u' = u - u_mean(time) over a given period (respectively for v') ; then summing/avering TKE over all domain. I'm also told to expect increasing TKE over time at the beginning of the run until it stabilizes, showing that model's spinned-up : but that is not what I observe with the quick script I wrote to do so.
Is there anyone here who'd be nice enough to help and provide some hints or methodological references ? Do I miss something ?
Thanks for your help
I read on the forum about spin-up time that WRF is quite rapid to stabilize : for mesoscale, 12h is commonly advised. My case runs nested domains (7km and 2.333 km) at time-step for d01 at 40 seconds (on climate modeling time-scale).
I've been advised to look at TKE in the first (few) month(s) of my simulation to determine ideal spin-up time. I wonder however how to proceed using U10 and V10 (and other Us and Vs at lower pressure levels for example).
My research weren't successful on that matter. I looked at Skamarock et al., 2004 (mentionned by Kelly somewhere on the forum) to plot log-log TKE vs. wavelength but it seems to focus on free troposphere and above ; hence I'm not quite sure that I can use it for surface winds.
I've been told to compute TKE as (u'² + v'²), with u' = u - u_mean(time) over a given period (respectively for v') ; then summing/avering TKE over all domain. I'm also told to expect increasing TKE over time at the beginning of the run until it stabilizes, showing that model's spinned-up : but that is not what I observe with the quick script I wrote to do so.
Is there anyone here who'd be nice enough to help and provide some hints or methodological references ? Do I miss something ?
Thanks for your help
