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