File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -244,6 +244,7 @@ void WindNinja::init(mesh& domain)
244244 H_forc = cfg.get (" H_forc" ,40.0 );
245245 Max_spdup = cfg.get (" Max_spdup" ,3 .);
246246 Min_spdup = cfg.get (" Min_spdup" ,0.1 );
247+ scale_factor = cfg.get (" scale_factor" ,1.0 );
247248 ninja_recirc = cfg.get (" ninja_recirc" ,false );
248249 Sx_crit = cfg.get (" Sx_crit" , 30 .);
249250
@@ -465,8 +466,8 @@ void WindNinja::run(mesh& domain)
465466 Atmosphere::Z_U_R , // UR is at our reference height
466467 0 ); // no canopy, no snow, but uses a snow roughness
467468
468- (*face)[ " U_R " _s]= W;
469-
469+ // scale_factor can be used to correct any low bias. Heavy handed approach
470+ (*face)[ " U_R " _s]= W * scale_factor;
470471
471472 (*face)[" zonal_u" _s]= U; // these are still H_forc
472473 (*face)[" zonal_v" _s]= V;
Original file line number Diff line number Diff line change @@ -185,6 +185,7 @@ REGISTER_MODULE_HPP(WindNinja);
185185 double H_forc; // Reference height for GEM forcing and WindNinja wind field library
186186 double Max_spdup; // Maximal value of crest speedup
187187 double Min_spdup; // Minimal value of crest speedup
188+ double scale_factor; // used to scale up the windspeeds in the case of a serious bias
188189 int L_avg; // Size of the square of averaging when computing the speed-up map
189190 // Not used by default and set to -1 if not specified in the config file.
190191 bool ninja_recirc; // Boolean to activate wind speed reduction on the leeside of mountainous terrain
You can’t perform that action at this time.
0 commit comments