CRCM5 - UQAM
CRCM5 Katja Winger, UQAM with contributions from Bernard Dugas, RPN, and Anne Frigon, Ouranos Ouranos, 14. Nov. 2008
Overview ●
What is CRCM5? ● GEM ● GEMCLIM ● GCM
●
Current and future projects for CRCM5 / GEMCLIM
●
Recent runs done with CRCM5 / GEMCLIM
Atmospheric models Dynamics
Physics
Dynamics
Physics
● grid setup ● parallelism (MPI and OpenMP) ● memory allocations ● advection ● horizontal diffusion ● lateral boundary conditions ● ●
● radiation ● convection ● condensation ● surface schemes ● boundary layer vertical diffusion ● specified surface forcings (when no ocean model) ● 3
Atmospheric models at Environment Canada Atmospheric Science and Technology Directorate
GCM Model for climate change studies CRD/CCCma
GEMCLIM Monthly and seasonal forecasting
(Victoria)
GEM
Model for weather forecasting MRD/RPN (Dorval)
GCM (General Circulation Model)
GEM (Global Environmental Multiscale)
CRD (Climate Research Division) CCCma (Canadian Centre for Climate Modelling and Analysis)
MRD (Meteorological Research Division) RPN (Recherche en prévision numérique) 4 4
GEMCLIM GEMCLIM is the climate configuration of Environment Canada's forecast model GEM (Global Environmental Multiscale) global uniform
global variable
limited area LAM
GEMCLIM was originally developed to study systematic errors of the forecast model GEM and is routinely run at RPN as a validation tool. Operationally it is used for monthly and seasonal forecasting.
5
GEM
GEMCLIM
The main differences between GEM and GEMCLIM are in the scripts around the model: ● the automatic re-launching of a multi-month job and ● the possibility of an automatic large amount of post processing (The differences in the scripts do not affect the model results!) Almost no differences in the dynamics and physics libraries.
6
Library differences GEM v_3.3.0
GEMCLIM v_3.3.0
Dynamic library ● handling of nesting information ● possible pilot of temperature at the top level ● possible output of CCC variables (MKS units) ● reset accumulators after output ● diagnostic accumulators in permanent bus ● time series (hours real*8)
Physics library ● diagnostic accumulators in permanent bus ● different ozone climatology file ● minimum wind speed (Force-Restore, glacier, sea ice): GEM: 0.01 m/s GEMCLIM: 2.5 m/s 7
Runtime differences GEM (currently operational at CMC)
GEMCLIM (current climate runs)
Dynamics ● different grids (region, resolution) ● different levels: GEM: 58 levels, top at 10 hPa (global, stretched) 28 levels, top at 10 hPa (EPS) eta levels (normalized sigma) GEMCLIM: 60 levels, top at 2 hPa (global, stretched) 53 levels, top at 10 hPa (LAM) hybrid levels (rcoef=1.6) ● GEMCLIM: pilot temperature at the top level (only LAM) ● GEMCLIM: reset accumulators after output 8
Runtime differences GEM (currently operational at CMC)
GEMCLIM (current climate runs)
Physics ● different radiation schemes: GEM: newrad (Fouquart-Bonnel + Garand, radfix) GEMCLIM: correlated-k radiation (Li and Barker) with trace gas “climatology” ● GEMCLIM: non-orographic gravity wave drag (Hines) ● GEMCLIM: daily update SST and sea ice fraction ● use different ozone climatology GEMCLIM: daily ozone updates
9
Dynamics GCM (Victoria)
● parallelized ● only global ● spectral ● options specified at compilation ● not aggregated (grid points either land or ocean) ● semi-implicit, eulerian, 3 time levels
GEMCLIM (Dorval)
● massively parallelized ● global and regional (LAM + variable) ● gridded (spectral version being worked on) ● namelists ● aggregated (land fractions) ● implicit, semi-lagrangian, 2 time levels
10
Physics configurations GCM
(Victoria)
GEMCLIM (Dorval)
SST & Sea Ice: ● coupled ocean model version (CGCM)
SST & Sea Ice: ● Prescribed SST and sea ice fraction ● soon possibility to coupled ocean model
Aerosols: ● aerosol physics (sources, sinks and transport of sea salt, greenhouse gases and sulfates)
Aerosols: ● Information in SST and sea ice fraction ● will soon have an aerosol climatology (build from GCM aerosols)
Radiation: ● correlated-K radiation (Li and Barker)
Radiation: => ● correlated-k radiation (Li and Barker) with trace gas “climatology” ● newrad (Fouquart-Bonnel (1980), short-wave Garand and Mailhot (1990), long-wave Fomichev, stratospheric long-wave) 11
Physics configurations GCM
(Victoria)
GEMCLIM (Dorval)
Surface scheme: ● CLASS 2.7 ● next version CLASS 3.4
Surface scheme: ● ISBA ● Force-Restore ● soon CLASS 3.4
Convection: ● Zang-McFarlane
Convection: ● deep: • Kain-Fritsch • classical Kuo ● shallow: (Girard and Pellerin)
Clouds: ● diagnostic cloud scheme ● next version prognostic cloud scheme
Clouds: ● prognostic cloud scheme (Sundquist)
Gravity wave drag: ● Scinocca
Gravity wave drag: ● McFarlane (orographic) ● Hines (non-orographic) 12
Physics configurations GCM
(Victoria) ● Mass fix (log of surface pressure) ● Closed internal water budget
GEMCLIM (Dorval)
● Mass fix (log of surface pressure) only global
● working on new carbon cycle with evolving vegetation
13
CRCM5 at UQAM / Ouranos Global and regional, current and future climate scenarios
GCM CRD/CCCma
CRCM5 UQAM / Ouranos (Montreal)
dynamics physics
MRD/RPN (Dorval)
(Victoria)
dynamics
GEMCLIM
?
physics physics
dynamics physics 14
GEMCLIM and CRCM5 Since almost 2 years GEMCLIM is installed on the Linux cluster of the CRCMD network as CRCM5.
UQAM / Ouranos
MRD/RPN GEMCLIM v_3.2.2 GEMCLIM v_3.3.0
=> =>
CRCM5.0.0 CRCM5.0.1
15
Outlook for GEMCLIM / CRCM5 ●
Vertical staggered grid (Charney-Phillips) (RPN)
●
Single column mode (RPN)
●
Coupling of ocean model(s) (RPN, UQAM)
●
CLASS 3.4 (UQAM, RPN)
●
Bechtold-Kain-Fritsch convection (UQAM)
●
GCM4 physics (CCCma, RPN)
16
Current projects at UQAM J.-P. Paquin and L. Duarte: Implementation & Validation of CLASS 3.4 in CRCM5 Y. Jiao: Implementation of an updated Bechtold-Kain-Fritsch mass flux convection scheme in CRCM5 and evaluation along a cross-section over the Pacific Ocean M. Qian: Coupling RCO (Rossby Centre Ocean) model to CRCM5 Z. Kothavala: Implementation & validation of CTEM (Canadian Terrestrial Ecosystem Model) in CRCM5
17
Recent runs done with GEMCLIM and CRCM5 All runs were done with the configuration of the GEMCLIM physics library as described above. Differences between model version: ● constantly developing radiation scheme correlated-k ● bugfixes • AMIP SST & sea ice (shift, lakes) • Correct call to lacs1 => -3% of the sea ice fraction OUT
18
global, 0.5°
global, 1.0°
global, 1.5°
GEMCLIM v_3.2.1 Period: 197801 - 200402 Output: full diagnostics (+ 3 hourly high res area) (+ pilot files from 1.5° global) SGMIP2 project: ● stretched grid, North America, 0.5°-1.5° ● stretched grid, Europe, 0.5°-1.5° ● global, 1.0° ● global, 0.5° (just Jan.+Jul.)
Also done: ● LAM over high res of NA stretched grid, 0.5° ● LAM over high res of EU stretched grid, 0.5° ● global, 1.5° (used as pilot for LAM's)
LAM, NA
stretched, NA
LAM, EU
stretched, EU 19
stretched grid over North America, 0.5°-1.5°
stretched grid over Europe, 0.5°-1.5°
20
GEMCLIM v_3.2.2 Global, 2.0° Period: 197801 - 200402 Output: full diagnostics LAM over North America, 0.5° piloted with ERA40 Period: 195709 - 200208 Output: full diagnostics + station time series starting 1998
CRCM5.0.0 Global, 2.0° Period: 197809 - 200611 Output: 3 hourly output + pilot files Stretched grid over hurricane tracks, 0.3°-2.0° Period: 197809 - 200611 Output: 3 hourly output 21
CRCM5.0.1 Stretched grid over hurricane tracks, 0.3°-2.0° Period: 197809 - 200611 Output: 3 hourly output 2 LAM grids over hurricane tracks, 0.3° Piloted with ERA40 resp. ECMWF analysis data and also with CRCM5 global 2.0°. Period: 197906 – 200611; June - November Output: 3 hourly output
Stretched grid over Indian Ocean and tropical Pacific, 0.5°-2.0° Period: 197801 - 200611 Output: 3 hourly output Stretched grid over North America, 0.5°-2.0° Period: 197801 - 200611 Output: 3 hourly output 22
stretched grid over Hurricane tracks, 0.3°-2.0°
stretched grid over Indian Ocean and tropical Pacific, 0.5°-2.0°
23
CRCM5.0.1 Global, 2.0° Period: 197801 - 200402 Output: full diagnostics LAM over North America, 0.5° Piloted with ERA40 Period: 195709 - 200208 Output: 3 hourly output LAM over North America, 0.25° Piloted with ERA40 Period: 195709 - 200208 Output: 3 hourly output
ENSEMBLES project LAM over Europe, 0.22° Piloted with ERA40 Period: 195709 - 200208 Output: full diagnostics + several 2D fields 3 hourly 24
GEMCLIM v_3.3.0 ICTS project (Inter-Continental Scale Experiments Transferability Study) LAM piloted with NCEP, 0.5°, 7 different regions: Canada, USA, South America, Europe, Africa, Asia, Australia Period: 199907 - 200411 Output: full diagnostics + time series
25
CRCM5.0.1 Arctic 1 stretched grid, 0.5°-2.0° Period: 197801 – 200611 7 LAM runs, 0.5° Piloted with ERA40 resp. CRCM5 global 2.0° Period: 197801 – 200611, resp. 197706 – max. 200208 Output: 3 hourly output
26
Arctic 1 stretched grid, 0.5°-2.0°
Largest LAM grid, 0.5°
27
Research activities North America D. Paquin-Ricard: A study on the representation of cloud microphysics and its interaction with radiation in the GEM model Reference: Evaluation of cloud representation in the Canadian GEM model using ARM data, D. Paquin-Ricard, C. G. Jones, P. Vaillancourt, Geophysical Research Abstracts, Vol. 9, 5-2-2007
M. Markovic: Evaluation of the Surface Radiation Budget over North America for a suite of Regional Climate Models Reference: The Surface Radiation Budget over North America: An assessment of Gridded Data Sets for Model Evaluation and the Evaluation of a Suite of Regional Climate Models, Markovic, M., C. Jones, K. Winger, D. Paquin, International Journal of Climatology, in print.
F. Dorsaz: Evaluation of regional climate model simulations of snow cover over Québec L. Separovic: Methodological approach to parameter perturbations in GEM-LAM seasonal simulations Y. He (UVic): Surface wind probability distributions over N. American regions: observations and RCM simulations E. P. Diaconescu: Analysis of Internal Variability of a regional Climate Model using Singular Vectors
28
North America & Europe M. Verville: Comparison of two regional climate modelling approaches using the GEM model, global variable-resolution versus one-way nested limited-area ICTS Z. Kothavala: The Transferability of Regional Climate Models to non-native domains Arctic M. Qian: The behavior of GEM-LAM over the Arctic using different simulation domains Hurricanes L.-P. Caron: A study on tropical cyclone activity using the GEM model Tropics - mid-latitudes M. Markovic: Tropical mid-latitude interactions Sheba P.-L. Carpentier: Evaluation of the stable boundary layer processes in GEM over the Arctic ocean during SHEBA
29
CRCM5.0.1 Africa A.-S. Daloz: Study of the benefits of increased resolution on the precipitation in Sub-Saharan Africa Divers C. Rosu and M. Radojevic (Philippe Gachon) P. Du, G. Dueymes, J. Dorais, and D. Simjanovski (Eric Girard) P. Dehasse and Y. Melin (Jean-Pierre Blanchet)
30
Merci... Check here fore recent runs: http://people.sca.uqam.ca/~winger/GEM/Version_3.3.0/Recent_runs.html
31
Overview ●
What is CRCM5? ● GEM ● GEMCLIM ● GCM
●
Current and future projects for CRCM5 / GEMCLIM
●
Recent runs done with CRCM5 / GEMCLIM
Atmospheric models Dynamics
Physics
Dynamics
Physics
● grid setup ● parallelism (MPI and OpenMP) ● memory allocations ● advection ● horizontal diffusion ● lateral boundary conditions ● ●
● radiation ● convection ● condensation ● surface schemes ● boundary layer vertical diffusion ● specified surface forcings (when no ocean model) ● 3
Atmospheric models at Environment Canada Atmospheric Science and Technology Directorate
GCM Model for climate change studies CRD/CCCma
GEMCLIM Monthly and seasonal forecasting
(Victoria)
GEM
Model for weather forecasting MRD/RPN (Dorval)
GCM (General Circulation Model)
GEM (Global Environmental Multiscale)
CRD (Climate Research Division) CCCma (Canadian Centre for Climate Modelling and Analysis)
MRD (Meteorological Research Division) RPN (Recherche en prévision numérique) 4 4
GEMCLIM GEMCLIM is the climate configuration of Environment Canada's forecast model GEM (Global Environmental Multiscale) global uniform
global variable
limited area LAM
GEMCLIM was originally developed to study systematic errors of the forecast model GEM and is routinely run at RPN as a validation tool. Operationally it is used for monthly and seasonal forecasting.
5
GEM
GEMCLIM
The main differences between GEM and GEMCLIM are in the scripts around the model: ● the automatic re-launching of a multi-month job and ● the possibility of an automatic large amount of post processing (The differences in the scripts do not affect the model results!) Almost no differences in the dynamics and physics libraries.
6
Library differences GEM v_3.3.0
GEMCLIM v_3.3.0
Dynamic library ● handling of nesting information ● possible pilot of temperature at the top level ● possible output of CCC variables (MKS units) ● reset accumulators after output ● diagnostic accumulators in permanent bus ● time series (hours real*8)
Physics library ● diagnostic accumulators in permanent bus ● different ozone climatology file ● minimum wind speed (Force-Restore, glacier, sea ice): GEM: 0.01 m/s GEMCLIM: 2.5 m/s 7
Runtime differences GEM (currently operational at CMC)
GEMCLIM (current climate runs)
Dynamics ● different grids (region, resolution) ● different levels: GEM: 58 levels, top at 10 hPa (global, stretched) 28 levels, top at 10 hPa (EPS) eta levels (normalized sigma) GEMCLIM: 60 levels, top at 2 hPa (global, stretched) 53 levels, top at 10 hPa (LAM) hybrid levels (rcoef=1.6) ● GEMCLIM: pilot temperature at the top level (only LAM) ● GEMCLIM: reset accumulators after output 8
Runtime differences GEM (currently operational at CMC)
GEMCLIM (current climate runs)
Physics ● different radiation schemes: GEM: newrad (Fouquart-Bonnel + Garand, radfix) GEMCLIM: correlated-k radiation (Li and Barker) with trace gas “climatology” ● GEMCLIM: non-orographic gravity wave drag (Hines) ● GEMCLIM: daily update SST and sea ice fraction ● use different ozone climatology GEMCLIM: daily ozone updates
9
Dynamics GCM (Victoria)
● parallelized ● only global ● spectral ● options specified at compilation ● not aggregated (grid points either land or ocean) ● semi-implicit, eulerian, 3 time levels
GEMCLIM (Dorval)
● massively parallelized ● global and regional (LAM + variable) ● gridded (spectral version being worked on) ● namelists ● aggregated (land fractions) ● implicit, semi-lagrangian, 2 time levels
10
Physics configurations GCM
(Victoria)
GEMCLIM (Dorval)
SST & Sea Ice: ● coupled ocean model version (CGCM)
SST & Sea Ice: ● Prescribed SST and sea ice fraction ● soon possibility to coupled ocean model
Aerosols: ● aerosol physics (sources, sinks and transport of sea salt, greenhouse gases and sulfates)
Aerosols: ● Information in SST and sea ice fraction ● will soon have an aerosol climatology (build from GCM aerosols)
Radiation: ● correlated-K radiation (Li and Barker)
Radiation: => ● correlated-k radiation (Li and Barker) with trace gas “climatology” ● newrad (Fouquart-Bonnel (1980), short-wave Garand and Mailhot (1990), long-wave Fomichev, stratospheric long-wave) 11
Physics configurations GCM
(Victoria)
GEMCLIM (Dorval)
Surface scheme: ● CLASS 2.7 ● next version CLASS 3.4
Surface scheme: ● ISBA ● Force-Restore ● soon CLASS 3.4
Convection: ● Zang-McFarlane
Convection: ● deep: • Kain-Fritsch • classical Kuo ● shallow: (Girard and Pellerin)
Clouds: ● diagnostic cloud scheme ● next version prognostic cloud scheme
Clouds: ● prognostic cloud scheme (Sundquist)
Gravity wave drag: ● Scinocca
Gravity wave drag: ● McFarlane (orographic) ● Hines (non-orographic) 12
Physics configurations GCM
(Victoria) ● Mass fix (log of surface pressure) ● Closed internal water budget
GEMCLIM (Dorval)
● Mass fix (log of surface pressure) only global
● working on new carbon cycle with evolving vegetation
13
CRCM5 at UQAM / Ouranos Global and regional, current and future climate scenarios
GCM CRD/CCCma
CRCM5 UQAM / Ouranos (Montreal)
dynamics physics
MRD/RPN (Dorval)
(Victoria)
dynamics
GEMCLIM
?
physics physics
dynamics physics 14
GEMCLIM and CRCM5 Since almost 2 years GEMCLIM is installed on the Linux cluster of the CRCMD network as CRCM5.
UQAM / Ouranos
MRD/RPN GEMCLIM v_3.2.2 GEMCLIM v_3.3.0
=> =>
CRCM5.0.0 CRCM5.0.1
15
Outlook for GEMCLIM / CRCM5 ●
Vertical staggered grid (Charney-Phillips) (RPN)
●
Single column mode (RPN)
●
Coupling of ocean model(s) (RPN, UQAM)
●
CLASS 3.4 (UQAM, RPN)
●
Bechtold-Kain-Fritsch convection (UQAM)
●
GCM4 physics (CCCma, RPN)
16
Current projects at UQAM J.-P. Paquin and L. Duarte: Implementation & Validation of CLASS 3.4 in CRCM5 Y. Jiao: Implementation of an updated Bechtold-Kain-Fritsch mass flux convection scheme in CRCM5 and evaluation along a cross-section over the Pacific Ocean M. Qian: Coupling RCO (Rossby Centre Ocean) model to CRCM5 Z. Kothavala: Implementation & validation of CTEM (Canadian Terrestrial Ecosystem Model) in CRCM5
17
Recent runs done with GEMCLIM and CRCM5 All runs were done with the configuration of the GEMCLIM physics library as described above. Differences between model version: ● constantly developing radiation scheme correlated-k ● bugfixes • AMIP SST & sea ice (shift, lakes) • Correct call to lacs1 => -3% of the sea ice fraction OUT
18
global, 0.5°
global, 1.0°
global, 1.5°
GEMCLIM v_3.2.1 Period: 197801 - 200402 Output: full diagnostics (+ 3 hourly high res area) (+ pilot files from 1.5° global) SGMIP2 project: ● stretched grid, North America, 0.5°-1.5° ● stretched grid, Europe, 0.5°-1.5° ● global, 1.0° ● global, 0.5° (just Jan.+Jul.)
Also done: ● LAM over high res of NA stretched grid, 0.5° ● LAM over high res of EU stretched grid, 0.5° ● global, 1.5° (used as pilot for LAM's)
LAM, NA
stretched, NA
LAM, EU
stretched, EU 19
stretched grid over North America, 0.5°-1.5°
stretched grid over Europe, 0.5°-1.5°
20
GEMCLIM v_3.2.2 Global, 2.0° Period: 197801 - 200402 Output: full diagnostics LAM over North America, 0.5° piloted with ERA40 Period: 195709 - 200208 Output: full diagnostics + station time series starting 1998
CRCM5.0.0 Global, 2.0° Period: 197809 - 200611 Output: 3 hourly output + pilot files Stretched grid over hurricane tracks, 0.3°-2.0° Period: 197809 - 200611 Output: 3 hourly output 21
CRCM5.0.1 Stretched grid over hurricane tracks, 0.3°-2.0° Period: 197809 - 200611 Output: 3 hourly output 2 LAM grids over hurricane tracks, 0.3° Piloted with ERA40 resp. ECMWF analysis data and also with CRCM5 global 2.0°. Period: 197906 – 200611; June - November Output: 3 hourly output
Stretched grid over Indian Ocean and tropical Pacific, 0.5°-2.0° Period: 197801 - 200611 Output: 3 hourly output Stretched grid over North America, 0.5°-2.0° Period: 197801 - 200611 Output: 3 hourly output 22
stretched grid over Hurricane tracks, 0.3°-2.0°
stretched grid over Indian Ocean and tropical Pacific, 0.5°-2.0°
23
CRCM5.0.1 Global, 2.0° Period: 197801 - 200402 Output: full diagnostics LAM over North America, 0.5° Piloted with ERA40 Period: 195709 - 200208 Output: 3 hourly output LAM over North America, 0.25° Piloted with ERA40 Period: 195709 - 200208 Output: 3 hourly output
ENSEMBLES project LAM over Europe, 0.22° Piloted with ERA40 Period: 195709 - 200208 Output: full diagnostics + several 2D fields 3 hourly 24
GEMCLIM v_3.3.0 ICTS project (Inter-Continental Scale Experiments Transferability Study) LAM piloted with NCEP, 0.5°, 7 different regions: Canada, USA, South America, Europe, Africa, Asia, Australia Period: 199907 - 200411 Output: full diagnostics + time series
25
CRCM5.0.1 Arctic 1 stretched grid, 0.5°-2.0° Period: 197801 – 200611 7 LAM runs, 0.5° Piloted with ERA40 resp. CRCM5 global 2.0° Period: 197801 – 200611, resp. 197706 – max. 200208 Output: 3 hourly output
26
Arctic 1 stretched grid, 0.5°-2.0°
Largest LAM grid, 0.5°
27
Research activities North America D. Paquin-Ricard: A study on the representation of cloud microphysics and its interaction with radiation in the GEM model Reference: Evaluation of cloud representation in the Canadian GEM model using ARM data, D. Paquin-Ricard, C. G. Jones, P. Vaillancourt, Geophysical Research Abstracts, Vol. 9, 5-2-2007
M. Markovic: Evaluation of the Surface Radiation Budget over North America for a suite of Regional Climate Models Reference: The Surface Radiation Budget over North America: An assessment of Gridded Data Sets for Model Evaluation and the Evaluation of a Suite of Regional Climate Models, Markovic, M., C. Jones, K. Winger, D. Paquin, International Journal of Climatology, in print.
F. Dorsaz: Evaluation of regional climate model simulations of snow cover over Québec L. Separovic: Methodological approach to parameter perturbations in GEM-LAM seasonal simulations Y. He (UVic): Surface wind probability distributions over N. American regions: observations and RCM simulations E. P. Diaconescu: Analysis of Internal Variability of a regional Climate Model using Singular Vectors
28
North America & Europe M. Verville: Comparison of two regional climate modelling approaches using the GEM model, global variable-resolution versus one-way nested limited-area ICTS Z. Kothavala: The Transferability of Regional Climate Models to non-native domains Arctic M. Qian: The behavior of GEM-LAM over the Arctic using different simulation domains Hurricanes L.-P. Caron: A study on tropical cyclone activity using the GEM model Tropics - mid-latitudes M. Markovic: Tropical mid-latitude interactions Sheba P.-L. Carpentier: Evaluation of the stable boundary layer processes in GEM over the Arctic ocean during SHEBA
29
CRCM5.0.1 Africa A.-S. Daloz: Study of the benefits of increased resolution on the precipitation in Sub-Saharan Africa Divers C. Rosu and M. Radojevic (Philippe Gachon) P. Du, G. Dueymes, J. Dorais, and D. Simjanovski (Eric Girard) P. Dehasse and Y. Melin (Jean-Pierre Blanchet)
30
Merci... Check here fore recent runs: http://people.sca.uqam.ca/~winger/GEM/Version_3.3.0/Recent_runs.html
31
