Ignore:
Timestamp:
Sep 29, 2016, 11:26:46 PM (8 years ago)
Author:
Laurent Fairhead
Message:

Merged trunk changes r2593:2640 into testing branch

Location:
LMDZ5/branches/testing
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/branches/testing

  • LMDZ5/branches/testing/libf/phylmd/cv3_routines.F90

    r2594 r2641  
    238238! gz = phi at the full levels (same as p).
    239239
    240   DO i = 1, len
    241     gz(i, 1) = 0.0
    242   END DO
     240!!  DO i = 1, len                    !jyg
     241!!    gz(i, 1) = 0.0                 !jyg
     242!!  END DO                           !jyg
     243    gz(:,:) = 0.                     !jyg: initialization of the whole array
    243244! ori      do 140 k=2,nlp
    244245  DO k = 2, nl ! convect3
     
    409410    IF (ok_new_feed) THEN
    410411      IF (iter==niter) THEN
    411         DO k = minorig, nl
     412        DO i = 1,len                         !jyg
     413          pfeedmin(i) = ph(i,minorig+1)      !jyg
     414        ENDDO                                !jyg
     415        DO k = minorig+1, nl                 !jyg
     416!!        DO k = minorig, nl                 !jyg
    412417          DO i = 1, len
    413418            IF (ph(i,k)>=plclfeed(i)) pfeedmin(i) = ph(i, k)
     
    10551060                         p, ph, h, tv, lv, lf, pbase, buoybase, plcl, &
    10561061                         inb, tp, tvp, clw, hp, ep, sigp, buoy, frac)
     1062  USE print_control_mod, ONLY: prt_level
    10571063  IMPLICIT NONE
    10581064
     
    11131119  REAL fracg
    11141120  REAL deltap
     1121
     1122  IF (prt_level >= 10) THEN
     1123    print *,'cv3_undilute2.0. t(1,k), q(1,k), qs(1,k) ', &
     1124                        (k, t(1,k), q(1,k), qs(1,k), k = 1,nl)
     1125  ENDIF
    11151126
    11161127! =====================================================================
     
    12991310    END DO
    13001311  END DO
     1312
     1313  IF (prt_level >= 10) THEN
     1314    print *,'cv3_undilute2.1. tp(1,k), tvp(1,k) ', &
     1315                          (k, tp(1,k), tvp(1,k), k = 1,nl)
     1316  ENDIF
    13011317
    13021318! =====================================================================
Note: See TracChangeset for help on using the changeset viewer.