Changeset 4674 for LMDZ6/trunk/libf/phylmd
- Timestamp:
- Sep 5, 2023, 12:13:50 PM (15 months ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/lmdz_cloudth.F90
r4664 r4674 602 602 INCLUDE "FCTTRE.h" 603 603 604 INTEGER itap,ind1,ind2 605 INTEGER ngrid,klev,klon,l,ig 604 integer, intent(in) :: ind2 605 integer, intent(in) :: ngrid,klev 606 607 real, dimension(ngrid,klev), intent(in) :: ztv 608 real, dimension(ngrid), intent(in) :: po 609 real, dimension(ngrid,klev), intent(in) :: zqta 610 real, dimension(ngrid,klev+1), intent(in) :: fraca 611 real, dimension(ngrid), intent(out) :: qcloud 612 real, dimension(ngrid,klev), intent(out) :: ctot 613 real, dimension(ngrid,klev), intent(out) :: ctot_vol 614 real, dimension(ngrid,klev), intent(in) :: zpspsk 615 real, dimension(ngrid,klev+1), intent(in) :: paprs 616 real, dimension(ngrid,klev), intent(in) :: pplay 617 real, dimension(ngrid,klev), intent(in) :: ztla 618 real, dimension(ngrid,klev), intent(inout) :: zthl 619 real, dimension(ngrid,klev), intent(in) :: ratqs 620 real, dimension(ngrid), intent(in) :: zqs 621 real, dimension(ngrid,klev), intent(in) :: t 606 622 real, dimension(ngrid,klev), intent(out) :: cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv 607 608 REAL ztv(ngrid,klev) 609 REAL po(ngrid) 623 624 610 625 REAL zqenv(ngrid) 611 REAL zqta(ngrid,klev)612 613 REAL fraca(ngrid,klev+1)614 REAL zpspsk(ngrid,klev)615 REAL paprs(ngrid,klev+1)616 REAL pplay(ngrid,klev)617 REAL ztla(ngrid,klev)618 REAL zthl(ngrid,klev)619 620 626 REAL zqsatth(ngrid,klev) 621 627 REAL zqsatenv(ngrid,klev) … … 628 634 REAL cth(ngrid,klev) 629 635 REAL cenv(ngrid,klev) 630 REAL ctot(ngrid,klev)631 636 REAL cth_vol(ngrid,klev) 632 637 REAL cenv_vol(ngrid,klev) 633 REAL ctot_vol(ngrid,klev)634 638 REAL rneb(ngrid,klev) 635 REAL t(ngrid,klev)636 639 REAL qsatmmussig1,qsatmmussig2,sqrt2pi,sqrt2,sqrtpi,pi 637 640 REAL rdd,cppd,Lv … … 641 644 REAL Tbef,zdelta,qsatbef,zcor 642 645 REAL qlbef 643 REAL ratqs(ngrid,klev) ! Determine the width of the vapour distribution644 646 REAL zpdf_sig(ngrid),zpdf_k(ngrid),zpdf_delta(ngrid) 645 647 REAL zpdf_a(ngrid),zpdf_b(ngrid),zpdf_e1(ngrid),zpdf_e2(ngrid) 646 REAL zqs(ngrid), qcloud(ngrid)647 648 REAL erf 648 649 650 651 INTEGER :: ind1,l, ig 649 652 650 653 IF (iflag_cloudth_vert.GE.1) THEN -
LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90
r4666 r4674 144 144 REAL, DIMENSION(klon,klev), INTENT(INOUT) :: zthl ! liquid potential temperature [K] 145 145 REAL, DIMENSION(klon,klev), INTENT(INOUT):: ratqs ! function of pressure that sets the large-scale 146 REAL, DIMENSION(klon,klev), INTENT( INOUT):: beta ! conversion rate of condensed water146 REAL, DIMENSION(klon,klev), INTENT(OUT):: beta ! conversion rate of condensed water 147 147 148 148 -
LMDZ6/trunk/libf/phylmd/lmdz_lscp_old.F90
r4670 r4674 76 76 REAL, DIMENSION(klon,klev), INTENT(IN) :: zqta, fraca 77 77 REAL, DIMENSION(klon,klev), INTENT(IN) :: zpspsk, ztla 78 REAL, DIMENSION(klon,klev), INTENT(IN ) :: zthl78 REAL, DIMENSION(klon,klev), INTENT(INOUT) :: zthl 79 79 ! 80 80 ! Input/output … … 173 173 !--------------------------------------------------------------- 174 174 ! 175 ! Fonctions en ligne: 176 ! 177 REAL :: fallvs,fallvc, zzz ! Vitesse de chute pour cristaux de glace 178 ! (Heymsfield & Donner, 1990) 179 fallvc (zzz) = 3.29/2.0 * ((zzz)**0.16) * ffallv_con 180 fallvs (zzz) = 3.29/2.0 * ((zzz)**0.16) * ffallv_lsc 181 ! 182 !--------------------------------------------------------------- 175 183 !AA Variables traceurs: 176 184 !AA Provisoire !!! Parametres alpha du lessivage … … 189 197 REAL, DIMENSION(klon) :: zlh_solid 190 198 REAL :: zm_solid 191 !--------------------------------------------------------------- 192 ! 193 ! Fonctions en ligne: 194 ! 195 REAL :: fallvs,fallvc,zzz ! Vitesse de chute pour cristaux de glace 196 ! (Heymsfield & Donner, 1990) 197 fallvc (zzz) = 3.29/2.0 * ((zzz)**0.16) * ffallv_con 198 fallvs (zzz) = 3.29/2.0 * ((zzz)**0.16) * ffallv_lsc 199 199 200 200 201 !
Note: See TracChangeset
for help on using the changeset viewer.