Changes between Version 13 and Version 14 of Documentation/UserGuide/VerifyClimate
- Timestamp:
- 2014-11-18T10:07:29+01:00 (10 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Documentation/UserGuide/VerifyClimate
v13 v14 1 1 Once you have run LMDZOR, the first question to ask yourself is: does the climate I'm getting make sense? There does not seem to be any straightforward way to answer this question, so this page is going to give ideas on how to procede. 2 2 3 == Setting the boundary conditions==3 == Setting the boundary conditions == 4 4 Rubbish in is rubbish out so first of all the boundary conditions of the coupled simulation should be acceptable: 5 5 * Run a spinup of ORCHIDEE offline to get reasonable biomass and soil water pools, which seems to be fairly important if you are using the 11-layer hydrology. The only way to do this that I know of is to regrid the driver files to the resolution you are interested in (such files now exist for the CRU dataset for 128x118) and use those to run offline for 50 years or so. Then you can use the resulting restart files to start the LMDZOR calculation. This is not a trivial step! You also need to make sure that the land mask you are using matches what is found in LMDZOR. The way we came up with to do this, was to run a one day simulation of LMDZOR without a restart and look at the resulting history files. There are variables inside the histday*nc files called CONTFRACOR and CONTFRACATM, which contain the fraction of each gridcell which have land. We made sure that all these pixels had climate variables in our regridded forcing files, even if it meant copying the forcing data from nearby pixels. 6 6 * It takes at least 10 simulation years to get a good idea of the climate. It also takes some time to equilibrate the climate after starting from a restart, so it seems that you can plan on a simulation of at least 15 years: 5 years to equilibrate, and 10 years over which to calculate the average climate. 7 7 8 == Top of the atmosphere radiation==8 == Top of the atmosphere radiation == 9 9 Prior to the industrial revolution the climate was stable. Hence, for a pre-industrial coupled simulation the net radiation at the top of the atmosphere should be close to zero. Even if your model experiment does not addresses pre-industrial conditions, you should launch a pre-industrial run and check for this condition. Then the same settings for the pre-industrial run and the experiments should be used. Remember to normalize the pixels for their surface area, remember to account for both land and ocean (especially if you are an ORCHIDEE-user) and use annual mean values as energy is stored and released over the season. The variables that should cancel each other out are TOPL and TOPS in ATM. 10 10 11 == Compare against other climate data==11 == Compare against other climate data == 12 12 Once you have the simulation, you need to do the averaging, find data to compare to, and create graphs. The tool ATLAS does exactly this. 13 13 … … 23 23 * Step 2 then compares this data to the reference simulation or observations, depending on what you are trying to do. 24 24 25 === Compare against other simulations===25 === Compare against other simulations === 26 26 As an example, let's say I have two simulations that I want to compare. 27 27 /ccc/store/cont003/dsm/p529grat/IGCM_OUT/LMDZOR/PROD/amip/MCOUP/ATM/Analyse/SE/MCOUP_SE_1980_1981_1M_histmth.nc … … 61 61 The output of this command is a series of folders all grouped into "atlas_model_vs_model/". In there you have pdf and html summaries of all the variables, plus individual directories for each variable if you just want to look at one of them. 62 62 63 === Zooming in on the area of interest ===63 === Zooming in on the area of interest === 64 64 65 65 Sometimes you are only interested in a certain region of the world. This is certainly the case when you are using a zoomed grid to run the simulation. Is there a way to get ATLAS to only print out the results over a single region? … … 95 95 NOTE: Do not use 20W; use -20E instead to mean 20 degrees west longitude. While it is true that 20W is the same as -20E, Ferret seems to have trouble with some of the commands if you use 20W, and the maps come out strange. I assume the same is true for 20S and -20N. 96 96 97 === Compare against observations ===97 === Compare against observations === 98 98 For pre-industrial runs this will be difficult but for post-1900 simulations this should be standard procedure. 99 99
