FUNCTION VOID  lwfac  (n_surf,n_all,surf_emit,fi,surf_area,
                       psi_lw)

/* Calculate absorption matrix for LW radiation
   [absorbed at surface, emitted at surface]
*/

LANGUAGE F77
INPUT
    INT n_surf,n_all;
    FLOAT surf_emit, fi, surf_area;
    
OUTPUT
    FLOAT psi_lw;


EXTERNAL  lwfac (n_surf,n_all,surf_emit,fi,surf_area,psi_lw)


END_FUNCTION