Changeset 740 for trunk/LMDZ.MARS/libf/phymars/simpleclouds.F

Timestamp:

Jul 25, 2012, 1:03:19 PM (12 years ago)

Author:

tnavarro

Message:

module for ice and radius radius computation

File:

• trunk/LMDZ.MARS/libf/phymars/simpleclouds.F (modified) (7 diffs)

trunk/LMDZ.MARS/libf/phymars/simpleclouds.F

@@ -3 +3 @@
      &             pq,pdq,pdqcloud,pdtcloud,
      &             nq,tau,rice)
+      USE updaterad
       implicit none
 c------------------------------------------------------------------
@@ -41 +42 @@
       integer nq                 ! nombre de traceurs
       REAL ptimestep             ! pas de temps physique (s)
-!      REAL pplev(ngrid,nlay+1)   ! pression aux inter-couches (Pa)
       REAL pplay(ngrid,nlay)     ! pression au milieu des couches (Pa)
       REAL pzlay(ngrid,nlay)     ! altitude at the middle of the layers
@@ -68 +68 @@
       SAVE firstcall
 
-      REAL CBRT
-      EXTERNAL CBRT
+            
+      REAL rhocloud(ngridmx,nlayermx)  ! Cloud density (kg.m-3)
 
       INTEGER ig,l
@@ -86 +86 @@
 c      REAL, PARAMETER :: ccn_factor = 4.5  !! comme TESTS_JB // 1. avant
       
-      REAL Mcon_out(ngridmx,nlayermx) ! mass to be condensed (not dMice !!)
 
 c-----------------------------------------------------------------------
@@ -110 +109 @@
       enddo
 
-      do l=1,nlay
-        do ig=1,ngrid
-
-c         Typical dust radius profile:
-          rdusttyp(ig,l)= max(.8e-6*exp(-pzlay(ig,l)/18000.),1.e-9)
-
-c         Typical CCN profile:
-c         Corrected equation, following Montmessin et al. 2004
-c           (N0=2e6 m-3 has been converted into N0=1.3e8 kg-1, otherwise
-c            the equation for rice is not homogeneous...)
-          ccntyp(ig,l)=
-     &       1.3e+8*max(tau(ig,1),0.001)/0.1*exp(-pzlay(ig,l)/10000.)
-c         The previously used profile was not correct:
-c         ccntyp(ig,l)=( epaisseur(ig,l)/masse(ig,l) ) *
-c    &       2.e+6/0.1*max(tau(ig,1),0.001)*exp(-pzlay(ig,l)/10000.)
-
-c         Reduce number of nuclei
-!         TEMPORAIRE : decrease the number of CCNs FF 04/200
-!         reduction facteur 3
-!         ccntyp(ig,l) = ccntyp(ig,l) / 27.
-!         reduction facteur 2
-!         ccntyp(ig,l) = ccntyp(ig,l) / 8.
-c -----------------------------------------------------------------
-          ccntyp(ig,l) = ccntyp(ig,l) / ccn_factor
-
-        enddo ! of do ig=1,ngrid
-      enddo ! of do l=1,nlay
 
       pdqcloud(1:ngrid,1:nlay,1:nq)=0
@@ -167 +139 @@
           zq(ig,l,igcm_h2o_vap)=zq(ig,l,igcm_h2o_vap)-dzq
           
-          Mcon_out(ig,l) = dzq
-
-c         Calcul du rayon moyen des particules de glace.
-c         Hypothese : Dans une couche, la glace presente se
-c         repartit uniformement autour du nbre de poussieres
-c         dont le rayon moyen est prescrit par rdusttyp.
-c         ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-          rice(ig,l)=max( CBRT ( (zq(ig,l,igcm_h2o_ice)/rho_ice
-     &      +ccntyp(ig,l)*(4./3.)*pi*rdusttyp(ig,l)**3.)
-     &      /(ccntyp(ig,l)*4./3.*pi) ), rdusttyp(ig,l)) 
-
 
         enddo ! of do ig=1,ngrid
@@ -195 +156 @@
       end do
 
-c------------------------------------------------------------------
-c     TEST_JBM
-!      IF (ngrid.eq.1) THEN
-!         call WRITEDIAGFI(ngrid,"mcond","h2o condensed mass","kg",1,
-!     &                    Mcon_out)
-!         call WRITEDIAGFI(ngrid,"rdusttyp","rdusttyp","m",1,
-!     &                    rdusttyp)
-!         call WRITEDIAGFI(ngrid,"ccntyp","ccntyp","kg-1",1,
-!     &                    ccntyp)
-!      ENDIF
+c     ice crystal radius
+      do l=1, nlay
+        do ig=1,ngridmx
+          call updaterice_typ(zq(ig,l,igcm_h2o_ice),
+     &       tau(ig,1),pzlay(ig,l),rice(ig,l))
+        end do
+      end do
+
 c------------------------------------------------------------------
       return