1. Canadian Network for Regional Climate and Weather Processes Practical Training: Plotting Model outputs using Matlab: Time series, vertical profiles, 2D plot Contact: oumarou.nikiema@uqam.ca nikiema@sca.uqam.ca Bureau: PK 6425

2. Model output Configuration • In practice, model data are in 5D format: • Variable of interest: temperature, Wind, … • Valid Time • Vertical position: pressure (hPa), height (m), … • Latitude position • Longitude position • Example: TT(2005010103,500hPa,30°N,60°W)

3. Login to the Computer … Access to Skynet3: ssh -Y wh01@skynet3.sca.uqam.ca Copy the directory in your home: cp –r /home/nikiema/Workshop_CNRCWP ./ Enter into the work directory: cdWorkshop_CNRCWP Check the list of files in your directory: ls –al Open Matlab: matlab&

4. Data files used by Matlab programs • In Directory: /skynet3_exec2/nikiema/Data_Workshop • 3 files: (Use VOIR to check these informations) • Model output: dp_... • 2005-01-01-03 to 2005-02-01-00 each 3h • Var.: TT , UU, VV and GZ • 26 pressure levels • Era-Interim data: era-… • 2005-01-01-03 to 2005-02-01-00 each 6h • Var.: TT , UU, VV and GZ • 34 pressure levels • LoLa: • latitude and longitude information

5. Data files used by Matlab programs voir -iment « Filename » -style | more Example: IP1: Pressure code CMC time code IP2: forecast hours ? For instance: IP2=1 for one hour prevision

6. Data files used by Matlab programs voir -iment « NOM DE FICHIER » –style Exemple: IP2: 30*24+3=723 Initial date: 2004-11-01-00h GZ data on: 2004-12-01-03h

7. Conversion: RPN to binary format Getfld6.m Matlab function used to convert RPN to binary fromat Getfld6.m: works only with E32 format It is possible to change data format by using: r.diag repack … Sometimes, you need to uncompress the data by using: fstuncompress -fstin « input file » -fstout « output file »

8. Conversion: RPN to binary format Syntax: Getfld6.m: Matlab function used to convert RPN to binary fromat

9. Conversion: RPN to binary format • TYPVAR: type de variable; 'A'(analysis),'P'(forecast),'X' • NOMVAR: variable name: TT, UU, VV, GZ, … • ETIKET: name of simulation • DATE: CMC date code • r.date -Sn 2005062600 (Code de la date: 324252800) • 1 janv. 1900 0h ===> 010100000 • 1 janv. 1900 1h ===> 010100010 • IP1: level code • IP2: prevision time • IP3: 0

10. Useful Fortran programs … r.date –Sn « DATE » options : -n to add end of line character, -S to get CMC time stamp (time code), -L to get date in this form YYYY MM DD HH MM SS -V to get in this form YYYYMMDDHHMMSS00 Example: Result: Result:

11. Useful Fortran programs … r.ip1 -kn « PRESSION » 2 options : -n to add end of line character, -k to get height code knowing the level in hPa, m, … -o to get the level value (in hPa, m, …) knowing code NB: k =====> kind kind : level type (the number express the level unity) ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 0 : m [metres] (height with respect to sea level) 1 : sg [sigma] (0.0 ----> 1.0) 2 : mb [mbars] (pressure in millibars) 3 : [others] (arbitrary code) 4 : M [metres] (height with respect to ground level) 5 : hy [hybrid] (0.0 ----> 1.0) 6 : th [theta]

12. Useful Fortran programs … r.ip1 -kn « PRESSION » 2 options : -n to add end of line character, -k to get height code knowing the level in hPa, m, … -o to get the level value (in hPa, m, …) knowing code Examples: Result: Result:

13. 1D plot of the Vertical profile Script: Vertical_Profile_1D.m

14. 1D plot of the time series Script: Time_series_1D.m

15. 2D plot: Latitude-height cross section Script: Vertical_Profile_2D_Along_long_v1.m

16. 2D plot: Latitude-height cross section Sometimes, it is better to plot the cross section With Log of the pressure. Script: Vertical_Profile_2D_Along_long_v2.m

17. 2D plot: Vertical profile of zonal-mean Temperature Script: Vertical_Profile_2D_Zavg.m

18. 2D plot: Longitude-height cross section Script: Vertical_Profile_2D_Along_lat_v2.m

19. 2D plot: horizontal map of temperature

20. 2D plot: horizontal map of temperature

21. quiverm ==> to seewindvectors

22. There are too much vectors, so you have to reduce the numbers of vectors

23. Ouvrir et exécuter: Visu_vent2D_2.m

24. Comparison: Model vs Era-Interim 1D vertical profile of the domaine-averaged

25. Comparison: Model vs Era-Interim Time series plot of the domaine-averaged

26. Comparison: Model vs Era-Interim