Changeset 735 for LMDZ4/branches/V3_test/libf/phylmd/physiq.F

Timestamp:

Nov 10, 2006, 3:50:16 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Modifications version parallele
YM/LF

File:

• LMDZ4/branches/V3_test/libf/phylmd/physiq.F (modified) (19 diffs)

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

@@ -251 +251 @@
       CHARACTER*3 ctetaSTD(nbteta)
       DATA ctetaSTD/'350','380','405'/
+      SAVE ctetaSTD
 c$OMP THREADPRIVATE(ctetaSTD)
       REAL rtetaSTD(nbteta)
       DATA rtetaSTD/350., 380., 405./
+      SAVE rtetaSTD
 c$OMP THREADPRIVATE(rtetaSTD)     
 c
@@ -299 +301 @@
      .60000., 50000., 40000., 30000., 25000., 20000.,
      .15000., 10000., 7000., 5000., 3000., 2000., 1000./
+      SAVE rlevstd
 c$OMP THREADPRIVATE(rlevSTD)
       CHARACTER*4 clevSTD(nlevSTD)
@@ -304 +307 @@
      .'500 ','400 ','300 ','250 ','200 ','150 ','100 ',
      .'70  ','50  ','30  ','20  ','10  '/
+      SAVE clevSTD
 c$OMP THREADPRIVATE(clevSTD)
 c
@@ -416 +420 @@
       DATA ifreq_isccp/3/
       SAVE ifreq_isccp
-c$OMP THREAPRIVATE(ifreq_isccp)
+c$OMP THREADPRIVATE(ifreq_isccp)
       CHARACTER*5 typinout(napisccp)
       DATA typinout/'i3od'/
-c$OMP THREAPRIVATE(typinout)
+      SAVE typinout
+c$OMP THREADPRIVATE(typinout)
 cIM verif boxptop BEG
       CHARACTER*1 verticaxe(napisccp)
       DATA verticaxe/'1'/ 
-c$OMP THREAPRIVATE(verticaxe)
+      SAVE verticaxe
+c$OMP THREADPRIVATE(verticaxe)
 cIM verif boxptop END
       INTEGER nvlev(napisccp)
@@ -446 +452 @@
 cIM verif boxptop END
 c
-      REAL tautab(0:255)
-      INTEGER invtau(-20:45000)
+      REAL,SAVE :: tautab_omp(0:255),tautab(0:255)
+      INTEGER,SAVE :: invtau_omp(-20:45000),invtau(-20:45000)
+c$OMP THREADPRIVATE(tautab,invtau)
       REAL emsfc_lw
       PARAMETER(emsfc_lw=0.99)
@@ -523 +530 @@
 c     save ok_isccp, ecrit_isccp, nid_isccp        
       save nid_isccp        
-c$OMP THREADPRIVATE(ok_isccp, ecrit_isccp, nid_isccp,nid_isccp)
+c$OMP THREADPRIVATE(nid_isccp)
+
 c     data ok_isccp,ecrit_isccp/.true.,0.125/      
 c     data ok_isccp,ecrit_isccp/.true.,1./      
@@ -533 +541 @@
       REAL zx_tau(kmaxm1), zx_pc(lmaxm1), zx_o500(iwmax)
       DATA zx_tau/0.0, 0.3, 1.3, 3.6, 9.4, 23., 60./
+      SAVE zx_tau
 cIM bad 151205     DATA zx_pc/50., 180., 310., 440., 560., 680., 800./
       DATA zx_pc/180., 310., 440., 560., 680., 800., 1000./
+      SAVE zx_pc
 c$OMP THREADPRIVATE(zx_tau,zx_pc)
 c cldtopres pression au sommet des nuages
@@ -541 +551 @@
       DATA cldtopres/180., 310., 440., 560., 680., 800., 1000./
       DATA cldtopres3/440., 680., 1000./
+      SAVE cldtopres,cldtopres3
 c$OMP THREADPRIVATE(cldtopres,cldtopres3)
 cIM 051005 BEG
@@ -556 +567 @@
       CHARACTER *3 pclev(lmaxm1)
       DATA pclev/'pc1','pc2','pc3','pc4','pc5','pc6','pc7'/
+      SAVE taulev,pclev
 c$OMP THREADPRIVATE(taulev,pclev)
 c
@@ -610 +622 @@
      .                'pc= 680-800hPa, tau> 60.',
      .                'pc= 800-1000hPa, tau> 60.'/
+       SAVE cnameisccp
 c$OMP THREADPRIVATE(cnameisccp)
 c
@@ -838 +851 @@
       REAL yv1(klon)            ! vents dans la premiere couche V
       REAL,SAVE,ALLOCATABLE :: ffonte(:,:)    !Flux thermique utilise pour fondre la neige
-c$OMP THREAPRIVATE(ffonte)
+c$OMP THREADPRIVATE(ffonte)
       REAL,SAVE,ALLOCATABLE :: fqcalving(:,:) !Flux d'eau "perdu" par la surface 
-c$OMP THREAPRIVATE(fqcalving)
+c$OMP THREADPRIVATE(fqcalving)
       REAL,SAVE,ALLOCATABLE :: fqfonte(:,:)  !Quantite d'eau de fonte des glaciers
-c$OMP THREAPRIVATE(fqcalving)
+c$OMP THREADPRIVATE(fqfonte)
 c                               !et necessaire pour limiter la
 c                               !hauteur de neige, en kg/m2/s
@@ -1074 +1087 @@
 c
       REAL,SAVE,ALLOCATABLE :: pblh(:, :)           ! Hauteur de couche limite
-c$OMP THREAPRIVATE(pblh)
+c$OMP THREADPRIVATE(pblh)
       REAL,SAVE,ALLOCATABLE :: plcl(:, :)           ! Niveau de condensation de la CLA
-c$OMP THREAPRIVATE(plcl)
+c$OMP THREADPRIVATE(plcl)
       REAL,SAVE,ALLOCATABLE :: capCL(:, :)          ! CAPE de couche limite
-c$OMP THREAPRIVATE(capCL)
+c$OMP THREADPRIVATE(capCL)
       REAL,SAVE,ALLOCATABLE :: oliqCL(:, :)          ! eau_liqu integree de couche limite
-c$OMP THREAPRIVATE(oliqCL)
+c$OMP THREADPRIVATE(oliqCL)
       REAL,SAVE,ALLOCATABLE :: cteiCL(:, :)          ! cloud top instab. crit. couche limite
-c$OMP THREAPRIVATE(cteiCL)
+c$OMP THREADPRIVATE(cteiCL)
       REAL,SAVE,ALLOCATABLE :: pblt(:, :)          ! T a la Hauteur de couche limite
-c$OMP THREAPRIVATE(pblt)
+c$OMP THREADPRIVATE(pblt)
       REAL,SAVE,ALLOCATABLE :: therm(:, :)
-c$OMP TREADPRIVATE(therm)
+c$OMP THREADPRIVATE(therm)
       REAL,SAVE,ALLOCATABLE :: trmb1(:, :)          ! deep_cape
-c$OMP TREADPRIVATE(trmb1)
+c$OMP THREADPRIVATE(trmb1)
       REAL,SAVE,ALLOCATABLE :: trmb2(:, :)          ! inhibition 
-c$OMP TREADPRIVATE(trmb2)
+c$OMP THREADPRIVATE(trmb2)
       REAL,SAVE,ALLOCATABLE :: trmb3(:, :)          ! Point Omega
-c$OMP TREADPRIVATE(trmb3)
+c$OMP THREADPRIVATE(trmb3)
 c Grdeurs de sorties
       REAL s_pblh(klon), s_lcl(klon), s_capCL(klon)
@@ -1289 +1302 @@
       REAL zx_tmp_fi2d(klon)      ! variable temporaire grille physique
       REAL zx_tmp_fi3d(klon,klev) ! variable temporaire pour champs 3D 
-#ifdef histmthNMC
+c#ifdef histmthNMC
 cym   A voir plus tard !!!!
 cym      REAL zx_tmp_NC(iim,jjmp1,nlevSTD)
       REAL zx_tmp_fiNC(klon,nlevSTD) 
-#endif
+c#endif
       REAL*8 zx_tmp2_fi3d(klon,klev) ! variable temporaire pour champs 3D 
       REAL zx_tmp_2d(iim,jjmp1), zx_tmp_3d(iim,jjmp1,klev)
@@ -1302 +1315 @@
       SAVE nid_day, nid_mth, nid_ins, nid_nmc, nid_day_seri
       SAVE nid_ctesGCM
-c$OMP THREADPRIVATE(nid_day, nid_mth, nid_ins, nid_nmc, nid_day_seri,nid_ctesGCM)
+c$OMP THREADPRIVATE(nid_day, nid_mth, nid_ins, nid_nmc)
+c$OMP THREADPRIVATE(nid_day_seri,nid_ctesGCM)
 c
 cIM 280405 BEG
@@ -1403 +1417 @@
 
       REAL,SAVE,ALLOCATABLE :: topswai(:), solswai(:) ! Aerosol indirect effect.
-c$OMP THREADPRIVATE(topswai(,solswai)
+c$OMP THREADPRIVATE(topswai,solswai)
       ! ok_aie=T ->
       !        ok_ade=T -AIE=topswai-topswad
@@ -1908 +1922 @@
 
 c#include "ini_histday_seri.h"
-
-c#include "ini_paramLMDZ_phy.h"
+#include "ini_histday_seri.h"
+
+#include "ini_paramLMDZ_phy.h"
 
 #endif
@@ -2857 +2872 @@
       ELSE IF (iflag_cldcon.eq.3) THEN
 c  On prend pour les nuages convectifs le max du calcul de la
-c  convection et du calcul du pas de temps pr��ent diminu�d'un facteur
+c  convection et du calcul du pas de temps precedent diminue d'un facteur
 c  facttemps
 c      facttemps=pdtphys/1.e4
@@ -2895 +2910 @@
 c
       IF (ok_isccp) THEN
-c#include "calcul_simulISCCP.h"
+#include "calcul_simulISCCP.h"
       ENDIF !ok_isccp
 
@@ -3581 +3596 @@
 #endif
 
-c#include "write_histday_seri.h"
-
-c#include "write_paramLMDZ_phy.h"
+#include "write_histday_seri.h"
+
+#include "write_paramLMDZ_phy.h"
 
 #endif