volcano data can not be loaded

[xiaohardy@163_com]

Hi, every one    I ran a case using cesm1_0_3 with -res T31_g37 -compset B1850. I want to do a only time-varying volcano experiment. So I modified the Buildconf/cam.buildnml.csh script, and added a piece of code in it as follows:&prescribed_volcaero_nl
prescribed_volcaero_datapath = '$DIN_LOC_ROOT/atm/cam/volc'
prescribed_volcaero_file = 'IVI2LoadingLatHeight501-2000_L18_c20100518.nc'
/
   After running, I found the history files are not forced by volcano data. My contents of atm_in file would be posted at the end.PS: Yes, I did the hybird run, set the RUN_STARTDATE "501-01-01" Any suggests must be appreciated. thanks in advancedxiao  atm_in:&aerodep_flx_nl
aerodep_flx_cycle_yr = 1850
aerodep_flx_datapath = '/gpfsES/geo/zywang/CESM/inputdata/atm/cam/chem/trop_mozart_aero/aero'
aerodep_flx_file = 'aerosoldep_monthly_1850_mean_1.9x2.5_c090421.nc'
aerodep_flx_type = 'CYCLICAL'
/
&cam_inparm
absems_data = '/gpfsES/geo/zywang/CESM/inputdata/atm/cam/rad/abs_ems_factors_fastvx.c030508.nc'
bnd_topo = '/gpfsES/geo/zywang/CESM/inputdata/atm/cam/topo/USGS-gtopo30_48x96_c050520.nc'
cam_branch_file = ' '
dtime = 1800
fincl1 = 'U'
fincl2 = 'TREFHTMN','TREFHTMX','TREFHT','PRECC','PRECL','PSL'
mfilt = 1,365,30,120,240
ncdata = 'xiao845.cam2.i.2404-01-01-00000.nc'
nhtfrq = 0,-24,-24,-6,-3
/
&chem_surfvals_nl
ch4vmr = 791.6e-9
co2vmr = 284.7e-6
f11vmr = 12.48e-12
f12vmr = 0.0
n2ovmr = 275.68e-9
/
&circ_diag_nl
do_circulation_diags = .false.
/
&cldfrc_nl
cldfrc_dp1 = 0.14D0
cldfrc_dp2 = 500.0D0
cldfrc_freeze_dry = .true.
cldfrc_ice = .false.
cldfrc_premit = 75000.0D0
cldfrc_rhminh = 0.800D0
cldfrc_rhminl = 0.845D0
cldfrc_sh1 = 0.07D0
cldfrc_sh2 = 500.0D0
/
&cldsed_nl
cldsed_ice_stokes_fac = 1.0D0
/
&cldwat_nl
cldwat_conke = 10.0e-6
cldwat_icritc = 1.0e-6
cldwat_icritw = 4.0e-4
cldwat_r3lcrit = 1.0e-6
/&cldwat_nl
cldwat_conke = 10.0e-6
cldwat_icritc = 1.0e-6
cldwat_icritw = 4.0e-4
cldwat_r3lcrit = 1.0e-6
/
&dyn_spectral_inparm
eul_nsplit = 1
/
&gw_drag_nl
fcrit2 = 1.0
/
&hkconv_nl
hkconv_c0 = 1.0e-4
hkconv_cmftau = 1800.0D0
/
&phys_ctl_nl
cam_chempkg = 'none'
cam_physpkg = 'cam4'
conv_water_in_rad = 0
deep_scheme = 'ZM'
do_iss = .false.
do_tms = .true.
eddy_scheme = 'HB'
history_microphysics = .false.
microp_scheme = 'RK'
shallow_scheme = 'Hack'
srf_flux_avg = 0
tms_orocnst = 1.0D0
tms_z0fac = 0.100D0
/
&prescribed_aero_nl
prescribed_aero_cycle_yr = 1850
prescribed_aero_datapath = '/gpfsES/geo/zywang/CESM/inputdata/atm/cam/chem/trop_mozart_aero/aero'
prescribed_aero_file = 'aero_1.9x2.5_L26_1850clim_c091112.nc'
prescribed_aero_specifier = 'sulf:SO4', 'bcar1:CB1', 'bcar2:CB2', 'ocar1:OC1', 'ocar2:OC2', 'sslt1:SSLT01', 'sslt2:SSLT02',
'sslt3:SSLT03', 'sslt4:SSLT04', 'dust1:DST01', 'dust2:DST02', 'dust3:DST03', 'dust4:DST04'
prescribed_aero_type = 'CYCLICAL'
/
&prescribed_ozone_nl
prescribed_ozone_cycle_yr = 1850
prescribed_ozone_datapath = '/gpfsES/geo/zywang/CESM/inputdata/atm/cam/ozone'
prescribed_ozone_file = 'ozone_1.9x2.5_L26_1850clim_c091112.nc'
prescribed_ozone_name = 'O3'
prescribed_ozone_type = 'CYCLICAL'/
&prescribed_volcaero_nl
prescribed_volcaero_datapath = '/gpfsES/geo/zywang/CESM/inputdata/atm/cam/volc'
prescribed_volcaero_file = 'IVI2LoadingLatHeight501-2000_L18_c20100518.nc'
/&rad_cnst_nl
rad_climate = 'P_Q:H2O', 'D_O2:O2', 'D_CO2:CO2', 'D_ozone:O3', 'D_N2O:N2O', 'D_CH4:CH4', 'D_CFC11:CFC11',
'D_CFC12:CFC12', 'D_sulf:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/sulfate_camrt_c080918.nc',
'D_dust1:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/dust1_camrt_c080918.nc',
'D_dust2:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/dust2_camrt_c080918.nc',
'D_dust3:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/dust3_camrt_c080918.nc',
'D_dust4:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/dust4_camrt_c080918.nc',
'D_bcar1:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/bcpho_camrt_c080918.nc',
'D_bcar2:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/bcphi_camrt_c080918.nc',
'D_ocar1:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/ocpho_camrt_c080918.nc',
'D_ocar2:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/ocphi_camrt_c080918.nc',
'D_SSLTA:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/ssam_camrt_c080918.nc',
'D_SSLTC:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/sscm_camrt_c080918.nc'
/
&solar_inparm
solar_const = 1360.89
solar_htng_spctrl_scl = .false.
/
&tropopause_nl
tropopause_climo_file = '/gpfsES/geo/zywang/CESM/inputdata/atm/cam/chem/trop_mozart/ub/clim_p_trop.nc'
/
&zmconv_nl
zmconv_c0_lnd = 0.0020D0
zmconv_c0_ocn = 0.0020D0
zmconv_ke = 3.0E-6
/

 

[eaton]

The volcano aerosols will not be felt by the radiation unless there is an entry for them in the rad_climate variable.  This would be added automatically by CAM's build-namelist utility if you used that to add the prescribed_volcaero_* variables.  If you edit the CAM namelist by hand, then you must also add the necessary modification to the rad_climate variable.  In this instance the string that should be appended to the list is
'D_VOLC_MMR:/gpfsES/geo/zywang/CESM/inputdata/atm/cam/physprops/sulfuricacid_cam3_c080918.nc'

 

[xiaohardy@163_com]

Yeah! Thanks eaton ! But how to use the CAM's build-namelist utility to add the prescribed_volcaero_* variabeles?xiao

 

[eaton]

Instead of editing the cam.buildnml.csh script you should put the same information into the $CASEROOT/user_nl_cam file as described in the user guide.  The namelist variables set in this file are used as input to cam's build-namelist utility.  It is the output from the build-namelist utility that is found in the cam.buildnml.csh script.see http://www.cesm.ucar.edu/models/cesm1.0/cesm/cesm_doc_1_0_4/x1062.html#customizing_cam 

 

[xiaohardy@163_com]

Thanks eaton, I got it

 

[wanlingfeng_123@163_com]

hi eaton    I want to add  prescribed volcano aerosols the same as a specifer volcanic episode, which give a daily forcing to the CESM1.0.4 model for last four months. Then study ocean response to this volcanic episode.I have used -res f19_g16 -compset B_1850 -mach generic_linux_intel  -max_tasks_per_node 12 to do this work. But i don't konw how to add the daily volcanic forcing to the model.Could  you give me some advices, thank you very much. 

 

[eaton]

First you should make sure you can run an existing compset which incorporates the volcanic forcings.  Next, copy the dataset which contains the forcings and use it as a template to create the forcing scenario you choose.  The only difficulty is learning how to manipulate netcdf datasets.  We don't provide tools for this, but there are many freely available general pupose tools for manipulating netcdf data, for example, the NCO operators.

 

[wanlingfeng_123@163_com]

thank you  for your advice.First  i don't konw whether B_1850 compset can do daily changed volcanic with flow namelist,&prescribed_volcaero_nl prescribed_volcaero_datapath           = '/NSFCGZ/nsfc2015_229/wlf/inputdata/atm/cam/volc' prescribed_volcaero_file               = 'IVI2LoadingLatHeight501-2000_L18_c20100518.nc'/'D_VOLC_MMR:$/inputdata/atm/cam/physprops/sulfuricacid_cam3_c080918.nc'Next the  'IVI2LoadingLatHeight501-2000_L18_c20100518.nc' data is monthly and years from 501 to 2000. what i want is  only specifier year with daily forcing and last four months.  

 

[eaton]

The B_1850 compset is designed to do a climatological run with 1850 forcings.  Either the forcing datasets contain a single year of data which is cycled year after year, or the dataset contains multiple years of data and a namelist option is available to specify a single year from the dataset to cycle over.  The B_1850 compset doesn't include volcanoe forcings because it doesn't make sense to cycle them year after year. 

 

[xuzq@mail_iap_ac_cn]

Hi Eaton,Based on the discussion above, it can be inferred that F_2000 (CESM1.0) compset doesn't include volcanic forcings. Am I right?Zhiqing

 

[kangwanying1992@gmail_com]

Hi Eaton,I noticed that the unit of the volcano forcing is kg/kg like a mixing ratio. Does this mean the areosol emitted by the volcano cannot be advected by the atmosphere? Is there a possibility to add areasol to the atmosphere and let the Brew-Dobson circulation advect it? Many thanks for any suggestion!Best,Wanying

 

[mmills]

Hi Eaton,I noticed that the unit of the volcano forcing is kg/kg like a mixing ratio. Does this mean the areosol emitted by the volcano cannot be advected by the atmosphere? Is there a possibility to add areasol to the atmosphere and let the Brew-Dobson circulation advect it? Many thanks for any suggestion!Best,Wanying

Prognostic stratospheric aerosol is a new feature in CESM2(WACCM6).