Skip to content

Radaer api ukca side - part two initialisation #45

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

Open
Alan J Hewitt (alanjhewitt) wants to merge 15 commits into
base: main
Choose a base branch
from
Open

Radaer api ukca side - part two initialisation #45

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
c013064
This could let radaer have a different i_mode_setup to ukca
alanjhewitt Aug 11, 2026
b55e95d
pipework to populate ukca_radaer structure
alanjhewitt Aug 11, 2026
b8840ee
Needed more UM repository functionality
alanjhewitt Aug 12, 2026
1054a38
Fixing build errors
alanjhewitt Aug 12, 2026
7e84066
More build errors
alanjhewitt Aug 12, 2026
d828455
didnt declare ereport
alanjhewitt Aug 14, 2026
7eca4d8
Needed by LFRic_apps
alanjhewitt Aug 14, 2026
5fed92b
Allocatable lists used by radaer kernel
alanjhewitt Aug 24, 2026
8e0f14c
call new ukca module
alanjhewitt Aug 24, 2026
d98da23
Replace hard coded settings with calculated values
alanjhewitt Aug 25, 2026
fae9033
tidy
alanjhewitt Aug 25, 2026
3d281a4
l_sustrat is a namelist option now
alanjhewitt Aug 25, 2026
d30ba8e
pointers were looking in wrong structure
alanjhewitt Aug 25, 2026
ff8947f
Dont need pointers here
alanjhewitt Aug 25, 2026
558a140
Applied code styling
alanjhewitt Aug 25, 2026
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
1 change: 1 addition & 0 deletions src/control/core/interface/ukca_config_specification_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -655,6 +655,7 @@ MODULE ukca_config_specification_mod

TYPE(glomap_variables_type), SAVE, TARGET :: glomap_variables
TYPE(glomap_variables_type), SAVE, TARGET :: glomap_variables_climatology
TYPE(glomap_variables_type), SAVE, TARGET :: glomap_variables_radaer

! ---------------------------------------------------------------------------
! -- Templates for parent callback procedures to be used in UKCA --
Expand Down
84 changes: 29 additions & 55 deletions src/control/core/top_level/ukca_radaer_lfric_api_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -55,6 +55,8 @@ SUBROUTINE ukca_radaer_lfric_interface( &
ukca_modal_wtv, &
! Logical to describe orientation
l_inverted, &
! Logical to account for optical properties of sulphate in atmosphere
l_sustrat, &
! Control option for prescribed single scattering albedo array
i_ukca_radaer_prescribe_ssa, &
! Model level of the tropopause (input)
Expand Down Expand Up @@ -97,6 +99,7 @@ SUBROUTINE ukca_radaer_lfric_interface( &
ip_ukca_mode_accum, &
ip_ukca_mode_coarse, &
i_ukca_bc_tuned, &
i_ukca_tune_bc_off, &
cp_su, cp_bc, cp_oc, &
cp_cl, cp_du, cp_so, &
cp_no3, cp_nn, cp_nh4, &
Expand All @@ -106,6 +109,8 @@ SUBROUTINE ukca_radaer_lfric_interface( &

USE ukca_radaer_compute_aod_mod, ONLY: ukca_radaer_compute_aod

USE ukca_radaer_lfric_saved_mod, ONLY: ukca_radaer_lfric

USE ukca_radaer_prepare_mod, ONLY: ukca_radaer_prepare

USE parkind1, ONLY: jpim, jprb
Expand Down Expand Up @@ -158,6 +163,9 @@ SUBROUTINE ukca_radaer_lfric_interface( &
! Logical to describe orientation
LOGICAL, INTENT(IN) :: l_inverted

! Logical to account for optical properties of sulphate in atmosphere
LOGICAL, INTENT(IN) :: l_sustrat

! When > 0, use a prescribed single scattering albedo field
INTEGER, INTENT(IN) :: i_ukca_radaer_prescribe_ssa

Expand Down Expand Up @@ -196,33 +204,33 @@ SUBROUTINE ukca_radaer_lfric_interface( &

! Modal mass-mixing ratios
REAL, INTENT(IN OUT) :: ukca_mode_mix_ratio( npd_profile, npd_layer, &
n_radaer_mode )
n_radaer_mode )

! Band-averaged modal optical properties
REAL, INTENT(IN OUT) :: aer_lw_absorption( npd_profile, npd_layer, &
n_radaer_mode, n_lw_band )
n_radaer_mode, n_lw_band )

REAL, INTENT(IN OUT) :: aer_sw_absorption( npd_profile, npd_layer, &
n_radaer_mode, n_sw_band )
n_radaer_mode, n_sw_band )

REAL, INTENT(IN OUT) :: aer_lw_scattering( npd_profile, npd_layer, &
n_radaer_mode, n_lw_band )
n_radaer_mode, n_lw_band )

REAL, INTENT(IN OUT) :: aer_sw_scattering( npd_profile, npd_layer, &
n_radaer_mode, n_sw_band )
n_radaer_mode, n_sw_band )

REAL, INTENT(IN OUT) :: aer_lw_asymmetry( npd_profile, npd_layer, &
n_radaer_mode, n_lw_band )
n_radaer_mode, n_lw_band )

REAL, INTENT(IN OUT) :: aer_sw_asymmetry( npd_profile, npd_layer, &
n_radaer_mode, n_sw_band )
n_radaer_mode, n_sw_band )

! Aerosol Optical Depth diagnostics
REAL, INTENT(IN OUT) :: aod_ukca_all_modes( npd_profile, npd_ukca_aod_wavel, &
n_ukca_mode )
n_ukca_mode )

REAL, INTENT(IN OUT) :: aaod_ukca_all_modes(npd_profile, npd_ukca_aod_wavel, &
n_ukca_mode )
n_ukca_mode )

! Local variables

Expand All @@ -239,18 +247,17 @@ SUBROUTINE ukca_radaer_lfric_interface( &

! -----------------------------------------------------------------

LOGICAL, PARAMETER :: l_nitrate = .FALSE. ! Make this a namelist option later
LOGICAL, PARAMETER :: l_sustrat = .TRUE. ! Make this a namelist option later
! l_sustrat=.true. for ga9
! Make this a namelist option later
LOGICAL, PARAMETER :: l_nitrate = .FALSE.

LOGICAL, PARAMETER :: l_cornarrow_ins = .FALSE.
! Make this a namelist option later
LOGICAL, PARAMETER :: l_cornarrow_ins = .FALSE.

! -----------------------------------------------------------------

! Maxwell-Garnett mixing approach logical control switches
INTEGER, PARAMETER :: i_ukca_tune_bc = i_ukca_bc_tuned
INTEGER, PARAMETER :: i_glomap_clim_tune_bc = 0 ! No tuning
INTEGER, PARAMETER :: i_glomap_clim_tune_bc = i_ukca_tune_bc_off

! Spectral information
INTEGER, PARAMETER :: ip_infra_red = 2
Expand All @@ -259,16 +266,13 @@ SUBROUTINE ukca_radaer_lfric_interface( &
LOGICAL, PARAMETER :: soluble_wanted = .TRUE.
LOGICAL, PARAMETER :: soluble_unwanted = .FALSE.

! Local copies from ukca_radaer_lfric
INTEGER :: i_cpnt_index( ncp_max, nmodes )

INTEGER :: i_cpnt_type( ncp_max_x_nmodes )
INTEGER :: i_mode_type( nmodes )

INTEGER :: n_cpnt_in_mode( nmodes )

LOGICAL :: l_soluble( nmodes )

INTEGER :: i_cpnt_type( ncp_max_x_nmodes )

! -----------------------------------------------------------------

INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0
Expand All @@ -280,42 +284,12 @@ SUBROUTINE ukca_radaer_lfric_interface( &

!-----------------------------------------------------------------------


! Note that these hard coded values are calculated in UM module
! ukca_radaer_init-ukca.F90
! Arrays with modes start on 2nd mode ( Aitkin Solvent )
! Radaer expects data to be structured in this way.

! No nucleation mode
l_soluble(1:nmodes) = [.TRUE., .TRUE., .TRUE., .FALSE., &
.FALSE.,.FALSE.,.FALSE.,.FALSE.]

! No nucleation mode
n_cpnt_in_mode(1:nmodes) = [ 3, 5, 5, 2, 1, 1, -1, -1 ]

! No nucleation mode
i_mode_type(1:nmodes) = [ 1, 2, 3, 1, 2, 3, -1, -1 ]

! No nucleation mode
i_cpnt_index(cp_su, 1:nmodes)=[ 1, 4, 9, 14, 16, 17, -1, -1 ]
i_cpnt_index(cp_bc, 1:nmodes)=[ 2, 5, 10, 15, -1, -1, -1, -1 ]
i_cpnt_index(cp_oc, 1:nmodes)=[ 3, 6, 11, -1, -1, -1, -1, -1 ]
i_cpnt_index(cp_cl, 1:nmodes)=[ -1, 7, 12, -1, -1, -1, -1, -1 ]
i_cpnt_index(cp_du, 1:nmodes)=[ -1, 8, 13, -1, -1, -1, -1, -1 ]
i_cpnt_index(cp_so, 1:nmodes)=[ -1, -1, -1, -1, -1, -1, -1, -1 ]
i_cpnt_index(cp_no3,1:nmodes)=[ -1, -1, -1, -1, -1, -1, -1, -1 ]
i_cpnt_index(cp_nn, 1:nmodes)=[ -1, -1, -1, -1, -1, -1, -1, -1 ]
i_cpnt_index(cp_nh4,1:nmodes)=[ -1, -1, -1, -1, -1, -1, -1, -1 ]

i_cpnt_type(1:ncp_max_x_nmodes) = &
[ 1, 2, 3, 1, 2, 3, 4, 5, 1, &
2, 3, 4, 5, 2, 3, 5, 5, -1, &
-1, -1, -1, -1, -1, -1, -1, -1, -1, &
-1, -1, -1, -1, -1, -1, -1, -1, -1, &
-1, -1, -1, -1, -1, -1, -1, -1, -1, &
-1, -1, -1, -1, -1, -1, -1, -1, -1, &
-1, -1, -1, -1, -1, -1, -1, -1, -1, &
-1, -1, -1, -1, -1, -1, -1, -1, -1 ]
! Local copies to make the code easier to read
i_cpnt_index = ukca_radaer_lfric%i_cpnt_index
i_cpnt_type = ukca_radaer_lfric%i_cpnt_type
i_mode_type = ukca_radaer_lfric%i_mode_type
n_cpnt_in_mode = ukca_radaer_lfric%n_cpnt_in_mode
l_soluble = ukca_radaer_lfric%l_soluble

!----------------------------------------------------------------------

Expand Down
Loading
Loading