...
Section | ||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
Initial conditions
Info |
---|
Note: The initial conditions data for the case studies described here are available for OpenIFS 40r1. Please contact OpenIFS Support if you require initial experiment data for more recent model version. |
Case study initial conditions for the Lothar storm are provided on the OpenIFS ftp site. The ftp site is password protected, only licensed institutes may be provided with the ftp password. Please contact: openifs-support@ecmwf.int.
The Lothar depression developed initially on 24th December off the North American east coast at 35N.
...
Info |
---|
As OpenIFS is a spectral model, the 'T' number refers to the triangular truncation in spectral space. Equivalent grid resolutions are: The number of vertical levels is given after the letter 'L' e.g. L62 means 62 vertical levels. Please note that higher resolutions progressively require more processors and computer memory to run. |
Panel | |||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| |||||||||||||||||||||||||||||||||||||||||||||||
|
...
Panel | |||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| |||||||||||||||||||||||||||||||||||
|
...
Download instructions
Code Block | ||
---|---|---|
| ||
% mkdir -p runs/lothar/t159 % cd runs % ftp ftp.ecmwf.int ftp> cd case_studies/lothar_storm ftp> binary ftp> get 1999122412_T159_fqar.tgz ftp> quit % tar zxf 1999122412_T159_fqar.tgz % ls 1999122412_T159.tgz ICMCLfqarINIT ICMGGfqarINIT ICMGGfqarINIUA ICMSHfqarINIT ecmwf % ls ecmwf NODE.001_01 namelistfc |
...
- What's the impact of the different 'lead times' on the forecast of the storm (i.e. starting from different dates)?
- What's the impact of resolution on the forecast of the storm: both for it's development and impact over areas worse hit in Europe?
- Does reducing the model timestep improve or worsen the forecast?
Reduce the gravity wave drag - how does this affect the forecast in the upper and lower levels?
Expand title How to change the code (click here to expand) Edit the source code to half the gravity wave drag coefficient
File: ifs/phys_ec/sugwd.F90, change:
Code Block Line 108: ! Revised gwd parameter values Line 109: GKDRAG =0.15_JPRB
to:
Code Block Line 108: ! Revised gwd parameter values Line 109: GKDRAG = 0.075_JPRB ! half GWD coefficient: 0.15_JPRB
Increase the precipitation auto conversion rate - what impact does this have?
Expand title How to change the code (click here to expand…) Edit the source code to increase the auto conversion rate by 20%
File: ifs/phys_ec/sucldp.F90, change:
Code Block line 123: RKCONV=1._JPRB/6000._JPRB ! 1/autoconversion time scale (s)
to:
Code Block line 123: ! RKCONV=1._JPRB/6000._JPRB ! 1/autoconversion time scale (s) line 124: RKCONV=1.2_JPRB/6000._JPRB ! default scaled by 20%: 1/autoconversion time scale (s)
Change the surface transfer coefficient in the turbulence scheme
Expand title How to change the model code (click here to expand) Reduce the coefficient by 20%.
Alter surf/module/surfexcdriver_ctl_mod.F90 from :
Code Block line 671: DO JL=KIDIA,KFDIA line 672: IF (JTILE == IFRMAX(JL)) THEN line 673: PKHLEV(JL)=ZKHLEV(JL) line 674: ENDIF line 675: ENDDO
to:
Code Block line 671: DO JL=KIDIA,KFDIA line 672: IF (JTILE == IFRMAX(JL)) THEN line 673: PKHLEV(JL)=ZKHLEV(JL) line 674: ENDIF line 675: ! reduce surface transfer coeff by 20% in turbulence scheme line 676: ZCFMTI(JL,JTILE)=0.8_JPRB*ZCFMTI(JL,JTILE) line 677: ENDDO
Reduce the asymptotic mixing length scale (K) - how does this affect surface & near-surface fields?
Expand title How to change the model code (click here to expand) For this change, two files need to be edited:
Code Block title ifs/phys_ec/suvdf.F90 line 53: RLAM = 75.0_jprb !! 150._JPRB: reduce to 75m
and:
Code Block title ifs/phys_ec/vdfexcu.F90 ZKLEN = 75.0_jprb !! 150.0_JPRB ! asymptotic K length scale troposphere - Reduce to 75m
- For these last 4 cases where the model's parametrizations have been altered, which make the biggest difference and why? Does any of the changes improve the forecast in any way?
- If you were providing forecasts for wind and precipitation to the general public based on these experiments, what could you say with certainty and what is less certain? How would this change over different countries?
...
Further reading
Ulbrich et al., 2001, Weather, 56, 70-80
...
A. Hannachi, J. Kjellsson, M. Tjernström, G. Carver, 2012, Teaching with the OpenIFS at Stockholm University, ECMWF Newsletter No. 134, Winter 2012/13.
Comments
...
We gratefully acknowledge: Dr Anton Beljaars (ECMWF) for suggestions and code changes for the parametrization changes in the list of sensitivity experiments; Prof Erland Kallen for reviewing & comments on the text.
HTML |
---|
<script type="text/javascript" src="https://softwarejira.ecmwf.int/issues/s/en_UKet2vtj/787/12/1.2.5/_/download/batch/com.atlassian.jira.collector.plugin.jira-issue-collector-plugin:issuecollector/com.atlassian.jira.collector.plugin.jira-issue-collector-plugin:issuecollector.js?collectorId=5fd84ec6"></script> |
Excerpt Include | ||||||
---|---|---|---|---|---|---|
|