<p><b>laura@ucar.edu</b> 2011-09-27 14:41:43 -0600 (Tue, 27 Sep 2011)</p><p>updated physics drivers<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/module_driver_convection_deep.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_convection_deep.F 2011-09-27 20:34:39 UTC (rev 1027)
+++ branches/atmos_physics/src/core_physics/module_driver_convection_deep.F 2011-09-27 20:41:43 UTC (rev 1028)
@@ -91,6 +91,8 @@
integer:: iCell
!---------------------------------------------------------------------------------------------
+ write(0,*)
+ write(0,*) '--- begin convection_deep initialization:'
convection_select: select case(conv_deep_scheme)
@@ -108,6 +110,8 @@
end select convection_select
+ write(0,*) '--- end convection_deep initialization'
+
end subroutine init_convection_deep
!=============================================================================================
@@ -192,19 +196,18 @@
ims = ims , ime = ime , jms = jms , jme = jme , kms = kds , kme = kme , &
its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte &
)
- write(0,*) 'max raincv=',maxval(raincv_p(its:ite,jts:jte))
201 format(i3,i6,1x,l1,5(1x,e15.8))
write(0,*) '--- end subroutine kf_eta_cps:'
! write(0,*)
- write(0,*) '--- deep convection:'
- do j = jts,jte
- do i = its,ite
- if(nca_p(i,j).gt.0. .and. raincv_p(i,j).gt.0.) then
- write(0,201) j,i,cu_act_flag(i,j),nca_p(i,j),raincv_p(i,j), &
- raincv_p(i,j)/dt_dyn,pratec_p(i,j)
- endif
- enddo
- enddo
+! write(0,*) '--- deep convection:'
+! do j = jts,jte
+! do i = its,ite
+! if(nca_p(i,j).gt.0. .and. raincv_p(i,j).gt.0.) then
+! write(0,201) j,i,cu_act_flag(i,j),nca_p(i,j),raincv_p(i,j), &
+! raincv_p(i,j)/dt_dyn,pratec_p(i,j)
+! endif
+! enddo
+! enddo
! write(0,*) '--- shallow convection:'
! do j = jts,jte
! do i = its,ite
Modified: branches/atmos_physics/src/core_physics/module_driver_lsm.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_lsm.F 2011-09-27 20:34:39 UTC (rev 1027)
+++ branches/atmos_physics/src/core_physics/module_driver_lsm.F 2011-09-27 20:41:43 UTC (rev 1028)
@@ -19,13 +19,13 @@
driver_lsm
logical,parameter:: rdxmaxalb = .true. !use snow albedo from geogrid;false use table values
- logical,parameter:: usemonalb = .false. !use climatological albedo instead of table values
logical,parameter:: myj = .false. !true if using Mellor-Yamada PBL scheme.
logical,parameter:: frpcpn = .false.
logical,parameter:: rdlai2d = .false.
- integer,parameter:: isice = 0
- integer,parameter:: isurban = 0
+!urban physics: since MPAS does not plan to run the urban physics option, the two options
+!below are defined locally:
+ integer,parameter,public:: isurban=1
integer,parameter:: sf_urban_physics = 0 !activate urban canopy model (=0: no urban canopy)
integer,private:: i,j,k,n
@@ -222,7 +222,6 @@
qsfc_p(i,j) = diag_physics % qsfc % array(i)
qz0_p(i,j) = diag_physics % qz0 % array(i)
rib_p(i,j) = diag_physics % rib % array(i)
- sfc_albbck_p(i,j) = diag_physics % sfc_albbck % array(i)
sfc_albedo_p(i,j) = diag_physics % sfc_albedo % array(i)
sfc_emibck_p(i,j) = diag_physics % sfc_emibck % array(i)
sfc_emiss_p(i,j) = diag_physics % sfc_emiss % array(i)
@@ -231,24 +230,25 @@
smstot_p(i,j) = diag_physics % smstot % array(i)
snotime_p(i,j) = diag_physics % snotime % array(i)
snopcx_p(i,j) = diag_physics % snopcx % array(i)
- snowh_p(i,j) = diag_physics % snowh % array(i)
sr_p(i,j) = diag_physics % sr % array(i)
udrunoff_p(i,j) = diag_physics % udrunoff % array(i)
- xland_p(i,j) = diag_physics % xland % array(i)
z0_p(i,j) = diag_physics % z0 % array(i)
znt_p(i,j) = diag_physics % znt % array(i)
- isltyp_p(i,j) = sfc_input % isltyp % array(i)
- ivgtyp_p(i,j) = sfc_input % ivgtyp % array(i)
- shdmin_p(i,j) = sfc_input % shdmin % array(i)
- shdmax_p(i,j) = sfc_input % shdmax % array(i)
- snoalb_p(i,j) = sfc_input % snoalb % array(i)
- snow_p(i,j) = sfc_input % snow % array(i)
- snowc_p(i,j) = sfc_input % snowc % array(i)
- tmn_p(i,j) = sfc_input % tmn % array(i)
- tsk_p(i,j) = sfc_input % skintemp % array(i)
- vegfra_p(i,j) = sfc_input % vegfra % array(i)
- xice_p(i,j) = sfc_input % xice % array(i)
+ isltyp_p(i,j) = sfc_input % isltyp % array(i)
+ ivgtyp_p(i,j) = sfc_input % ivgtyp % array(i)
+ shdmin_p(i,j) = sfc_input % shdmin % array(i)
+ shdmax_p(i,j) = sfc_input % shdmax % array(i)
+ snoalb_p(i,j) = sfc_input % snoalb % array(i)
+ sfc_albbck_p(i,j) = sfc_input % sfc_albbck % array(i)
+ snow_p(i,j) = sfc_input % snow % array(i)
+ snowc_p(i,j) = sfc_input % snowc % array(i)
+ snowh_p(i,j) = sfc_input % snowh % array(i)
+ tmn_p(i,j) = sfc_input % tmn % array(i)
+ tsk_p(i,j) = sfc_input % skintemp % array(i)
+ vegfra_p(i,j) = sfc_input % vegfra % array(i)
+ xice_p(i,j) = sfc_input % xice % array(i)
+ xland_p(i,j) = sfc_input % xland % array(i)
enddo
enddo
@@ -260,8 +260,6 @@
/ (1. - diag_physics % sfc_albedo % array(i))
enddo
enddo
- write(0,*) 'max rainbl_p=',maxval(rainbl_p(:,:))
-
end subroutine lsm_from_MPAS
@@ -312,7 +310,6 @@
diag_physics % qsfc % array(i) = qsfc_p(i,j)
diag_physics % qz0 % array(i) = qz0_p(i,j)
diag_physics % rib % array(i) = rib_p(i,j)
- diag_physics % sfc_albbck % array(i) = sfc_albbck_p(i,j)
diag_physics % sfc_albedo % array(i) = sfc_albedo_p(i,j)
diag_physics % sfc_emibck % array(i) = sfc_emibck_p(i,j)
diag_physics % sfc_emiss % array(i) = sfc_emiss_p(i,j)
@@ -321,29 +318,29 @@
diag_physics % smstot % array(i) = smstot_p(i,j)
diag_physics % snotime % array(i) = snotime_p(i,j)
diag_physics % snopcx % array(i) = snopcx_p(i,j)
- diag_physics % snowh % array(i) = snowh_p(i,j)
diag_physics % sr % array(i) = sr_p(i,j)
diag_physics % udrunoff % array(i) = udrunoff_p(i,j)
- diag_physics % xland % array(i) = xland_p(i,j)
diag_physics % z0 % array(i) = z0_p(i,j)
diag_physics % znt % array(i) = znt_p(i,j)
!not needed ?:
- sfc_input % isltyp % array(i) = isltyp_p(i,j)
- sfc_input % ivgtyp % array(i) = ivgtyp_p(i,j)
- sfc_input % shdmin % array(i) = shdmin_p(i,j)
- sfc_input % shdmax % array(i) = shdmax_p(i,j)
- sfc_input % snoalb % array(i) = snoalb_p(i,j)
- sfc_input % snow % array(i) = snow_p(i,j)
- sfc_input % snowc % array(i) = snowc_p(i,j)
- sfc_input % skintemp % array(i) = tsk_p(i,j)
- sfc_input % tmn % array(i) = tmn_p(i,j)
- sfc_input % vegfra % array(i) = vegfra_p(i,j)
- sfc_input % xice % array(i) = xice_p(i,j)
+ sfc_input % isltyp % array(i) = isltyp_p(i,j)
+ sfc_input % ivgtyp % array(i) = ivgtyp_p(i,j)
+ sfc_input % shdmin % array(i) = shdmin_p(i,j)
+ sfc_input % shdmax % array(i) = shdmax_p(i,j)
+ sfc_input % snoalb % array(i) = snoalb_p(i,j)
+ sfc_input % sfc_albbck % array(i) = sfc_albbck_p(i,j)
+ sfc_input % snow % array(i) = snow_p(i,j)
+ sfc_input % snowc % array(i) = snowc_p(i,j)
+ sfc_input % snowh % array(i) = snowh_p(i,j)
+ sfc_input % skintemp % array(i) = tsk_p(i,j)
+ sfc_input % tmn % array(i) = tmn_p(i,j)
+ sfc_input % vegfra % array(i) = vegfra_p(i,j)
+ sfc_input % xice % array(i) = xice_p(i,j)
+ sfc_input % xland % array(i) = xland_p(i,j)
enddo
enddo
- write(0,*) '--- end subroutine lsm_to_MPAS:'
end subroutine lsm_to_MPAS
@@ -360,15 +357,15 @@
type(sfc_input_type),intent(inout):: sfc_input
!---------------------------------------------------------------------------------------------
-
+ write(0,*)
write(0,*) '--- begin land surface model initialization:'
lsm_select: select case (trim(lsm_scheme))
case ("noah")
- write(0,*) ' enter subroutine noah_init_forMPAS:'
+ write(0,*) ' enter subroutine noah_init_forMPAS:'
call noah_init_forMPAS(dminfo,mesh,diag_physics,sfc_input)
- write(0,*) ' end subroutine noah_init_forMPAS'
+ write(0,*) ' end subroutine noah_init_forMPAS'
case default
@@ -391,7 +388,11 @@
type(sfc_input_type),intent(inout) :: sfc_input
!---------------------------------------------------------------------------------------------
+ write(0,*)
write(0,*) '--- enter subroutine driver_lsm:'
+ write(0,*) '--- isice = ',isice
+ write(0,*) '--- iswater = ', iswater
+ write(0,*) '--- isurban = ', isurban
!formats:
101 format(2i6,8(1x,e15.8))
@@ -428,7 +429,7 @@
rovcp = rcp , dt = dt_pbl , myj = myj , &
itimestep = itimestep , frpcpn = frpcpn , rdlai2d = rdlai2d , &
xice_threshold = xice_threshold , &
- usemonalb = input_sfc_albedo , &
+ usemonalb = config_sfc_albedo , &
mminlu = input_landuse_data , &
num_soil_layers = num_soil_layers , &
num_roof_layers = num_soil_layers , &
@@ -447,7 +448,7 @@
!copy all arrays back to MPAS geodesic grid:
call lsm_to_MPAS(mesh,diag_physics,sfc_input)
- write(0,*) '--- end subroutine lsm_to_MPAS:'
+ write(0,*) '--- end subroutine driver_lsm'
end subroutine driver_lsm
Modified: branches/atmos_physics/src/core_physics/module_driver_microphysics.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_microphysics.F 2011-09-27 20:34:39 UTC (rev 1027)
+++ branches/atmos_physics/src/core_physics/module_driver_microphysics.F 2011-09-27 20:41:43 UTC (rev 1028)
@@ -146,6 +146,8 @@
!=============================================================================================
subroutine microphysics_init
!=============================================================================================
+ write(0,*)
+ write(0,*) '--- begin cloud microphysics initialization:'
microp_select: select case(microp_scheme)
@@ -163,6 +165,8 @@
end select microp_select
+ write(0,*) '--- end cloud microphysics initialization'
+
end subroutine microphysics_init
!=============================================================================================
@@ -296,7 +300,7 @@
! dynamics grid:
#ifdef non_hydrostatic_core
- call microphysics_to_MPAS(mesh,state,diag,tend)
+ call microphysics_to_MPAS(mesh,state,diag,tend,itimestep)
#elif hydrostatic_core
call microphysics_to_MPAS(state,tend)
#endif
Modified: branches/atmos_physics/src/core_physics/module_driver_pbl.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_pbl.F 2011-09-27 20:34:39 UTC (rev 1027)
+++ branches/atmos_physics/src/core_physics/module_driver_pbl.F 2011-09-27 20:41:43 UTC (rev 1028)
@@ -81,11 +81,12 @@
end subroutine deallocate_pbl
!=============================================================================================
- subroutine pbl_from_MPAS(diag_physics)
+ subroutine pbl_from_MPAS(sfc_input,diag_physics)
!=============================================================================================
!input arguments:
type(diag_physics_type),intent(in):: diag_physics
+ type(sfc_input_type),intent(in) :: sfc_input
!---------------------------------------------------------------------------------------------
@@ -104,7 +105,8 @@
v10_p(i,j) = diag_physics % v10 % array(i)
wspd_p(i,j) = diag_physics % wspd % array(i)
znt_p(i,j) = diag_physics % znt % array(i)
- xland_p(i,j) = diag_physics % xland % array(i)
+
+ xland_p(i,j) = sfc_input % xland % array(i)
!initialization for YSU PBL scheme:
kpbl_p(i,j) = 1
enddo
@@ -161,20 +163,21 @@
end subroutine pbl_to_MPAS
!=============================================================================================
- subroutine driver_pbl(diag_physics,tend_physics)
+ subroutine driver_pbl(sfc_input,diag_physics,tend_physics)
!=============================================================================================
!input and output arguments:
!---------------------------
+ type(sfc_input_type),intent(inout) :: sfc_input
type(diag_physics_type),intent(inout):: diag_physics
type(tend_physics_type),intent(inout):: tend_physics
!---------------------------------------------------------------------------------------------
write(0,*)
- write(0,*) '--- enter pbl driver: dt_pbl=',dt_pbl
+ write(0,*) '--- enter subroutine driver_pbl: dt_pbl=',dt_pbl
!copy all MPAS arrays to rectanguler grid arrays:
- call pbl_from_MPAS(diag_physics)
+ call pbl_from_MPAS(sfc_input,diag_physics)
pbl_select: select case (trim(pbl_scheme))
@@ -206,6 +209,7 @@
!copy all arrays back to the MPAS grid:
call pbl_to_MPAS(diag_physics,tend_physics)
+ write(0,*) '--- enter subroutine driver_pbl'
end subroutine driver_pbl
Modified: branches/atmos_physics/src/core_physics/module_driver_radiation_lw.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_radiation_lw.F 2011-09-27 20:34:39 UTC (rev 1027)
+++ branches/atmos_physics/src/core_physics/module_driver_radiation_lw.F 2011-09-27 20:41:43 UTC (rev 1028)
@@ -98,10 +98,10 @@
do j = jts,jte
do i = its,ite
sfc_emiss_p(i,j) = diag_physics % sfc_emiss % array(i)
- xland_p(i,j) = diag_physics % xland % array(i)
tsk_p(i,j) = sfc_input % skintemp % array(i)
snow_p(i,j) = sfc_input % snow % array(i)
xice_p(i,j) = sfc_input % xice % array(i)
+ xland_p(i,j) = sfc_input % xland % array(i)
enddo
enddo
do j = jts,jte
@@ -198,13 +198,14 @@
!---------------------------------------------------------------------------------------------
+ write(0,*)
write(0,*) '--- enter radiation_lw initialization:'
radiation_lw_select: select case (trim(radt_lw_scheme))
case ("rrtmg_lw")
- write(0,*) ' enter subroutine rrtmg_lwinit:'
+ write(0,*) ' enter subroutine rrtmg_lwinit:'
call rrtmg_initlw_forMPAS(dminfo)
- write(0,*) ' end subroutine rrtmg_lwinit'
+ write(0,*) ' end subroutine rrtmg_lwinit'
case default
@@ -229,7 +230,7 @@
write(0,100) itimestep
!formats:
- 100 format(/,'--- enter subroutine driver_radiation_lw: ',i6)
+ 100 format(/,' --- enter subroutine driver_radiation_lw: ',i6)
101 format(i8,12(1x,e15.8))
!copy all MPAS arrays to rectangular grid:
@@ -275,6 +276,8 @@
!copy all arrays back to MPAS geodesic grid:
call radiation_lw_to_MPAS(diag_physics,tend_physics)
+ write(0,*) '--- end subroutine driver_radiation_lw:'
+
end subroutine driver_radiation_lw
!=============================================================================================
Modified: branches/atmos_physics/src/core_physics/module_driver_radiation_sw.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_radiation_sw.F 2011-09-27 20:34:39 UTC (rev 1027)
+++ branches/atmos_physics/src/core_physics/module_driver_radiation_sw.F 2011-09-27 20:41:43 UTC (rev 1028)
@@ -2,7 +2,6 @@
module module_driver_radiation_sw
use grid_types
- !use module_physics_aquaplanet !temporary until we have a time manager and sfc conditions.
use module_physics_manager
use module_physics_constants
use module_physics_rrtmg_swinit
@@ -111,10 +110,10 @@
xlon_p(i,j) = (mesh % lonCell % array(i)) / degrad
sfc_albedo_p(i,j) = diag_physics % sfc_albedo % array(i)
- xland_p(i,j) = diag_physics % xland % array(i)
snow_p(i,j) = sfc_input % snow % array(i)
tsk_p(i,j) = sfc_input % skintemp % array(i)
xice_p(i,j) = sfc_input % xice % array(i)
+ xland_p(i,j) = sfc_input % xland % array(i)
enddo
enddo
do j = jts,jte
@@ -211,16 +210,16 @@
type(dm_info), intent(in):: dminfo
!---------------------------------------------------------------------------------------------
-
+ write(0,*)
write(0,*) '--- begin radiation_sw initialization:'
!call to shortwave radiation scheme:
radiation_sw_select: select case (trim(radt_sw_scheme))
case ("rrtmg_sw")
- write(0,*) ' enter subroutine rrtmg_swinit:'
+ write(0,*) ' enter subroutine rrtmg_swinit:'
call rrtmg_initsw_forMPAS(dminfo)
- write(0,*) ' end subroutine rrtmg_swinit'
+ write(0,*) ' end subroutine rrtmg_swinit'
case default
@@ -246,7 +245,6 @@
type(sfc_input_type) ,intent(inout):: sfc_input
type(tend_physics_type),intent(inout):: tend_physics
-
!local variables:
real(kind=RKIND):: radt,xtime_m
@@ -256,18 +254,28 @@
write(0,100) itimestep
!formats:
- 100 format(/,'--- enter subroutine driver_radiation_sw: ',i6)
- 101 format(2i6,8(1x,e15.8))
- 102 format(3i6,8(1x,e15.8))
+ 100 format(/,' --- enter subroutine driver_radiation_sw: ',i6)
+ 101 format(' ITIMESTEP =', i6 ,/, &
+ ' YEAR =', i6 ,/, &
+ ' JULDAY =', i6 ,/, &
+ ' GMT =', f16.9,/, &
+ ' XTIME_M =', f16.9,/, &
+ ' CURR_JULDAY =', f16.9,/, &
+ ' SOLCON =', f16.9,/, &
+ ' DECLIN =', f16.9)
!xtime => state % xtime % scalar
- xtime_m = xtime_s/60.
+ !ldf (2011-09-19): change needed until AdvanceClock is corrected by Conrad.
+ xtime_m = xtime_s/60.
+ !xtime_m = xtime_s/60. - dt_dyn/60.
!copy all MPAS arrays to rectangular grid:
call radiation_sw_from_MPAS(mesh,diag_physics,sfc_input)
!... calculates solar declination:
- call radconst(declin,solcon,julday,degrad,dpd)
+!call radconst(declin,solcon,julday,degrad,dpd)
+ call radconst(declin,solcon,curr_julday,degrad,dpd)
+ write(0,101) itimestep,year,julday,gmt,xtime_m,curr_julday,solcon,declin
!... convert the radiation time_step to minutes:
radt = dt_radtsw/60.
@@ -276,33 +284,10 @@
radiation_sw_select: select case (trim(radt_sw_scheme))
case ("rrtmg_sw")
-! write(0,*) '--- enter subroutine rrtmg_sw:'
-! write(0,*) ' gmt = ', gmt
-! write(0,*) ' dt_radtsw = ', dt_radtsw
-! write(0,*)
-! do j = jts,jte
-! do i = its,its+20
-! write(0,101) j,i,xlon_p(i,j),xlat_p(i,j),xland_p(i,j),xice_p(i,j), &
-! sfc_albedo_p(i,j),tsk_p(i,j)
-! enddo
-! enddo
-! write(0,*)
-! do j=jts,jte
-! do i=its,its
-! do k=kte,kts,-1
-! write(0,102) j,i,k,pres_p(i,k,j),pi_p(i,k,j),dz_p(i,k,j),rho_p(i,k,j),t_p(i,k,j)
-! enddo
-! write(0,*)
-! do k=kte,kts,-1
-! write(0,102) j,i,k,cldfrac_p(i,k,j),qv_p(i,k,j),qc_p(i,k,j),qr_p(i,k,j), &
-! qi_p(i,k,j),qs_p(i,k,j),qg_p(i,k,j)
-! enddo
-! enddo
-! enddo
-
f_ice(:,:) = 0.
f_rain(:,:) = 0.
+ write(0,*) '--- enter subroutine rrtmg_swrad:'
call rrtmg_swrad( &
rthratensw = rthratensw_p , swupt = swupt_p , swuptc = swuptc_p , &
swdnt = swdnt_p , swdntc = swdntc_p , swupb = swupb_p , &
@@ -330,6 +315,7 @@
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,*) '--- exit subroutine rrtmg_swrad'
case default
@@ -339,7 +325,6 @@
call radiation_sw_to_MPAS(diag_physics,tend_physics)
write(0,*) '--- end subroutine driver_radiation_sw:'
- write(0,*)
end subroutine driver_radiation_sw
@@ -348,7 +333,8 @@
!=============================================================================================
!input arguments:
- integer,intent(in):: julian
+!integer,intent(in):: julian
+ real(kind=RKIND),intent(in):: julian
real(kind=RKIND),intent(in):: degrad,dpd
!output arguments:
Modified: branches/atmos_physics/src/core_physics/module_driver_sfclayer.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_sfclayer.F 2011-09-27 20:34:39 UTC (rev 1027)
+++ branches/atmos_physics/src/core_physics/module_driver_sfclayer.F 2011-09-27 20:41:43 UTC (rev 1028)
@@ -128,8 +128,8 @@
area_p(i,j) = mesh % areaCell % array(i)
hpbl_p(i,j) = diag_physics % hpbl % array(i)
mavail_p(i,j) = diag_physics % mavail % array(i)
- xland_p(i,j) = diag_physics % xland % array(i)
tsk_p(i,j) = sfc_input % skintemp % array(i)
+ xland_p(i,j) = sfc_input % xland % array(i)
!inout variables:
hfx_p(i,j) = diag_physics % hfx % array(i)
qfx_p(i,j) = diag_physics % qfx % array(i)
@@ -226,7 +226,7 @@
logical:: allowed_to_read
!---------------------------------------------------------------------------------------------
-
+ write(0,*)
write(0,*) '--- enter sfclayer_initialization:'
sfclayer_select: select case (trim(sfclayer_scheme))
@@ -256,7 +256,7 @@
!---------------------------------------------------------------------------------------------
write(0,*)
- write(0,*) '--- enter sfclayer_driver:'
+ write(0,*) '--- enter subroutine driver_sfclayer:'
!copy all MPAS arrays to rectanguler grid:
call sfclayer_from_MPAS(mesh,diag_physics,sfc_input)
@@ -299,6 +299,8 @@
!copy all arrays back to the MPAS grid:
call sfclayer_to_MPAS(diag_physics)
+ write(0,*) '--- end subroutine driver_sfclayer'
+
end subroutine driver_sfclayer
!=============================================================================================
</font>
</pre>