Changeset 4072 for LMDZ6/trunk/libf/phylmd/calcratqs_multi_mod.F90

Timestamp:

Feb 1, 2022, 6:34:34 PM (2 years ago)

Author:

evignon

Message:

Mise à jour de la routine de nuages lscp
les principaux changements consistent en:

• des corrections de bug (déclaration et division

par zero) pour que la routine compile avec les modifications d'OB

• des restructurations pour que les fonctions de lscp_tools_mod

travaillent sur des tableaux et non des scalaires

• une séparation des coefficients d'évaporation de la neige et de la pluie

File:

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

LMDZ6/trunk/libf/phylmd/calcratqs_multi_mod.F90

@@ -230 +230 @@
 
 INTEGER :: i,k,nsrf
-REAL :: ratiom, qsat2m, dqsatdT
-REAL, DIMENSION(klon) :: xsi0
+REAL, DIMENSION(klon) :: xsi0, ratiom, qsat2m, dqsatdT
 REAL, DIMENSION (klon,klev) :: zlay
 
@@ -239 +238 @@
 !-------------------------------------------
 
-DO i=1,klon
     
-    ratiom=0.
-    xsi0(i)=0.
+    ratiom(:)=0.
+    xsi0(:)=0.
     
     DO nsrf=1,nbsrf
-    CALL CALC_QSAT_ECMWF(t2m(i,nsrf),q2m(i,nsrf),paprs(i,1),RTT,0,.false.,qsat2m,dqsatdT)
-    ratiom=ratiom+pctsrf(i,nsrf)*(q2m(i,nsrf)/qsat2m)
-    xsi0(i)=xsi0(i)+pctsrf(i,nsrf)*((q2m(i,nsrf)/qsat2m-ratiom)**2)
+    CALL CALC_QSAT_ECMWF(klon,t2m(:,nsrf),q2m(:,nsrf),paprs(:,1),RTT,0,.false.,qsat2m,dqsatdT)
+    ratiom(:)=ratiom(:)+pctsrf(:,nsrf)*(q2m(:,nsrf)/qsat2m(:))
+    xsi0(:)=xsi0(:)+pctsrf(:,nsrf)*((q2m(:,nsrf)/qsat2m(:)-ratiom(:))**2)
     END DO
     
-    xsi0(i)=sqrt(xsi0(i))/(ratiom+1E-6)
-END DO
+    xsi0(:)=sqrt(xsi0(:))/(ratiom(:)+1E-6)