Changeset 1991 for trunk

Timestamp:

Aug 28, 2018, 5:35:44 PM (6 years ago)

Author:

jleconte

Message:

forgot to commit the file gfluxv

File:

• trunk/LMDZ.GENERIC/libf/phystd/gfluxv.F (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/phystd/gfluxv.F

@@ -45 +45 @@
 !!      PARAMETER (NLP=101) ! MUST BE LARGER THAN NLEVEL
 
-      REAL*8 EM, EP
+      REAL*8 EM, EP, EXPTRM
       REAL*8 W0(L_NLAYRAD), COSBAR(L_NLAYRAD), DTAU(L_NLAYRAD)
       REAL*8 TAU(L_NLEVRAD), WDEL(L_NLAYRAD), CDEL(L_NLAYRAD)
@@ -54 +54 @@
       REAL*8 G3(L_NLAYRAD), GAMA(L_NLAYRAD),CP(L_NLAYRAD),CM(L_NLAYRAD)
       REAL*8 CPM1(L_NLAYRAD),CMM1(L_NLAYRAD), E1(L_NLAYRAD)
-      REAL*8 E2(L_NLAYRAD),E3(L_NLAYRAD),E4(L_NLAYRAD),EXPTRM(L_NLAYRAD)
+      REAL*8 E2(L_NLAYRAD),E3(L_NLAYRAD),E4(L_NLAYRAD)
       REAL*8 FLUXUP, FLUXDN
       REAL*8 FACTOR, TAUCUMIN(L_LEVELS), TAUCUM(L_LEVELS)
@@ -184 +184 @@
 
       DO L=1,L_NLAYRAD
-        EXPTRM(L) = MIN(TAUMAX,LAMDA(L)*DTAU(L))  ! CLIPPED EXPONENTIAL
-        EP = EXP(EXPTRM(L))
+        EXPTRM = MIN(TAUMAX,LAMDA(L)*DTAU(L))  ! CLIPPED EXPONENTIAL
+        EP = EXP(EXPTRM)
 
         EM        = 1.0/EP
@@ -200 +200 @@
  
       DO L=1,L_NLAYRAD-1
-        EXPTRM(L) = MIN(TAUMAX,LAMDA(L)*(TAUCUM(2*L+1)-TAUCUM(2*L)))
- 
-        EP = EXP(EXPTRM(L))
+        EXPTRM = MIN(TAUMAX,LAMDA(L)*(TAUCUM(2*L+1)-TAUCUM(2*L)))
+ 
+        EP = EXP(EXPTRM)
 
         EM    = 1.0/EP
@@ -243 +243 @@
 !      EM    = 1.0
 C JL18 correction to account for the fact that the radiative top is not at zero optical depth.
-      EXPTRM(L) = MIN(TAUMAX,LAMDA(L)*(TAUCUM(2)))
-      EP = EXP(EXPTRM(L))
+      EXPTRM = MIN(TAUMAX,LAMDA(L)*(TAUCUM(2)))
+      EP = EXP(EXPTRM)
       EM    = 1.0/EP
       G4    = 1.0-G3(1)
@@ -286 +286 @@
 
       L     = L_NLAYRAD 
-      EXPTRM(L) = MIN(TAUMAX,LAMDA(L)*(TAUCUM(L_LEVELS)-
+      EXPTRM = MIN(TAUMAX,LAMDA(L)*(TAUCUM(L_LEVELS)-
      *                                 TAUCUM(L_LEVELS-1)))
 
-      EP    = EXP(EXPTRM(L))
+      EP    = EXP(EXPTRM)
       EM    = 1.0/EP
       G4    = 1.0-G3(L)