A ROUTINE FOR THE CALCULATION OF THE TIME-DEPENDENT HEIGHT OF THE ATMOSPHERIC BOUNDARY LAYER FROM SURFACE-LAYER PARAMETERS