<p><b>duda</b> 2012-02-06 11:57:28 -0700 (Mon, 06 Feb 2012)</p><p>BRANCH COMMIT<br>
<br>
Change &quot;config_sst_prefix&quot; to &quot;config_sfc_prefix&quot;, since the file referred to<br>
by this variable contains more than just SST.<br>
<br>
Also, add code to handle polar stereographic and Gaussian grids in physics_initialize_real.<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>
M    src/core_atmos_physics/mpas_atmphys_initialize_real.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/namelist.input.init_nhyd_atmos
===================================================================
--- branches/atmos_physics/namelist.input.init_nhyd_atmos        2012-02-06 17:35:05 UTC (rev 1466)
+++ branches/atmos_physics/namelist.input.init_nhyd_atmos        2012-02-06 18:57:28 UTC (rev 1467)
@@ -15,7 +15,7 @@
 &amp;data_sources
    config_geog_data_path  = '/mmm/users/wrfhelp/WPS_GEOG/'
    config_met_prefix      = 'CFSR'
-   config_sst_prefix      = 'SST'
+   config_sfc_prefix      = 'SST'
    config_fg_interval     = 21600
 /
 

Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_initialize_real.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_initialize_real.F        2012-02-06 17:35:05 UTC (rev 1466)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_initialize_real.F        2012-02-06 18:57:28 UTC (rev 1467)
@@ -6,7 +6,7 @@
                            config_input_sst,   &amp;
                            config_nsoillevels, &amp;
                            config_start_time,  &amp;
-                           config_sst_prefix
+                           config_sfc_prefix
  use mpas_grid_types
  use init_atm_hinterp
  use init_atm_llxy
@@ -52,10 +52,10 @@
  interp_list(3) = 0
 
 !open intermediate file:
- call read_met_init(trim(config_sst_prefix),.false.,config_start_time(1:13),istatus)
+ call read_met_init(trim(config_sfc_prefix),.false.,config_start_time(1:13),istatus)
  if(istatus /= 0) &amp;
-    write(0,*) 'Error reading ',trim(config_sst_prefix)//':'//config_start_time(1:13)
- write(0,*) 'Processing ',trim(config_sst_prefix)//':'//config_start_time(1:13)
+    write(0,*) 'Error reading ',trim(config_sfc_prefix)//':'//config_start_time(1:13)
+ write(0,*) 'Processing ',trim(config_sfc_prefix)//':'//config_start_time(1:13)
 
 !scan through all the fields in the file:
  call read_next_met_field(field,istatus)
@@ -89,6 +89,22 @@
                        knownj = 1.0_RKIND, &amp;
                        lat1 = real(field % startlat,RKIND), &amp;
                        lon1 = real(field % startlon,RKIND))
+       else if (field % iproj == PROJ_GAUSS) then
+          call map_set(PROJ_GAUSS, proj, &amp;
+                       nlat = nint(field % deltalat), &amp;
+                       loninc = real(field % deltalon,RKIND), &amp;
+                       lat1 = real(field % startlat,RKIND), &amp;
+                       lon1 = real(field % startlon,RKIND))
+!                       nxmax = nint(360.0 / field % deltalon), &amp;
+       else if (field % iproj == PROJ_PS) then
+          call map_set(PROJ_PS, proj, &amp;
+                       dx = real(field % dx,RKIND), &amp;
+                       truelat1 = real(field % truelat1,RKIND), &amp;
+                       stdlon = real(field % xlonc,RKIND), &amp;
+                       knowni = real(field % nx / 2.0,RKIND), &amp;
+                       knownj = real(field % ny / 2.0,RKIND), &amp;
+                       lat1 = real(field % startlat,RKIND), &amp;
+                       lon1 = real(field % startlon,RKIND))
        end if
    
        !Interpolate field to each MPAS grid cell:

Modified: branches/atmos_physics/src/core_init_nhyd_atmos/Registry
===================================================================
--- branches/atmos_physics/src/core_init_nhyd_atmos/Registry        2012-02-06 17:35:05 UTC (rev 1466)
+++ branches/atmos_physics/src/core_init_nhyd_atmos/Registry        2012-02-06 18:57:28 UTC (rev 1467)
@@ -14,7 +14,7 @@
 namelist integer   dimensions config_months               12
 namelist character data_sources config_geog_data_path     /data3/mp/wrfhelp/WPS_GEOG/
 namelist character data_sources config_met_prefix         FILE
-namelist character data_sources config_sst_prefix         FILE
+namelist character data_sources config_sfc_prefix         FILE
 namelist integer   data_sources config_fg_interval        21600
 namelist real      vertical_grid  config_ztop             24000.0
 namelist integer   vertical_grid  config_nsmterrain       2

Modified: branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-02-06 17:35:05 UTC (rev 1466)
+++ branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-02-06 18:57:28 UTC (rev 1467)
@@ -4346,12 +4346,12 @@
       curr_time = mpas_get_clock_time(fg_clock, MPAS_NOW) 
       do while (curr_time &lt;= stop_time)
          call mpas_get_time(curr_time, dateTimeString=timeString)
-         write(0,*) 'Processing ',trim(config_sst_prefix)//':'//timeString(1:13)
+         write(0,*) 'Processing ',trim(config_sfc_prefix)//':'//timeString(1:13)
 
          ! Open intermediate file
-         call read_met_init(trim(config_sst_prefix), .false., timeString(1:13), istatus)
+         call read_met_init(trim(config_sfc_prefix), .false., timeString(1:13), istatus)
          if (istatus /= 0) then
-            write(0,*) 'Error reading ',trim(config_sst_prefix)//':'//timeString(1:13)
+            write(0,*) 'Error reading ',trim(config_sfc_prefix)//':'//timeString(1:13)
             exit
          end if
 

</font>
</pre>