Skip to content

Ufs/dev ugwp improvements #22

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

Merged
ChunxiZhang-NOAA merged 15 commits into from
Dec 22, 2022
Merged

Ufs/dev ugwp improvements #22

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
15 commits
Select commit Hold shift + click to select a range
fa06f28
Enhanced stability of drag_suite module and added diagnostics
mdtoyNOAA Nov 8, 2022
3ad8050
Merge branch 'ufs/dev' of https://github.com/ufs-community/ccpp-physi…
mdtoyNOAA Nov 8, 2022
7e84121
Merge branch 'ufs/dev' of https://github.com/ufs-community/ccpp-physi…
mdtoyNOAA Nov 14, 2022
0af5474
Merge branch 'ufs/dev' of https://github.com/ufs-community/ccpp-physi…
mdtoyNOAA Nov 21, 2022
1a9fcfe
Updated standard_names for do_gsl_drag_ls_bl, do_gsl_drag_ss and do_g…
mdtoyNOAA Nov 21, 2022
253fcea
Incorporated dustinswales suggestions in PR#22
mdtoyNOAA Nov 23, 2022
b1b2692
Merge branch 'ufs/dev' of https://github.com/ufs-community/ccpp-physi…
mdtoyNOAA Nov 23, 2022
86ff927
Merge branch 'ufs/dev' of https://github.com/ufs-community/ccpp-physi…
mdtoyNOAA Nov 29, 2022
0d488db
Make input assumed-shape
Dec 1, 2022
b82f0d3
Merge pull request #1 from dustinswales/hotfix_uwm1488
mdtoyNOAA Dec 2, 2022
4b572f6
Merge branch 'ufs/dev' of https://github.com/ufs-community/ccpp-physi…
mdtoyNOAA Dec 5, 2022
a23b876
Modified calculation of dtfac in drag_suite.F90 and revised various l…
mdtoyNOAA Dec 7, 2022
9277a29
Removed passage of variables 'ak' and 'bk' to drag_suite -- these are…
mdtoyNOAA Dec 7, 2022
5a091c8
Added clarifying remarks in drag_suite and unified_ugwp modules
mdtoyNOAA Dec 8, 2022
712bca2
Merge branch 'ufs/dev' of https://github.com/ufs-community/ccpp-physi…
mdtoyNOAA Dec 15, 2022
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
17 changes: 12 additions & 5 deletions physics/cires_ugwp.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -195,11 +195,13 @@ end subroutine cires_ugwp_finalize
! \section det_cires_ugwp CIRES UGWP V0 Scheme Detailed Algorithm
subroutine cires_ugwp_run(do_ugwp, me, master, im, levs, ntrac, dtp, kdt, lonr, &
oro, oro_uf, hprime, nmtvr, oc, theta, sigma, gamma, elvmax, clx, oa4, &
do_tofd, ldiag_ugwp, cdmbgwd, xlat, xlat_d, sinlat, coslat, area, &
ugrs, vgrs, tgrs, qgrs, prsi, prsl, prslk, phii, phil, &
do_tofd, ldiag_ugwp, cdmbgwd, xlat, xlat_d, sinlat, coslat, &
area, ugrs, vgrs, tgrs, qgrs, prsi, prsl, prslk, phii, phil, &
del, kpbl, dusfcg, dvsfcg, gw_dudt, gw_dvdt, gw_dtdt, gw_kdis, &
tau_tofd, tau_mtb, tau_ogw, tau_ngw, zmtb, zlwb, zogw, &
dudt_mtb,dudt_ogw, dudt_tms, du3dt_mtb, du3dt_ogw, du3dt_tms, &
dusfc_ms, dvsfc_ms, dusfc_bl, dvsfc_bl, &
dudt_ogw, dtauy2d_ms, dtaux2d_bl, dtauy2d_bl, &
dudt_mtb, dudt_tms, du3dt_mtb, du3dt_ogw, du3dt_tms, &
dudt, dvdt, dtdt, rdxzb, con_g, con_pi, con_cp, con_rd, con_rv, con_fvirt, &
con_omega, rain, ntke, q_tke, dqdt_tke, lprnt, ipr, &
dtend, dtidx, index_of_x_wind, index_of_y_wind, index_of_temperature, &
Expand Down Expand Up @@ -228,6 +230,9 @@ subroutine cires_ugwp_run(do_ugwp, me, master, im, levs, ntrac, dtp, kdt, lonr
real(kind=kind_phys), intent(out), dimension(:) :: tau_mtb, tau_ogw, tau_tofd, tau_ngw
real(kind=kind_phys), intent(out), dimension(:, :):: gw_dudt, gw_dvdt, gw_dtdt, gw_kdis
real(kind=kind_phys), intent(out), dimension(:, :):: dudt_mtb, dudt_ogw, dudt_tms
real(kind=kind_phys), intent(out), dimension(:) :: dusfc_ms, dvsfc_ms, dusfc_bl, dvsfc_bl
real(kind=kind_phys), intent(out), dimension(:, :) :: dtauy2d_ms
real(kind=kind_phys), intent(out), dimension(:, :) :: dtaux2d_bl, dtauy2d_bl

! dtend is only allocated if ldiag=.true.
real(kind=kind_phys), optional, intent(inout) :: dtend(:,:,:)
Expand Down Expand Up @@ -314,9 +319,11 @@ subroutine cires_ugwp_run(do_ugwp, me, master, im, levs, ntrac, dtp, kdt, lonr
ugrs, vgrs, tgrs, qgrs(:,:,1), &
kpbl, prsi, del, prsl, prslk, phii, phil, dtp, kdt, &
hprime, oc, oa4, clx, theta, sigma, gamma, &
elvmax, dusfcg, dvsfcg, &
elvmax, dusfcg, dvsfcg, dudt_ogw, dtauy2d_ms, &
dtaux2d_bl, dtauy2d_bl, dusfc_ms, dvsfc_ms, &
dusfc_bl, dvsfc_bl, &
con_g, con_cp, con_rd, con_rv, lonr, &
nmtvr, cdmbgwd, me, lprnt, ipr, rdxzb, &
nmtvr, cdmbgwd, me, lprnt, ipr, rdxzb, ldiag_ugwp, &
errmsg, errflg)
if (errflg/=0) return
endif
Expand Down
66 changes: 61 additions & 5 deletions physics/cires_ugwp.meta
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -596,11 +596,35 @@
type = real
kind = kind_phys
intent = out
[dudt_mtb]
standard_name = instantaneous_change_in_x_wind_due_to_mountain_blocking_drag
long_name = instantaneous change in x wind due to mountain blocking drag
units = m s-2
dimensions = (horizontal_loop_extent,vertical_layer_dimension)
[dusfc_ms]
standard_name = vertically_integrated_x_momentum_flux_due_to_mesoscale_orographic_gravity_wave_drag
long_name = integrated x momentum flux from mesoscale gwd
units = Pa
dimensions = (horizontal_loop_extent)
type = real
kind = kind_phys
intent = out
[dvsfc_ms]
standard_name = vertically_integrated_y_momentum_flux_due_to_mesoscale_orographic_gravity_wave_drag
long_name = integrated y momentum flux from mesoscale gwd
units = Pa
dimensions = (horizontal_loop_extent)
type = real
kind = kind_phys
intent = out
[dusfc_bl]
standard_name = vertically_integrated_x_momentum_flux_due_to_blocking_drag
long_name = integrated x momentum flux from blocking drag
units = Pa
dimensions = (horizontal_loop_extent)
type = real
kind = kind_phys
intent = out
[dvsfc_bl]
standard_name = vertically_integrated_y_momentum_flux_due_to_blocking_drag
long_name = integrated y momentum flux from blocking drag
units = Pa
dimensions = (horizontal_loop_extent)
type = real
kind = kind_phys
intent = out
Expand All @@ -612,6 +636,38 @@
type = real
kind = kind_phys
intent = out
[dtauy2d_ms]
Comment thread
ChunxiZhang-NOAA marked this conversation as resolved.
standard_name = tendency_of_y_wind_due_to_mesoscale_orographic_gravity_wave_drag
long_name = instantaneous change in y wind due to orographic gw drag
units = m s-2
dimensions = (horizontal_loop_extent,vertical_layer_dimension)
type = real
kind = kind_phys
intent = out
[dtaux2d_bl]
standard_name = tendency_of_x_wind_due_to_blocking_drag
long_name = x wind tendency from blocking drag
units = m s-2
dimensions = (horizontal_loop_extent,vertical_layer_dimension)
type = real
kind = kind_phys
intent = out
[dtauy2d_bl]
standard_name = tendency_of_y_wind_due_to_blocking_drag
long_name = y wind tendency from blocking drag
units = m s-2
dimensions = (horizontal_loop_extent,vertical_layer_dimension)
type = real
kind = kind_phys
intent = out
[dudt_mtb]
standard_name = instantaneous_change_in_x_wind_due_to_mountain_blocking_drag
long_name = instantaneous change in x wind due to mountain blocking drag
units = m s-2
dimensions = (horizontal_loop_extent,vertical_layer_dimension)
type = real
kind = kind_phys
intent = out
[dudt_tms]
standard_name = tendency_of_x_wind_due_to_turbulent_orographic_form_drag
long_name = instantaneous change in x wind due to TOFD
Expand Down
Loading