module_damping_em.F

1 !WRF:MODEL_LAYER:DYNAMICS
2 !
3 MODULE module_damping_em
4 
5   USE module_wrf_error
6 
7 CONTAINS
8 
9 !------------------------------------------------------------------------------
10 
11   SUBROUTINE held_suarez_damp( ru_tend, rv_tend, ru, rv, p, pb,  &
12                                ids,ide, jds,jde, kds,kde, &
13                                ims,ime, jms,jme, kms,kme, &
14                                its,ite, jts,jte, kts,kte )
15 
16     IMPLICIT NONE
17 
18     INTEGER,      INTENT(IN   )    :: ids,ide, jds,jde, kds,kde
19     INTEGER,      INTENT(IN   )    :: ims,ime, jms,jme, kms,kme
20     INTEGER,      INTENT(IN   )    :: its,ite, jts,jte, kts,kte
21 
22     REAL, DIMENSION( ims:ime , kms:kme, jms:jme ),          &
23           INTENT(INOUT) ::                         ru_tend, &
24                                                    rv_tend
25 
26     REAL, DIMENSION( ims:ime , kms:kme, jms:jme ),          &
27           INTENT(IN) ::                            ru, rv, p, pb
28 
29     integer :: i,j,k
30 
31     REAL :: delty,delthez,sigb,kka,kkf
32     REAL :: sig,sigterm,kkt,kkv,daylensec
33 
34     sigb=0.7
35     daylensec=60.0*60.0*24.0
36     kkf=1.0/daylensec
37 
38 !  fixed limits so no divide by zero, WCS 070509
39 
40     DO j=max(jds+1,jts),min(jde-1,jte)
41     DO k=kts,MIN(kte,kde-1)
42     DO i=its,ite
43 
44        sig=    (p(i,k,j-1)+pb(i,k,j-1)+p(i,k,j)+pb(i,k,j))/     &
45                (p(i,1,j-1)+pb(i,1,j-1)+p(i,1,j)+pb(i,1,j))
46        sigterm=max(0.0,(sig-sigb)/(1.0-sigb))
47        kkv=kkf*sigterm
48        rv_tend(i,k,j)=rv_tend(i,k,j)-kkv*rv(i,k,j)
49  
50     END DO
51     END DO
52     END DO
53 
54     DO j=jts,min(jde-1,jte)
55     DO k=kts,MIN(kte,kde-1)
56     DO i=its,ite
57 
58        sig=    (p(i-1,k,j)+pb(i-1,k,j)+p(i,k,j)+pb(i,k,j))/     &
59                (p(i-1,1,j)+pb(i-1,1,j)+p(i,1,j)+pb(i,1,j))
60        sigterm=max(0.0,(sig-sigb)/(1.0-sigb))
61        kkv=kkf*sigterm
62        ru_tend(i,k,j)=ru_tend(i,k,j)-kkv*ru(i,k,j)
63  
64     END DO
65     END DO
66     END DO
67 
68   END SUBROUTINE held_suarez_damp
69 
70 !------------------------------------------------------------------------------
71 
72 END MODULE module_damping_em