- Timestamp:
- Jul 10, 2023, 12:24:05 AM (17 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/cloudth_mod.F90
r4613 r4617 30 30 INTEGER itap,ind1,ind2 31 31 INTEGER ngrid,klev,klon,l,ig 32 INTEGER iflag_ratqs32 !ym INTEGER iflag_ratqs 33 33 34 34 REAL ztv(ngrid,klev) … … 275 275 INTEGER itap,ind1,ind2 276 276 INTEGER ngrid,klev,klon,l,ig 277 INTEGER iflag_ratqs277 !ym INTEGER iflag_ratqs 278 278 279 279 REAL ztv(ngrid,klev) … … 615 615 INTEGER itap,ind1,ind2 616 616 INTEGER ngrid,klev,klon,l,ig 617 INTEGER iflag_ratqs617 !ym INTEGER iflag_ratqs 618 618 619 619 REAL ztv(ngrid,klev) … … 838 838 INTEGER itap,ind1,ind2 839 839 INTEGER ngrid,klev,klon,l,ig 840 INTEGER iflag_ratqs840 !ym INTEGER iflag_ratqs 841 841 842 842 REAL ztv(ngrid,klev) … … 1026 1026 ENDIF 1027 1027 sigma1s = sigma1s_fraca + sigma1s_ratqs 1028 IF (iflag_ratqs.eq.11) then1029 sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv1030 ENDIF1028 !ym IF (iflag_ratqs.eq.11) then 1029 !ym sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv 1030 !ym ENDIF 1031 1031 sigma2s=(sigma2s_factor*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**sigma2s_power))+0.002*zqta(ind1,ind2) 1032 1032 ! tests … … 1607 1607 1608 1608 INTEGER itap,ind1,l,ig,iter,k 1609 INTEGER iflag_ratqs1609 ! INTEGER iflag_ratqs 1610 1610 INTEGER iflag_topthermals, niter 1611 1611 LOGICAL falseklon(klon) … … 1902 1902 1903 1903 sigma1s = sigma1s_fraca + sigma1s_ratqs 1904 IF (iflag_ratqs.eq.11) then1905 sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv1906 ENDIF1904 !ym IF (iflag_ratqs.eq.11) then 1905 !ym sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv 1906 !ym ENDIF 1907 1907 sigma2s=(sigma2s_factor*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**sigma2s_power))+0.002*zqta(ind1,ind2) 1908 1908 … … 2147 2147 2148 2148 sigma1s = sigma1s_fraca + sigma1s_ratqs 2149 IF (iflag_ratqs.eq.11) then2150 sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv2151 ENDIF2152 IF (iflag_ratqs.eq.11) then2153 2154 ENDIF2149 !ym IF (iflag_ratqs.eq.11) then 2150 !ym sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv 2151 !ym ENDIF 2152 !ym IF (iflag_ratqs.eq.11) then 2153 !ym sigma1s = ratqs(ind1,ind2)*po(ind1)*aenvl 2154 !ym ENDIF 2155 2155 deltasenv=aenvl*vert_alpha*sigma1s 2156 2156 xenv1=-(senvl+deltasenv)/(sqrt(2.)*sigma1s)
Note: See TracChangeset
for help on using the changeset viewer.