<p><b>laura@ucar.edu</b> 2011-02-09 09:58:18 -0700 (Wed, 09 Feb 2011)</p><p>added statements for initialization of Thompson cloud microphysics<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/physics_wrf/module_mp_thompson.F
===================================================================
--- branches/atmos_physics/src/core_physics/physics_wrf/module_mp_thompson.F        2011-02-09 16:54:00 UTC (rev 730)
+++ branches/atmos_physics/src/core_physics/physics_wrf/module_mp_thompson.F        2011-02-09 16:58:18 UTC (rev 731)
@@ -728,6 +728,26 @@
          read(11) tnr_racg
          read(11) tnr_gacr
          close(unit=11)
+         write(0,*) '       end read table_qr_acr_qg.dat'
+         write(0,*)
+         write(0,*) '--- end subroutine qr_acr_qg:'
+         write(0,*) 'max tcg_racg =',maxval(tcg_racg)
+         write(0,*) 'min tcg_racg =',minval(tcg_racg)
+      
+         write(0,*) 'max tmr_racg =',maxval(tmr_racg)
+         write(0,*) 'min tmr_racg =',minval(tmr_racg)
+
+         write(0,*) 'max tcr_gacr =',maxval(tcr_gacr)
+         write(0,*) 'min tcr_gacr =',minval(tcr_gacr)
+
+         write(0,*) 'max tmg_gacr =',maxval(tmg_gacr)
+         write(0,*) 'min tmg_gacr =',minval(tmg_gacr)
+
+         write(0,*) 'max tnr_racg =',maxval(tnr_racg)
+         write(0,*) 'min tnr_racg =',minval(tnr_racg)
+
+         write(0,*) 'max tnr_gacr =',maxval(tnr_gacr)
+         write(0,*) 'min tnr_gacr =',minval(tnr_gacr)
 !        DO ig1 = 1, ntb_g1
 !        DO ig  = 1, ntb_g
 !        DO ir1 = 1, ntb_r1
@@ -776,7 +796,38 @@
          read(11) tnr_sacr2
          close(unit=11)
          write(0,*) '       end read table_qr_acr_qs.dat'
+         write(0,*)
+         write(0,*) '--- end subroutine qr_acr_qs:'
+         write(0,*) 'max tcs_racs1 =',maxval(tcs_racs1)
+         write(0,*) 'min tcs_racs1 =',minval(tcs_racs1)      
 
+         write(0,*) 'max tmr_racs1 =',maxval(tmr_racs1)
+         write(0,*) 'min tmr_racs1 =',minval(tmr_racs1)
+
+         write(0,*) 'max tcr_sacr1 =',maxval(tcr_sacr1)
+         write(0,*) 'min tcr_sacr1 =',minval(tcr_sacr1)
+
+         write(0,*) 'max tms_sacr1 =',maxval(tms_sacr1)      
+         write(0,*) 'min tms_sacr1 =',minval(tms_sacr1)      
+
+         write(0,*) 'max tcr_sacr2 =',maxval(tcr_sacr2)
+         write(0,*) 'min tcr_sacr2 =',minval(tcr_sacr2)
+
+         write(0,*) 'max tcr_sacr2 =',maxval(tcr_sacr2)
+         write(0,*) 'min tcr_sacr2 =',minval(tcr_sacr2)
+
+         write(0,*) 'max tnr_racs1 =',maxval(tnr_racs1)
+         write(0,*) 'min tnr_racs1 =',minval(tnr_racs1)
+
+         write(0,*) 'max tnr_racs2 =',maxval(tnr_racs2)      
+         write(0,*) 'min tnr_racs2 =',minval(tnr_racs2)      
+
+         write(0,*) 'max tnr_sacr1 =',maxval(tnr_sacr1)
+         write(0,*) 'min tnr_sacr1 =',minval(tnr_sacr1)      
+
+         write(0,*) 'max tnr_sacr2 =',maxval(tnr_sacr2)
+         write(0,*) 'min tnr_sacr2 =',minval(tnr_sacr2)      
+
 !..Cloud water and rain freezing (Bigg, 1953).
          open(unit=11,file='./LOOKUP_TABLES/table_freezeH2O.dat', &amp;
               form='unformatted',status='old',action='read')
@@ -788,7 +839,26 @@
          read(11) tni_qcfz
          close(unit=11)
          write(0,*) '       end read table_freezeH2O.dat'
+         write(0,*)
+         write(0,*) '--- end subroutine freezeH2O:'
+         write(0,*) 'max tpi_qrfz =',maxval(tpi_qrfz)
+         write(0,*) 'min tpi_qrfz =',minval(tpi_qrfz)      
 
+         write(0,*) 'max tni_qrfz =',maxval(tni_qrfz)
+         write(0,*) 'min tni_qrfz =',minval(tni_qrfz)      
+
+         write(0,*) 'max tpg_qrfz =',maxval(tpg_qrfz)
+         write(0,*) 'min tpg_qrfz =',minval(tpg_qrfz)      
+
+         write(0,*) 'max tnr_qrfz =',maxval(tnr_qrfz)
+         write(0,*) 'min tnr_qrfz =',minval(tnr_qrfz)      
+
+         write(0,*) 'max tpi_qcfz =',maxval(tpi_qcfz)
+         write(0,*) 'min tpi_qcfz =',minval(tpi_qcfz)      
+
+         write(0,*) 'max tni_qcfz =',maxval(tni_qcfz)
+         write(0,*) 'min tni_qcfz =',minval(tni_qcfz)      
+
 !..Conversion of some ice mass into snow category.
          open(unit=11,file='./LOOKUP_TABLES/table_qi_aut_qs.dat', &amp;
               form='unformatted',status='old',action='read')
@@ -797,6 +867,12 @@
          read(11) tni_iaus
          close(unit=11)
          write(0,*) '       end read table_qi_aut_qs.dat'
+         write(0,*) '--- end subroutine qi_aut_qs:'
+         write(0,*) 'max tps_iaus =',maxval(tps_iaus)
+         write(0,*) 'min tps_iaus =',minval(tps_iaus)      
+      
+         write(0,*) 'max tni_iaus =',maxval(tni_iaus)
+         write(0,*) 'min tni_iaus =',minval(tni_iaus)      
 
          write(0,*) '    end pre-calculated look-up tables'
          iiwarm = .false.

</font>
</pre>