Changeset 227 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/condense_cloud.F90

Timestamp:

07/16/14 18:05:01 (10 years ago)

Author:

milmd

Message:

Last LMDZ version (1315) with OpenMP directives and other stuff

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/condense_cloud.F90 (modified) (9 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/condense_cloud.F90

@@ -32 +32 @@
 !     ptsrf(ngrid)          Surface temperature
 !     
-!     pdt(ngrid,nlayermx)   Time derivative before condensation/sublimation of pt
+!     pdt(ngrid,nlayer)   Time derivative before condensation/sublimation of pt
 !     pdtsrf(ngrid)         Time derivative before condensation/sublimation of ptsrf
 !     pqsurf(ngrid,nq)      Sedimentation flux at the surface (kg.m-2.s-1)
@@ -38 +38 @@
 !     Outputs
 !     -------
-!     pdpsrf(ngrid)         \  Contribution of condensation/sublimation
-!     pdtc(ngrid,nlayermx)  /  to the time derivatives of Ps, pt, and ptsrf
-!     pdtsrfc(ngrid)       /
+!     pdpsrf(ngrid)       \  Contribution of condensation/sublimation
+!     pdtc(ngrid,nlayer)  /  to the time derivatives of Ps, pt, and ptsrf
+!     pdtsrfc(ngrid)     /
 !     
 !     Both
@@ -56 +56 @@
 !==================================================================
 
-#include "dimensions.h"
-#include "dimphys.h"
+!#include "dimensions.h"
+!#include "dimphys.h"
 #include "comcstfi.h"
-#include "comvert.h"
+!#include "comvert.h"
 #include "callkeys.h"
 
@@ -99 +99 @@
 
       REAL reffrad(ngrid,nlayer) ! radius (m) of the co2 ice particles
-      REAL*8 zt(ngrid,nlayermx)
-      REAL zq(ngrid,nlayermx,nq)
+      REAL*8 zt(ngrid,nlayer)
+      REAL zq(ngrid,nlayer,nq)
       REAL zcpi
-      REAL ztcond (ngrid,nlayermx)
-      REAL ztnuc (ngrid,nlayermx)
+      REAL ztcond (ngrid,nlayer)
+      REAL ztnuc (ngrid,nlayer)
       REAL ztcondsol(ngrid) 
       REAL zdiceco2(ngrid)
-      REAL zcondicea(ngrid,nlayermx), zcondices(ngrid)
+      REAL zcondicea(ngrid,nlayer), zcondices(ngrid)
       REAL zfallice(ngrid), Mfallice(ngrid) 
-      REAL zmflux(nlayermx+1)
-      REAL zu(nlayermx),zv(nlayermx)
+      REAL zmflux(nlayer+1)
+      REAL zu(nlayer),zv(nlayer)
       REAL ztsrf(ngrid)
-      REAL ztc(nlayermx), ztm(nlayermx+1) 
-      REAL zum(nlayermx+1) , zvm(nlayermx+1)
+      REAL ztc(nlayer), ztm(nlayer+1) 
+      REAL zum(nlayer+1) , zvm(nlayer+1)
       LOGICAL condsub(ngrid)
       REAL subptimestep
       Integer Ntime
-      real masse (ngrid,nlayermx), w(ngrid,nlayermx,nq)
-      real wq(ngrid,nlayermx+1)
+      real masse (ngrid,nlayer), w(ngrid,nlayer,nq)
+      real wq(ngrid,nlayer+1)
       real vstokes,reff
 
 !     Special diagnostic variables
-      real tconda1(ngrid,nlayermx)
-      real tconda2(ngrid,nlayermx)
-      real zdtsig (ngrid,nlayermx)
-      real zdt (ngrid,nlayermx)
+      real tconda1(ngrid,nlayer)
+      real tconda2(ngrid,nlayer)
+      real zdtsig (ngrid,nlayer)
+      real zdt (ngrid,nlayer)
 
 !-----------------------------------------------------------------------
@@ -133 +133 @@
       REAL,SAVE :: cpice=1000.
       REAL,SAVE,ALLOCATABLE,DIMENSION(:) :: emisref
+!$OMP THREADPRIVATE(latcond,ccond,cpice,emisref)
 
       LOGICAL,SAVE :: firstcall=.true.
+!$OMP THREADPRIVATE(firstcall)
       REAL,EXTERNAL :: SSUM
 
@@ -140 +142 @@
 
       INTEGER,SAVE :: i_co2ice=0      ! co2 ice
+!$OMP THREADPRIVATE(i_co2ice)
       CHARACTER(LEN=20) :: tracername ! to temporarily store text
 
@@ -205 +208 @@
 !     zcondices(ngrid)     condensation rate on the ground  (kg/m2/s)
 !     zfallice(ngrid)      flux of ice falling on surface   (kg/m2/s)
-!     pdtc(ngrid,nlayermx) dT/dt due to phase changes       (K/s)
+!     pdtc(ngrid,nlayer) dT/dt due to phase changes       (K/s)
      
 
@@ -301 +304 @@
             
 !     sedimentation computed from radius computed from q in module radii_mod
-         call co2_reffrad(ngrid,nq,zq,reffrad)
+         call co2_reffrad(ngrid,nlayer,nq,zq,reffrad)
          
          do  ilay=1,nlayer
@@ -321 +324 @@
 !     Computing q after sedimentation
 
-         call vlz_fi(ngrid,zq(1,1,i_co2ice),2.,masse,w(1,1,i_co2ice),wq)
+         call vlz_fi(ngrid,nlayer,zq(1,1,i_co2ice),2.,masse,w(1,1,i_co2ice),wq)