source: LMDZ4/trunk/libf/dyn3dpar/qminimum_p.F @ 2629

Last change on this file since 2629 was 1146, checked in by Laurent Fairhead, 16 years ago

Réintegration dans le tronc des modifications issues de la branche LMDZ-dev
comprises entre la révision 1074 et 1145
Validation: une simulation de 1 jour en séquentiel sur PC donne les mêmes
résultats entre la trunk et la dev
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.0 KB
RevLine 
[630]1      SUBROUTINE qminimum_p( q,nq,deltap )
2      USE parallel
3      IMPLICIT none
4c
5c  -- Objet : Traiter les valeurs trop petites (meme negatives)
6c             pour l'eau vapeur et l'eau liquide
7c
8#include "dimensions.h"
9#include "paramet.h"
10#include "comvert.h"
11c
12      INTEGER nq
13      REAL q(ip1jmp1,llm,nq), deltap(ip1jmp1,llm)
14c
15      INTEGER iq_vap, iq_liq
16      PARAMETER ( iq_vap = 1 ) ! indice pour l'eau vapeur
17      PARAMETER ( iq_liq = 2 ) ! indice pour l'eau liquide
18      REAL seuil_vap, seuil_liq
19      PARAMETER ( seuil_vap = 1.0e-10 ) ! seuil pour l'eau vapeur
20      PARAMETER ( seuil_liq = 1.0e-11 ) ! seuil pour l'eau liquide
21c
22c  NB. ....( Il est souhaitable mais non obligatoire que les valeurs des
23c            parametres seuil_vap, seuil_liq soient pareilles a celles
24c            qui  sont utilisees dans la routine    ADDFI       )
25c     .................................................................
26c
27      INTEGER i, k, iq
28      REAL zx_defau, zx_abc, zx_pump(ip1jmp1), pompe
29c
30      REAL SSUM
31      EXTERNAL SSUM
32c
33      INTEGER imprim
34      SAVE imprim
35      DATA imprim /0/
[985]36c$OMP THREADPRIVATE(imprim)
[630]37      INTEGER ijb,ije
[985]38      INTEGER Index_pump(ip1jmp1)
39      INTEGER nb_pump
[630]40c
41c Quand l'eau liquide est trop petite (ou negative), on prend
42c l'eau vapeur de la meme couche et la convertit en eau liquide
43c (sans changer la temperature !)
44c
[985]45
[630]46      ijb=ij_begin
47      ije=ij_end
[985]48
49c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)       
[630]50      DO 1000 k = 1, llm
51      DO 1040 i = ijb, ije
[1146]52            if (seuil_liq - q(i,k,iq_liq) .gt. 0.d0 ) then
53               q(i,k,iq_vap) = q(i,k,iq_vap) + q(i,k,iq_liq) - seuil_liq
54               q(i,k,iq_liq) = seuil_liq
55            endif
[630]56 1040 CONTINUE
57 1000 CONTINUE
[985]58c$OMP END DO NOWAIT
59c$OMP BARRIER
60c --->  SYNCHRO OPENMP ICI
61
[630]62c
63c Quand l'eau vapeur est trop faible (ou negative), on complete
64c le defaut en prennant de l'eau vapeur de la couche au-dessous.
65c
66      iq = iq_vap
67c
68      DO k = llm, 2, -1
69ccc      zx_abc = dpres(k) / dpres(k-1)
[985]70c$OMP DO SCHEDULE(STATIC)
[630]71      DO i = ijb, ije
[1146]72         if ( seuil_vap - q(i,k,iq) .gt. 0.d0 ) then
73            q(i,k-1,iq) =  q(i,k-1,iq) - ( seuil_vap - q(i,k,iq) ) *
74     &           deltap(i,k) / deltap(i,k-1)
75            q(i,k,iq)   =  seuil_vap 
76         endif
[630]77      ENDDO
[985]78c$OMP END DO NOWAIT
[630]79      ENDDO
[985]80c$OMP BARRIER
[630]81c
82c Quand il s'agit de la premiere couche au-dessus du sol, on
83c doit imprimer un message d'avertissement (saturation possible).
84c
[985]85      nb_pump=0
86c$OMP DO SCHEDULE(STATIC)
[630]87      DO i = ijb, ije
88         zx_pump(i) = AMAX1( 0.0, seuil_vap - q(i,1,iq) )
89         q(i,1,iq)  = AMAX1( q(i,1,iq), seuil_vap )
[985]90         IF (zx_pump(i) > 0.0) THEN
91            nb_pump = nb_pump+1
92            Index_pump(nb_pump)=i
93         ENDIF
[630]94      ENDDO
[985]95c$OMP END DO 
96!      pompe = SSUM(ije-ijb+1,zx_pump(ijb),1)
97
98      IF (imprim.LE.100 .AND. nb_pump .GT. 0 ) THEN
99         PRINT *, 'ATT!:on pompe de l eau au sol'
100         DO i = 1, nb_pump
[630]101               imprim = imprim + 1
[985]102               PRINT*,'  en ',index_pump(i),zx_pump(index_pump(i))
[630]103         ENDDO
104      ENDIF
105c
106      RETURN
107      END
Note: See TracBrowser for help on using the repository browser.