Question about DCOMPUTETD in wrf_user.f

The NCL graphics package.

Question about DCOMPUTETD in wrf_user.f

Postby Julius » Wed Nov 08, 2017 3:52 pm

Hello there,

I have a question about DCOMPUTETD.

Code: Select all
      SUBROUTINE DCOMPUTETD(TD,PRESSURE,QV_IN,NX)
      IMPLICIT NONE
      INTEGER NX
      DOUBLE PRECISION PRESSURE(NX)
      DOUBLE PRECISION QV_IN(NX)
      DOUBLE PRECISION TD(NX)
C NCLEND
      DOUBLE PRECISION QV,TDC

      INTEGER I

      DO I = 1,NX
          QV = DMAX1(QV_IN(I),0.D0)
c vapor pressure
          TDC = QV*PRESSURE(I)/ (.622D0+QV)

c avoid problems near zero
          TDC = DMAX1(TDC,0.001D0)
          TD(I) = (243.5D0*LOG(TDC)-440.8D0)/ (19.48D0-LOG(TDC))
      END DO

      RETURN
      END


What are the QV_IN and PRESSURE units? kg/kg and Pa?
Julius
 
Posts: 3
Joined: Thu Oct 31, 2013 12:06 pm

Return to NCL

Who is online

Users browsing this forum: No registered users and 2 guests