[WIP] Add mixed-phase partition plot

[chengzhuzhang]

Adding mixed-phase partition plot form E3SM Diags.

@golaz @wlin7 @shaochengx Just a heads-up that this is an example plot of the mixed phase partition plot with the first historical ensemble:
https://web.lcrc.anl.gov/public/e3sm/diagnostic_output/ac.zhang40/E3SMv3_mp/v3.LR.historical_0051/e3sm_diags/atm_monthly_180x360_aave/model_vs_obs_1985-2014/viewer/mp_partition/variable/lcf/plot.html
The deviation of v3 from v2 is obvious. It is expected because the cloud ice (CLDICE) used in this diagnostics has a changed definition that both ice and snow are in this same category, thus results lower liquid condensation fraction.

@zyuying and I planned to use CLD_CAL_TMPLIQ/ICE from simulator to replace CLDICE to make more meaningful comparison. It needs additional development in E3SM Diags, therefore, we won't have this set ready for next zppy release.

Also tagging @susburrows because this feature in e3sm diags was also requested by Aerosol Working Group.

[chengzhuzhang]

Update my run script here:

[default]
#input = /lcrc/group/e3sm2/ac.golaz/E3SMv3/v3.LR.piControl
input = /lcrc/group/e3sm2/ac.wlin/E3SMv3/v3.LR.historical_0051
output = /lcrc/group/e3sm2/ac.zhang40/E3SMv3/v3.LR.historical_trop_mp
case = v3.LR.historical_0051
www = /lcrc/group/e3sm/public_html/diagnostic_output/ac.zhang40/E3SMv3_mp
partition = compute
#environment_commands = "source /lcrc/soft/climate/e3sm-unified/load_latest_e3sm_unified_chrysalis.sh"
environment_commands = "source /lcrc/soft/climate/e3sm-unified/test_e3sm_unified_1.9.3rc4_chrysalis.sh"
campaign = "water_cycle"
debug = True

[ts]
active = True
years = "1985:2014:30"
walltime = "00:50:00"

  [[ atm_monthly_180x360_aave ]]
  input_subdir = "archive/atm/hist"
  input_files = "eam.h0"
  frequency = "monthly"
  mapping_file = /home/ac.zender/data/maps/map_ne30pg2_to_cmip6_180x360_aave.20200201.nc
  vars = "FSNTOA,FLUT,FSNT,FLNT,FSNS,FLNS,SHFLX,QFLX,TAUX,TAUY,PRECC,PRECL,PRECSC,PRECSL,TS,TREFHT,CLDTOT,CLDHGH,CLDMED,CLDLOW,U,ICEFRAC,LANDFRAC,OCNFRAC,PS,CLDICE,CLDLIQ,T"

[e3sm_diags]
active = True
walltime = "4:00:00"
#years = "0001:0100:50", "0001:0100:100"
years = "1985:2014:30",
ts_num_years = 30
ref_start_yr = 1979
ref_final_yr = 2016
multiprocessing = True
num_workers = 8

  [[ atm_monthly_180x360_aave ]]
  environment_commands = "source /home/ac.zhang40/y/etc/profile.d/conda.sh; conda activate edv290"
  short_name = 'v3.LR.historical_0051'
  grid = '180x360_aave'
#  sets="lat_lon","zonal_mean_xy","zonal_mean_2d","polar","cosp_histogram","meridional_mean_2d","annual_cycle_zonal_mean","qbo","diurnal_cycle","zonal_mean_2d_stratosphere","aerosol_budget","tropical_subseasonal",
  sets="mp_partition",
  output_format_subplot = "pdf",

[susburrows]

Thanks, @chengzhuzhang, for the heads-up.