| 1 | Update 05/02/2018 |
|---|
| 2 | Thibaut Lurton |
|---|
| 3 | Scripts now have paths directed to version 3.2 of the data. |
|---|
| 4 | /prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/SOLARIS-HEPPA/SOLARIS-HEPPA-3-2/atmos/ |
|---|
| 5 | |
|---|
| 6 | |
|---|
| 7 | Solar forcing dataset from CMIP6 |
|---|
| 8 | See http://www.wcrp-climate.org/wgcm-cmip/wgcm-cmip6 |
|---|
| 9 | and http://solarisheppa.geomar.de/cmip6 |
|---|
| 10 | |
|---|
| 11 | Version 3.1 |
|---|
| 12 | Not distributed through ESGF yet |
|---|
| 13 | The data is available on ciclad in /prodigfs/project/input4MIPs/SOLAR/ |
|---|
| 14 | |
|---|
| 15 | Data was processed by Rong-Ming in June 2016 and reprocessed in October 2016 by Olivier Boucher. |
|---|
| 16 | The final versions of the script are CMIP6_historical.m and CMIP6_piControl.m |
|---|
| 17 | it can be run in matlab on ciclad only on a line command |
|---|
| 18 | > run CMIP6_historical.m |
|---|
| 19 | A script process.sh runs all scripts and post-processes the files for use in IPSL-CM6. |
|---|
| 20 | |
|---|
| 21 | The data are used to prescribe the solar constant (tsi) and its distribution (ssi_frac) in the new |
|---|
| 22 | (RRTM) SW scheme of the LMDz model. This is activated with the ok_suntime_rrtm=y flag. |
|---|
| 23 | The solar constant and RSUN are output in the paramphy netcdf output file. |
|---|
| 24 | |
|---|
| 25 | Daily data are used. The number of timesteps in the nc file has to match the number of days |
|---|
| 26 | of the year in the model (ie 360, 365 or 366 days). Otherwise the model stops. |
|---|
| 27 | A tolerance was added so that a 366-day file can be used in a 365-day model year. |
|---|
| 28 | |
|---|
| 29 | Files are prepared for two bands (NSW=2) or 6 bands (NSW=6) option. Default NSW=6 is used ! |
|---|
| 30 | |
|---|
| 31 | For piControl and pdControl, a single value (repeated over the number of days) is used: |
|---|
| 32 | solarforcing_v3.1_XXControl_ave_Ybands_ZZZdays.nc |
|---|
| 33 | XX=pi or pd for pre-industrial or present-day control |
|---|
| 34 | Y=2 or 6 depending on whether NSW=2 or 6 |
|---|
| 35 | ZZZ=360, 365 or 366 days according to which calendar (and number of days in the year) is chosen. |
|---|
| 36 | |
|---|
| 37 | The average for piControl consists of the first two solar cycles after 1850 (as per CMIP6 protocol) |
|---|
| 38 | The average for pdControl consists of a 1991-2010 average (not an official CMIP6 run). |
|---|
| 39 | |
|---|
| 40 | Annual files are available for the historical and future period: |
|---|
| 41 | solarforcing_v3.1_daily_6bands_XXXX.nc |
|---|
| 42 | where XXXX is the year. The files are for a Gregorian calendar (365 and 366 days according to the year). |
|---|
| 43 | |
|---|
| 44 | There is also a scenario with a lower solar constant in 21st century for sensitivity study. |
|---|
| 45 | This is not an official CMIP6 scenario though. |
|---|
| 46 | |
|---|
| 47 | The solar forcing file is to be renamed solarforcing.nc in the running directory. |
|---|
| 48 | |
|---|
| 49 | Olivier Boucher |
|---|
| 50 | 11/10/2016 |
|---|
| 51 | revised 27/09/2017 |
|---|
