CALL WRITE_0D_R8( temp_EvPrRn, INDEX_NONE,
IF (temp_EvPrRn.NE.UNSET_RL) THEN
& *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
& temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,
& convertFW2Salt, temp_EvPrRn, salt_EvPrRn,
_RL temp_EvPrRn
temp_EvPrRn = UNSET_RL