/[lmdze]/trunk/phylmd/Interface_surf/coef_diff_turb.f
ViewVC logotype

Diff of /trunk/phylmd/Interface_surf/coef_diff_turb.f

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

revision 298 by guez, Thu Jul 26 16:45:51 2018 UTC revision 303 by guez, Thu Sep 6 14:25:07 2018 UTC
# Line 12  contains Line 12  contains
12      use nr_util, only: assert      use nr_util, only: assert
13    
14      USE clesphys, ONLY: ok_kzmin      USE clesphys, ONLY: ok_kzmin
     use comconst, only: dtphys  
15      use coefkz_m, only: coefkz      use coefkz_m, only: coefkz
16      use coefkzmin_m, only: coefkzmin      use coefkzmin_m, only: coefkzmin
17      use coefkz2_m, only: coefkz2      use coefkz2_m, only: coefkz2
# Line 64  contains Line 63  contains
63              * (paprs(:, 1) / pplay(:, k))**rkappa * (1. + 0.61 * q(:, k))              * (paprs(:, 1) / pplay(:, k))**rkappa * (1. + 0.61 * q(:, k))
64    
65         zlev(:, 1) = 0.         zlev(:, 1) = 0.
        zlev(:, klev + 1) = 2. * zlay(:, klev) - zlay(:, klev - 1)  
66         forall (k = 2:klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1))         forall (k = 2:klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1))
67           zlev(:, klev + 1) = 2. * zlay(:, klev) - zlev(:, klev)
68    
69         CALL yamada4(dtphys, zlev, zlay, u, v, teta, q2, coefm, coefh, &         CALL yamada4(zlev, zlay, u, v, teta, q2, coefm, coefh, &
70              ustarhb(u(:, 1), v(:, 1), cdragm))              ustarhb(u(:, 1), v(:, 1), cdragm))
71      else      else
72         CALL coefkz(nsrf, paprs, pplay, ts, u, v, t, q, zgeop, coefm, coefh)         CALL coefkz(nsrf, paprs, pplay, ts, u, v, t, q, zgeop, coefm, coefh)
Legend:
Removed from v.298  
changed lines
  Added in v.303
  ViewVC Help
Powered by ViewVC 1.1.21