Changeset 704 for LMDZ4/branches/V3_test/libf/phylmd/thermcell.F

Timestamp:

Aug 17, 2006, 5:41:51 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

File:

• LMDZ4/branches/V3_test/libf/phylmd/thermcell.F (modified) (23 diffs)

LMDZ4/branches/V3_test/libf/phylmd/thermcell.F

@@ -6 +6 @@
 c    s                  ,pu_therm,pv_therm
      s                  ,r_aspect,l_mix,w2di,tho)
-
+      USE dimphy
       IMPLICIT NONE
 
@@ -14 +14 @@
 c   de "thermiques" explicitement representes
 c
-c   Réécriture à partir d'un listing papier à Habas, le 14/02/00
-c
-c   le thermique est supposé homogène et dissipé par mélange avec
-c   son environnement. la longueur l_mix contrôle l'efficacité du
-c   mélange
-c
-c   Le calcul du transport des différentes espèces se fait en prenant
+c   R�criture �partir d'un listing papier �Habas, le 14/02/00
+c
+c   le thermique est suppos�homog�e et dissip�par m�ange avec
+c   son environnement. la longueur l_mix contr�e l'efficacit�du
+c   m�ange
+c
+c   Le calcul du transport des diff�entes esp�es se fait en prenant
 c   en compte:
 c     1. un flux de masse montant
@@ -33 +33 @@
 c   -------------
 
-#include "dimensions.h"
-#include "dimphy.h"
+cym#include "dimensions.h"
+cym#include "dimphy.h"
 #include "YOMCST.h"
 
@@ -52 +52 @@
       save idetr
       data idetr/3/
-
+c$OMP THREADPRIVATE(idetr)
 c   local:
 c   ------
@@ -85 +85 @@
       real fracc(klon,klev+1)
       real zf,zf2
-      real thetath2(klon,klev),wth2(klon,klev)
-      common/comtherm/thetath2,wth2
+      real,allocatable,save :: thetath2(:,:),wth2(:,:)
+c$OMP THREADPRIVATE(thetath2,wth2)
+cym      common/comtherm/thetath2,wth2
 
       real count_time
@@ -93 +94 @@
       data isplit/0/
       save isplit
-
+c$OMP THREADPRIVATE(isplit)
       logical sorties
       real rho(klon,klev),rhobarz(klon,klev+1),masse(klon,klev)
@@ -120 +121 @@
       data first /.false./
       save first
+c$OMP THREADPRIVATE(first)
 cRC
 
@@ -132 +134 @@
       save ncorrec
       data ncorrec/0/
-      
+c$OMP THREADPRIVATE(ncorrec)
+      logical,save :: firstCall=.true.
+c$OMP THREADPRIVATE(firstCall)
 c
 c-----------------------------------------------------------------------
@@ -138 +142 @@
 c   ---------------
 c
+      if (firstcall) then
+        allocate(thetath2(klon,klev),wth2(klon,klev))
+        thetath2(:,:)=0.
+        wth2(:,:)=0.
+        firstcall=.false.
+      endif
+      
        sorties=.true.
       IF(ngrid.NE.klon) THEN
@@ -501 +512 @@
 c   calcul de la largeur de chaque ascendance dans le cas conservatif.
 c   dans ce cas simple, on suppose que la largeur de l'ascendance provenant
-c   d'une couche est égale à la hauteur de la couche alimentante.
+c   d'une couche est �ale �la hauteur de la couche alimentante.
 c   La vitesse maximale dans l'ascendance est aussi prise comme estimation
 c   de la vitesse d'entrainement horizontal dans la couche alimentante.
@@ -540 +551 @@
 c      print*,'10 OK convect8'
 c     print*,'WA2 ',wa_moy
-c   calcul de la fraction de la maille concernée par l'ascendance en tenant
+c   calcul de la fraction de la maille concern� par l'ascendance en tenant
 c   compte de l'epluchage du thermique.
 c
@@ -979 +990 @@
       subroutine dqthermcell(ngrid,nlay,ptimestep,fm,entr,masse
      .    ,q,dq,qa)
+      use dimphy
       implicit none
 
@@ -989 +1001 @@
 c=======================================================================
 
-#include "dimensions.h"
-#include "dimphy.h"
+cym#include "dimensions.h"
+cym#include "dimphy.h"
 
       integer ngrid,nlay
@@ -1053 +1065 @@
      .    ,fraca,larga
      .    ,u,v,du,dv,ua,va)
+       use dimphy
       implicit none
 
@@ -1063 +1076 @@
 c=======================================================================
 
-#include "dimensions.h"
-#include "dimphy.h"
+cym#include "dimensions.h"
+cym#include "dimphy.h"
 
       integer ngrid,nlay
@@ -1106 +1119 @@
             if ((fm(ig,k+1)+detr(ig,k))*ptimestep.gt.
      s         1.e-5*masse(ig,k)) then
-c   On itère sur la valeur du coeff de freinage.
+c   On it�e sur la valeur du coeff de freinage.
 c              gamma0=rho(ig,k)*(zlev(ig,k+1)-zlev(ig,k))
                gamma0=masse(ig,k)
@@ -1112 +1125 @@
      s         *0.5/larga(ig)
 c              gamma0=0.
-c   la première fois on multiplie le coefficient de freinage
+c   la premi�e fois on multiplie le coefficient de freinage
 c   par le module du vent dans la couche en dessous.
                dua=ua(ig,k-1)-u(ig,k-1)
@@ -1166 +1179 @@
       subroutine dqthermcell2(ngrid,nlay,ptimestep,fm,entr,masse,frac
      .    ,q,dq,qa)
+      use dimphy
       implicit none
 
@@ -1176 +1190 @@
 c=======================================================================
 
-#include "dimensions.h"
-#include "dimphy.h"
+cym#include "dimensions.h"
+cym#include "dimphy.h"
 
       integer ngrid,nlay
@@ -1246 +1260 @@
      .    ,fraca,larga
      .    ,u,v,du,dv,ua,va)
+      use dimphy
       implicit none
 
@@ -1256 +1271 @@
 c=======================================================================
 
-#include "dimensions.h"
-#include "dimphy.h"
+cym#include "dimensions.h"
+cym#include "dimphy.h"
 
       integer ngrid,nlay
@@ -1302 +1317 @@
             if ((fm(ig,k+1)+detr(ig,k))*ptimestep.gt.
      s         1.e-5*masse(ig,k)) then
-c   On itère sur la valeur du coeff de freinage.
+c   On it�e sur la valeur du coeff de freinage.
 c              gamma0=rho(ig,k)*(zlev(ig,k+1)-zlev(ig,k))
                gamma0=masse(ig,k)
@@ -1313 +1328 @@
                zf=0.
                zf2=1./(1.-zf)
-c   la première fois on multiplie le coefficient de freinage
+c   la premi�e fois on multiplie le coefficient de freinage
 c   par le module du vent dans la couche en dessous.
                dua=ua(ig,k-1)-u(ig,k-1)