Skip to content

Commit 1caa539

Browse files
mgdudadouglasjacobsen
mgduda
authored and
douglasjacobsen
committed
Minor updates to MPAS-A:
- modify settings in default namelists - add code to try acquiring physics lookup tables using 'curl' if 'git' and 'svn' aren't available - remove unused variables relhum and qsat from Registry - halt initialization if required meteorological data are not available
1 parent 1b40b47 commit 1caa539

File tree

6 files changed

+59
-18
lines changed

6 files changed

+59
-18
lines changed

‎namelist.input.atmosphere

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,7 @@
11
&nhyd_model
22
config_dt = 450.0
33
config_start_time = '2010-10-23_00:00:00'
4-
config_run_duration = '12:00:00'
4+
config_run_duration = '5_00:00:00'
55
config_number_of_sub_steps = 6
66
config_h_mom_eddy_visc2 = 0.0
77
config_h_mom_eddy_visc4 = 0.0
@@ -25,7 +25,7 @@
2525
config_coef_3rd_order = 0.25
2626
config_epssm = 0.1
2727
config_smdiv = 0.1
28-
config_h_ScaleWithMesh = .false.
28+
config_h_ScaleWithMesh = .true.
2929
config_newpx = .false.
3030
/
3131
config_stop_time = '0000-01-16_00:00:00'
@@ -43,6 +43,7 @@
4343
config_restart_interval = '1_00:00:00'
4444
config_sfc_update_name = 'x1.40962.sfc_update.nc'
4545
config_sfc_update_interval = 'none'
46+
config_hifreq_output_interval = 'none'
4647
config_frames_per_outfile = 1
4748
config_pio_num_iotasks = 0
4849
config_pio_stride = 1

‎namelist.input.init_atmosphere

Lines changed: 6 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -13,15 +13,15 @@
1313
/
1414

1515
&data_sources
16-
config_geog_data_path = '/mmm/users/wrfhelp/WPS_GEOG/'
16+
config_geog_data_path = '/glade/p/work/wrfhelp/WPS_GEOG/'
1717
config_met_prefix = 'CFSR'
1818
config_sfc_prefix = 'SST'
1919
config_fg_interval = 21600
2020
/
2121

2222
&vertical_grid
2323
config_ztop = 30000.0
24-
config_nsmterrain = 2
24+
config_nsmterrain = 1
2525
config_smooth_surfaces = .true.
2626
/
2727

@@ -33,14 +33,15 @@
3333
/
3434

3535
&io
36-
config_input_name = 'x1.40962.static.nc'
37-
config_output_name = 'x1.40962.init.nc'
36+
config_input_name = 'grid.nc'
37+
config_output_name = 'init.nc'
38+
config_sfc_update_name = 'sfc_update.nc'
3839
config_pio_num_iotasks = 0
3940
config_pio_stride = 1
4041
/
4142

4243
&decomposition
43-
config_block_decomp_file_prefix = 'x1.40962.graph.info.part.'
44+
config_block_decomp_file_prefix = 'graph.info.part.'
4445
/
4546

4647
&restart

‎src/core_atmos_physics/checkout_data_files.sh

Lines changed: 26 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -56,6 +56,32 @@ else
5656
fi
5757

5858

59+
#
60+
# Try using 'curl'
61+
#
62+
which curl
63+
if [ $? == 0 ]; then
64+
echo "*** trying curl to obtain WRF physics tables ***"
65+
curl -o master.zip https://codeload.github.com/MPAS-Dev/MPAS-Data/zip/master
66+
if [ $? == 0 ]; then
67+
which unzip
68+
if [ $? == 0 ]; then
69+
unzip master.zip
70+
mv MPAS-Data-master/atmosphere/physics_wrf/files physics_wrf/
71+
rm -rf master.zip MPAS-Data-master
72+
exit 0
73+
else
74+
echo "*** unzip not in path -- unable to unzip WRF physics tables"
75+
rm -f master.zip
76+
fi
77+
else
78+
echo "*** failed to obtain WRF physics tables using curl ***"
79+
fi
80+
else
81+
echo "*** curl not in path ***"
82+
fi
83+
84+
5985
echo "***************************************************"
6086
echo "Unable to obtain WRF physics tables by any means"
6187
echo "***************************************************"

‎src/core_atmosphere/Registry.xml

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -477,8 +477,6 @@
477477
<var name="graupelnc" type="real" dimensions="nCells Time" streams="ro"/>
478478

479479
<var name="precipw" type="real" dimensions="nCells Time" streams="o"/>
480-
<var name="qsat" type="real" dimensions="nVertLevels nCells Time" streams="o"/>
481-
<var name="relhum" type="real" dimensions="nVertLevels nCells Time" streams="o"/>
482480

483481
<!-- ================================================================================================== -->
484482
<!-- ... PARAMETERIZATION OF CONVECTION: -->

‎src/core_init_atmosphere/mpas_init_atm_cases.F

Lines changed: 16 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -2390,8 +2390,11 @@ subroutine init_atm_case_gfs(grid, fg, state, diag, diag_physics, init_case)
23902390
call read_met_init(trim(config_met_prefix), .false., config_start_time(1:13), istatus)
23912391
23922392
if (istatus /= 0) then
2393-
write(0,*) 'Error reading initial met data'
2394-
return
2393+
write(0,*) '********************************************************************************'
2394+
write(0,*) 'Error opening initial meteorological data file '// &
2395+
trim(config_met_prefix)//':'//config_start_time(1:13)
2396+
write(0,*) '********************************************************************************'
2397+
call mpas_dmpar_abort(grid % block % domain % dminfo)
23952398
end if
23962399
23972400
call read_next_met_field(field, istatus)
@@ -2772,8 +2775,11 @@ subroutine init_atm_case_gfs(grid, fg, state, diag, diag_physics, init_case)
27722775
call read_met_init(trim(config_met_prefix), .false., config_start_time(1:13), istatus)
27732776

27742777
if (istatus /= 0) then
2775-
write(0,*) 'Error reading initial met data'
2776-
return
2778+
write(0,*) '********************************************************************************'
2779+
write(0,*) 'Error opening initial meteorological data file '// &
2780+
trim(config_met_prefix)//':'//config_start_time(1:13)
2781+
write(0,*) '********************************************************************************'
2782+
call mpas_dmpar_abort(grid % block % domain % dminfo)
27772783
end if
27782784

27792785
call read_next_met_field(field, istatus)
@@ -2811,8 +2817,11 @@ subroutine init_atm_case_gfs(grid, fg, state, diag, diag_physics, init_case)
28112817
call read_met_init(trim(config_met_prefix), .false., config_start_time(1:13), istatus)
28122818

28132819
if (istatus /= 0) then
2814-
write(0,*) 'Error reading initial met data'
2815-
return
2820+
write(0,*) '********************************************************************************'
2821+
write(0,*) 'Error opening initial meteorological data file '// &
2822+
trim(config_met_prefix)//':'//config_start_time(1:13)
2823+
write(0,*) '********************************************************************************'
2824+
call mpas_dmpar_abort(grid % block % domain % dminfo)
28162825
end if
28172826

28182827
call read_next_met_field(field, istatus)
@@ -3318,7 +3327,7 @@ subroutine init_atm_case_gfs(grid, fg, state, diag, diag_physics, init_case)
33183327
call read_met_init('SEAICE_FRACTIONAL', .true., config_start_time(1:13), istatus)
33193328
33203329
if (istatus /= 0) then
3321-
write(0,*) 'Error reading SEAICE_FRACTIONAL data'
3330+
write(0,*) 'SEAICE_FRACTIONAL file not found...'
33223331
end if
33233332
33243333
if (istatus == 0) then

‎src/core_init_atmosphere/mpas_init_atm_surface.F

Lines changed: 8 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -115,7 +115,10 @@ subroutine interp_sfc_to_MPAS(timeString,mesh,fg,dminfo)
115115
!open intermediate file:
116116
call read_met_init(trim(config_sfc_prefix),.false.,timeString,istatus)
117117
if(istatus /= 0) then
118-
write(0,*) 'Error reading ',trim(config_sfc_prefix)//':'//timeString(1:13)
118+
write(0,*) '********************************************************************************'
119+
write(0,*) 'Error opening surface file '// &
120+
trim(config_sfc_prefix)//':'//timeString(1:13)
121+
write(0,*) '********************************************************************************'
119122
call mpas_dmpar_abort(dminfo)
120123
else
121124
write(0,*) 'Processing file ',trim(config_sfc_prefix)//':'//timeString(1:13)
@@ -143,7 +146,10 @@ subroutine interp_sfc_to_MPAS(timeString,mesh,fg,dminfo)
143146
!read sea-surface temperatures and seaice data. open intermediate file:
144147
call read_met_init(trim(config_sfc_prefix),.false.,timeString(1:13),istatus)
145148
if(istatus /= 0) then
146-
write(0,*) 'Error reading ',trim(config_sfc_prefix)//':'//timeString(1:13)
149+
write(0,*) '********************************************************************************'
150+
write(0,*) 'Error opening surface file '// &
151+
trim(config_sfc_prefix)//':'//timeString(1:13)
152+
write(0,*) '********************************************************************************'
147153
call mpas_dmpar_abort(dminfo)
148154
endif
149155

0 commit comments

Comments
 (0)