<p><b>laura@ucar.edu</b> 2010-12-22 12:28:30 -0700 (Wed, 22 Dec 2010)</p><p>updated physics for long- and short-wave radiation codes<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/Makefile
===================================================================
--- branches/atmos_physics/src/core_physics/Makefile 2010-12-22 19:25:55 UTC (rev 662)
+++ branches/atmos_physics/src/core_physics/Makefile 2010-12-22 19:28:30 UTC (rev 663)
@@ -17,6 +17,8 @@
module_driver_cloudiness.o \
module_driver_convection_deep.o \
module_driver_microphysics.o \
+ module_driver_radiation_lw.o \
+ module_driver_radiation_sw.o \
module_driver_sfclayer.o \
module_driver_pbl.o \
module_physics_aquaplanet.o \
@@ -67,6 +69,17 @@
module_physics_vars.o \
./physics_wrf/module_bl_ysu.o
+module_driver_radiation_lw.o: \
+ module_physics_constants.o \
+ module_physics_rrtmg_lwinit.o \
+ module_physics_vars.o
+
+module_driver_radiation_sw.o: \
+ module_physics_aquaplanet.o \
+ module_physics_constants.o \
+ module_physics_rrtmg_swinit.o \
+ module_physics_vars.o
+
module_driver_sfclayer.o: \
module_physics_aquaplanet.o \
module_physics_constants.o \
@@ -80,9 +93,9 @@
module_physics_init.o: \
module_driver_convection_deep.o \
module_driver_microphysics.o \
- module_driver_sfclayer.o \
- module_physics_rrtmg_lwinit.o \
- module_physics_rrtmg_swinit.o
+ module_driver_radiation_lw.o \
+ module_driver_radiation_sw.o \
+ module_driver_sfclayer.o
module_physics_manager.o: \
module_driver_convection_deep.o \
@@ -142,6 +155,8 @@
module_physics_driver.o: \
module_driver_convection_deep.o \
module_driver_pbl.o \
+ module_driver_radiation_lw.o \
+ module_driver_radiation_sw.o \
module_driver_sfclayer.o \
module_physics_constants.o \
module_physics_interface_nhyd.o \
Modified: branches/atmos_physics/src/core_physics/module_driver_cloudiness.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_cloudiness.F 2010-12-22 19:25:55 UTC (rev 662)
+++ branches/atmos_physics/src/core_physics/module_driver_cloudiness.F 2010-12-22 19:28:30 UTC (rev 663)
@@ -42,7 +42,7 @@
do j = jts,jte
do k = kts,kte
do i = its,ite
- cldfrac_p(i,k,j) = diag_physics % cldfrac % array(i,k)
+ cldfrac_p(i,k,j) = diag_physics % cldfrac % array(k,i)
enddo
enddo
enddo
@@ -61,7 +61,7 @@
do j = jts,jte
do k = kts,kte
do i = its,ite
- diag_physics % cldfrac % array(i,k) = cldfrac_p(i,k,j)
+ diag_physics % cldfrac % array(k,i) = cldfrac_p(i,k,j)
enddo
enddo
enddo
Modified: branches/atmos_physics/src/core_physics/module_physics_control.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_control.F 2010-12-22 19:25:55 UTC (rev 662)
+++ branches/atmos_physics/src/core_physics/module_physics_control.F 2010-12-22 19:28:30 UTC (rev 663)
@@ -135,7 +135,6 @@
config_sfclayer_scheme .eq. 'off') moist_physics = .false.
write(0,*) '--- end subroutine physics_namelist_check:'
- write(0,*)
end subroutine physics_namelist_check
@@ -230,46 +229,6 @@
diag_physics % zol % array(icell) = 0.
diag_physics % znt % array(icell) = 0.
- !.. shortwave radiation:
- diag_physics % gsw % array(icell) = 0.
- diag_physics % swcf % array(icell) = 0.
- diag_physics % swdnb % array(icell) = 0.
- diag_physics % swdnbc % array(icell) = 0.
- diag_physics % swdnt % array(icell) = 0.
- diag_physics % swdntc % array(icell) = 0.
- diag_physics % swupb % array(icell) = 0.
- diag_physics % swupbc % array(icell) = 0.
- diag_physics % swupt % array(icell) = 0.
- diag_physics % swuptc % array(icell) = 0.
- diag_physics % swdnflx % array(icell) = 0.
- diag_physics % swdnflxc % array(icell) = 0.
- diag_physics % swupflx % array(icell) = 0.
- diag_physics % swupflxc % array(icell) = 0.
-
- !.. longwave radiation:
- diag_physics % glw % array(icell) = 0.
- diag_physics % lwcf % array(icell) = 0.
- diag_physics % lwdnb % array(icell) = 0.
- diag_physics % lwdnbc % array(icell) = 0.
- diag_physics % lwdnt % array(icell) = 0.
- diag_physics % lwdntc % array(icell) = 0.
- diag_physics % lwupb % array(icell) = 0.
- diag_physics % lwupbc % array(icell) = 0.
- diag_physics % lwupt % array(icell) = 0.
- diag_physics % lwuptc % array(icell) = 0.
- diag_physics % lwdnflx % array(icell) = 0.
- diag_physics % lwdnflxc % array(icell) = 0.
- diag_physics % lwupflx % array(icell) = 0.
- diag_physics % lwupflxc % array(icell) = 0.
- diag_physics % olrtoa % array(icell) = 0.
-
- !.. surface characteristics:
- diag_physics % emiss % array(icell) = 0.
- diag_physics % snow % array(icell) = 0.
- diag_physics % tsk % array(icell) = 0.
- diag_physics % xice % array(icell) = 0.
- diag_physics % xland % array(icell) = 0.
-
enddo
!2-d variables:
Modified: branches/atmos_physics/src/core_physics/module_physics_init.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_init.F 2010-12-22 19:25:55 UTC (rev 662)
+++ branches/atmos_physics/src/core_physics/module_physics_init.F 2010-12-22 19:28:30 UTC (rev 663)
@@ -10,9 +10,9 @@
use module_driver_convection_deep
use module_driver_microphysics
+ use module_driver_radiation_lw, only: init_radiation_lw
+ use module_driver_radiation_sw, only: init_radiation_sw
use module_driver_sfclayer
- use module_physics_rrtmg_lwinit
- use module_physics_rrtmg_swinit
implicit none
private
@@ -50,13 +50,12 @@
!if(config_radt_sw_scheme.ne.'off') call init_radiation_sw(dminfo)
!initialization of longwave radiation processes:
- if(config_radt_lw_scheme.ne.'off') call rrtmg_lwinit_forMPAS(dminfo)
+ if(config_radt_lw_scheme.ne.'off') call init_radiation_lw(dminfo)
!initialization of shortwave radiation processes:
- if(config_radt_sw_scheme.ne.'off') call rrtmg_swinit_forMPAS(dminfo)
+ if(config_radt_sw_scheme.ne.'off') call init_radiation_sw(dminfo)
- write(0,*) '--- end subroutine physics_init:'
- write(0,*)
+ write(0,*) '--- end subroutine physics_init'
end subroutine physics_init
Modified: branches/atmos_physics/src/core_physics/module_physics_manager.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_manager.F 2010-12-22 19:25:55 UTC (rev 662)
+++ branches/atmos_physics/src/core_physics/module_physics_manager.F 2010-12-22 19:28:30 UTC (rev 663)
@@ -97,6 +97,11 @@
dt_cu = dt_physics !for now.
dt_pbl = dt_physics !for now.
+!longwave radiation scheme:
+ radt_lw_scheme = trim(config_radt_lw_scheme)
+
+!shortwave radiation scheme:
+ radt_sw_scheme = trim(config_radt_sw_scheme)
dt_radtsw = (dt_dyn * config_n_radt_sw) / 60.
!allocation of all physics arrays:
Modified: branches/atmos_physics/src/core_physics/module_physics_rrtmg_lwinit.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_rrtmg_lwinit.F 2010-12-22 19:25:55 UTC (rev 662)
+++ branches/atmos_physics/src/core_physics/module_physics_rrtmg_lwinit.F 2010-12-22 19:28:30 UTC (rev 663)
@@ -12,12 +12,12 @@
implicit none
private
- public:: rrtmg_lwinit_forMPAS
+ public:: rrtmg_initlw_forMPAS
contains
!=============================================================================================
- subroutine rrtmg_lwinit_forMPAS(dminfo)
+ subroutine rrtmg_initlw_forMPAS(dminfo)
!=============================================================================================
!input arguments:
@@ -32,7 +32,7 @@
!flux to heating rate conversion factor.
call rrtmg_lw_ini(cp)
- end subroutine rrtmg_lwinit_forMPAS
+ end subroutine rrtmg_initlw_forMPAS
!=============================================================================================
subroutine rrtmg_lwlookuptable(dminfo)
</font>
</pre>