Changeset 2596

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

Location:

LMDZ5/trunk/libf/phylmd

Files:

• newmicro.F90 (modified) (7 diffs)
• rrtm/swr.F90 (modified) (2 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

LMDZ5/trunk/libf/phylmd/rrtm/swr.F90

@@ -192 +192 @@
     ZCLOUD(JL) = ZSS1(JL)
   ELSEIF (NOVLP == 2) THEN
-   stop 'provisoire pour verifier option novlp=1b'
+!IM150716  stop 'provisoire pour verifier option novlp=1b'
+   print*,'rrtm provisoire pour verifier option novlp=2 maximum'
 !* maximum
     ZCLOUD(JL) = MAX( ZSS1(JL) , ZCLOUD(JL) )
@@ -198 +199 @@
 !++MODIFCODE
   ELSEIF ((NOVLP == 3).OR.((NOVLP  >=  5).AND.(NOVLP /= 8))) THEN
-   stop 'provisoire pour verifier option novlp=1c'
+!IM150716  stop 'provisoire pour verifier option novlp=1c'
+    print*,'rrtm provisoire pour verifier option novlp=3 random'
 !--MODIFCODE
 !* random