Skip to content

Commit a45e1b5

Browse files
mdtoyNOAAmdtoyNOAA
authored
Develop ugwp improvements (#600)
* Enhanced stability of drag_suite module and added diagnostics * Updated standard_names for do_gsl_drag_ls_bl, do_gsl_drag_ss, do_gsl_drag_tofd and ugwp_seq_update * Assumed-shape in gwdps.f * Added logical flag 'ugwp_seq_update' to gfs_physics_nml namelist * Modified calculation of dtfac in drag_suite.F90 and revised various long_name's of diagnostic variables * Removed passage of variables 'ak' and 'bk' to drag_suite -- these are no longer used * Removed unused variables from GFS_typedefs.F90 and added clarifying comments in ccpp-physics
1 parent c6646f3 commit a45e1b5

7 files changed

+693
-93
lines changed

‎ccpp/data/CCPP_typedefs.F90

+4-2
Original file line numberDiff line numberDiff line change
@@ -922,7 +922,8 @@ subroutine gfs_interstitial_create (Interstitial, IM, Model)
922922
allocate (Interstitial%zngw (IM) )
923923

924924
! CIRES UGWP v1
925-
if (Model%do_ugwp_v1) then
925+
if (Model%ldiag_ugwp .or. Model%do_ugwp_v0 .or. Model%do_ugwp_v0_nst_only &
926+
.or. Model%do_ugwp_v1) then
926927
allocate (Interstitial%dudt_ngw (IM,Model%levs))
927928
allocate (Interstitial%dvdt_ngw (IM,Model%levs))
928929
allocate (Interstitial%dtdt_ngw (IM,Model%levs))
@@ -1550,7 +1551,8 @@ subroutine gfs_interstitial_phys_reset (Interstitial, Model)
15501551
Interstitial%zngw = clear_val
15511552

15521553
! CIRES UGWP v1
1553-
if (Model%do_ugwp_v1) then
1554+
if (Model%ldiag_ugwp .or. Model%do_ugwp_v0 .or. Model%do_ugwp_v0_nst_only &
1555+
.or. Model%do_ugwp_v1) then
15541556
Interstitial%dudt_ngw = clear_val
15551557
Interstitial%dvdt_ngw = clear_val
15561558
Interstitial%dtdt_ngw = clear_val

‎ccpp/data/CCPP_typedefs.meta

+4-4
Original file line numberDiff line numberDiff line change
@@ -2324,31 +2324,31 @@
23242324
dimensions = (horizontal_loop_extent,vertical_layer_dimension)
23252325
type = real
23262326
kind = kind_phys
2327-
active = (flag_for_ugwp_version_1 .or. control_for_drag_suite_gravity_wave_drag==33 .or. control_for_drag_suite_gravity_wave_drag==22 .or. control_for_drag_suite_gravity_wave_drag==3 .or. control_for_drag_suite_gravity_wave_drag==2)
2327+
active = (flag_for_unified_gravity_wave_physics_diagnostics .or. flag_for_ugwp_version_0 .or. flag_for_ugwp_version_0_nonorographic_gwd .or. flag_for_ugwp_version_1)
23282328
[dvdt_ngw]
23292329
standard_name = tendency_of_y_wind_due_to_nonorographic_gravity_wave_drag
23302330
long_name = meridional wind tendency due to non-stationary GWs
23312331
units = m s-2
23322332
dimensions = (horizontal_loop_extent,vertical_layer_dimension)
23332333
type = real
23342334
kind = kind_phys
2335-
active = (flag_for_ugwp_version_1 .or. control_for_drag_suite_gravity_wave_drag==33 .or. control_for_drag_suite_gravity_wave_drag==22 .or. control_for_drag_suite_gravity_wave_drag==3 .or. control_for_drag_suite_gravity_wave_drag==2)
2335+
active = (flag_for_unified_gravity_wave_physics_diagnostics .or. flag_for_ugwp_version_0 .or. flag_for_ugwp_version_0_nonorographic_gwd .or. flag_for_ugwp_version_1)
23362336
[dtdt_ngw]
23372337
standard_name = tendency_of_air_temperature_due_to_nonorographic_gravity_wave_drag
23382338
long_name = air temperature tendency due to non-stationary GWs
23392339
units = K s-1
23402340
dimensions = (horizontal_loop_extent,vertical_layer_dimension)
23412341
type = real
23422342
kind = kind_phys
2343-
active = (flag_for_ugwp_version_1 .or. control_for_drag_suite_gravity_wave_drag==33 .or. control_for_drag_suite_gravity_wave_drag==22 .or. control_for_drag_suite_gravity_wave_drag==3 .or. control_for_drag_suite_gravity_wave_drag==2)
2343+
active = (flag_for_unified_gravity_wave_physics_diagnostics .or. flag_for_ugwp_version_0 .or. flag_for_ugwp_version_0_nonorographic_gwd .or. flag_for_ugwp_version_1)
23442344
[kdis_ngw]
23452345
standard_name = atmosphere_momentum_diffusivity_due_to_nonorographic_gravity_wave_drag
23462346
long_name = eddy mixing due to non-stationary GWs
23472347
units = m2 s-1
23482348
dimensions = (horizontal_loop_extent,vertical_layer_dimension)
23492349
type = real
23502350
kind = kind_phys
2351-
active = (flag_for_ugwp_version_1 .or. control_for_drag_suite_gravity_wave_drag==33 .or. control_for_drag_suite_gravity_wave_drag==22 .or. control_for_drag_suite_gravity_wave_drag==3 .or. control_for_drag_suite_gravity_wave_drag==2)
2351+
active = (flag_for_unified_gravity_wave_physics_diagnostics .or. flag_for_ugwp_version_0 .or. flag_for_ugwp_version_0_nonorographic_gwd .or. flag_for_ugwp_version_1)
23522352
[zlwb]
23532353
standard_name = height_of_low_level_wave_breaking
23542354
long_name = height of low level wave breaking

‎ccpp/data/GFS_typedefs.F90

+77-49
Original file line numberDiff line numberDiff line change
@@ -688,6 +688,7 @@ module GFS_typedefs
688688
!vay 2018 GW physics switches
689689

690690
logical :: ldiag_ugwp
691+
logical :: ugwp_seq_update ! flag to update winds between UGWP steps
691692
logical :: do_ugwp ! do mesoscale UGWP + TOFD + RF
692693
logical :: do_tofd ! tofd flag in gwdps.f
693694
logical :: do_gwd ! logical for gravity wave drag (gwd)
@@ -978,7 +979,7 @@ module GFS_typedefs
978979
logical :: do_ugwp_v0 !< flag for version 0 ugwp GWD
979980
logical :: do_ugwp_v0_orog_only !< flag for version 0 ugwp GWD (orographic drag only)
980981
logical :: do_ugwp_v0_nst_only !< flag for version 0 ugwp GWD (non-stationary GWD only)
981-
logical :: do_gsl_drag_ls_bl !< flag for GSL drag (large-scale GWD and blocking only)
982+
logical :: do_gsl_drag_ls_bl !< flag for GSL drag (mesoscale GWD and blocking only)
982983
logical :: do_gsl_drag_ss !< flag for GSL drag (small-scale GWD only)
983984
logical :: do_gsl_drag_tofd !< flag for GSL drag (turbulent orog form drag only)
984985
logical :: do_ugwp_v1 !< flag for version 1 ugwp GWD
@@ -1836,14 +1837,22 @@ module GFS_typedefs
18361837
real (kind=kind_phys), pointer :: dudt_ofd(:,:) => null() !<
18371838
real (kind=kind_phys), pointer :: dvdt_ofd(:,:) => null() !<
18381839

1839-
real (kind=kind_phys), pointer :: du_ogwcol(:) => null() !<
1840-
real (kind=kind_phys), pointer :: dv_ogwcol(:) => null() !<
1841-
real (kind=kind_phys), pointer :: du_oblcol(:) => null() !<
1842-
real (kind=kind_phys), pointer :: dv_oblcol(:) => null() !<
1843-
real (kind=kind_phys), pointer :: du_osscol(:) => null() !<
1844-
real (kind=kind_phys), pointer :: dv_osscol(:) => null() !<
1845-
real (kind=kind_phys), pointer :: du_ofdcol(:) => null() !<
1846-
real (kind=kind_phys), pointer :: dv_ofdcol(:) => null() !<
1840+
real (kind=kind_phys), pointer :: du_ogwcol(:) => null() !< instantaneous sfc u-momentum flux from OGW
1841+
real (kind=kind_phys), pointer :: dv_ogwcol(:) => null() !< instantaneous sfc v-momentum flux from OGW
1842+
real (kind=kind_phys), pointer :: du_oblcol(:) => null() !< instantaneous sfc u-momentum flux from blocking
1843+
real (kind=kind_phys), pointer :: dv_oblcol(:) => null() !< instantaneous sfc v-momentum flux from blocking
1844+
real (kind=kind_phys), pointer :: du_osscol(:) => null() !< instantaneous sfc u-momentum flux from SSGWD
1845+
real (kind=kind_phys), pointer :: dv_osscol(:) => null() !< instantaneous sfc v-momentum flux from SSGWD
1846+
real (kind=kind_phys), pointer :: du_ofdcol(:) => null() !< instantaneous sfc u-momentum flux from TOFD
1847+
real (kind=kind_phys), pointer :: dv_ofdcol(:) => null() !< instantaneous sfc v-momentum flux from TOFD
1848+
real (kind=kind_phys), pointer :: du3_ogwcol(:) => null() !< time-averaged sfc u-momentum flux from OGW
1849+
real (kind=kind_phys), pointer :: dv3_ogwcol(:) => null() !< time-averaged sfc v-momentum flux from OGW
1850+
real (kind=kind_phys), pointer :: du3_oblcol(:) => null() !< time-averaged sfc u-momentum flux from blocking
1851+
real (kind=kind_phys), pointer :: dv3_oblcol(:) => null() !< time-averaged sfc v-momentum flux from blocking
1852+
real (kind=kind_phys), pointer :: du3_osscol(:) => null() !< time-averaged sfc u-momentum flux from SSGWD
1853+
real (kind=kind_phys), pointer :: dv3_osscol(:) => null() !< time-averaged sfc v-momentum flux from SSGWD
1854+
real (kind=kind_phys), pointer :: du3_ofdcol(:) => null() !< time-averaged sfc u-momentum flux from TOFD
1855+
real (kind=kind_phys), pointer :: dv3_ofdcol(:) => null() !< time-averaged sfc v-momentum flux from TOFD
18471856
!
18481857
!---vay-2018 UGWP-diagnostics daily mean
18491858
!
@@ -1856,28 +1865,27 @@ module GFS_typedefs
18561865
real (kind=kind_phys), pointer :: dt3dt_pbl(:,:) => null() !< daily aver GFS_phys tend for Temp pbl
18571866
!
18581867
real (kind=kind_phys), pointer :: du3dt_ogw(:,:) => null() !< daily aver GFS_phys tend for WE-U OGW
1859-
real (kind=kind_phys), pointer :: dv3dt_ogw(:,:) => null() !< daily aver GFS_phys tend for SN-V OGW
1860-
real (kind=kind_phys), pointer :: dt3dt_ogw(:,:) => null() !< daily aver GFS_phys tend for Temp OGW
1868+
!
1869+
real (kind=kind_phys), pointer :: ldu3dt_ogw(:,:) => null() !< time aver GFS_phys tend for WE-U OGW
1870+
real (kind=kind_phys), pointer :: ldu3dt_obl(:,:) => null() !< time aver GFS_phys tend for WE-U OBL
1871+
real (kind=kind_phys), pointer :: ldu3dt_oss(:,:) => null() !< time aver GFS_phys tend for WE-U OSS
1872+
real (kind=kind_phys), pointer :: ldu3dt_ofd(:,:) => null() !< time aver GFS_phys tend for WE-U OFD
18611873
!
18621874
real (kind=kind_phys), pointer :: du3dt_mtb(:,:) => null() !< daily aver GFS_phys tend for WE-U MTB
1863-
real (kind=kind_phys), pointer :: dv3dt_mtb(:,:) => null() !< daily aver GFS_phys tend for SN-V MTB
1864-
real (kind=kind_phys), pointer :: dt3dt_mtb(:,:) => null() !< daily aver GFS_phys tend for Temp MTB
18651875
!
18661876
real (kind=kind_phys), pointer :: du3dt_tms(:,:) => null() !< daily aver GFS_phys tend for WE-U TMS
1867-
real (kind=kind_phys), pointer :: dv3dt_tms(:,:) => null() !< daily aver GFS_phys tend for SN-V TMS
1868-
real (kind=kind_phys), pointer :: dt3dt_tms(:,:) => null() !< daily aver GFS_phys tend for Temp TMS
18691877
!
18701878
real (kind=kind_phys), pointer :: du3dt_ngw(:,:) => null() !< daily aver GFS_phys tend for WE-U NGW
18711879
real (kind=kind_phys), pointer :: dv3dt_ngw(:,:) => null() !< daily aver GFS_phys tend for SN-V NGW
1872-
real (kind=kind_phys), pointer :: dt3dt_ngw(:,:) => null() !< daily aver GFS_phys tend for Temp NGW
18731880
!
1874-
real (kind=kind_phys), pointer :: du3dt_cgw(:,:) => null() !< daily aver GFS_phys tend for WE-U NGW
1875-
real (kind=kind_phys), pointer :: dv3dt_cgw(:,:) => null() !< daily aver GFS_phys tend for SN-V NGW
1876-
real (kind=kind_phys), pointer :: dt3dt_cgw(:,:) => null() !< daily aver GFS_phys tend for Temp NGW
1881+
real (kind=kind_phys), pointer :: dws3dt_ogw(:,:) => null() !< time aver GFS_phys tend for windspeed OGW
1882+
real (kind=kind_phys), pointer :: dws3dt_obl(:,:) => null() !< time aver GFS_phys tend for windspeed OBL
1883+
real (kind=kind_phys), pointer :: dws3dt_oss(:,:) => null() !< time aver GFS_phys tend for windspeed OSS
1884+
real (kind=kind_phys), pointer :: dws3dt_ofd(:,:) => null() !< time aver GFS_phys tend for windspeed OFD
18771885
!
1878-
real (kind=kind_phys), pointer :: du3dt_moist(:,:) => null() !< daily aver GFS_phys tend for WE-U MOIST
1879-
real (kind=kind_phys), pointer :: dv3dt_moist(:,:) => null() !< daily aver GFS_phys tend for SN-V MOIST
1880-
real (kind=kind_phys), pointer :: dt3dt_moist(:,:) => null() !< daily aver GFS_phys tend for Temp MOIST
1886+
real (kind=kind_phys), pointer :: ldu3dt_ngw(:,:) => null() !< time aver GFS_phys tend for u wind NGW
1887+
real (kind=kind_phys), pointer :: ldv3dt_ngw(:,:) => null() !< time aver GFS_phys tend for v wind NGW
1888+
real (kind=kind_phys), pointer :: ldt3dt_ngw(:,:) => null() !< time aver GFS_phys tend for temperature NGW
18811889
!
18821890
!--- Instantaneous UGWP-diagnostics 16-variables
18831891
! Diag%gwp_ax, Diag%gwp_axo, Diag%gwp_axc, Diag%gwp_axf, &
@@ -3165,16 +3173,17 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
31653173
logical :: do_ugwp_v0 = .true. !< flag for version 0 ugwp GWD
31663174
logical :: do_ugwp_v0_orog_only = .false. !< flag for version 0 ugwp GWD (orographic drag only)
31673175
logical :: do_ugwp_v0_nst_only = .false. !< flag for version 0 ugwp GWD (non-stationary GWD only)
3168-
logical :: do_gsl_drag_ls_bl = .false. !< flag for GSL drag (large-scale GWD and blocking only)
3176+
logical :: do_gsl_drag_ls_bl = .false. !< flag for GSL drag (mesoscale GWD and blocking only)
31693177
logical :: do_gsl_drag_ss = .false. !< flag for GSL drag (small-scale GWD only)
31703178
logical :: do_gsl_drag_tofd = .false. !< flag for GSL drag (turbulent orog form drag only)
31713179
logical :: do_ugwp_v1 = .false. !< flag for version 1 ugwp GWD
31723180
logical :: do_ugwp_v1_orog_only = .false. !< flag for version 1 ugwp GWD (orographic drag only)
31733181
logical :: do_ugwp_v1_w_gsldrag = .false. !< flag for version 1 ugwp GWD (orographic drag only)
31743182
!--- vay-2018
3175-
logical :: ldiag_ugwp = .false. !< flag for UGWP diag fields
3176-
logical :: do_ugwp = .false. !< flag do UGWP+RF
3177-
logical :: do_tofd = .false. !< flag do Turb oro Form Drag
3183+
logical :: ldiag_ugwp = .false. !< flag for UGWP diag fields
3184+
logical :: ugwp_seq_update = .false. !< flag for updating winds between UGWP steps
3185+
logical :: do_ugwp = .false. !< flag do UGWP+RF
3186+
logical :: do_tofd = .false. !< flag do Turb oro Form Drag
31783187

31793188
logical :: do_gwd = .false. !< flag for running gravity wave drag
31803189
logical :: do_cnvgwd = .false. !< flag for running conv gravity wave drag
@@ -3535,7 +3544,7 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
35353544
do_ugwp_v0_nst_only, &
35363545
do_gsl_drag_ls_bl, do_gsl_drag_ss, do_gsl_drag_tofd, &
35373546
do_ugwp_v1, do_ugwp_v1_orog_only, do_ugwp_v1_w_gsldrag, &
3538-
var_ric, coef_ric_l, coef_ric_s, hurr_pbl, &
3547+
ugwp_seq_update, var_ric, coef_ric_l, coef_ric_s, hurr_pbl, &
35393548
do_myjsfc, do_myjpbl, &
35403549
hwrf_samfdeep, hwrf_samfshal,progsigma, &
35413550
h2o_phys, pdfcld, shcnvcw, redrag, hybedmf, satmedmf, &
@@ -3723,6 +3732,7 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
37233732
!VAY-ugwp --- set some GW-related switches
37243733
!
37253734
Model%ldiag_ugwp = ldiag_ugwp
3735+
Model%ugwp_seq_update = ugwp_seq_update
37263736
Model%do_ugwp = do_ugwp
37273737
Model%do_tofd = do_tofd
37283738

@@ -6924,28 +6934,29 @@ subroutine diag_create (Diag, IM, Model)
69246934
allocate (Diag%dv3dt_pbl (IM,Model%levs) )
69256935
allocate (Diag%dt3dt_pbl (IM,Model%levs) )
69266936
allocate (Diag%du3dt_ogw (IM,Model%levs) )
6927-
allocate (Diag%dv3dt_ogw (IM,Model%levs) )
6928-
allocate (Diag%dt3dt_ogw (IM,Model%levs) )
69296937
allocate (Diag%du3dt_mtb (IM,Model%levs) )
6930-
allocate (Diag%dv3dt_mtb (IM,Model%levs) )
6931-
allocate (Diag%dt3dt_mtb (IM,Model%levs) )
69326938
allocate (Diag%du3dt_tms (IM,Model%levs) )
6933-
allocate (Diag%dv3dt_tms (IM,Model%levs) )
6934-
allocate (Diag%dt3dt_tms (IM,Model%levs) )
69356939
allocate (Diag%du3dt_ngw (IM,Model%levs) )
69366940
allocate (Diag%dv3dt_ngw (IM,Model%levs) )
6937-
allocate (Diag%dt3dt_ngw (IM,Model%levs) )
6938-
allocate (Diag%du3dt_cgw (IM,Model%levs) )
6939-
allocate (Diag%dv3dt_cgw (IM,Model%levs) )
6940-
allocate (Diag%dt3dt_moist (IM,Model%levs))
69416941
allocate (Diag%dudt_tot (IM,Model%levs) )
69426942
allocate (Diag%dvdt_tot (IM,Model%levs) )
69436943
allocate (Diag%dtdt_tot (IM,Model%levs) )
69446944
allocate (Diag%uav_ugwp (IM,Model%levs) )
69456945
allocate (Diag%tav_ugwp (IM,Model%levs) )
6946+
allocate (Diag%dws3dt_ogw (IM,Model%levs) )
6947+
allocate (Diag%dws3dt_obl (IM,Model%levs) )
6948+
allocate (Diag%dws3dt_oss (IM,Model%levs) )
6949+
allocate (Diag%dws3dt_ofd (IM,Model%levs) )
6950+
allocate (Diag%ldu3dt_ogw (IM,Model%levs) )
6951+
allocate (Diag%ldu3dt_obl (IM,Model%levs) )
6952+
allocate (Diag%ldu3dt_oss (IM,Model%levs) )
6953+
allocate (Diag%ldu3dt_ofd (IM,Model%levs) )
6954+
allocate (Diag%ldu3dt_ngw (IM,Model%levs) )
6955+
allocate (Diag%ldv3dt_ngw (IM,Model%levs) )
6956+
allocate (Diag%ldt3dt_ngw (IM,Model%levs) )
69466957
endif
69476958

6948-
if (Model%do_ugwp_v1 .or. Model%gwd_opt==33 .or. Model%gwd_opt==22) then
6959+
if (Model%do_ugwp_v1 .or. Model%ldiag_ugwp) then
69496960
allocate (Diag%dudt_ogw (IM,Model%levs))
69506961
allocate (Diag%dvdt_ogw (IM,Model%levs))
69516962
allocate (Diag%dudt_obl (IM,Model%levs))
@@ -6962,6 +6973,14 @@ subroutine diag_create (Diag, IM, Model)
69626973
allocate (Diag%dv_osscol (IM) )
69636974
allocate (Diag%du_ofdcol (IM) )
69646975
allocate (Diag%dv_ofdcol (IM) )
6976+
allocate (Diag%du3_ogwcol (IM) )
6977+
allocate (Diag%dv3_ogwcol (IM) )
6978+
allocate (Diag%du3_oblcol (IM) )
6979+
allocate (Diag%dv3_oblcol (IM) )
6980+
allocate (Diag%du3_osscol (IM) )
6981+
allocate (Diag%dv3_osscol (IM) )
6982+
allocate (Diag%du3_ofdcol (IM) )
6983+
allocate (Diag%dv3_ofdcol (IM) )
69656984
else
69666985
allocate (Diag%dudt_ogw (IM,Model%levs))
69676986
endif
@@ -7239,7 +7258,7 @@ subroutine diag_phys_zero (Diag, Model, linit, iauwindow_center)
72397258
Diag%dtdt_gw = zero
72407259
Diag%kdis_gw = zero
72417260

7242-
if (Model%do_ugwp_v1 .or. Model%gwd_opt==33 .or. Model%gwd_opt==22) then
7261+
if (Model%do_ugwp_v1 .or. Model%ldiag_ugwp) then
72437262
Diag%dudt_ogw = zero
72447263
Diag%dvdt_ogw = zero
72457264
Diag%dudt_obl = zero
@@ -7256,6 +7275,14 @@ subroutine diag_phys_zero (Diag, Model, linit, iauwindow_center)
72567275
Diag%dv_osscol = zero
72577276
Diag%du_ofdcol = zero
72587277
Diag%dv_ofdcol = zero
7278+
Diag%du3_ogwcol = zero
7279+
Diag%dv3_ogwcol = zero
7280+
Diag%du3_oblcol = zero
7281+
Diag%dv3_oblcol = zero
7282+
Diag%du3_osscol = zero
7283+
Diag%dv3_osscol = zero
7284+
Diag%du3_ofdcol = zero
7285+
Diag%dv3_ofdcol = zero
72597286
else
72607287
Diag%dudt_ogw = zero
72617288
end if
@@ -7265,25 +7292,26 @@ subroutine diag_phys_zero (Diag, Model, linit, iauwindow_center)
72657292
Diag%dv3dt_pbl = zero
72667293
Diag%dt3dt_pbl = zero
72677294
Diag%du3dt_ogw = zero
7268-
Diag%dv3dt_ogw = zero
7269-
Diag%dt3dt_ogw = zero
72707295
Diag%du3dt_mtb = zero
7271-
Diag%dv3dt_mtb = zero
7272-
Diag%dt3dt_mtb = zero
72737296
Diag%du3dt_tms = zero
7274-
Diag%dv3dt_tms = zero
7275-
Diag%dt3dt_tms = zero
72767297
Diag%du3dt_ngw = zero
72777298
Diag%dv3dt_ngw = zero
7278-
Diag%dt3dt_ngw = zero
7279-
Diag%du3dt_moist = zero
7280-
Diag%dv3dt_moist = zero
7281-
Diag%dt3dt_moist = zero
72827299
Diag%dudt_tot = zero
72837300
Diag%dvdt_tot = zero
72847301
Diag%dtdt_tot = zero
72857302
Diag%uav_ugwp = zero
72867303
Diag%tav_ugwp = zero
7304+
Diag%dws3dt_ogw = zero
7305+
Diag%dws3dt_obl = zero
7306+
Diag%dws3dt_oss = zero
7307+
Diag%dws3dt_ofd = zero
7308+
Diag%ldu3dt_ogw = zero
7309+
Diag%ldu3dt_obl = zero
7310+
Diag%ldu3dt_oss = zero
7311+
Diag%ldu3dt_ofd = zero
7312+
Diag%ldu3dt_ngw = zero
7313+
Diag%ldv3dt_ngw = zero
7314+
Diag%ldt3dt_ngw = zero
72877315
!COORDE
72887316
Diag%du3dt_dyn = zero
72897317
endif

0 commit comments

Comments
 (0)