Changeset 5144 for LMDZ6/branches/Amaury_dev/libf/phylmd/evappot.F90

Timestamp:

Jul 29, 2024, 11:01:04 PM (8 weeks ago)

Author:

abarral

Message:

Put YOMCST.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/evappot.F90 (modified) (1 diff)

LMDZ6/branches/Amaury_dev/libf/phylmd/evappot.F90

@@ -1 @@
-SUBROUTINE evappot(klon,nbsrf,ftsol,pplay,cdragh,  &
-      t_seri,q_seri,u_seri,v_seri,evap_pot)
-USE lmdz_YOETHF
-USE lmdz_fcttre, ONLY: foeew, foede, qsats, qsatl, dqsats, dqsatl, thermcep
+SUBROUTINE evappot(klon, nbsrf, ftsol, pplay, cdragh, &
+        t_seri, q_seri, u_seri, v_seri, evap_pot)
+  USE lmdz_yoethf
+  USE lmdz_fcttre, ONLY: foeew, foede, qsats, qsatl, dqsats, dqsatl, thermcep
+  USE lmdz_yomcst
 
-IMPLICIT NONE
+  IMPLICIT NONE
 
-INCLUDE "YOMCST.h"
+  INTEGER :: klon, nbsrf
+  REAL, DIMENSION(klon, nbsrf) :: ftsol, evap_pot
+  REAL, DIMENSION(klon) :: pplay, t_seri, wind, q_seri, u_seri, v_seri, cdragh
 
+  INTEGER :: nsrf, i
+  REAL, DIMENSION(klon, nbsrf) :: qsat_ftsol
+  REAL, DIMENSION(klon) :: rhos, norme_u
+  REAL :: t_coup
 
-INTEGER :: klon, nbsrf
-REAL, DIMENSION(klon,nbsrf) :: ftsol,evap_pot
-REAL, DIMENSION(klon) :: pplay,t_seri,wind,q_seri,u_seri,v_seri,cdragh
+  t_coup = 234.   ! Quelle horreur !!!!!
 
-INTEGER :: nsrf,i
-REAL, DIMENSION(klon,nbsrf) :: qsat_ftsol
-REAL, DIMENSION(klon) :: rhos, norme_u
-REAL :: t_coup
+  DO nsrf = 1, nbsrf
+    DO i = 1, klon
+      IF (ftsol(i, nsrf)<t_coup) THEN
+        qsat_ftsol(i, nsrf) = qsats(ftsol(i, nsrf)) / pplay(i)
+      ELSE
+        qsat_ftsol(i, nsrf) = qsatl(ftsol(i, nsrf)) / pplay(i)
+      ENDIF
+    ENDDO
+  ENDDO
+  ! ========================================================== c
+  ! Calcul de l'evaporation Potentielle
 
-      t_coup=234.   ! Quelle horreur !!!!!
-
-DO nsrf = 1, nbsrf
-   DO i = 1, klon
-      IF (ftsol(i,nsrf)<t_coup) THEN
-         qsat_ftsol(i,nsrf) = qsats(ftsol(i,nsrf))/pplay(i)
-      ELSE
-         qsat_ftsol(i,nsrf) = qsatl(ftsol(i,nsrf))/pplay(i)
-      ENDIF
-   ENDDO
-ENDDO
-! ========================================================== c
-! Calcul de l'evaporation Potentielle
-
-
-rhos(:) = pplay(:)/(RD*t_seri(:))
-norme_u(:)=1.+sqrt(u_seri(:)*u_seri(:)+v_seri(:)*v_seri(:))
-DO nsrf = 1, nbsrf
-  evap_pot(:,nsrf)=rhos(:)*cdragh(:)*norme_u(:)*(qsat_ftsol(:,nsrf)-q_seri(:))
-ENDDO
-RETURN
+  rhos(:) = pplay(:) / (RD * t_seri(:))
+  norme_u(:) = 1. + sqrt(u_seri(:) * u_seri(:) + v_seri(:) * v_seri(:))
+  DO nsrf = 1, nbsrf
+    evap_pot(:, nsrf) = rhos(:) * cdragh(:) * norme_u(:) * (qsat_ftsol(:, nsrf) - q_seri(:))
+  ENDDO
+  RETURN
 
 END