_RL tmpbulk
tmpbulk= exf_BulkqSat(ssttmp + cen2kel)
ssq = saltsat*tmpbulk/atmrho
tmpbulk= exf_BulkCdn(sh)
rdn = sqrt(tmpbulk)
tmpbulk= exf_BulkRhn(stable)
tstar = tmpbulk*deltap
tmpbulk= exf_BulkCdn(uzn)
rh = tmpbulk/( exf_one +
& tmpbulk/karman*(aln - psixh) )
ustress(i,j,bi,bj) = atmrho*tmpbulk*us*
vstress(i,j,bi,bj) = atmrho*tmpbulk*us*