Changeset 2596 for LMDZ5/trunk/libf/phylmd/newmicro.F90

Timestamp:

Jul 21, 2016, 5:46:23 PM (8 years ago)

Author:

musat

Message:

In newmicro allow use of cloud overlap hypothesis defined
in radopt.h as for the radiation code.
In rrtm allow max and random cloud overlap hypothesis.
IM

File:

• LMDZ5/trunk/libf/phylmd/newmicro.F90 (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/newmicro.F90

@@ -17 +17 @@
   ! Auteur(s): Z.X. Li (LMD/CNRS) date: 19930910
   ! O.   Boucher (LMD/CNRS) mise a jour en 201212
+  ! I. Musat (LMD/CNRS) : prise en compte de la meme hypothese de recouvrement
+  !                       pour les nuages que pour le rayonnement rrtm via
+  !                       le parametre novlp de radopt.h : 20160721
   ! Objet: Calculer epaisseur optique et emmissivite des nuages
   ! ======================================================================
@@ -53 +56 @@
   include "radopt.h"
 
-  ! choix de l'hypothese de recouvrememnt nuageuse
-  LOGICAL random, maximum_random, maximum
-  PARAMETER (random=.FALSE., maximum_random=.TRUE., maximum=.FALSE.)
+  ! choix de l'hypothese de recouvrement nuageuse via radopt.h (IM, 19.07.2016)
+  ! !novlp=1: max-random
+  ! !novlp=2: maximum
+  ! !novlp=3: random
+! LOGICAL random, maximum_random, maximum
+! PARAMETER (random=.FALSE., maximum_random=.TRUE., maximum=.FALSE.)
 
   LOGICAL, SAVE :: first = .TRUE.
@@ -512 +518 @@
       reffclwtop(i) = 0.
       cldncl(i) = 0.
-      IF (random .OR. maximum_random) tcc(i) = 1.
-      IF (maximum) tcc(i) = 0.
+      IF (novlp.EQ.3 .OR. novlp.EQ.1) tcc(i) = 1.
+      IF (novlp.EQ.2) tcc(i) = 0.
     END DO
 
@@ -524 +530 @@
         IF (pcltau(i,k)>thres_tau .AND. pclc(i,k)>thres_neb) THEN
 
-          IF (maximum) THEN
+          IF (novlp.EQ.2) THEN
             IF (first) THEN
               WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM'
@@ -533 +539 @@
           END IF
 
-          IF (random) THEN
+          IF (novlp.EQ.3) THEN
             IF (first) THEN
               WRITE (*, *) 'Hypothese de recouvrement: RANDOM'
@@ -542 +548 @@
           END IF
 
-          IF (maximum_random) THEN
+          IF (novlp.EQ.1) THEN
             IF (first) THEN
               WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM_ &
@@ -567 +573 @@
       END DO ! loop over k
 
-      IF (random .OR. maximum_random) tcc(i) = 1. - tcc(i)
+      IF (novlp.EQ.3 .OR. novlp.EQ.1) tcc(i) = 1. - tcc(i)
 
     END DO ! loop over i