Changeset 253 for trunk/LMDZ.GENERIC/libf/phystd/newsedim.F

Timestamp:

Aug 2, 2011, 11:13:07 AM (13 years ago)

Author:

emillour

Message:

Generic GCM

• Massive update to version 0.7

EM+RW

File:

• trunk/LMDZ.GENERIC/libf/phystd/newsedim.F (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/phystd/newsedim.F

@@ -12 +12 @@
 !==================================================================
 
-c-----------------------------------------------------------------------
-c   declarations:
-c   -------------
+!-----------------------------------------------------------------------
+!   declarations
+!   ------------
 
 #include "dimensions.h"
 #include "dimphys.h"
 #include "comcstfi.h"
-c
-c   arguments:
-c   ----------
+
+!   arguments
+!   ---------
 
       INTEGER ngrid,nlay,naersize
@@ -57 +57 @@
 c    Physical constant
 c    ~~~~~~~~~~~~~~~~~
-!      REAL visc, molrad
+      REAL visc, molrad
 c     Gas molecular viscosity (N.s.m-2)
-!      data visc/1.e-5/       ! CO2
+      data visc/1.e-5/       ! CO2
 c     Effective gas molecular radius (m)
-!      data molrad/2.2e-10/   ! CO2
+      data molrad/2.2e-10/   ! CO2
 
 c     local and saved variable
@@ -69 +69 @@
 c    ** un petit test de coherence
 c       --------------------------
+
+
+      !print*,'temporary fixed particle rad in newsedim!!'
 
       IF (firstcall) THEN
@@ -81 +84 @@
 
 
+
 !=======================================================================
 !     Preliminary calculations for sedimenation velocity
@@ -95 +99 @@
 c   (correction = 0.85 for irregular particles, 0.5 for disk shaped particles)
 c        a = a    * 0.85
+
+
       ENDIF
-
-
-
-
 
 c-----------------------------------------------------------------------
@@ -116 +118 @@
             else
               i=ngrid*(l-1)+ig
-              rfall=rd(i)
+              rfall=rd(i) ! how can this be correct!!?
             endif  
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! TEMPORARY MODIF BY RDW !!!!
+            !rfall=5.e-6
+            if(rfall.lt.1.e-7)then
+               rfall=1.e-7
+            endif
+            !if(rfall.gt.5.e-5)then
+            !   rfall=5.e-5
+            !endif
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
             vstokes(ig,l) = b * rho * rfall**2 *
      &      (1 + 1.333* ( a*pt(ig,l)/pplev(ig,l) )/rfall)