Changeset 4446 for LMDZ6/branches/Portage_acc/libf/phylmd/fisrtilp.F90

Timestamp:

Mar 1, 2023, 6:22:39 PM (16 months ago)

Author:

Laurent Fairhead

Message:

Merged trunk revisions from 4127 to 4443 (HEAD) into branch

File:

• LMDZ6/branches/Portage_acc/libf/phylmd/fisrtilp.F90 (modified) (3 diffs)

LMDZ6/branches/Portage_acc/libf/phylmd/fisrtilp.F90

@@ -89 +89 @@
   ! Coeffients de fraction lessivee : pour OFF-LINE
   !
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: pfrac_nucl
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: pfrac_1nucl
-  REAL, DIMENSION(klon,klev),      INTENT(OUT)  :: pfrac_impa
+  REAL, DIMENSION(klon,klev),      INTENT(INOUT)  :: pfrac_nucl
+  REAL, DIMENSION(klon,klev),      INTENT(INOUT)  :: pfrac_1nucl
+  REAL, DIMENSION(klon,klev),      INTENT(INOUT)  :: pfrac_impa
   !
   ! Fraction d'aerosols lessivee par impaction et par nucleation
@@ -1565 +1565 @@
         DO i=1, klon
             IF (zrflclr(i) + ziflclr(i) .GT. 0 ) THEN
-                znebprecipclr(i) = min(znebprecipclr(i), max(zrflclr(i)/(znebprecipclr(i)*rain_int_min), ziflclr(i)/(znebprecipclr(i)*rain_int_min)))
+               znebprecipclr(i) = min(znebprecipclr(i), max(zrflclr(i) &
+                    /(znebprecipclr(i)*rain_int_min), &
+                    ziflclr(i)/(znebprecipclr(i)*rain_int_min)))
             ELSE
                 znebprecipclr(i)=0.
@@ -1571 +1573 @@
 
             IF (zrflcld(i) + ziflcld(i) .GT. 0 ) THEN
-                znebprecipcld(i) = min(znebprecipcld(i), max(zrflcld(i)/(znebprecipcld(i)*rain_int_min), ziflcld(i)/(znebprecipcld(i)*rain_int_min)))
+               znebprecipcld(i) = min(znebprecipcld(i), &
+                    max(zrflcld(i)/(znebprecipcld(i)*rain_int_min), &
+                    ziflcld(i)/(znebprecipcld(i)*rain_int_min)))
             ELSE
                 znebprecipcld(i)=0.