<p><b>duda</b> 2012-07-26 12:06:33 -0600 (Thu, 26 Jul 2012)</p><p>BRANCH COMMIT<br>
<br>
Add new namelist option config_rotatation_rate_scale to scale the <br>
planet rotation rate for the J&W baroclinic wave case.<br>
<br>
<br>
M namelist.input.init_nhyd_atmos<br>
M src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F<br>
M src/core_init_nhyd_atmos/Registry<br>
</p><hr noshade><pre><font color="gray">Modified: branches/dcmip/namelist.input.init_nhyd_atmos
===================================================================
--- branches/dcmip/namelist.input.init_nhyd_atmos        2012-07-26 17:36:53 UTC (rev 2059)
+++ branches/dcmip/namelist.input.init_nhyd_atmos        2012-07-26 18:06:33 UTC (rev 2060)
@@ -1,15 +1,17 @@
&nhyd_model
+ config_test_case = 9
config_start_time = '0000-01-01_00:00:00'
config_theta_adv_order = 3
/
&dcmip
- config_dcmip_case = '3-1'
- config_planet_scale = 500.0
+ config_dcmip_case = '3-1'
+ config_planet_scale = 500.0
+ config_rotation_rate_scale = 500.0
/
&dimensions
- config_nvertlevels = 41
+ config_nvertlevels = 30
/
&data_sources
Modified: branches/dcmip/src/core_init_nhyd_atmos/Registry
===================================================================
--- branches/dcmip/src/core_init_nhyd_atmos/Registry        2012-07-26 17:36:53 UTC (rev 2059)
+++ branches/dcmip/src/core_init_nhyd_atmos/Registry        2012-07-26 18:06:33 UTC (rev 2060)
@@ -9,6 +9,7 @@
namelist real nhyd_model config_coef_3rd_order 0.25
namelist character dcmip config_dcmip_case 2-0-0
namelist real dcmip config_planet_scale 1.0
+namelist real dcmip config_rotation_rate_scale 1.0
namelist integer dimensions config_nvertlevels 26
namelist integer dimensions config_nsoillevels 4
namelist integer dimensions config_nfglevels 27
Modified: branches/dcmip/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F
===================================================================
--- branches/dcmip/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-07-26 17:36:53 UTC (rev 2059)
+++ branches/dcmip/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-07-26 18:06:33 UTC (rev 2060)
@@ -362,7 +362,7 @@
etavs = (1.-0.252)*pii/2.
r_earth = grid % sphere_radius
- omega_e = omega
+ omega_e = omega * config_rotation_rate_scale
p0 = 1.e+05
write(0,*) ' point 1 in test case setup '
@@ -806,14 +806,14 @@
! Generate rotated Coriolis field
!
- grid % fEdge % array(iEdge) = 2.0 * omega * &
+ grid % fEdge % array(iEdge) = 2.0 * omega_e * &
( -cos(grid%lonEdge%array(iEdge)) * cos(grid%latEdge%array(iEdge)) * sin(alpha_grid) + &
sin(grid%latEdge%array(iEdge)) * cos(alpha_grid) &
)
end do
do iVtx=1,grid % nVertices
- grid % fVertex % array(iVtx) = 2.0 * omega * &
+ grid % fVertex % array(iVtx) = 2.0 * omega_e * &
(-cos(grid%lonVertex%array(iVtx)) * cos(grid%latVertex%array(iVtx)) * sin(alpha_grid) + &
sin(grid%latVertex%array(iVtx)) * cos(alpha_grid) &
)
@@ -1017,7 +1017,7 @@
integer :: k,i, itr
r_earth = rad
- omega_e = omega
+ omega_e = omega * config_rotation_rate_scale
rdx = 1./(dlat*r_earth)
do i=1,nlat-1
</font>
</pre>