Changeset 4717 for LMDZ6

Timestamp:

Oct 5, 2023, 5:33:27 PM (15 months ago)

Author:

Laurent Fairhead

Message:

There's a bug in icefrac_lsc_mod that did not appear previously because of the
conditions in which it was called. Modifying the lmdz_call_cloud_optics_prop routine waiting for the bug to be corrected

File:

• LMDZ6/trunk/libf/phylmd/lmdz_call_cloud_optics_prop.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_call_cloud_optics_prop.F90

@@ -20 +20 @@
   USE lmdz_cloud_optics_prop_ini , ONLY : ok_cdnc
   USE lmdz_cloud_optics_prop_ini , ONLY : bl95_b0, bl95_b1
-  USE lmdz_cloud_optics_prop_ini , ONLY : iflag_ice_thermo, ok_new_lscp
+  USE lmdz_cloud_optics_prop_ini , ONLY : iflag_ice_thermo, ok_new_lscp, iflag_t_glace
   USE lmdz_cloud_optics_prop, ONLY : cloud_optics_prop
   
@@ -95 +95 @@
   REAL :: dzfice(klon, klev)
 
-  DO k = 1, klev
-    IF (ok_new_lscp) THEN
-      CALL icefrac_lscp(klon,temp(:,k),iflag_ice_thermo,distcltop(:,k),temp_cltop(:,k), &
+  IF (iflag_t_glace .EQ. 0) THEN
+     icefrac_optics(:, :) = 0.
+  ELSE
+     DO k = 1, klev
+        IF (ok_new_lscp) THEN
+          CALL icefrac_lscp(klon,temp(:,k),iflag_ice_thermo,distcltop(:,k),temp_cltop(:,k), &
    &         icefrac_optics(:,k),dzfice(:,k))
-    ELSE
-      CALL icefrac_lsc(klon,temp(:,k),pplay(1:klon,k)/paprs(1:klon,1),icefrac_optics(:,k))
-    ENDIF
-  ENDDO
+        ELSE
+          CALL icefrac_lsc(klon,temp(:,k),pplay(1:klon,k)/paprs(1:klon,1),icefrac_optics(:,k))
+        ENDIF
+     ENDDO
+  ENDIF
+