Skip to content

Hybrid vertical coordinate for ARW #5

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
davegill wants to merge 9 commits into from
Closed

Hybrid vertical coordinate for ARW #5

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
dbb4ebb
Branch for ARW Hybrid Coordinate, root=9429
davegill Jul 5, 2016
2410470
Branch for ARW Hybrid Coordinate, root=9429
davegill Jul 5, 2016
9232e6c
Initial branch for the ARW hybrid coordinate modifications.
davegill Jul 5, 2016
232dbf4
Add the hybrid vertical coordinate to the DFI option. This
davegill Jul 19, 2016
dc407db
Opps, forgot to swap a "k" for "kk". This is in the total
davegill Jul 19, 2016
e386867
Geez, how many of these (k) arrays will I mess up? Also,
davegill Jul 19, 2016
9e9d5bd
Diffusion mods + MP Morrison + tweak DFI
davegill Sep 15, 2016
fbc5b48
Fix the big step calc_ww_cp (k) that should be (k-1)
davegill Sep 15, 2016
23d7b5e
Throw away set of merges for hybrid testing
davegill Sep 16, 2016
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions Registry/Registry.EM
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -18,6 +18,7 @@ include registry.afwa
include registry.sbm
include registry.elec
include registry.bdy_perturb
include registry.hyb_coord

# added to output 5 for ESMF
state real landmask ij misc 1 - i0125rh05d=(interp_fcnm_imask)u=(copy_fcnm) "LANDMASK" "LAND MASK (1 FOR LAND, 0 FOR WATER)" ""
Expand Down
56 changes: 56 additions & 0 deletions Registry/registry.hyb_coord
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,56 @@
# Dry pressure, Pd
# Dry surface pressure = Pds
# Model top pressure = Pt
# Mass in column, Pc = Pds - Pt

# Total dry pressure
# Pd = BF ( Pds - Pt ) + ( eta - BF ) ( P0 - Pt ) + Pt

# Hybrid coordinate: mu is still d(Pd)/d(eta)
# new MUT = d Pd / d eta = d BF / d eta * ( Pcb + Pc ) + ( 1 - d BF / d eta ) * ( P0 - Pt )

# Define two columnar constants, function only of eta and other constants, specifically for the "mu" replacements.
# C1 = d BF / d eta
# C2 = ( 1 - d BF / d eta ) * ( P0 - Pt )

# Total field
# new MUT(i,k,j) = C1(k) * ( Pcb(i,j) + Pc(i,j) ) + C2(k)

# Base-state, background field
# new MUB(i,k,j) = C1(k) * Pcb(i,j) + C2(k)

# Perturbation field
# new MU(i,k,j) = C1(k) * Pc(i,j)

# Define two columnar constants, function only of eta and other constants, specifically for the "pressure" replacements.
# C3 = BF
# C4 = ( eta - BF ) * ( P0 - Pt )

# new dry pressure Pd(i,k,j) = C3(k) * ( Pds(i,j) - Pt ) + C4(k) + Pt


#<Table> <Type> <Sym> <Dims> <Use> <NumTLev> <Stagger> <IO> <DNAME> <DESCRIP> <UNITS>

state real bf k misc 1 Z i02rh0{22}{23}{24} "BF" "full levels, bf=0 => isobaric; bf=znw => sigma" "Dimensionless"
state real c1h k misc 1 - i02rh0{22}{23}{24} "C1H" "half levels, c1h = d bf / d eta, using znw" "Dimensionless"
state real c2h k misc 1 - i02rh0{22}{23}{24} "C2H" "half levels, c2h = (1-c1h)*(p0-pt)" "Pa"

state real bh k misc 1 - i02rh0{22}{23}{24} "BH" "half levels, bh=0 => isobaric; bh=znu => sigma" "Dimensionless"
state real c1f k misc 1 Z i02rh0{22}{23}{24} "C1F" "full levels, c1f = d bf / d eta, using znu" "Dimensionless"
state real c2f k misc 1 Z i02rh0{22}{23}{24} "C2F" "full levels, c2f = (1-c1f)*(p0-pt)" "Pa"

state real c3h k misc 1 - i02rh0{22}{23}{24} "C3H" "half levels, c3h = bh" "Dimensionless"
state real c4h k misc 1 - i02rh0{22}{23}{24} "C4H" "half levels, c4h = (eta-bh)*(p0-pt)+pt, using znu" "Pa"

state real c3f k misc 1 Z i02rh0{22}{23}{24} "C3F" "full levels, c3f = bf" "Dimensionless"
state real c4f k misc 1 Z i02rh0{22}{23}{24} "C4F" "full levels, c4f = (eta-bf)*(p0-pt)+pt, using znw" "Pa"

state real pcb ij dyn_em 1 - irhdus "PCB" "base state dry air mass in column" "Pa"
state real pc ijb dyn_em 2 - irhusdf=(bdy_interp:dt) "PC" "perturbation dry air mass in column" "Pa"




#<Table> <Type> <Sym> <Nml> <NumEntry> <Default> <DESCRIP> <UNITS>
rconfig integer hybrid_opt namelist,dynamics 1 0 i0 "HYBRID_OPT" "0=Original WRF coordinate, 1=Terrain Following using hybrid formulation, 2=Klemp cubic form with etac" "Flag"
rconfig real etac namelist,dynamics 1 -1. i0 "ETAC" "znw(k) < etac, eta surfaces are isobaric" "Dimensionless"
Loading