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

Diff of /trunk/phylmd/hbtm.f

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

revision 38 by guez, Thu Jan 6 17:52:19 2011 UTC revision 49 by guez, Wed Aug 24 11:43:14 2011 UTC
# Line 8  contains Line 8  contains
8         flux_q, u, v, t, q, pblh, cape, EauLiq, ctei, pblT, therm, trmb1, &         flux_q, u, v, t, q, pblh, cape, EauLiq, ctei, pblT, therm, trmb1, &
9         trmb2, trmb3, plcl)         trmb2, trmb3, plcl)
10    
11      use dimens_m      ! D'après Holstag et Boville et Troen et Mahrt
     use dimphy  
     use SUPHEC_M  
     use yoethf_m  
     use fcttre  
   
     ! D'apres Holstag & Boville et Troen & Mahrt  
12      ! JAS 47 BLM      ! JAS 47 BLM
13      ! Algorithme thèse Anne Mathieu      ! Algorithme thèse Anne Mathieu
14      ! Critère d'entraînement Peter Duynkerke (JAS 50)      ! Critère d'entraînement Peter Duynkerke (JAS 50)
# Line 35  contains Line 29  contains
29      ! mais on garde la possibilité de changer si besoin est (jusqu'à présent      ! mais on garde la possibilité de changer si besoin est (jusqu'à présent
30      ! la forme de HB avec le 1er niveau modele etait conservee)      ! la forme de HB avec le 1er niveau modele etait conservee)
31    
32        USE dimphy, ONLY: klev, klon, max
33        USE suphec_m, ONLY: rcpd, rd, retv, rg, rkappa, rlvtt, rtt, rv
34        USE yoethf_m, ONLY: r2es, rvtmp2
35        USE fcttre, ONLY: foeew
36    
37      REAL RLvCp, REPS      REAL RLvCp, REPS
38      ! Arguments:      ! Arguments:
39    
# Line 447  contains Line 446  contains
446    
447      ! Main level loop to compute the diffusivities and      ! Main level loop to compute the diffusivities and
448      ! counter-gradient terms:      ! counter-gradient terms:
449      DO k = 2, isommet      loop_level: DO k = 2, isommet
450         ! Find levels within boundary layer:         ! Find levels within boundary layer:
451         DO i = 1, knon         DO i = 1, knon
452            unslev(i) = .FALSE.            unslev(i) = .FALSE.
# Line 517  contains Line 516  contains
516    
517         ! For all layers, compute integral info and CTEI         ! For all layers, compute integral info and CTEI
518         DO i = 1, knon         DO i = 1, knon
519            if (check(i).or.omegafl(i)) then            if (check(i) .or. omegafl(i)) then
520               if (.not.Zsat(i)) then               if (.not. Zsat(i)) then
521                  T2 = T2m(i) * s(i, k)                  T2 = T2m(i) * s(i, k)
522                  ! thermodyn functions                  ! thermodyn functions
523                  zdelta=MAX(0., SIGN(1., RTT - T2))                  zdelta=MAX(0., SIGN(1., RTT - T2))
# Line 543  contains Line 542  contains
542               ! cette ligne a deja ete faite normalement ?               ! cette ligne a deja ete faite normalement ?
543            endif            endif
544         ENDDO         ENDDO
545      end DO      end DO loop_level
546    
547    END SUBROUTINE HBTM    END SUBROUTINE HBTM
548    
Legend:
Removed from v.38  
changed lines
  Added in v.49
  ViewVC Help
Powered by ViewVC 1.1.21