Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp_old.F90
- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp_old.F90
r5153 r5158 757 757 IF (iflag_pdf==0) THEN 758 758 ! version creneau de (Li, 1998) 759 doi=1,klon759 DO i=1,klon 760 760 zdelq = min(ratqs(i,k),0.99) * zq(i) 761 761 rneb(i,k) = (zq(i)+zdelq-zqs(i)) / (2.0*zdelq) … … 768 768 ! les valeurs de T et Q initiales 769 769 ! ---------------------------------------------------------------- 770 doi=1,klon770 DO i=1,klon 771 771 IF(zq(i)<1.e-15) THEN 772 772 ncoreczq=ncoreczq+1 … … 800 800 801 801 endif 802 doi=1,klon802 DO i=1,klon 803 803 rneb(i,k)=ctot(i,k) 804 804 rneblsvol(i,k)=ctot_vol(i,k) … … 821 821 !CR: variation de qsat avec T en presence de glace ou non 822 822 !initialisations 823 doi=1,klon823 DO i=1,klon 824 824 DT(i) = 0. 825 825 n_i(i)=0 … … 846 846 ! Iteration pour condensation avec variation de qsat(T) 847 847 ! ----------------------------------------------------- 848 doiter=1,iflag_fisrtilp_qsat+1848 DO iter=1,iflag_fisrtilp_qsat+1 849 849 850 doi=1,klon850 DO i=1,klon 851 851 ! do while ((abs(DT(i)).gt.DDT0).OR.(n_i(i).EQ.0)) 852 852 ! !! convergence = .TRUE. tant que l'on n'a pas converge !! … … 929 929 IF (.NOT. ice_thermo) THEN 930 930 ! -------------------------- 931 doi=1,klon931 DO i=1,klon 932 932 IF ((convergence(i).OR.(n_i(i)==0)).AND.lognormale(i)) THEN 933 933 qlbef(i)=max(0.,zqn(i)-zqs(i)) … … 949 949 endif 950 950 951 doi=1,klon951 DO i=1,klon 952 952 IF ((convergence(i).OR.(n_i(i)==0)).AND.lognormale(i)) THEN 953 953 IF (iflag_t_glace==0) THEN
Note: See TracChangeset
for help on using the changeset viewer.