test/em_scm_xy/README.scm

   1 This is a single column model test case.
   2
   3 README updated 9 Feb. 2009.
   4
   5 -------------------------------------------------------------------------
   6
   7 Author: Josh Hacker, NCAR
   8
   9 -------------------------------------------------------------------------
  10
  11 Please communicate your changes so that we can all benefit from model
  12 improvements.
  13
  14 -------------------------------------------------------------------------
  15
  16 Description:  The SCM-WRF runs on a 3x3 stencil with periodic lateral
  17 boundary conditions in X and Y.  No horizontal gradients are supported
  18 in this configuration so advection is absent unless explicitly imposed.
  19 The formulation for advection in the current version, coded in
  20 module_scm_force.F, is an upstream relaxation following Ghan et al. (1999).
  21
  22 -------------------------------------------------------------------------
  23
  24 Initialization: As with the other test cases, initialization is handled
  25 with text files input to ideal.exe (via module_initialize_scm_xy.F).
  26 Here we use both files input_sounding and input_soil.  The format for
  27 input_sounding is the surface followed by the sounding:
  28
  29 z_terrain u_10 v_10 t_2 q_2 psfc
  30 z u v theta qv
  31 . . .   .    .
  32 . . .   .    .
  33 . . .   .    .
  34
  35 Here z is in meters, u and v are in m/s, theta is in K, and qv is in
  36 kg^3/kg^3.  The number of levels does not matter as long as z_top <=
  37 the top level in input_sounding.  The format for input soil is the skin
  38 and deep soil temperatures followed by the soil level data:
  39
  40 0.0  TSK   TMN
  41 z    SOILT SOILM
  42 .    .     .
  43 .    .     .
  44 .    .     .
  45
  46 Here z is positive downward, in meters.  TSK, SOILT, and TMN are in K,
  47 and SOILM is the volumetric soil moisture fraction.  Initialization of
  48 the soil uses the modules for real-data cases with your specified LSM,
  49 so use wisely.
  50
  51 -------------------------------------------------------------------------
  52
  53 Forcing:  Solar forcing will follow the dates as specified in the
  54 namelist.  All other forcing (U_g, V_g, W, and advection) are specified
  55 in a file that comes in via an auxilliary input.  For example:
  56
  57 &time_control
  58 auxinput3_inname                    = "force_ideal.nc"
  59 auxinput3_interval_h                = 59
  60
  61 This file looks very much like wrfbdy_d01 and in fact is used in the
  62 same way, with values specified at the intervals and tendencies on
  63 those values valid for the period within the corresponding interval.
  64 Thus time-dependent forcing on time-dependent height coordinates is
  65 possible; one simply needs to use multiple intervals.
  66
  67 A script to produce a two-time level (i.e. constant) forcing is
  68 available, and is called make_scm_forcing.ncl.  See the documentation
  69 at the beginning of this.  The example that comes with the release is
  70 documented below.  I have a much more sophisticated script to create
  71 forcing from a series of met_em files.  Please contact me if you would
  72 like to work with me on that (hacker at ucar.edu).
  73
  74 -------------------------------------------------------------------------
  75
  76 Namelist: For namelist parameters, see README.namelist.  Non-SCM namelist
  77 parameters that must be set to get the expected behavior are:
  78
  79 &dynamics
  80 pert_coriolis                       = .true.,
  81
  82 &bdy_control
  83 periodic_x                          = .true.,
  84 periodic_y                          = .true.,
  85 (all others false)
  86
  87 -------------------------------------------------------------------------
  88
  89 Example:  This example is based on the GABLS II intercomparison case.
  90 It it not the GABLS II intercomparison case for a number of reasons,
  91 but the two most important are (1) no advection, and (2) it does not
  92 specify surface fluxes for forcing.  The advection in GABLES is specified
  93 differently and this approach was not implemented in this version.
  94 Also, code to force with surface fluxes still needs to be developed.
  95 Again, please communicate with me if you would like to do this.
  96