scapps@uci_edu
Member
I am adding some additional surface variables to CAM such as different wind speed bins, etc.. These fields will be used throughout CAM/CLM and will need to be output into the history file. Is it safe to append these variables to the comsrf.F90..surface_state type, or should I create a new type?
type surface_state
real(r8) :: tbot(pcols) ! bot level temperature
real(r8) :: zbot(pcols) ! bot level height above surface
real(r8) :: ubot(pcols) ! bot level u wind
real(r8) :: vbot(pcols) ! bot level v wind
real(r8) :: qbot(pcols) ! bot level specific humidity
real(r8) :: pbot(pcols) ! bot level pressure
real(r8) :: flwds(pcols) !
real(r8) :: precsc(pcols) !
real(r8) :: precsl(pcols) !
real(r8) :: precc(pcols) !
real(r8) :: precl(pcols) !
real(r8) :: soll(pcols) !
real(r8) :: sols(pcols) !
real(r8) :: solld(pcols) !
real(r8) :: solsd(pcols) !
real(r8) :: srfrad(pcols) !
real(r8) :: thbot(pcols) !
real(r8) :: tssub(pcols,plevmx) ! cam surface/subsurface temperatures
! ADDITIONAL VARIABLES
real(r8) :: wnd_spd_bins(pcols,bin_nbr) ! ADDED
real(r8) :: dnmss(pcols) ! Downdraft mass flux at surface
real(r8) :: tke(pcols) ! Turbulence kinetic energy at surface
end type surface_state
And, add the outfld calls to the diagnostics.F90 module:
call outfld('DNMSS', surface%dnmss, pcols, lchnk)
.....
Thank you for your time,
Scott
type surface_state
real(r8) :: tbot(pcols) ! bot level temperature
real(r8) :: zbot(pcols) ! bot level height above surface
real(r8) :: ubot(pcols) ! bot level u wind
real(r8) :: vbot(pcols) ! bot level v wind
real(r8) :: qbot(pcols) ! bot level specific humidity
real(r8) :: pbot(pcols) ! bot level pressure
real(r8) :: flwds(pcols) !
real(r8) :: precsc(pcols) !
real(r8) :: precsl(pcols) !
real(r8) :: precc(pcols) !
real(r8) :: precl(pcols) !
real(r8) :: soll(pcols) !
real(r8) :: sols(pcols) !
real(r8) :: solld(pcols) !
real(r8) :: solsd(pcols) !
real(r8) :: srfrad(pcols) !
real(r8) :: thbot(pcols) !
real(r8) :: tssub(pcols,plevmx) ! cam surface/subsurface temperatures
! ADDITIONAL VARIABLES
real(r8) :: wnd_spd_bins(pcols,bin_nbr) ! ADDED
real(r8) :: dnmss(pcols) ! Downdraft mass flux at surface
real(r8) :: tke(pcols) ! Turbulence kinetic energy at surface
end type surface_state
And, add the outfld calls to the diagnostics.F90 module:
call outfld('DNMSS', surface%dnmss, pcols, lchnk)
.....
Thank you for your time,
Scott