Skip to content

Commit

Permalink
Based on suggestions from code reviewers of PR NOAA-EMC#730,
Browse files Browse the repository at this point in the history 
    adding more comments to explain the usage of
    corp_gust/hwllp_gust/oerr_gust for analysis of wind gust in 3drtma.
  • Loading branch information
@GangZhao-NOAA
GangZhao-NOAA committed Apr 5, 2024
1 parent b18a972 commit e1512de
Showing 1 changed file with 32 additions and 3 deletions.
35 changes: 32 additions & 3 deletions src/gsi/rapidrefresh_cldsurf_mod.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -198,8 +198,25 @@ module rapidrefresh_cldsurf_mod
! in hyrbid run is as similar as the analysis of howv in pure 3dvar run,
! the static BE of howv used in hybrid run needs to be tuned (inflated actually).
! corp_gust - namelist real, static BE of gust (standard error deviation)
! note: 1. initialised to be an arbitary negative value, in order to skip this
! negative value, instead to use value (3.0 m/s) set in subroutine
! berror_read_wgt_reg as default.
! 2. (3drtma only) if a user-specified value (e.g., 2.0 m/s) is preferred
! for corp_gust, in GSI namelist session "rapidrefresh_cldsurf",
! set "corp_gust=2.0,"
! hwllp_gust - namelist real, static BE de-correlation length scale of gust
! note: 1. initialised to be an arbitary negative value, in order to skip this
! negative value, instead to use value (same value for q) set in
! subroutine berror_read_wgt_reg as default
! 2. (3drtma only) if a user-specified value (e.g., 100 km) is preferred
! for hwllp_gust, in GSI namelist session "rapidrefresh_cldsurf",
! set "hwllp_gust=100000.0,"
! oerr_gust - namelist real, observation error of gust
! note: 1. initialised to be an arbitary negative value, in order to skip this
! negative value, instead to use value (1.0 m/s) set in read_prepbufr.f90
! 2. (3drtma only) if a user-specified value (e.g., 1.5 m/s ) is preferred
! for oerr_gust, in GSI namelist session "rapidrefresh_cldsurf",
! set "oerr_gust=1.5,"
! i_gust_3dda - integer, control the analysis of gust in 3D analysis (either var or hybrid)
! = 0 (gust-off: default) : no analysis of gust in 3D analysis.
! = 1 (gust-on) : if variable name "gust" is found in anavinfo,
Expand Down Expand Up @@ -458,9 +475,21 @@ subroutine init_rapidrefresh_cldsurf
corp_howv = 0.42_r_kind ! 0.42 meters (default)
hwllp_howv = 170000.0_r_kind ! 170,000.0 meters (170km as default for 3DRTMA, 50km is used in 2DRTMA)
i_howv_3dda = 0 ! no analysis of significant wave height (howv) in 3D analysis (default)
corp_gust = -3.00_r_kind ! default to use the value (3.0 m/s) set in m_berror_stats_reg.f90
hwllp_gust = -170000.0_r_kind ! default to use the value (same as q) set in m_berror_stats_reg.f90
oerr_gust = -1.0_r_kind ! default to use the value (1.0 m/s) set in read_prepbufr.f90
corp_gust = -1.50_r_kind ! initialised as negative & void to be skipped, in order to use
! the value (3.0 m/s) set in sub berror_read_wgt_reg (as default).
! If user-specified value is preferred, set it in session
! "rapidrefresh_cldsurf" of GSI namelist file

hwllp_gust = -90000.0_r_kind ! initialised as a value, in order to skip this negative value
! and to use the value (used for q) set in sub berror_read_wgt_reg.
! If user-specified value is preferred, set it in session
! "rapidrefresh_cldsurf" of GSI namelist file

oerr_gust = -2.5_r_kind ! initialised as a negative value, in order to skip this negative value
! and to use the value (1.0 m/s) set in read_prepbufr.f90
! If user-specified value is preferred, set it in session
! "rapidrefresh_cldsurf" of GSI namelist file

i_gust_3dda = 0 ! no analysis of wind gust (gust) in 3D analysis (default)

!-- searching for specific variable in state variable list (reading from anavinfo)
Expand Down

0 comments on commit e1512de

Please sign in to comment.