From d1c1dcb68febe22f42c9a10a37ab05a099234e4a Mon Sep 17 00:00:00 2001 From: Sven Karsten Date: Thu, 16 Feb 2023 13:33:38 +0100 Subject: [PATCH] Apply albedo for direct and diffusive radiation in MOM code directly --- src/oasis_interface/oas_exchange_fields.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/oasis_interface/oas_exchange_fields.F90 b/src/oasis_interface/oas_exchange_fields.F90 index ae90b7f..3bd216c 100644 --- a/src/oasis_interface/oas_exchange_fields.F90 +++ b/src/oasis_interface/oas_exchange_fields.F90 @@ -196,7 +196,7 @@ INTEGER, INTENT(IN) :: dt_cpld jn = 31 !IF( nrcvinfo(jn) == OASIS_Rcv ) THEN ! shortwave radiation direct 31..36 do k = 1, size(Ice_boundary%sw_flux_vis_dir,3) - WHERE (maskt == 0) Ice_boundary%sw_flux_vis_dir(isc:iec,jsc:jec,k) = -frcv(isc:iec,jsc:jec,jn) + WHERE (maskt == 0) Ice_boundary%sw_flux_vis_dir(isc:iec,jsc:jec,k) = -frcv(isc:iec,jsc:jec,jn)*(1.0-Ice%albedo_vis_dir(isc:iec,jsc:jec,k)) jn = jn + 1 enddo !ENDIF @@ -204,7 +204,7 @@ INTEGER, INTENT(IN) :: dt_cpld jn = 37 !IF( nrcvinfo(jn) == OASIS_Rcv ) THEN ! shortware radiation diffusive 37) do k = 1, size(Ice_boundary%sw_flux_vis_dif,3) - WHERE (maskt == 0) Ice_boundary%sw_flux_vis_dif(isc:iec,jsc:jec,k) = -frcv(isc:iec,jsc:jec,jn) + WHERE (maskt == 0) Ice_boundary%sw_flux_vis_dif(isc:iec,jsc:jec,k) = -frcv(isc:iec,jsc:jec,jn)*(1.0-Ice%albedo_vis_dif(isc:iec,jsc:jec,k)) enddo !ENDIF