Changeset 1992 for LMDZ5/trunk/libf/phylmd/tilft43.F90

Timestamp:

Mar 5, 2014, 2:19:12 PM (11 years ago)

Author:

lguez

Message:

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.

File:

• LMDZ5/trunk/libf/phylmd/tilft43.F90 (moved) (moved from LMDZ5/trunk/libf/phylmd/tilft43.F) (1 diff)

LMDZ5/trunk/libf/phylmd/tilft43.F90

@@ -1 @@
-!
+
 ! $Header$
-!
-        SUBROUTINE TLIFT43(P,T,Q,QS,GZ,ICB,NK,TVP,TPK,CLW,ND,NL,KK)
-        REAL GZ(ND),TPK(ND),CLW(ND),P(ND)
-        REAL T(ND),Q(ND),QS(ND),TVP(ND),LV0
-C
-C   ***   ASSIGN VALUES OF THERMODYNAMIC CONSTANTS     ***
-C
-c -- sb:
-c!      CPD=1005.7
-c!      CPV=1870.0
-c!      CL=4190.0
-c!      RV=461.5
-c!      RD=287.04
-c!      LV0=2.501E6
-c!      G=9.8
-c!      ROWL=1000.0
-c ajouts:
-#include "YOMCST.h"
-        CPD = RCPD
-        CPV = RCPV
-        CL = RCW
-        LV0 = RLVTT
-        G = RG
-        ROWL= RATM/100.
-        GRAVITY = RG !sb: Pr que gravite ne devienne pas humidite!
-C sb --
-C
-        CPVMCL=CL-CPV
-        EPS=RD/RV
-        EPSI=1./EPS
-C
-C   ***  CALCULATE CERTAIN PARCEL QUANTITIES, INCLUDING STATIC ENERGY   ***
-C
-        AH0=(CPD*(1.-Q(NK))+CL*Q(NK))*T(NK)+Q(NK)*(LV0-CPVMCL*(
-     1   T(NK)-273.15))+GZ(NK)
-        CPP=CPD*(1.-Q(NK))+Q(NK)*CPV
-        CPINV=1./CPP
-C
-        IF(KK.EQ.1)THEN
-C
-C   ***   CALCULATE LIFTED PARCEL QUANTITIES BELOW CLOUD BASE   ***
-C
-        DO 50 I=1,ICB-1
-         CLW(I)=0.0
-   50   CONTINUE
-        DO 100 I=NK,ICB-1
-         TPK(I)=T(NK)-(GZ(I)-GZ(NK))*CPINV
-         TVP(I)=TPK(I)*(1.+Q(NK)*EPSI)
-  100   CONTINUE
-        END IF
-C
-C    ***  FIND LIFTED PARCEL QUANTITIES ABOVE CLOUD BASE    ***
-C
-        NST=ICB
-        NSB=ICB
-        IF(KK.EQ.2)THEN  
-         NST=NL
-         NSB=ICB+1
-        END IF
-        DO 300 I=NSB,NST
-         TG=T(I)
-         QG=QS(I)
-         ALV=LV0-CPVMCL*(T(I)-273.15)
-         DO 200 J=1,2
-          S=CPD+ALV*ALV*QG/(RV*T(I)*T(I))
-          S=1./S
-          AHG=CPD*TG+(CL-CPD)*Q(NK)*T(I)+ALV*QG+GZ(I)
-          TG=TG+S*(AH0-AHG)
-          TG=MAX(TG,35.0)
-          TC=TG-273.15
-          DENOM=243.5+TC
-          IF(TC.GE.0.0)THEN  
-           ES=6.112*EXP(17.67*TC/DENOM)
-          ELSE  
-           ES=EXP(23.33086-6111.72784/TG+0.15215*LOG(TG))
-          END IF  
-          QG=EPS*ES/(P(I)-ES*(1.-EPS))
-  200    CONTINUE
-         ALV=LV0-CPVMCL*(T(I)-273.15)
-         TPK(I)=(AH0-(CL-CPD)*Q(NK)*T(I)-GZ(I)-ALV*QG)/CPD
-         CLW(I)=Q(NK)-QG
-         CLW(I)=MAX(0.0,CLW(I))
-         RG=QG/(1.-Q(NK))
-         TVP(I)=TPK(I)*(1.+RG*EPSI)
-  300   CONTINUE
 
-c -- sb:
-        RG = GRAVITY  ! RG redevient la gravite de YOMCST (sb)
-c sb --
+SUBROUTINE tlift43(p, t, q, qs, gz, icb, nk, tvp, tpk, clw, nd, nl, kk)
+  REAL gz(nd), tpk(nd), clw(nd), p(nd)
+  REAL t(nd), q(nd), qs(nd), tvp(nd), lv0
 
-        RETURN
-        END
+  ! ***   ASSIGN VALUES OF THERMODYNAMIC CONSTANTS     ***
 
+  ! -- sb:
+  ! !      CPD=1005.7
+  ! !      CPV=1870.0
+  ! !      CL=4190.0
+  ! !      RV=461.5
+  ! !      RD=287.04
+  ! !      LV0=2.501E6
+  ! !      G=9.8
+  ! !      ROWL=1000.0
+  ! ajouts:
+  include "YOMCST.h"
+  cpd = rcpd
+  cpv = rcpv
+  cl = rcw
+  lv0 = rlvtt
+  g = rg
+  rowl = ratm/100.
+  gravity = rg !sb: Pr que gravite ne devienne pas humidite!
+  ! sb --
+
+  cpvmcl = cl - cpv
+  eps = rd/rv
+  epsi = 1./eps
+
+  ! ***  CALCULATE CERTAIN PARCEL QUANTITIES, INCLUDING STATIC ENERGY   ***
+
+  ah0 = (cpd*(1.-q(nk))+cl*q(nk))*t(nk) + q(nk)*(lv0-cpvmcl*(t(nk)-273.15)) + &
+    gz(nk)
+  cpp = cpd*(1.-q(nk)) + q(nk)*cpv
+  cpinv = 1./cpp
+
+  IF (kk==1) THEN
+
+    ! ***   CALCULATE LIFTED PARCEL QUANTITIES BELOW CLOUD BASE   ***
+
+    DO i = 1, icb - 1
+      clw(i) = 0.0
+    END DO
+    DO i = nk, icb - 1
+      tpk(i) = t(nk) - (gz(i)-gz(nk))*cpinv
+      tvp(i) = tpk(i)*(1.+q(nk)*epsi)
+    END DO
+  END IF
+
+  ! ***  FIND LIFTED PARCEL QUANTITIES ABOVE CLOUD BASE    ***
+
+  nst = icb
+  nsb = icb
+  IF (kk==2) THEN
+    nst = nl
+    nsb = icb + 1
+  END IF
+  DO i = nsb, nst
+    tg = t(i)
+    qg = qs(i)
+    alv = lv0 - cpvmcl*(t(i)-273.15)
+    DO j = 1, 2
+      s = cpd + alv*alv*qg/(rv*t(i)*t(i))
+      s = 1./s
+      ahg = cpd*tg + (cl-cpd)*q(nk)*t(i) + alv*qg + gz(i)
+      tg = tg + s*(ah0-ahg)
+      tg = max(tg, 35.0)
+      tc = tg - 273.15
+      denom = 243.5 + tc
+      IF (tc>=0.0) THEN
+        es = 6.112*exp(17.67*tc/denom)
+      ELSE
+        es = exp(23.33086-6111.72784/tg+0.15215*log(tg))
+      END IF
+      qg = eps*es/(p(i)-es*(1.-eps))
+    END DO
+    alv = lv0 - cpvmcl*(t(i)-273.15)
+    tpk(i) = (ah0-(cl-cpd)*q(nk)*t(i)-gz(i)-alv*qg)/cpd
+    clw(i) = q(nk) - qg
+    clw(i) = max(0.0, clw(i))
+    rg = qg/(1.-q(nk))
+    tvp(i) = tpk(i)*(1.+rg*epsi)
+  END DO
+
+  ! -- sb:
+  rg = gravity ! RG redevient la gravite de YOMCST (sb)
+  ! sb --
+
+  RETURN
+END SUBROUTINE tlift43
+