O t3d, q3d, u3d, v3d,
I deltaT, t3d, q3d, pFull3d, pHalf3d, coldT,
q3d_tmp = q3d + cond_dt_qg
q3d_tmp = q3d
I sNx,sNy, myTime, lat2d, pHalf3d, t3d, q3d,
I t3d(:,:,Nr), q3d(:,:,Nr), u3d(:,:,Nr), v3d(:,:,Nr),
I u3d, v3d, t3d, q3d,
_RL q3d (sNx,sNy,Nr)
q3d(:,:,k) = MAX( salt(1:sNx,1:sNy,kc,bi,bj), 0. _d 0 )