Changeset 4444 for LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_delta_eddington.h

Timestamp:

Feb 21, 2023, 3:26:41 PM (17 months ago)

Author:

idelkadi

Message:

Update of the ECRAD radiative code version implemented in the LMDZ model.
Upgrade to the :
​https://github.com/ecmwf/ecrad/trunk
Version svn : 749
UUID du dépôt : 44b0ca93-0ed8-356e-d663-ce57b7db7bff

File:

• LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_delta_eddington.h (modified) (2 diffs)

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_delta_eddington.h

@@ -1 @@
-! radiation_delta_eddington.h - Delta-Eddington scaling
+! radiation_delta_eddington.h - Delta-Eddington scaling -*- f90 -*-
 !
 ! (C) Copyright 2015- ECMWF.
@@ -92 +92 @@
 end subroutine delta_eddington_scat_od
 
+
+!---------------------------------------------------------------------
+! Revert delta-Eddington-scaled quantities in-place, back to their
+! original state
+elemental subroutine revert_delta_eddington(od, ssa, g)
+
+  use parkind1, only : jprb
+  
+  ! Total optical depth, single scattering albedo and asymmetry
+  ! factor
+  real(jprb), intent(inout) :: od, ssa, g
+  
+  ! Fraction of the phase function deemed to be in the forward lobe
+  ! and therefore treated as if it is not scattered at all
+  real(jprb) :: f
+  
+  g   = g / (1.0_jprb - g)
+  f   = g*g
+  ssa = ssa / (1.0_jprb - f + f*ssa);
+  od  = od / (1.0_jprb - ssa*f)
+  
+end subroutine revert_delta_eddington