C Theta-e calculcation
        e=prs(i,j,k)*q(i,j,k)/(0.622+q(i,j,k))/100.
        tlcl=55.0+2840.0/(3.5*alog(t(i,j,k))-alog(e)-4.805)
        tlcl=amin1(tlcl,t(i,j,k))
        theta=t(i,j,k)*(100000./prs(i,j,k))**(0.2854*
     &        (1.0-0.28*q(i,j,k)))
        thetae(i,j,k)=theta*exp(((3376./tlcl)-2.54)
     &                    *q(i,j,k)*(1.0+0.81*q(i,j,k))  )
 


C RH calculation
        if(t(i,j,k).le.273.15 .and. iice.ne.0)then
          es=6.11*exp(22.514-(6150./t(i,j,k)))
        else
          es=6.112*exp(17.67*((t(i,j,k)-273.15)/(t(i,j,k)-29.65)))
        endif
        qs=0.622*es/((prs(I,J,K)/100.)-es)
        rh(i,j,k)=100.*q(I,J,K)/qs