phylmd/CV_routines/cv_prelim.f90


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

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

! inputs:
      integer len, nd, ndp1
      real, intent(in):: t(len,nd)
      real 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
2	SUBROUTINE cv_prelim(len,nd,ndp1,t,q,p,ph &
3	,lv,cpn,tv,gz,h,hm)
4	use cvthermo
5	use cvparam
6	implicit none
7
8	!=====================================================================
9	! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
10	!=====================================================================
11
12	! inputs:
13	integer len, nd, ndp1
14	real, intent(in):: t(len,nd)
15	real q(len,nd), p(len,nd), ph(len,ndp1)
16
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