/[lmdze]/trunk/phylmd/cdrag.f
ViewVC logotype

Diff of /trunk/phylmd/cdrag.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/phylmd/clcdrag.f revision 274 by guez, Wed Jul 11 16:50:27 2018 UTC trunk/phylmd/cdrag.f revision 275 by guez, Wed Jul 11 17:06:09 2018 UTC
# Line 1  Line 1 
1  module clcdrag_m  module cdrag_m
2    
3    IMPLICIT NONE    IMPLICIT NONE
4    
5  contains  contains
6    
7    SUBROUTINE clcdrag(nsrf, speed, t, q, zgeop, psol, ts, qsurf, rugos, pcfm, &    SUBROUTINE cdrag(nsrf, speed, t, q, zgeop, psol, ts, qsurf, rugos, pcfm, &
8         pcfh, pref)         pcfh, pref)
9    
10      ! From LMDZ4/libf/phylmd/clcdrag.F90, version 1.1.1.1, 2004/05/19 12:53:07      ! From LMDZ4/libf/phylmd/clcdrag.F90 and
11        ! LMDZ4/libf/phylmd/coefcdrag.F90, version 1.1.1.1, 2004/05/19
12        ! 12:53:07
13    
14      ! Objet : calcul des drag coefficients au sol pour le moment et      ! Objet : calcul des drag coefficients au sol pour le moment et
15      ! les flux de chaleur sensible et latente et calcul de la pression      ! les flux de chaleur sensible et latente et calcul de la pression
# Line 60  contains Line 62  contains
62    
63      knon = assert_eq([size(speed), size(t), size(q), size(zgeop), size(ts), &      knon = assert_eq([size(speed), size(t), size(q), size(zgeop), size(ts), &
64           size(qsurf), size(rugos), size(pcfm), size(pcfh), size(pcfm)], &           size(qsurf), size(rugos), size(pcfm), size(pcfh), size(pcfm)], &
65           "clcdrag knon")           "cdrag knon")
66            
67      DO i = 1, knon      DO i = 1, knon
68         zdu2 = max(cepdu2, speed(i)**2)         zdu2 = max(cepdu2, speed(i)**2)
# Line 100  contains Line 102  contains
102      if (present(pref)) &      if (present(pref)) &
103           pref = exp(log(psol) - zgeop / (RD * t * (1. + RETV * max(q, 0.))))           pref = exp(log(psol) - zgeop / (RD * t * (1. + RETV * max(q, 0.))))
104    
105    END SUBROUTINE clcdrag    END SUBROUTINE cdrag
106    
107  end module clcdrag_m  end module cdrag_m
Legend:
Removed from v.274  
changed lines
  Added in v.275
  ViewVC Help
Powered by ViewVC 1.1.21