Difference between revisions of "GENIE ExampleGENIE1expt"

From SourceWiki
Jump to navigation Jump to search
Line 35: Line 35:
 
** e.g. /genie-output/genie_eb_go_gs_expt2/goldstein/spn.opsit - this file lists maximum and minimum meridional overturning streamfunctions in the Pacific and Atlantic and their locations.
 
** e.g. /genie-output/genie_eb_go_gs_expt2/goldstein/spn.opsit - this file lists maximum and minimum meridional overturning streamfunctions in the Pacific and Atlantic and their locations.
 
*** compare the values between expt1 and expt2 to see the effect of doubling CO2  
 
*** compare the values between expt1 and expt2 to see the effect of doubling CO2  
* You can look at the netCDF results using Panoply or Matlab scripts.
+
* You can look at the netCDF results using e.g. Panoply.

Revision as of 11:38, 24 July 2007

A simple experiment examining the effects of atmosphere parameters on climate and its sensitivity

Aim: To examine the spatial pattern of temperature for genie_eb_go_gs (default 36x36x8 grid), examine its sensitivity to varying parameters in the simple atmosphere model (EMBM), and evaluate the model climate sensitivity

Setting up the model

  • You will be working with genie_eb_go_gs configuration
  • Navigate to genie-main and open runtime_defaults.sh in order to see all model parameters
  • Go to the EMBM (ea_...) parameters and note values of ea_12 to ea_15 which are all the parameters controlling the diffusion of heat and moisture
  • Go to /configs and copy genie_eb_go_gs.config to genie_eb_go_gs_expt1.config (or similar)
  • In the new config file change the EXPID to genie_eb_go_gs_expt1
  • Set the run length to 2000 years by changing ma_koverall_total=1000000
    • (Note: this is the total number of atmospheric timesteps and by default the EMBM takes 5 timesteps per ocean timestep and the ocean takes 100 timesteps per year, therefore 500 timesteps per year for the EMBM and 500*2000=1000000)
  • You need to control model data output. We recommend the following choices:
    • changing the frequency of ‘health checks’, by setting ea_3, go_3, gs_3 = 10000 (every 100 years, in goldstein timesteps)
    • changing the frequency of ‘full’ (or restart) model datatsets, by setting ea_4, go_4, gs_4 = 20000 (every 200 years, in goldstein timesteps – to give 10 datasets)
    • changing the frequency of timeseries output, by setting ea_5, go_5, gs_5 = 100 (to give annual 1 January snapshots)
  • First we are going to run at pre-industrial (278ppmv) CO2 level, which requires no change to ea_20
  • However, you can generate a unique instance of the model by varying key EMBM parameters ea_12 to ea_15. Some suggested reasonable ranges for these parameters are as follows:
    • ea_12 = 1.0e6 to 8.0e6 (default 5.0e6 already quite high)
    • ea_13 = 1.0e5 to 1.0e7 (default 1e6)
    • ea_14 = 0.5 to 1.5 (default 1.0 radians)
    • ea_15 = 0.0 to 0.2 (default 0.1)

Running the model

  • Log into a node by using
    • qrsh –q interactive.q
  • Start the model experiment from /genie-main by using
    • ./genie_example.job –f configs/ genie_eb_go_gs_expt1.config
  • You will now get output to the screen and the experiment will take around an hour
  • Meanwhile output datasets will be written to /genie-output/genie_eb_go_gs_expt1

Further work

  • Repeat the above but create genie_eb_go_gs_expt2.config in which you change EXPID to genie_eb_go_gs_expt2, and set ea_20 = 2.0 (for 2x pre-industrial CO2)
  • You can view the ascii output files
    • e.g. /genie-output/genie_eb_go_gs_expt2/embm/spn.airt – this file lists hemispheric mean air temperatures and global maximum and minimum air temperatures and their locations.
      • compare the values between expt1 and expt2 to see the effect of doubling CO2
    • e.g. /genie-output/genie_eb_go_gs_expt2/goldstein/spn.opsit - this file lists maximum and minimum meridional overturning streamfunctions in the Pacific and Atlantic and their locations.
      • compare the values between expt1 and expt2 to see the effect of doubling CO2
  • You can look at the netCDF results using e.g. Panoply.