Changeset 2450
- Timestamp:
- Feb 24, 2016, 3:53:09 PM (9 years ago)
- Location:
- LMDZ5/trunk/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90
r2410 r2450 305 305 306 306 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012 307 !! REAL, DIMENSION(klon,klev), INTENT(IN) :: t_x ! Temp érature hors poche froide308 !! REAL, DIMENSION(klon,klev), INTENT(IN) :: t_w ! Temp érature dans la poches froide307 !! REAL, DIMENSION(klon,klev), INTENT(IN) :: t_x ! Temp\'erature hors poche froide 308 !! REAL, DIMENSION(klon,klev), INTENT(IN) :: t_w ! Temp\'erature dans la poches froide 309 309 !! REAL, DIMENSION(klon,klev), INTENT(IN) :: q_x ! 310 !! REAL, DIMENSION(klon,klev), INTENT(IN) :: q_w ! Pareil pour l'humidit é310 !! REAL, DIMENSION(klon,klev), INTENT(IN) :: q_w ! Pareil pour l'humidit\'e 311 311 REAL, DIMENSION(klon,klev), INTENT(IN) :: wake_dlt !temperature difference between (w) and (x) (K) 312 312 REAL, DIMENSION(klon,klev), INTENT(IN) :: wake_dlq !humidity difference between (w) and (x) (K) … … 679 679 680 680 !!! jyg le 25/03/2013 681 !! Variables intermediaires pour le raccord des deux colonnes àla surface681 !! Variables intermediaires pour le raccord des deux colonnes \`a la surface 682 682 REAL :: dd_Ch 683 683 REAL :: dd_Cm … … 1572 1572 1573 1573 !!! nrlmd le 13/06/2011 1574 !----- On finit le calcul des coefficients d' échange:on multiplie le cdrag par le module du vent et la densité dans la première couche1574 !----- On finit le calcul des coefficients d'\'echange:on multiplie le cdrag par le module du vent et la densit\'e dans la premi\`ere couche 1575 1575 ! Kech_h_x(j) = ycdragh_x(j) * & 1576 1576 ! (1.0+SQRT(yu_x(j,1)**2+yv_x(j,1)**2)) * & … … 1664 1664 ENDIF 1665 1665 ! 1666 ! Calcul des coef A, B équivalents dans la couche 11666 ! Calcul des coef A, B \'equivalents dans la couche 1 1667 1667 ! 1668 1668 AcoefH(j) = AcoefH_x(j) + ywake_s(j)*(Kech_H_wp(j)/Kech_Hp(j))*dd_AH … … 1684 1684 1685 1685 ! 1686 ! Calcul des cdrag équivalents dans la couche1686 ! Calcul des cdrag \'equivalents dans la couche 1687 1687 ! 1688 1688 ycdragm(j) = ycdragm_x(j) + ywake_s(j)*dd_CM 1689 1689 ycdragh(j) = ycdragh_x(j) + ywake_s(j)*dd_CH 1690 1690 ! 1691 ! Calcul de T, q, u et v équivalents dans la couche 11691 ! Calcul de T, q, u et v \'equivalents dans la couche 1 1692 1692 yt(j,1) = yt_x(j,1) + ywake_s(j)*(Kech_h_w(j)/Kech_h(j))*dd_t 1693 1693 yq(j,1) = yq_x(j,1) + ywake_s(j)*(Kech_h_w(j)/Kech_h(j))*dd_q … … 1993 1993 !!jyg!! ENDIF 1994 1994 !!jyg!! 1995 !!jyg!!!-----Calcul de ybeta (evap_r éelle/evap_potentielle)1995 !!jyg!!!-----Calcul de ybeta (evap_r\'eelle/evap_potentielle) 1996 1996 !!jyg!!!!!!! jyg le 23/02/2012 1997 1997 !!jyg!!!!!!! -
LMDZ5/trunk/libf/phylmd/physiq_mod.F90
r2448 r2450 451 451 real w0(klon) ! Vitesse des thermiques au LCL 452 452 real w_conv(klon) ! Vitesse verticale de grande \'echelle au LCL 453 real tke0(klon,klev+1) ! TKE au d ébut du pas de temps453 real tke0(klon,klev+1) ! TKE au début du pas de temps 454 454 real therm_tke_max0(klon) ! TKE dans les thermiques au LCL 455 455 real env_tke_max0(klon) ! TKE dans l'environnement au LCL … … 1768 1768 else 1769 1769 1770 !CR: on r é-évapore eau liquide et glace1770 !CR: on ré-évapore eau liquide et glace 1771 1771 1772 1772 ! zdelta = MAX(0.,SIGN(1.,RTT-t_seri(i,k))) … … 1780 1780 q_seri(i,k) = q_seri(i,k) + zb 1781 1781 ql_seri(i,k) = 0.0 1782 !on évapore la glace1782 !on évapore la glace 1783 1783 qs_seri(i,k) = 0.0 1784 1784 d_t_eva(i,k) = za … … 2591 2591 ! 2592 2592 !!! nrlmd le 22/03/2011---Si on met les poches hors des thermiques il faut rajouter cette 2593 !------------------------- tendance calcul ée hors des poches froides2593 !------------------------- tendance calcul\'ee hors des poches froides 2594 2594 ! 2595 2595 if (iflag_wake>=1) then … … 2655 2655 DO i=1,klon 2656 2656 IF (rneb(i,k)==0.) THEN 2657 ! On ne tient compte des tendances qu'en dehors des nuages (c'est �|dire2657 ! On ne tient compte des tendances qu'en dehors des nuages (c'est-\`a-dire 2658 2658 ! a priri dans une region ou l'eau se reevapore). 2659 2659 dt_dwn(i,k)= dt_dwn(i,k)+ & … … 2806 2806 IF (mod(iflag_pbl_split/2,2) .EQ. 1) THEN 2807 2807 ! Si les thermiques ne sont presents que hors des poches, la tendance moyenne 2808 ! associ ée doit etre multipliee par la fraction surfacique qu'ils couvrent.2808 ! associ\'ee doit etre multipliee par la fraction surfacique qu'ils couvrent. 2809 2809 DO k=1,klev 2810 2810 DO i=1,klon
Note: See TracChangeset
for help on using the changeset viewer.