Changeset 2605 for LMDZ5/trunk/libf


Ignore:
Timestamp:
Jul 28, 2016, 11:46:18 AM (8 years ago)
Author:
jyg
Message:

Adding some initialization in cv3_routines.F90.
Should help guarantee 1+1=2.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/cv3_routines.F90

    r2579 r2605  
    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)
Note: See TracChangeset for help on using the changeset viewer.