Registry/registry.hyb_coord

   1 #       Dry pressure, Pd
   2 #       Dry surface pressure = Pds
   3 #       Model top pressure = Pt
   4 #       Dry mass in column (base + perturbation), Pcb + Pc = Pds - Pt
   5 #       1d column weighting term, B: BF is full levels, BH is half levels
   6
   7 #       Total dry pressure
   8 #       Pd = BF ( Pds - Pt ) + ( eta - BF ) ( P0 - Pt ) + Pt
   9
  10 #       Note that when B (full levels, BF) is identically eta (full levels, ZNW),
  11 #       the dry pressure, Pd = the terrain following coordinate.
  12 #       Pd = eta ( Pds - Pt ) + Pt
  13
  14 #       Note that when B (full levels, BF) is zero,
  15 #       the dry pressure, Pd = isobaric.
  16 #       Pd = eta ( P0 - Pt ) + Pt
  17
  18 #       Hybrid coordinate: mu is still d(Pd)/d(eta)
  19 #       new MUT = d Pd / d eta = d BF / d eta * ( Pcb + Pc ) + ( 1 - d BF / d eta ) * ( P0 - Pt )
  20
  21 #       Define two columnar constants, function only of eta and other constants, specifically for the "mu" replacements.
  22 #       C1 = d BF / d eta
  23 #       C2 = ( 1 - d BF / d eta ) * ( P0 - Pt )
  24
  25 #       Total field
  26 #       new              MUT(i,k,j) = C1(k) * ( Pcb(i,j) + Pc(i,j) ) + C2(k)
  27
  28 #       Base-state, background field
  29 #       new              MUB(i,k,j) = C1(k) *   Pcb(i,j)             + C2(k)
  30
  31 #       Perturbation field
  32 #       new              MU(i,k,j)  = C1(k) *              Pc(i,j)
  33
  34 #       Define two columnar constants, function only of eta and other constants, specifically for the "pressure" replacements.
  35 #       C3 = BF
  36 #       C4 = ( eta - BF ) * ( P0 - Pt )
  37
  38 #       new dry pressure Pd(i,k,j)  = C3(k) * ( Pds(i,j) - Pt      ) + C4(k) + Pt
  39
  40
  41 #<Table> <Type>   <Sym>         <Dims>   <Use>           <NumTLev> <Stagger>        <IO>            <DNAME>                        <DESCRIP>                            <UNITS>
  42
  43 state    real      c1h             k     misc                1         -     i02rh0{22}{23}{24}      "C1H"       "half levels, c1h = d bf / d eta, using znw"        "Dimensionless"
  44 state    real      c2h             k     misc                1         -     i02rh0{22}{23}{24}      "C2H"       "half levels, c2h = (1-c1h)*(p0-pt)"                "Pa"
  45
  46 state    real      c1f             k     misc                1         Z     i02rh0{22}{23}{24}      "C1F"       "full levels, c1f = d bf / d eta, using znu"        "Dimensionless"
  47 state    real      c2f             k     misc                1         Z     i02rh0{22}{23}{24}      "C2F"       "full levels, c2f = (1-c1f)*(p0-pt)"                "Pa"
  48
  49 state    real      c3h             k     misc                1         -     i02rh0{22}{23}{24}      "C3H"       "half levels, c3h = bh"                             "Dimensionless"
  50 state    real      c4h             k     misc                1         -     i02rh0{22}{23}{24}      "C4H"       "half levels, c4h = (eta-bh)*(p0-pt), using znu"    "Pa"
  51
  52 state    real      c3f             k     misc                1         Z     i02rh0{22}{23}{24}      "C3F"       "full levels, c3f = bf"                             "Dimensionless"
  53 state    real      c4f             k     misc                1         Z     i02rh0{22}{23}{24}      "C4F"       "full levels, c4f = (eta-bf)*(p0-pt), using znw"    "Pa"
  54
  55 state    real      pcb            ij    dyn_em               1         -     irhdus                  "PCB"       "base state dry air mass in column"                 "Pa"
  56 state    real      pc             ijb   dyn_em               2         -     irhusdf=(bdy_interp:dt) "PC"        "perturbation dry air mass in column"               "Pa"
  57
  58
  59
  60
  61 #<Table> <Type>   <Sym>         <Nml>                    <NumEntry> <Default>                                                     <DESCRIP>                             <UNITS>
  62 rconfig  integer  hybrid_opt    namelist,dynamics            1          2     i0                  "HYBRID_OPT"   "0=Original WRF coordinate, 1=Terrain Following using hybrid formulation, 2=Klemp cubic form with etac"         "Flag"
  63 rconfig  real     etac          namelist,dynamics            1         0.2    i0                  "ETAC"         "znw(k) < etac, eta surfaces are isobaric, 0.2 is a good default"          "Pa/Pa"