Skip to content

Run-time option for climate GHG for radiation #1625

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
davegill merged 18 commits into from
Jan 19, 2022
Merged

Run-time option for climate GHG for radiation #1625

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
18 commits
Select commit Hold shift + click to select a range
79cb321
Run-time option for climate GHG for radiation
davegill Jan 6, 2022
a93ce74
RRTM no CLWRFGHG ifdef
davegill Jan 6, 2022
2f5b931
RRTMG LW and SW are b4b serial and OpenMP
davegill Jan 6, 2022
378fc95
Halfway with RRTMG fast
davegill Jan 6, 2022
15cf21b
Now CAM radiation
davegill Jan 6, 2022
aa60c19
Provide users with initial approx value of GHGs used
davegill Jan 7, 2022
b28695b
Remove ghg_input args from get_aerosol and vert_interpolate
davegill Jan 7, 2022
cb27dde
Add SSP to other GHG RCPs
davegill Jan 7, 2022
d13e52b
Per recommendation, SSP245 is default for GHG
davegill Jan 7, 2022
a57f0a0
Include ghg_input in the README.namelist file
davegill Jan 8, 2022
71fef8d
RRTM uses CO2 = 379E-6 22:56
davegill Jan 8, 2022
f708aa4
RRTM: make CO2 a time dependent function
davegill Jan 8, 2022
fa3a8f9
Put RRTM description back to the time-dependent CO2
davegill Jan 9, 2022
0f11722
ghg description 10:04
davegill Jan 12, 2022
bf9047f
ghg CO2 function only in RRTMs 11:11
davegill Jan 12, 2022
d467800
GHG error messages 12:56
davegill Jan 12, 2022
0c7e9ee
ghg allowable combinations 13:39
davegill Jan 12, 2022
6a357e3
Add ghg_input to the metadata output
davegill Jan 13, 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
16 changes: 14 additions & 2 deletions Makefile
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -546,9 +546,15 @@ em_real : wrf
ln -sf ../../run/CAM_AEROPT_DATA . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.RCP4.5 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.RCP6 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.RCP8.5 CAMtr_volume_mixing_ratio ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.RCP8.5 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.A1B . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.A2 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP119 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP126 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP245 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP245 CAMtr_volume_mixing_ratio ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP370 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP585 . ; \
ln -sf ../../run/CLM_ALB_ICE_DFS_DATA . ; \
ln -sf ../../run/CLM_ALB_ICE_DRC_DATA . ; \
ln -sf ../../run/CLM_ASM_ICE_DFS_DATA . ; \
Expand Down Expand Up @@ -611,9 +617,15 @@ em_real : wrf
ln -sf ../../run/CAM_AEROPT_DATA . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.RCP4.5 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.RCP6 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.RCP8.5 CAMtr_volume_mixing_ratio ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.RCP8.5 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.A1B . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.A2 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP119 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP126 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP245 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP245 CAMtr_volume_mixing_ratio ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP370 . ; \
ln -sf ../../run/CAMtr_volume_mixing_ratio.SSP585 . ; \
ln -sf ../../run/CLM_ALB_ICE_DFS_DATA . ; \
ln -sf ../../run/CLM_ALB_ICE_DRC_DATA . ; \
ln -sf ../../run/CLM_ASM_ICE_DFS_DATA . ; \
Expand Down
1 change: 1 addition & 0 deletions Registry/Registry.EM_COMMON
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2397,6 +2397,7 @@ rconfig integer compute_radar_ref derived 1 0
rconfig integer ra_lw_physics namelist,physics max_domains -1 rh "ra_lw_physics" "" ""
rconfig integer ra_sw_physics namelist,physics max_domains -1 rh "ra_sw_physics" "" ""
rconfig integer ra_sw_eclipse namelist,physics 1 0 rh "ra_sw_eclipse" "0/1 flag: 1=turn eclipse on" ""
rconfig integer ghg_input namelist,physics 1 1 rh "ghg_input" "for CAM, RRTM, RRTMG, RRTMG_fast: 0/1 flag: 0=constant (CO2 is a function of year for RRTM*); 1=time-varying GHG from CAMtr climate file"
rconfig real radt namelist,physics max_domains 0 h "RADT" "" ""
rconfig real naer namelist,physics max_domains 1e9 rh "NAER" "" ""
rconfig integer sf_sfclay_physics namelist,physics max_domains -1 rh "sf_sfclay_physics" "" ""
Expand Down
3 changes: 3 additions & 0 deletions clean
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -57,6 +57,9 @@ if ( "$arg" == '-a' || "$arg" == '-aa' ) then
*/CCN_ACTIVATE.BIN \
*/CAMtr_volume_mixing_ratio.RCP4.5 */CAMtr_volume_mixing_ratio.RCP6 */CAMtr_volume_mixing_ratio.RCP8.5 \
*/CAMtr_volume_mixing_ratio.A1B */CAMtr_volume_mixing_ratio.A2 */CAMtr_volume_mixing_ratio \
*/CAMtr_volume_mixing_ratio.SSP119 */CAMtr_volume_mixing_ratio.SSP126 \
*/CAMtr_volume_mixing_ratio.SSP245 */CAMtr_volume_mixing_ratio.SSP370 \
*/CAMtr_volume_mixing_ratio.SSP585 \
*/CLM_*DATA */RRTMG_LW_DATA */RRTMG_SW_DATA \
*/p3_lookup* */BROADBAND_CLOUD_GODDARD.bin \
*/ozone.formatted */ozone_lat.formatted */ozone_plev.formatted \
Expand Down
2 changes: 1 addition & 1 deletion dyn_em/module_first_rk_step_part1.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -438,7 +438,7 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
& ,progn=config_flags%progn &
#endif
& ,slope_rad=config_flags%slope_rad,topo_shading=config_flags%topo_shading &
& ,shadowmask=grid%shadowmask &
& ,shadowmask=grid%shadowmask,ghg_input=config_flags%ghg_input &
& ,ht=grid%ht,dx=grid%dx,dy=grid%dy,dx2d=grid%dx2d,area2d=grid%area2d &
& ,diffuse_frac=grid%diffuse_frac &
& ,obscur=grid%ECOBSC, mask=grid%ECMASK &
Expand Down
1 change: 1 addition & 0 deletions dyn_em/start_em.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1025,6 +1025,7 @@ SUBROUTINE start_domain_em ( grid, allowed_to_read &
grid%rublten,grid%rvblten,grid%rthblten, &
grid%rqvblten,grid%rqcblten,grid%rqiblten, &
grid%rthraten,grid%rthratenlw,grid%rthratensw, &
grid%this_is_an_ideal_run, &
!BSINGH - For WRFCuP scheme(11/12/2013)
grid%cupflag,grid%cldfra_cup,grid%cldfratend_cup, & !wig, 18-Sep-2006
grid%shall, & !wig, 18-Sep-2006
Expand Down
1 change: 1 addition & 0 deletions main/depend.common
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -516,6 +516,7 @@ module_physics_init.o : \
module_ra_cam.o \
$(PHYS_CU) $(PHYS_BL) \
module_ra_cam_support.o \
module_ra_clWRF_support.o \
module_ra_sw.o \
module_ra_gsfcsw.o \
module_ra_gfdleta.o \
Expand Down
38 changes: 37 additions & 1 deletion phys/module_physics_init.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,7 @@ MODULE module_physics_init
#if ( defined( DM_PARALLEL ) && ( ! defined( STUBMPI ) ) )
USE module_dm, ONLY : wrf_dm_max_real
#endif
USE module_ra_clWRF_support

! USE module_ssib_veg , ONLY : init_module_ssib_veg !fds (SSiB constants)
!Local data for CAM's MG MP scheme
Expand All @@ -39,6 +40,7 @@ SUBROUTINE phy_init ( id, config_flags, DT, restart, zfull, zhalf, &
RUBLTEN,RVBLTEN,RTHBLTEN, &
RQVBLTEN,RQCBLTEN,RQIBLTEN, &
RTHRATEN,RTHRATENLW,RTHRATENSW, &
this_is_an_ideal_run, &
#if ( EM_CORE == 1 )
!BSINGH - For WRFCuP scheme(11/12/2013)
cupflag,cldfra_cup,cldfratend_cup, & !wig, 18-Sep-2006
Expand Down Expand Up @@ -308,7 +310,7 @@ SUBROUTINE phy_init ( id, config_flags, DT, restart, zfull, zhalf, &
INTEGER , INTENT(OUT) ,OPTIONAL :: nyear
REAL , INTENT(OUT) ,OPTIONAL :: nday

LOGICAL, INTENT(IN) :: start_of_simulation
LOGICAL, INTENT(IN) :: start_of_simulation, this_is_an_ideal_run
REAL, INTENT(IN) :: DT, p_top, DX, DY
REAL, DIMENSION(ims:ime,jms:jme), INTENT(INOUT), OPTIONAL :: DX2D, AREA2D
LOGICAL, INTENT(IN) :: restart
Expand Down Expand Up @@ -910,6 +912,11 @@ SUBROUTINE phy_init ( id, config_flags, DT, restart, zfull, zhalf, &
INTEGER, OPTIONAL :: irr_ph,irr_freq
!-----------------------------------------------------------------

! Climate GHG file read for radiation
REAL(KIND=8) :: co2dum,n2odum,ch4dum,f11dum,f12dum
CHARACTER(LEN=8) :: name
!-----------------------------------------------------------------

#if ( EM_CORE == 1 )

! Compute 2d grid distance and 2d grid cell area. For use with
Expand All @@ -926,6 +933,35 @@ SUBROUTINE phy_init ( id, config_flags, DT, restart, zfull, zhalf, &
end if
#endif

name = " "
IF ( ( config_flags%ghg_input .EQ. 1 ) .AND. &
( .NOT. this_is_an_ideal_run ) ) THEN
IF ( config_flags%ra_lw_physics .EQ. RRTMSCHEME ) THEN
name = "RRTM"
ELSE IF ( config_flags%ra_lw_physics .EQ. CAMLWSCHEME ) THEN
name = "CAM"
ELSE IF ( config_flags%ra_lw_physics .EQ. RRTMG_LWSCHEME ) THEN
name = "RRTMG"
ELSE IF ( config_flags%ra_lw_physics .EQ. RRTMG_LWSCHEME_FAST ) THEN
name = "RRTMG"
END IF
CALL read_CAMgases(julyr,float(julday),.true.,TRIM(name),co2dum,n2odum,ch4dum,f11dum,f12dum)
WRITE(message,*) 'GHG annual values from CAM trace gas file'
CALL wrf_message(TRIM(message))
WRITE(message,*) 'Year = ',julyr,', Julian day = ',julday
CALL wrf_message(TRIM(message))
WRITE(message,*) 'CO2 = ',co2dum,' volume mixing ratio'
CALL wrf_message(TRIM(message))
WRITE(message,*) 'N2O = ',n2odum,' volume mixing ratio'
CALL wrf_message(TRIM(message))
WRITE(message,*) 'CH4 = ',ch4dum,' volume mixing ratio'
CALL wrf_message(TRIM(message))
WRITE(message,*) 'CFC11 = ',f11dum,' volume mixing ratio'
CALL wrf_message(TRIM(message))
WRITE(message,*) 'CFC12 = ',f12dum,' volume mixing ratio'
CALL wrf_message(TRIM(message))
END IF

aercu_opt=config_flags%aercu_opt !PSH/TWG 06/10/16
aercu_fct=config_flags%aercu_fct !PSH/TWG 06/10/16
sf_urban_physics=config_flags%sf_urban_physics
Expand Down
Loading