Changeset 1740 for LMDZ5/trunk
- Timestamp:
- Apr 2, 2013, 3:46:34 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/physiq.F
r1738 r1740 546 546 c================================================================================================= 547 547 cCR04.12.07: on ajoute les nouvelles variables du nouveau schema de convection avec poches froides 548 c Variables li ées àla poche froide (jyg)548 c Variables li\'ees \`a la poche froide (jyg) 549 549 550 550 REAL mip(klon,klev) ! mass flux shed by the adiab ascent at each level … … 579 579 580 580 cRC 581 c Variables li ées àla poche froide (jyg et rr)582 c Version diagnostique pour l'instant : pas de r étroaction sur la convection581 c Variables li\'ees \`a la poche froide (jyg et rr) 582 c Version diagnostique pour l'instant : pas de r\'etroaction sur la convection 583 583 584 584 REAL t_wake(klon,klev),q_wake(klon,klev) ! wake pour la convection … … 640 640 641 641 c--------Stochastic Boundary Layer Triggering: ALE_BL-------- 642 c---Propri étés du thermiques au LCL643 real zlcl_th(klon) ! Altitude du LCL calcul écontinument (pcon dans thermcell_main.F90)642 c---Propri\'et\'es du thermiques au LCL 643 real zlcl_th(klon) ! Altitude du LCL calcul\'e continument (pcon dans thermcell_main.F90) 644 644 real fraca0(klon) ! Fraction des thermiques au LCL 645 645 real w0(klon) ! Vitesse des thermiques au LCL 646 real w_conv(klon) ! Vitesse verticale de grande échelle au LCL646 real w_conv(klon) ! Vitesse verticale de grande \'echelle au LCL 647 647 real therm_tke_max0(klon) ! TKE dans les thermiques au LCL 648 648 real env_tke_max0(klon) ! TKE dans l'environnement au LCL … … 652 652 real ale_bl_stat(klon) 653 653 654 c---D éclenchement stochastique654 c---D\'eclenchement stochastique 655 655 integer :: tau_trig(klon) 656 656 real proba_notrig(klon) … … 664 664 665 665 c---Fermeture statistique 666 real alp_bl_det(klon) ! ALP d éterministe du thermique unique667 real alp_bl_fluct_m(klon) ! ALP li ée aux fluctuations de flux de masse sous-nuageux668 real alp_bl_fluct_tke(klon) ! ALP li ée aux fluctuations d'énergie cinétique sous-nuageuse669 real alp_bl_conv(klon) ! ALP li ée à grande échelle666 real alp_bl_det(klon) ! ALP d\'terministe du thermique unique 667 real alp_bl_fluct_m(klon) ! ALP li\'ee aux fluctuations de flux de masse sous-nuageux 668 real alp_bl_fluct_tke(klon) ! ALP li\'ee aux fluctuations d'\'energie cin\'etique sous-nuageuse 669 real alp_bl_conv(klon) ! ALP li\'ee \`a grande \'echelle 670 670 real alp_bl_stat(klon) ! ALP totale 671 671 … … 1363 1363 1364 1364 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1365 !! Un petit travail àfaire ici.1365 !! Un petit travail \`a faire ici. 1366 1366 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1367 1367 … … 1513 1513 ALLOCATE(iGCM(0), jGCM(0)) 1514 1514 end if 1515 ENDIF !debut 1515 else 1516 ALLOCATE(tabijGCM(0)) 1517 ALLOCATE(lonGCM(0), latGCM(0)) 1518 ALLOCATE(iGCM(0), jGCM(0)) 1519 ENDIF 1516 1520 1517 1521 DO i=1,klon … … 1526 1530 ! justement quand ok_orodr = false. 1527 1531 ! ce rugoro est utilise par la couche limite et fait double emploi 1528 ! avec les param étrisations spécifiques de Francois Lott.1532 ! avec les param\'etrisations sp\'ecifiques de Francois Lott. 1529 1533 ! DO i=1,klon 1530 1534 ! rugoro(i) = MAX(1.0e-05, zstd(i)*zsig(i)/2.0) … … 1932 1936 ! Calculs de l'orbite. 1933 1937 ! Necessaires pour le rayonnement et la surface (calcul de l'albedo). 1934 ! doit donc etre plac éavant radlwsw et pbl_surface1938 ! doit donc etre plac\'e avant radlwsw et pbl_surface 1935 1939 1936 1940 !!! jyg 17 Sep 2010 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 1961 1965 ! Pour une solarlong0=1000., on calcule un ensoleillement moyen sur 1962 1966 ! l'annee a partir d'une formule analytique. 1963 ! Cet ensoleillement est sym métrique autour de l'équateur et1967 ! Cet ensoleillement est sym\'etrique autour de l'\'equateur et 1964 1968 ! non nul aux poles. 1965 1969 IF (abs(solarlong0-1000.)<1.e-4) then … … 2232 2236 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 2233 2237 ! Modif FH 2010/04/27. Sans doute temporaire. 2234 ! Deux options pour le alp_offset : constant si > à 0 ou proportionnel Ãa2238 ! Deux options pour le alp_offset : constant si >?? 0 ou proportionnel ??a 2235 2239 ! w si <0 2236 2240 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 2628 2632 ENDIF 2629 2633 2630 c----Tirage al éatoire et calcul de ale_bl_trig2634 c----Tirage al\'eatoire et calcul de ale_bl_trig 2631 2635 do i=1,klon 2632 2636 if ( (ale_bl_stat(i) .gt. abs(cin(i))+1.e-10) ) then … … 2701 2705 c ============== 2702 2706 2703 ! Dans le cas o ùon active les thermiques, on fait partir l'ajustement2707 ! Dans le cas o\`u on active les thermiques, on fait partir l'ajustement 2704 2708 ! a partir du sommet des thermiques. 2705 2709 ! Dans le cas contraire, on demarre au niveau 1.
Note: See TracChangeset
for help on using the changeset viewer.