Skip to content

Commit

Permalink
remove alpha from buoy grad
Browse files Browse the repository at this point in the history
  • Loading branch information
yairchn committed May 19, 2022
1 parent ba29a2b commit d3cf1a7
Show file tree
Hide file tree
Showing 7 changed files with 148 additions and 565 deletions.
4 changes: 0 additions & 4 deletions .buildkite/pipeline.yml
Original file line number Diff line number Diff line change
Expand Up @@ -153,10 +153,6 @@ steps:
command: "julia --color=yes --project=integration_tests integration_tests/3dBomex.jl"
artifact_paths: "integration_tests/output/3dBomex/*"

- label: ":earth_americas: 3D sphere bariclinic wave"
command: "julia --color=yes --project=integration_tests integration_tests/baroclinic_wave.jl"
artifact_paths: "integration_tests/output/baroclinic_wave/*"

- group: "Performance monitoring"
steps:

Expand Down
8 changes: 4 additions & 4 deletions driver/Cases.jl
Original file line number Diff line number Diff line change
Expand Up @@ -565,7 +565,7 @@ ForcingBase(case::TRMM_LBA, param_set::APS; kwargs...) =
function surface_ref_state(::TRMM_LBA, param_set::APS, namelist)
molmass_ratio = CPP.molmass_ratio(param_set)
Pg = 991.3 * 100 #Pressure at ground
Tg = 296.85 # surface values for reference state (RS) which outputs p rho0 alpha0
Tg = 296.85 # surface values for reference state (RS) which outputs p, ρ
pvg = TD.saturation_vapor_pressure(param_set, Tg, TD.Liquid())
qtg = (1 / molmass_ratio) * pvg / (Pg - pvg) #Total water mixing ratio at surface
return TD.PhaseEquil_pTq(param_set, Pg, Tg, qtg)
Expand Down Expand Up @@ -647,7 +647,7 @@ ForcingBase(case::ARM_SGP, param_set::APS; kwargs...) =

function surface_ref_state(::ARM_SGP, param_set::APS, namelist)
Pg = 970.0 * 100 #Pressure at ground
Tg = 299.0 # surface values for reference state (RS) which outputs p rho0 alpha0
Tg = 299.0 # surface values for reference state (RS) which outputs p, ρ
qtg = 15.2 / 1000 #Total water mixing ratio at surface
return TD.PhaseEquil_pTq(param_set, Pg, Tg, qtg)
end
Expand Down Expand Up @@ -731,7 +731,7 @@ ForcingBase(case::GATE_III, param_set::APS; kwargs...) = ForcingBase(get_forcing

function surface_ref_state(::GATE_III, param_set::APS, namelist)
Pg = 1013.0 * 100 #Pressure at ground
Tg = 299.184 # surface values for reference state (RS) which outputs p rho0 alpha0
Tg = 299.184 # surface values for reference state (RS) which outputs p, ρ
qtg = 16.5 / 1000 #Total water mixing ratio at surface
return TD.PhaseEquil_pTq(param_set, Pg, Tg, qtg)
end
Expand Down Expand Up @@ -1028,7 +1028,7 @@ end
function surface_ref_state(::SP, param_set::APS, namelist)
Pg = 1.0e5 #Pressure at ground
Tg = 300.0 #Temperature at ground
qtg = 1.0e-4 #Total water mixing ratio at TC. if set to 0, alpha0, rho0, p are NaN.
qtg = 1.0e-4 #Total water mixing ratio at TC. if set to 0, p, ρ are NaN.
return TD.PhaseEquil_pTq(param_set, Pg, Tg, qtg)
end

Expand Down
Loading

0 comments on commit d3cf1a7

Please sign in to comment.