Running MPAS

Everything you need

Contents for today

• Generate IC from LENS2
• Get PR Mesh
• Static Field init
• Interpolate init
• Run
• Post-processing

Get IC from LENS2 – process_ihesp_data.csh

• /glade/u/home/kenyan/mpas_lens/debug_ihesp_all_hourly.csh
• Choose among
• /glade/u/home/kenyan/mpas_lens/edit_nc_files.py
• /glade/u/home/kenyan/mpas_lens/convert_cesm_hybrid_nc_to_pressure_int_1month.ncl
• /glade/u/home/kenyan/mpas_lens/output_ihesp_intmd.ncl

For py file

For ncl file

• SST…

Comments – h2 files

• hour_6 near-surface: U10 / V10 / T2 / Q2 / PSFC  — from d651 cam.h2
• V10 missing!
• But the py file calculates u10 v10 t2 q2 instead

​

• Why use monthly surface pressure instead of 6 hourly? – 1b vs 4

​

• Conflicting Snow Water Equivalent file – monthly atm: SNOWHLND; monthly land: SNOWLIQ ?

​

• Fix sea ice?

​

• SST remapping….

Results!

Get Mesh – 60km-3km

• ln -s $...meshes/x20.835586.grid.nc .
• /…/grid_rotate x20.835586.nc PR.grid.nc

​

​

Get init 1 – static (sec 3.2)

• In streams.init_atmosphere:

<immutable_stream name="input"

type="input"

filename_template=“PR.grid.nc"

input_interval="initial_only"/>

​

<immutable_stream name="output"

type="output"

filename_template=“PR.static.nc"

packages="initial_conds"

output_interval="initial_only" />

​

​

• In streams.init_atmosphere:

/

&io

config_pio_num_iotasks = 0

config_pio_stride = 1

/

&decomposition

config_block_decomp_file_prefix = 'x20.835586.graph.info.part.'

/

​

​

Results! – plot terrain (sec 3.4 - PR.static.nc)

Get init 2 – interpolation (sec 3.4)

ln -s ${MPAS_TUTORIAL}/met_data/GFS:2014-09-10_00 .

Edit namelist.init_atmosphere

Results!

• TODO: use pywinter to read intermediate files..?

Run

$ ln -s ${HOME}/MPAS-Model/atmosphere_model .

$ ln -s ${HOME}/MPAS-Model/src/core_atmosphere/physics/physics_wrf/files/* .

​

$ cp ${HOME}/MPAS-Model/namelist.atmosphere .

$ cp ${HOME}/MPAS-Model/streams.atmosphere .

$ cp ${HOME}/MPAS-Model/stream_list.atmosphere.* .

​

ln -s ${MPAS_TUTORIAL}/meshes/ x20.835586.graph.info.part.64 .

​

Edit namelist.atmosphere, streams.atmosphere

​

$ mpiexec -n 64 ./atmosphere_model

​

PIO stride?