_RL frace, rs, hq
hq = hnew(k)*qicen(k)-dhi*Lfresh
qicen(k)=hq/hnew(k)
hqol = huol*hq/hu
hq = 2. _d 0
_RL ustar,tstar,qstar,ht,hu,hq
_RL hq
& hq,
hq = ht
CALL WRITE_0D_RL( hq, INDEX_NONE, 'hq =',