<p><b>duda</b> 2011-12-07 14:42:34 -0700 (Wed, 07 Dec 2011)</p><p>Fix restartability with physics, plus other minor fixes, in non-hydrostatic physics.<br>
(From branches/atmos_physics r1206, r1207, and r1208)<br>
<br>
<br>
M src/core_atmos_physics/mpas_atmphys_date_time.F<br>
M src/core_atmos_physics/mpas_atmphys_init.F<br>
M src/core_atmos_physics/mpas_atmphys_driver_lsm.F<br>
M src/core_atmos_physics/mpas_atmphys_landuse.F<br>
</p><hr noshade><pre><font color="gray">Modified: trunk/mpas/src/core_atmos_physics/mpas_atmphys_date_time.F
===================================================================
--- trunk/mpas/src/core_atmos_physics/mpas_atmphys_date_time.F        2011-12-07 21:35:44 UTC (rev 1244)
+++ trunk/mpas/src/core_atmos_physics/mpas_atmphys_date_time.F        2011-12-07 21:42:34 UTC (rev 1245)
@@ -3,8 +3,7 @@
implicit none
private
- public:: get_current_date, &
- get_julgmt, &
+ public:: get_julgmt, &
monthly_interp_to_date, &
monthly_min_max
@@ -13,17 +12,6 @@
contains
!=============================================================================================
- subroutine get_current_date
-!=============================================================================================
-
-!temporary subroutine to provide the date at which the forecast begins.right now, that date is
-!hard-wired. This will change when the time manager is in place.
-
- current_date = '2010-10-23_00'
-
- end subroutine get_current_date
-
-!=============================================================================================
subroutine get_julgmt(date_str,julyr,julday,gmt)
!=============================================================================================
Modified: trunk/mpas/src/core_atmos_physics/mpas_atmphys_driver_lsm.F
===================================================================
--- trunk/mpas/src/core_atmos_physics/mpas_atmphys_driver_lsm.F        2011-12-07 21:35:44 UTC (rev 1244)
+++ trunk/mpas/src/core_atmos_physics/mpas_atmphys_driver_lsm.F        2011-12-07 21:42:34 UTC (rev 1245)
@@ -18,10 +18,10 @@
deallocate_lsm, &
driver_lsm
- logical,parameter:: rdxmaxalb = .true. !use snow albedo from geogrid;false use table values
- logical,parameter:: myj = .false. !true if using Mellor-Yamada PBL scheme.
- logical,parameter:: frpcpn = .false.
- logical,parameter:: rdlai2d = .false.
+ logical,parameter:: rdmaxalb = .false. !use snow albedo from geogrid;false use table values
+ logical,parameter:: myj = .false. !true if using Mellor-Yamada PBL scheme.
+ logical,parameter:: frpcpn = .false.
+ logical,parameter:: rdlai2d = .false.
!urban physics: since MPAS does not plan to run the urban physics option, the two options
!below are defined locally:
Modified: trunk/mpas/src/core_atmos_physics/mpas_atmphys_init.F
===================================================================
--- trunk/mpas/src/core_atmos_physics/mpas_atmphys_init.F        2011-12-07 21:35:44 UTC (rev 1244)
+++ trunk/mpas/src/core_atmos_physics/mpas_atmphys_init.F        2011-12-07 21:42:34 UTC (rev 1245)
@@ -73,8 +73,7 @@
!initialization of global surface properties. set here for now, but may be moved when time
!manager is implemented:
- if(.not. config_do_restart) &
- call landuse_init_forMPAS(dminfo,julday,mesh,diag_physics,sfc_input)
+ call landuse_init_forMPAS(dminfo,julday,mesh,diag_physics,sfc_input)
!initialization of parameterized deep convective processes:
if(config_conv_deep_scheme .ne. 'off') &
Modified: trunk/mpas/src/core_atmos_physics/mpas_atmphys_landuse.F
===================================================================
--- trunk/mpas/src/core_atmos_physics/mpas_atmphys_landuse.F        2011-12-07 21:35:44 UTC (rev 1244)
+++ trunk/mpas/src/core_atmos_physics/mpas_atmphys_landuse.F        2011-12-07 21:42:34 UTC (rev 1245)
@@ -6,7 +6,8 @@
module mpas_atmphys_landuse
use mpas_configure,only: input_landuse_data, &
config_sfc_albedo, &
- config_frac_seaice
+ config_frac_seaice, &
+ config_do_restart
use mpas_dmpar
use mpas_grid_types
@@ -184,6 +185,11 @@
DM_BCAST_MACRO(scfx)
write(0,*) '--- isice =',isice
write(0,*) '--- iswater =',iswater
+ if(config_do_restart) then
+ write(0,*) '--- config_do_restart =', config_do_restart
+ write(0,*) '--- skip the end of landuse_init_forMPAS'
+ return
+ endif
!defines the surface properties over the entire domain:
do iCell = 1, nCells
</font>
</pre>