radiative forcing calculation using rad_diag

[chris_fletcher@utoronto_ca]

I am trying to compute the radiative forcing of aerosols via the diagnostic radiative heating rate calculation using rad_diag_1 in the namelist. I am running cesm1_0_4.

For instance, to diagnose the radiative forcing of sulfate, I am assuming I need to insert the line:
rad_diag_1 = 'D_sulf:$DIN_LOC_ROOT/atm/cam/physprops/sulfate_camrt_c080918.nc'
into the namelist.

This appears to work correctly, since when I run the model I see the following in atm.log:

init_lists: namelist input for diagnostic list:01
   rad namelist spec: D_sulf:/scratch/p/pjk/cgf/cesm1_0_4/inputdata/atm/cam/physprops/sulfate_camrt_c080918.nc
 
 gas list for diag01 calculations
 
 aerosol list for diag01 calculations
   D_sulf optics and phys props in :/scratch/p/pjk/cgf/cesm1_0_4/inputdata/atm/cam/physprops/sulfate_camrt_c080918.nc

However, my question is: how do I see the output from the diagnostic calculation? I was hoping to be able to save FLNT, FSNT etc. from the diagnostic calculation so that I could evaluate the change in TOA net radiation due to sulfate. But I am not sure which fields I need to add to fincl, or if the output is switched on in another way?

Thanks in advance.

 

[eaton]

The diagnostic radiation fields will appear in the cam logfile. Search for ******* MASTER FIELD LIST *******. The diagnostic fields will have '_d1' appended.

Note that the diagnostic calculations are only done for the rrtmg radiation code which is part of the cam5 physics package. I'm wondering whether you have run with rrtmg because I see I see the physprop file specified as sulfate_camrt_c080918.nc which can only be used with the older camrt radiation code which is part of the cam4 physics.

 

[chris_fletcher@utoronto_ca]

Thanks for the response.

Yes, sorry I forgot to say that I am running CAM4 physics, not CAM5. So it sounds like this isn't going to work.