/[lmdze]/trunk/Sources/phylmd/CV30_routines/cv30_unsat.f

Diff of /trunk/Sources/phylmd/CV30_routines/cv30_unsat.f

Parent Directory | Revision Log | View Patch Patch

-revision 195 by guez,
Wed May 18 17:56:44 2016 UTC
+revision 198 by guez,
Tue May 31 16:17:35 2016 UTC
 Line 10 
 contains
      ! Unsaturated (precipitating) downdrafts
      use cv30_param_m, only: nl, sigd
-     use cv_thermo_m, only: cpd, ginv, grav
+     use cv_thermo_m, only: cpd, ginv
+     use SUPHEC_M, only: rg
      integer, intent(in):: icb(:) ! (ncum)
+     ! {2 <= icb <= nl - 3}
      integer, intent(in):: inb(:) ! (ncum)
      ! first model level above the level of neutral buoyancy of the
-Line 20 
 contains
+Line 22 
 contains
      real, intent(in):: t(:, :), q(:, :), qs(:, :) ! (ncum, nl)
      real, intent(in):: gz(:, :) ! (klon, klev)
-     real, intent(in):: u(:, :), v(:, :) ! (klon, klev)
+     real, intent(in):: u(:, :), v(:, :) ! (ncum, nl)
      real, intent(in):: p(:, :) ! (klon, klev) pressure at full level, in hPa
-     real, intent(in):: ph(:, :) ! (klon, klev + 1)
+     real, intent(in):: ph(:, :) ! (ncum, klev + 1)
      real, intent(in):: th(:, :) ! (ncum, nl - 1)
      real, intent(in):: tv(:, :) ! (klon, klev)
-     real, intent(in):: lv(:, :) ! (klon, klev)
+     real, intent(in):: lv(:, :) ! (ncum, nl)
      real, intent(in):: cpn(:, :) ! (klon, klev)
      real, intent(in):: ep(:, :) ! (ncum, klev)
      real, intent(in):: clw(:, :) ! (ncum, klev)
-Line 33 
 contains
+Line 35 
 contains
      real, intent(in):: ment(:, :, :) ! (ncum, klev, klev)
      real, intent(in):: elij(:, :, :) ! (ncum, klev, klev)
      real, intent(in):: delt
-     real, intent(in):: plcl(:) ! (klon)
+     real, intent(in):: plcl(:) ! (ncum)
      real, intent(out):: mp(:, :) ! (klon, klev)
+     ! mass flux of the unsaturated downdraft, defined positive downward
+     ! M_p in Emanuel (1991 928)
      real, intent(out):: qp(:, :), up(:, :), vp(:, :) ! (ncum, nl)
      real, intent(out):: wt(:, :) ! (ncum, nl)
-     real, intent(out):: water(:, :), evap(:, :) ! (ncum, nl)
+     real, intent(out):: water(:, :) ! (ncum, nl)
+     ! precipitation mixing ratio, l_p in Emanuel (1991 928)
+     real, intent(out):: evap(:, :) ! (ncum, nl)
+     ! sigt * rate of evaporation of precipitation, in s-1
+     ! \sigma_s E in Emanuel (1991 928)
      real, intent(out):: b(:, :) ! (ncum, nl - 1)
      ! Local:
-Line 51 
 contains
+Line 63 
 contains
      integer i, il, imax
      real tinv, delti
      real afac, afac1, afac2, bfac
-     real pr1, pr2, sigt, b6, c6, revap, tevap, delth
+     real pr1, sigt, b6, c6, revap, tevap, delth
      real amfac, amp2, xf, tf, fac2, ur, sru, fac, d, af, bf
      real ampmax
      real lvcp(size(icb), nl) ! (ncum, nl)
-Line 93 
 contains
+Line 105 
 contains
         ! and condensed water flux
         ! Calculate detrained precipitation
-        forall (il = 1:ncum, inb(il) >= i .and. lwork(il)) wdtrain(il) = grav &
+        forall (il = 1:ncum, inb(il) >= i .and. lwork(il)) wdtrain(il) = rg &
              * (ep(il, i) * m(il, i) * clw(il, i) &
              + sum(max(elij(il, :i - 1, i) - (1. - ep(il, i)) * clw(il, i), 0.) &
              * ment(il, :i - 1, i)))
-Line 135 
 contains
+Line 147 
 contains
               afac = max(afac, 0.)
               bfac = 1. / (sigd * wt(il, i))
-              ! Prise en compte de la variation progressive de sigt dans
+              if (i <= icb(il)) then
-              ! les couches icb et icb - 1:
+                 ! Prise en compte de la variation progressive de sigt dans
-              ! pour plcl <= ph(i + 1), pr1 = 0 et pr2 = 1
+                 ! les couches icb et icb - 1 :
-              ! pour plcl >= ph(i), pr1 = 1 et pr2 = 0
+                 ! pour plcl <= ph(i + 1), pr1 = 0
-              ! pour ph(i + 1) < plcl < ph(i), pr1 est la proportion \`a cheval
+                 ! pour plcl >= ph(i), pr1 = 1
-              ! sur le nuage, et pr2 est la proportion sous la base du
+                 ! pour ph(i + 1) < plcl < ph(i), pr1 est la proportion
-              ! nuage.
+                 ! \`a cheval sur le nuage.
-              pr1 = max(0., min(1., &
+                 pr1 = max(0., min(1., &
-                   (plcl(il) - ph(il, i + 1)) / (ph(il, i) - ph(il, i + 1))))
+                      (plcl(il) - ph(il, i + 1)) / (ph(il, i) - ph(il, i + 1))))
-              pr2 = max(0., min(1., &
+                 sigt = sigp * pr1 + 1. - pr1
-                   (ph(il, i) - plcl(il)) / (ph(il, i) - ph(il, i + 1))))
+              else
-              sigt = sigp * pr1 + pr2
+                 ! {i >= icb(il) + 1}
+                 sigt = sigp
+              end if
               b6 = bfac * 50. * sigd * (ph(il, i) - ph(il, i + 1)) * sigt * afac
               c6 = water(il, i + 1) + bfac * wdtrain(il) - 50. * sigd * bfac &
-Line 209 
 contains
+Line 223 
 contains
                     mp(il, i) = max(0., mp(il, i))
-                    ! Il y a vraisemblablement une erreur dans la
+                    ! Il y a vraisemblablement une erreur dans la ligne
-                    ! ligne suivante : il faut diviser par (mp(il,
+                    ! suivante : il faut diviser par (mp(il, i) * sigd
-                    ! i) * sigd * grav) et non par (mp(il, i) + sigd
+                    ! * rg) et non par (mp(il, i) + sigd * 0.1).  Et il
-                    ! * 0.1).  Et il faut bien revoir les facteurs
+                    ! faut bien revoir les facteurs 100.
-                    ! 100.
                     b(il, i - 1) = b(il, i) + 100. * (p(il, i - 1) - p(il, i)) &
                          * tevap / (mp(il, i) + sigd * 0.1) - 10. * (th(il, i) &
                          - th(il, i - 1)) * t(il, i) / (lvcp(il, i) * sigd &

 Legend:



Removed from v.195
 


changed lines


 
Added in v.198
 Legend:



Removed from v.195
 


changed lines


 
Added in v.198
-Removed from v.195
+Added in v.198

	ViewVC Help
Powered by ViewVC 1.1.21