Changeset 5105 for LMDZ6/branches/Amaury_dev/libf/misc/q_sat.f90

Timestamp:

Jul 23, 2024, 7:14:34 PM (8 weeks ago)

Author:

abarral

Message:

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F to *.f90
Fix gradsdef.h formatting
Remove unnecessary "RETURN" at the end of functions/subroutines

File:

• LMDZ6/branches/Amaury_dev/libf/misc/q_sat.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/misc/q_sat.F) (1 diff)

LMDZ6/branches/Amaury_dev/libf/misc/q_sat.f90

@@ -4 +4 @@
 
 
-      SUBROUTINE q_sat(np,temp,pres,qsat)
+SUBROUTINE q_sat(np,temp,pres,qsat)
 
-      IMPLICIT none
-!======================================================================
-! Autheur(s): Z.X. Li (LMD/CNRS)
-!  reecriture vectorisee par F. Hourdin.
-! Objet: calculer la vapeur d'eau saturante (formule Centre Euro.)
-!======================================================================
-! Arguments:
-! kelvin---input-R: temperature en Kelvin
-! millibar--input-R: pression en mb
+  IMPLICIT none
+  !======================================================================
+  ! Autheur(s): Z.X. Li (LMD/CNRS)
+  !  reecriture vectorisee par F. Hourdin.
+  ! Objet: calculer la vapeur d'eau saturante (formule Centre Euro.)
+  !======================================================================
+  ! Arguments:
+  ! kelvin---input-R: temperature en Kelvin
+  ! millibar--input-R: pression en mb
 
-! q_sat----output-R: vapeur d'eau saturante en kg/kg
-!======================================================================
+  ! q_sat----output-R: vapeur d'eau saturante en kg/kg
+  !======================================================================
 
-      integer np
-      REAL temp(np),pres(np),qsat(np)
+  integer :: np
+  REAL :: temp(np),pres(np),qsat(np)
 
-      REAL r2es
-      PARAMETER (r2es=611.14 *18.0153/28.9644)
+  REAL :: r2es
+  PARAMETER (r2es=611.14 *18.0153/28.9644)
 
-      REAL r3les, r3ies, r3es
-      PARAMETER (R3LES=17.269)
-      PARAMETER (R3IES=21.875)
+  REAL :: r3les, r3ies, r3es
+  PARAMETER (R3LES=17.269)
+  PARAMETER (R3IES=21.875)
 
-      REAL r4les, r4ies, r4es
-      PARAMETER (R4LES=35.86)
-      PARAMETER (R4IES=7.66)
+  REAL :: r4les, r4ies, r4es
+  PARAMETER (R4LES=35.86)
+  PARAMETER (R4IES=7.66)
 
-      REAL rtt
-      PARAMETER (rtt=273.16)
+  REAL :: rtt
+  PARAMETER (rtt=273.16)
 
-      REAL retv
-      PARAMETER (retv=28.9644/18.0153 - 1.0)
+  REAL :: retv
+  PARAMETER (retv=28.9644/18.0153 - 1.0)
 
-      real zqsat
-      integer ip
+  real :: zqsat
+  integer :: ip
 
-!     ------------------------------------------------------------------
+  ! ------------------------------------------------------------------
 
 
-      do ip=1,np
+  do ip=1,np
 
-!      write(*,*)'kelvin,millibar=',kelvin,millibar
-!       write(*,*)'temp,pres=',temp(ip),pres(ip)
+   ! write(*,*)'kelvin,millibar=',kelvin,millibar
+   !  write(*,*)'temp,pres=',temp(ip),pres(ip)
 
-         IF (temp(ip) <= rtt) THEN
-            r3es = r3ies
-            r4es = r4ies
-         ELSE
-            r3es = r3les
-            r4es = r4les
-         ENDIF
+     IF (temp(ip) <= rtt) THEN
+        r3es = r3ies
+        r4es = r4ies
+     ELSE
+        r3es = r3les
+        r4es = r4les
+     ENDIF
 
-         zqsat=r2es/pres(ip)*EXP(r3es*(temp(ip)-rtt)/(temp(ip)-r4es))
-         zqsat=MIN(0.5,ZQSAT)
-         zqsat=zqsat/(1.-retv *zqsat)
+     zqsat=r2es/pres(ip)*EXP(r3es*(temp(ip)-rtt)/(temp(ip)-r4es))
+     zqsat=MIN(0.5,ZQSAT)
+     zqsat=zqsat/(1.-retv *zqsat)
 
-         qsat(ip)= zqsat
-!      write(*,*)'qsat=',qsat(ip)
+     qsat(ip)= zqsat
+   ! write(*,*)'qsat=',qsat(ip)
 
-      enddo
+  enddo
 
-      RETURN
-      END
+
+END SUBROUTINE q_sat