/[lmdze]/trunk/Sources/phylmd/yamada4.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/yamada4.f

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

revision 250 by guez, Wed Nov 15 13:56:45 2017 UTC revision 251 by guez, Mon Jan 8 14:12:02 2018 UTC
# Line 8  module yamada4_m Line 8  module yamada4_m
8    
9  contains  contains
10    
11    SUBROUTINE yamada4(dt, g, zlev, zlay, u, v, teta, q2, km, kn, ustar)    SUBROUTINE yamada4(dt, zlev, zlay, u, v, teta, q2, km, kn, ustar)
12    
13      ! From LMDZ4/libf/phylmd/yamada4.F, version 1.1 2004/06/22 11:45:36      ! From LMDZ4/libf/phylmd/yamada4.F, version 1.1 2004/06/22 11:45:36
14    
15      USE conf_phys_m, ONLY: iflag_pbl      USE conf_phys_m, ONLY: iflag_pbl
16      USE dimphy, ONLY: klev      USE dimphy, ONLY: klev
17      use nr_util, only: assert, assert_eq      use nr_util, only: assert, assert_eq
18        USE suphec_m, ONLY: rg
19    
20      REAL, intent(in):: dt ! pas de temps      REAL, intent(in):: dt ! pas de temps
     real, intent(in):: g  
21    
22      REAL zlev(:, :) ! (knon, klev + 1)      REAL zlev(:, :) ! (knon, klev + 1)
23      ! altitude \`a chaque niveau (interface inf\'erieure de la couche de      ! altitude \`a chaque niveau (interface inf\'erieure de la couche de
# Line 122  contains Line 122  contains
122            m2(ig, k) = ((u(ig, k)-u(ig, k-1))**2 + (v(ig, k)-v(ig, k-1))**2) &            m2(ig, k) = ((u(ig, k)-u(ig, k-1))**2 + (v(ig, k)-v(ig, k-1))**2) &
123                 /(dz(ig, k)*dz(ig, k))                 /(dz(ig, k)*dz(ig, k))
124            dtetadz(ig, k) = (teta(ig, k)-teta(ig, k-1))/dz(ig, k)            dtetadz(ig, k) = (teta(ig, k)-teta(ig, k-1))/dz(ig, k)
125            n2(ig, k) = g*2.*dtetadz(ig, k)/(teta(ig, k-1) + teta(ig, k))            n2(ig, k) = rg*2.*dtetadz(ig, k)/(teta(ig, k-1) + teta(ig, k))
126            ri = n2(ig, k)/max(m2(ig, k), 1.e-10)            ri = n2(ig, k)/max(m2(ig, k), 1.e-10)
127            if (ri.lt.ric) then            if (ri.lt.ric) then
128               rif(ig, k) = frif(ri)               rif(ig, k) = frif(ri)
Legend:
Removed from v.250  
changed lines
  Added in v.251
  ViewVC Help
Powered by ViewVC 1.1.21