search

daliagachc / wrf_management

package created to manage the run of wrf at uhel supercomputer (taito)
0 stars 0 forks source link

PBL height #10

Open Rulock opened 5 years ago

Rulock commented 5 years ago

I've made some gif animations for time evolution of PBL height (every hour for 4 days starting on 5th day of the run): https://imgur.com/a/7lJiExL

d02 domain - shallow stable PBL over colder-than-air lakes during daytime and convective PBL over Titicaca Lake at night as we wanted it to be: ezgif-2-8675db167365

I did it in Matlab since I had all the scripts ready. Should it rewrite it in Python? Any advices how to improve the plots? Not sure if isolines of terrain height illustrate topography well.

Rulock commented 5 years ago

In our run we use Mellor–Yamada–Janjic (MYJ) PBL scheme which is a 1.5-order prognostic turbulent kinetic energy (TKE) scheme with local formulation to determine the eddy diffusion coefficients and the vertical mixing. The PBL height is determined from the value of TKE as the level at which the TKE profile decreases to a small prescribed value ( 0.2 [m/s]^2 ).

Rulock commented 5 years ago

Same animations but with histograms:

PBLH_d01_histogram_opt PBLH_d02_histogram_opt PBLH_d03_histogram_opt PBLH_d04_histogram_opt

Rulock commented 5 years ago

More draft animations: d01 and inner domains: d01domain-min d04 domain (Lake mask looks ugly, I will try to make it prettier): d04domain-min d04 T2m plots over 3 days (same problem with Lake mask): d04T2-min d04 10m wind (still some work to be done): d04wind-min I aslo will do temperature and wind vertical cross-sections along La Paz - Chacaltaya station line. Any ideas for more plots? PBL height?

Rulock commented 5 years ago

Cross sections: d04domaincross-min Potential temperature: d04Thetacross-min TKE and PBL height: d04TKEcross-min Winds to be added.