Changeset 2529 for LMDZ5/trunk/libf/phylmd/rrtm/recmwf_aero.F90

Timestamp:

May 30, 2016, 10:46:40 PM (8 years ago)

Author:

oboucher

Message:

Reviewing the way swaero_diag is computed to correct an old bug that
makes topswad and co diagnostics wrong for the first timestep...
And introducing swaero_diag in rrtm/recmwf_aero so that double and
triple radiation calls when diagnostics are not needed to save time

File:

• LMDZ5/trunk/libf/phylmd/rrtm/recmwf_aero.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/rrtm/recmwf_aero.F90

@@ -161 +161 @@
 IMPLICIT NONE
 INCLUDE "clesphys.h"
-
 
 INTEGER(KIND=JPIM),INTENT(IN)    :: KPROMA 
@@ -306 +305 @@
 ! not ok_ade and not ok_aie = 1 to proceed
 ! therefore the cases have the following corresponding switches
-! 1 = not ok_ade or not ok_aie
-! 2 = ok_aie
-! 3 = ok_ade
+! 1 = not ok_ade and not ok_aie OR not ok_ade and ok_aie and swaero_diag OR ok_ade and not ok_aie and swaero_diag
+! 2 = not ok_ade and ok_aie OR ok_aie and ok_ade and swaero_diag
+! 3 = ok_ade and not ok_aie OR ok_aie and ok_ade and swaero_diag
 ! 4 = ok_ade and ok_aie
 ! 5 = no aerosol feedback wanted or no aerosol at all 
@@ -445 +444 @@
 IF (flag_aerosol .GT. 0 .OR. flag_aerosol_strat) THEN
 
-IF ( .not. ok_ade .or. .not. ok_aie ) THEN
+!--Case 1
+IF ( ( .not. ok_ade .AND. .not. ok_aie ) .OR.             & 
+   & ( .not. ok_ade .AND. ok_aie .AND. swaero_diag ) .OR. & 
+   & ( ok_ade .AND. .not. ok_aie .AND. swaero_diag ) ) THEN
 
 ! natural aerosols for direct and indirect effect 
@@ -483 +485 @@
 ENDIF 
 
-IF (ok_aie) THEN
+!--Case 2
+IF ( ( .not. ok_ade .AND. ok_aie ) .OR. & 
+   & ( ok_ade .AND. ok_aie .AND. swaero_diag ) ) THEN
 
 ! natural aerosols for direct indirect effect 
@@ -521 +525 @@
 ENDIF ! ok_aie      
 
-IF (ok_ade) THEN
+!--Case 3
+IF ( ( ok_ade .AND. .not. ok_aie ) .OR. &
+   & ( ok_ade .AND. ok_aie .AND. swaero_diag ) ) THEN
 
 ! direct effect of total aerosol activated
@@ -559 +565 @@
 ENDIF !-end ok_ade
 
+!--Case 4
 IF (ok_ade .and. ok_aie) THEN