<p><b>laura@ucar.edu</b> 2011-06-09 10:31:23 -0600 (Thu, 09 Jun 2011)</p><p>uncommented call to kessler cloud microphysics<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/module_driver_microphysics.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_microphysics.F        2011-06-09 16:21:12 UTC (rev 885)
+++ branches/atmos_physics/src/core_physics/module_driver_microphysics.F        2011-06-09 16:31:23 UTC (rev 886)
@@ -214,7 +214,7 @@
!... initialization of soundings for non-hydrostatic or hydrostatic dynamical cores.
#ifdef non_hydrostatic_core
- call microphysics_from_MPAS(mesh,state,diag)
+ call microphysics_from_MPAS(mesh,state,tend,diag)
#elif hydrostatic_core
call microphysics_from_MPAS(state)
#endif
@@ -233,14 +233,18 @@
case ("kessler")
-! call kessler( &
-! th_p , qv_p , qc_p , qr_p , rho_p , pi_p , &
-! dt_microp , z_phy , xlv , cp , ep_2 , svp1 , &
-! svp2 , svp3 , svpt0 , rho_w , dz_phy , rainnc_phy , &
-! rainncv_phy , &
-! ids,ide,jds,jde,kds,kde, & ! domain dimensions
-! ims,ime,jms,jme,kms,kme, & ! memory dimensions
-! its,ite,jts,jte,kts,kte) ! tile dimensions
+ call kessler( &
+ t = th_p , qv = qv_p , qc = qc_p , &
+ qr = qr_p , rho = rho_p , pii = pi_p , &
+ dt_in = dt_microp , z = z_p , xlv = xlv , &
+ cp = cp , ep2 = ep_2 , svp1 = svp1 , &
+ svp2 = svp2 , svp3 = svp3 , svpt0 = svpt0 , &
+ rhowater = rho_w , dz8w = dz_p , rainnc = rainnc_p , &
+ rainncv = rainncv_p , &
+ ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde , &
+ ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme , &
+ its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte &
+ )
write(0,*) '--- end kessler:',istep
case ("thompson")
@@ -257,12 +261,6 @@
write(0,*) '--- end thompson:',istep
case ("wsm6")
-! write(0,*) '--- max qv=',maxval(qv_p(:,:,:))
-! write(0,*) '--- max qc=',maxval(qc_p(:,:,:))
-! write(0,*) '--- max qr=',maxval(qr_p(:,:,:))
-! write(0,*) '--- max qi=',maxval(qi_p(:,:,:))
-! write(0,*) '--- max qs=',maxval(qs_p(:,:,:))
-! write(0,*) '--- max qg=',maxval(qg_p(:,:,:))
call wsm6( &
th = th_p , q = qv_p , qc = qc_p , &
qr = qr_p , qi = qi_p , qs = qs_p , &
@@ -280,11 +278,6 @@
ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme , &
its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte &
)
-! do j = jts,jte
-! do i = its,ite
-! if(rainncv_p(i,j).gt.0.) write(0,200) j,i,rainncv_p(i,j),rainnc_p(i,j)
-! enddo
-! enddo
write(0,*) '--- end wsm6:',istep
case default
</font>
</pre>