Changeset 227 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/callsedim.F

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/callsedim.F (modified) (3 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/callsedim.F

@@ -28 +28 @@
 c   -------------
 
-#include "dimensions.h"
-#include "dimphys.h"
+!#include "dimensions.h"
+!#include "dimphys.h"
 #include "comcstfi.h"
 #include "callkeys.h"
@@ -65 +65 @@
       real epaisseur (ngrid,nlay) ! Layer thickness (m)
       real wq(ngrid,nlay+1) ! displaced tracer mass (kg.m-2)
-c      real dens(ngrid,nlayermx) ! Mean density of the ice part. accounting for dust core
+c      real dens(ngrid,nlay) ! Mean density of the ice part. accounting for dust core
 
 
       LOGICAL,SAVE :: firstcall=.true.
+!$OMP THREADPRIVATE(firstcall)
 
 c    ** un petit test de coherence
@@ -125 +126 @@
           if (water.and.(iq.eq.igcm_h2o_ice)) then
             ! compute radii for h2o_ice 
-             call h2o_reffrad(ngrid,zqi(1,1,igcm_h2o_ice),zt,
+             call h2o_reffrad(ngrid,nlay,zqi(1,1,igcm_h2o_ice),zt,
      &                reffrad(1,1,iaero_h2o),nueffrad(1,1,iaero_h2o))
             ! call sedimentation for h2o_ice