<div>Dear Colleges,</div><div> </div><div>My aim is to plot the WRF simulated vertical cross-sections as shown in the attached example which I downloaded from <a href="https://gradsaddict.blogspot.com/2013/05/tutorial-vertical-cross-section-with.html">https://gradsaddict.blogspot.com/2013/05/tutorial-vertical-cross-section-with.html</a>. However, they used pressure levels instead of model (sigma) level heights. In my ctl file I have 64 fixed model level heights relative to ground surface, also I have the model topography field (HGT=>h 1 t,y,x  compare). I am interested in 39 S (set lat -39) west-east vertical cross-section for the QRAIN variable (QRAIN=>qr 64 t,z,y,x  compare). The question is if it is possible to reset my fixed 64 vertical  levels and plot the cross-section considering terrain height (HGT variable) after opening the ctl file and plot absolute heights (levs + HGT) in my vertical cross-section, then shade terrain height as it was presented in the attached Figure?</div><div> </div><div>Kind regards,</div><div>Artur</div><div> </div>