phylmd/CV_routines/cv_prelim.f


      SUBROUTINE cv_prelim(len,nd,ndp1,t,q,p,ph &
                          ,lv,cpn,tv,gz,h,hm)
            use cvthermo
            use cv_param
      implicit none

!=====================================================================
! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
!=====================================================================

! inputs:
      integer, intent(in):: len, nd, ndp1
      real, intent(in):: t(len,nd)
      real, intent(in):: q(len,nd), p(len,nd), ph(len,ndp1)

! outputs:
      real lv(len,nd), cpn(len,nd), tv(len,nd)
      real gz(len,nd), h(len,nd), hm(len,nd)

! local variables:
      integer k, i
      real cpx(len,nd)


      do 110 k=1,nlp
        do 100 i=1,len
          lv(i,k)= lv0-clmcpv*(t(i,k)-t0)
          cpn(i,k)=cpd*(1.0-q(i,k))+cpv*q(i,k)
          cpx(i,k)=cpd*(1.0-q(i,k))+cl*q(i,k)
          tv(i,k)=t(i,k)*(1.0+q(i,k)*epsim1)
 100    continue
 110  continue
!
! gz = phi at the full levels (same as p).
!
      do 120 i=1,len
        gz(i,1)=0.0
 120  continue
      do 140 k=2,nlp
        do 130 i=1,len
          gz(i,k)=gz(i,k-1)+hrd*(tv(i,k-1)+tv(i,k)) &
               *(p(i,k-1)-p(i,k))/ph(i,k)
 130    continue
 140  continue
!
! h  = phi + cpT (dry static energy).
! hm = phi + cp(T-Tbase)+Lq
!
      do 170 k=1,nlp
        do 160 i=1,len
          h(i,k)=gz(i,k)+cpn(i,k)*t(i,k)
          hm(i,k)=gz(i,k)+cpx(i,k)*(t(i,k)-t(i,1))+lv(i,k)*q(i,k)
 160    continue
 170  continue

      return
      end
1	guez	52
2			SUBROUTINE cv_prelim(len,nd,ndp1,t,q,p,ph &
3			,lv,cpn,tv,gz,h,hm)
4			use cvthermo
5	guez	103	use cv_param
6	guez	52	implicit none
7
8			!=====================================================================
9			! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
10			!=====================================================================
11
12			! inputs:
13	guez	97	integer, intent(in):: len, nd, ndp1
14	guez	52	real, intent(in):: t(len,nd)
15	guez	103	real, intent(in):: q(len,nd), p(len,nd), ph(len,ndp1)
16	guez	52
17			! outputs:
18			real lv(len,nd), cpn(len,nd), tv(len,nd)
19			real gz(len,nd), h(len,nd), hm(len,nd)
20
21			! local variables:
22			integer k, i
23			real cpx(len,nd)
24
25
26
27			do 110 k=1,nlp
28			do 100 i=1,len
29			lv(i,k)= lv0-clmcpv*(t(i,k)-t0)
30			cpn(i,k)=cpd(1.0-q(i,k))+cpvq(i,k)
31			cpx(i,k)=cpd(1.0-q(i,k))+clq(i,k)
32			tv(i,k)=t(i,k)(1.0+q(i,k)epsim1)
33			100 continue
34			110 continue
35			!
36			! gz = phi at the full levels (same as p).
37			!
38			do 120 i=1,len
39			gz(i,1)=0.0
40			120 continue
41			do 140 k=2,nlp
42			do 130 i=1,len
43			gz(i,k)=gz(i,k-1)+hrd*(tv(i,k-1)+tv(i,k)) &
44			*(p(i,k-1)-p(i,k))/ph(i,k)
45			130 continue
46			140 continue
47			!
48			! h = phi + cpT (dry static energy).
49			! hm = phi + cp(T-Tbase)+Lq
50			!
51			do 170 k=1,nlp
52			do 160 i=1,len
53			h(i,k)=gz(i,k)+cpn(i,k)*t(i,k)
54			hm(i,k)=gz(i,k)+cpx(i,k)(t(i,k)-t(i,1))+lv(i,k)q(i,k)
55			160 continue
56			170 continue
57
58			return
59			end