<p><b>laura@ucar.edu</b> 2010-10-14 15:39:11 -0600 (Thu, 14 Oct 2010)</p><p>miscellaneous updates with revised dynamical core<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/module_physics_control.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_control.F 2010-10-14 18:25:59 UTC (rev 555)
+++ branches/atmos_physics/src/core_physics/module_physics_control.F 2010-10-14 21:39:11 UTC (rev 556)
@@ -135,6 +135,7 @@
config_sfclayer_scheme .eq. 'off') moist_physics = .false.
write(0,*) '--- end subroutine physics_namelist_check:'
+ write(0,*)
end subroutine physics_namelist_check
Modified: branches/atmos_physics/src/core_physics/module_physics_init.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_init.F 2010-10-14 18:25:59 UTC (rev 555)
+++ branches/atmos_physics/src/core_physics/module_physics_init.F 2010-10-14 21:39:11 UTC (rev 556)
@@ -26,6 +26,8 @@
type(state_type),intent(inout):: s
!=============================================================================================
+ write(0,*)
+ write(0,*) '--- enter subroutine physics_init:'
!initialization of parameterized deep convective processes:
if(config_conv_deep_scheme .ne. 'off') then
@@ -42,6 +44,7 @@
write(0,*) '--- end subroutine sfclayer_interface_init:'
endif
+ write(0,*) '--- end subroutine physics_init:'
write(0,*)
end subroutine physics_init
Modified: branches/atmos_physics/src/core_physics/module_physics_interface_hyd.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_interface_hyd.F 2010-10-14 18:25:59 UTC (rev 555)
+++ branches/atmos_physics/src/core_physics/module_physics_interface_hyd.F 2010-10-14 21:39:11 UTC (rev 556)
@@ -151,7 +151,7 @@
!=============================================================================================
!output variables:
- type(state_type),intent(out):: state
+ type(state_type),intent(inout):: state
!---------------------------------------------------------------------------------------------
Modified: branches/atmos_physics/src/core_physics/module_physics_interface_nhyd.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_interface_nhyd.F 2010-10-14 18:25:59 UTC (rev 555)
+++ branches/atmos_physics/src/core_physics/module_physics_interface_nhyd.F 2010-10-14 21:39:11 UTC (rev 556)
@@ -206,22 +206,16 @@
z_p(i,k,j) = 0.5*(zgrid(k+1,i) + zgrid(k,i))
dz_p(i,k,j) = zgrid(k+1,i) - zgrid(k,i)
- qv_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qv,k,i))
- qc_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qc,k,i))
- qr_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qr,k,i))
+! qv_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qv,k,i))
+! qc_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qc,k,i))
+! qr_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qr,k,i))
+ qv_p(i,k,j) = vars % scalars % array(vars%index_qv,k,i)
+ qc_p(i,k,j) = vars % scalars % array(vars%index_qc,k,i)
+ qr_p(i,k,j) = vars % scalars % array(vars%index_qr,k,i)
enddo
enddo
enddo
-!j=1
-!do i = 5000,5000
-!do k = kme, kms, -1
-! write(0,201) j,i,k,pres_p(i,k,j),pi_p(i,k,j),rho_p(i,k,j),th_p(i,k,j), &
-! qv_p(i,k,j),qc_p(i,k,j),qr_p(i,k,j)
-!enddo
-!enddo
-!write(0,*)
-
!additional initialization as function of cloud microphysics scheme:
microp_select_init: select case(microp_scheme)
@@ -232,13 +226,18 @@
do i = its, itf
!mass mixing ratios:
- qi_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qi,k,i))
- qs_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qs,k,i))
- qg_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qg,k,i))
+! qi_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qi,k,i))
+! qs_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qs,k,i))
+! qg_p(i,k,j) = max(0.,vars % scalars % array(vars%index_qg,k,i))
+ qi_p(i,k,j) = vars % scalars % array(vars%index_qi,k,i)
+ qs_p(i,k,j) = vars % scalars % array(vars%index_qs,k,i)
+ qg_p(i,k,j) = vars % scalars % array(vars%index_qg,k,i)
!number concentrations:
- qnr_p(i,k,j) = max(0., vars % scalars % array(vars%index_qnr,k,i))
- qni_p(i,k,j) = max(0., vars % scalars % array(vars%index_qni,k,i))
+! qnr_p(i,k,j) = max(0., vars % scalars % array(vars%index_qnr,k,i))
+! qni_p(i,k,j) = max(0., vars % scalars % array(vars%index_qni,k,i))
+ qnr_p(i,k,j) = vars % scalars % array(vars%index_qnr,k,i)
+ qni_p(i,k,j) = vars % scalars % array(vars%index_qni,k,i)
enddo
enddo
@@ -248,62 +247,6 @@
end select microp_select_init
- max_theta = maxval(vars%theta%array(kts:ktf,its:itf))
- min_theta = minval(vars%theta%array(kts:ktf,its:itf))
- write(0,*) ' max_theta = ',max_theta,' min_theta = ',min_theta
-
- max_theta = maxval(vars%theta%array(:,:))
- min_theta = minval(vars%theta%array(:,:))
- write(0,*) ' max_theta = ',max_theta,' min_theta = ',min_theta
-
- max_thp = maxval(th_p(its:itf,kts:ktf,jts:jtf))
- min_thp = minval(th_p(its:itf,kts:ktf,jts:jtf))
- write(0,*) ' max_thp = ',max_thp ,' min_thp = ',min_thp
-
- max_thp = maxval(th_p(:,:,:))
- min_thp = minval(th_p(:,:,:))
- write(0,*) ' max_thp = ',max_thp ,' min_thp = ',min_thp
-
- max_tp = maxval(t_p(its:itf,kts:ktf,jts:jtf))
- min_tp = minval(t_p(its:itf,kts:ktf,jts:jtf))
- write(0,*) ' max_tp = ',max_tp ,' min_tp = ',min_tp
-
- max_tp = maxval(t_p(:,:,:))
- min_tp = minval(t_p(:,:,:))
- write(0,*) ' max_tp = ',max_tp ,' min_tp = ',min_tp
-
- max_qv = maxval(qv_p(:,:,:))
- min_qv = minval(qv_p(:,:,:))
- write(0,*) ' max_qv = ',max_qv,' min_qv = ',min_qv
-
- max_qc = maxval(qc_p(:,:,:))
- min_qc = minval(qc_p(:,:,:))
- write(0,*) ' max_qc = ',max_qc,' min_qc = ',min_qc
-
- max_qr = maxval(qr_p(:,:,:))
- min_qr = minval(qr_p(:,:,:))
- write(0,*) ' max_qr = ',max_qr,' min_qr = ',min_qr
-
- max_qi = maxval(qi_p(:,:,:))
- min_qi = minval(qi_p(:,:,:))
- write(0,*) ' max_qi = ',max_qi,' min_qi = ',min_qi
-
- max_qs = maxval(qs_p(:,:,:))
- min_qs = minval(qs_p(:,:,:))
- write(0,*) ' max_qs = ',max_qs,' min_qs = ',min_qs
-
- max_qg = maxval(qg_p(:,:,:))
- min_qg = minval(qg_p(:,:,:))
- write(0,*) ' max_qg = ',max_qg,' min_qg = ',min_qg
-
- max_qnr = maxval(qnr_p(:,:,:))
- min_qnr = minval(qnr_p(:,:,:))
- write(0,*) ' max_qnr = ',max_qnr,' min_qnr = ',min_qnr
-
- max_qni = maxval(qni_p(:,:,:))
- min_qni = minval(qni_p(:,:,:))
- write(0,*) ' max_qni = ',max_qni,' min_qni = ',min_qni
-
!formats:
201 format(3i6,10(1x,e15.8))
@@ -318,7 +261,7 @@
type(mesh_type),intent(in):: grid
!output variables:
- type(state_type),intent(out):: vars
+ type(state_type),intent(inout):: vars
real(kind=RKIND):: min_theta,min_thp,min_tp
real(kind=RKIND):: max_theta,max_thp,max_tp
@@ -409,62 +352,6 @@
end select microp_select_init
- max_theta = maxval(vars%theta%array(kts:ktf,its:itf))
- min_theta = minval(vars%theta%array(kts:ktf,its:itf))
- write(0,*) ' max_theta = ',max_theta,' min_theta = ',min_theta
-
- max_theta = maxval(vars%theta%array(:,:))
- min_theta = minval(vars%theta%array(:,:))
- write(0,*) ' max_theta = ',max_theta,' min_theta = ',min_theta
-
- max_thp = maxval(th_p(its:itf,kts:ktf,jts:jtf))
- min_thp = minval(th_p(its:itf,kts:ktf,jts:jtf))
- write(0,*) ' max_thp = ',max_thp ,' min_thp = ',min_thp
-
- max_thp = maxval(th_p(:,:,:))
- min_thp = minval(th_p(:,:,:))
- write(0,*) ' max_thp = ',max_thp ,' min_thp = ',min_thp
-
- max_tp = maxval(t_p(its:itf,kts:ktf,jts:jtf))
- min_tp = minval(t_p(its:itf,kts:ktf,jts:jtf))
- write(0,*) ' max_tp = ',max_tp ,' min_tp = ',min_tp
-
- max_tp = maxval(t_p(:,:,:))
- min_tp = minval(t_p(:,:,:))
- write(0,*) ' max_tp = ',max_tp ,' min_tp = ',min_tp
-
- max_qv = maxval(vars%scalars%array(vars%index_qv,:,:))
- min_qv = minval(vars%scalars%array(vars%index_qv,:,:))
- write(0,*) ' max_qv = ',max_qv,' min_qv = ',min_qv
-
- max_qc = maxval(vars%scalars%array(vars%index_qc,:,:))
- min_qc = minval(vars%scalars%array(vars%index_qc,:,:))
- write(0,*) ' max_qc = ',max_qc,' min_qc = ',min_qc
-
- max_qr = maxval(vars%scalars%array(vars%index_qr,:,:))
- min_qr = minval(vars%scalars%array(vars%index_qr,:,:))
- write(0,*) ' max_qr = ',max_qr,' min_qr = ',min_qr
-
- max_qi = maxval(vars%scalars%array(vars%index_qi,:,:))
- min_qi = minval(vars%scalars%array(vars%index_qi,:,:))
- write(0,*) ' max_qi = ',max_qi,' min_qi = ',min_qi
-
- max_qs = maxval(vars%scalars%array(vars%index_qs,:,:))
- min_qs = minval(vars%scalars%array(vars%index_qs,:,:))
- write(0,*) ' max_qs = ',max_qs,' min_qs = ',min_qs
-
- max_qg = maxval(vars%scalars%array(vars%index_qg,:,:))
- min_qg = minval(vars%scalars%array(vars%index_qg,:,:))
- write(0,*) ' max_qg = ',max_qg,' min_qg = ',min_qg
-
- max_qnr = maxval(vars%scalars%array(vars%index_qnr,:,:))
- min_qnr = minval(vars%scalars%array(vars%index_qnr,:,:))
- write(0,*) ' max_qnr = ',max_qnr,' min_qnr = ',min_qnr
-
- max_qni = maxval(vars%scalars%array(vars%index_qni,:,:))
- min_qni = minval(vars%scalars%array(vars%index_qni,:,:))
- write(0,*) ' max_qni = ',max_qni,' min_qni = ',min_qni
-
!formats:
201 format(3i6,10(1x,e15.8))
</font>
</pre>