<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&amp;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 @@
 &amp;nhyd_model
+   config_test_case       = 9
    config_start_time      = '0000-01-01_00:00:00'
    config_theta_adv_order = 3
 /
 
 &amp;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
 /
 
 &amp;dimensions
-   config_nvertlevels     = 41
+   config_nvertlevels     = 30
 /
 
 &amp;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 * &amp;
+         grid % fEdge % array(iEdge) = 2.0 * omega_e * &amp;
                                        ( -cos(grid%lonEdge%array(iEdge)) * cos(grid%latEdge%array(iEdge)) * sin(alpha_grid) + &amp;
                                          sin(grid%latEdge%array(iEdge)) * cos(alpha_grid) &amp;
                                        )
       end do
 
       do iVtx=1,grid % nVertices
-         grid % fVertex % array(iVtx) = 2.0 * omega * &amp;
+         grid % fVertex % array(iVtx) = 2.0 * omega_e * &amp;
                                          (-cos(grid%lonVertex%array(iVtx)) * cos(grid%latVertex%array(iVtx)) * sin(alpha_grid) + &amp;
                                           sin(grid%latVertex%array(iVtx)) * cos(alpha_grid) &amp;
                                          )
@@ -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>