/[lmdze]/trunk/phylmd/physiq.f90
ViewVC logotype

Diff of /trunk/phylmd/physiq.f90

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

revision 339 by guez, Thu Sep 26 17:08:42 2019 UTC revision 340 by guez, Thu Sep 26 17:29:50 2019 UTC
# Line 19  contains Line 19  contains
19      use ajsec_m, only: ajsec      use ajsec_m, only: ajsec
20      use calltherm_m, only: calltherm      use calltherm_m, only: calltherm
21      USE clesphys, ONLY: cdhmax, cdmmax, ecrit_ins, ok_instan      USE clesphys, ONLY: cdhmax, cdmmax, ecrit_ins, ok_instan
22      USE clesphys2, ONLY: conv_emanuel, nbapp_rad, new_oliq, ok_orodr, ok_orolf      USE clesphys2, ONLY: conv_emanuel, nbapp_rad, ok_orodr, ok_orolf
23      USE conf_interface_m, ONLY: conf_interface      USE conf_interface_m, ONLY: conf_interface
24      USE pbl_surface_m, ONLY: pbl_surface      USE pbl_surface_m, ONLY: pbl_surface
25      use clouds_gno_m, only: clouds_gno      use clouds_gno_m, only: clouds_gno
# Line 103  contains Line 103  contains
103    
104      LOGICAL:: firstcal = .true.      LOGICAL:: firstcal = .true.
105    
     LOGICAL, PARAMETER:: ok_stratus = .FALSE.  
     ! Ajouter artificiellement les stratus  
   
106      ! pour phystoke avec thermiques      ! pour phystoke avec thermiques
107      REAL fm_therm(klon, llm + 1)      REAL fm_therm(klon, llm + 1)
108      REAL entr_therm(klon, llm)      REAL entr_therm(klon, llm)
# Line 727  contains Line 724  contains
724            q_seri(i, k) = q_seri(i, k) + d_q_lsc(i, k)            q_seri(i, k) = q_seri(i, k) + d_q_lsc(i, k)
725            ql_seri(i, k) = ql_seri(i, k) + d_ql_lsc(i, k)            ql_seri(i, k) = ql_seri(i, k) + d_ql_lsc(i, k)
726            cldfra(i, k) = rneb(i, k)            cldfra(i, k) = rneb(i, k)
           IF (.NOT. new_oliq) cldliq(i, k) = ql_seri(i, k)  
727         ENDDO         ENDDO
728      ENDDO      ENDDO
729    
# Line 784  contains Line 780  contains
780         cldliq = cldliq + rnebcon * clwcon         cldliq = cldliq + rnebcon * clwcon
781      ENDIF      ENDIF
782    
     ! 2. Nuages stratiformes  
   
     IF (ok_stratus) THEN  
        CALL diagcld2(paprs, play, t_seri, q_seri, diafra, dialiq)  
        DO k = 1, llm  
           DO i = 1, klon  
              IF (diafra(i, k) > cldfra(i, k)) THEN  
                 cldliq(i, k) = dialiq(i, k)  
                 cldfra(i, k) = diafra(i, k)  
              ENDIF  
           ENDDO  
        ENDDO  
     ENDIF  
   
783      ! Precipitation totale      ! Precipitation totale
784      DO i = 1, klon      DO i = 1, klon
785         rain_fall(i) = rain_con(i) + rain_lsc(i)         rain_fall(i) = rain_con(i) + rain_lsc(i)
Legend:
Removed from v.339  
changed lines
  Added in v.340
  ViewVC Help
Powered by ViewVC 1.1.21