Units problem in ANA_SRFLUX

[LauraB]

Hi all,

The description in ANA_SRFLUX.h reads:

Code: Select all

!  This subroutine sets kinematic surface solar shortwave radiation    !
!  flux "srflx" (degC m/s) using an analytical expression.             !

However, when the short wave input has to be given, the title is:

Code: Select all

!-----------------------------------------------------------------------
!  Set incoming solar shortwave radiation (W/m2).
!-----------------------------------------------------------------------

After getting crazy temperatures with srflx given in W/m2 (and with SOLAR_SOURCE defined), I think this last title is incorrect, and that srflx should be given in units of degC m/s (i.e., W/m2 /(rho0*cp) )...

Thanks and best regards,
Laura

[wilkin]

I believe you are correct. If you look at the preceeding code segment where srflux is computed from an analytical expression considering local time of day and latitude it documents the steps:

Code: Select all

120 	!-----------------------------------------------------------------------
121 	!  Compute shortwave radiation (degC m/s):
122 	!
123 	!  ALBEDO option: Compute shortwave radiation flux using the Laevastu
124 	!                 cloud correction to the Zillman equation for cloudless
125 	!  radiation (Parkinson and Washington 1979, JGR, 84, 311-337).  Notice
126 	!  that flux is scaled from W/m2 to degC m/s by dividing by (rho0*Cp).