Ignore:
Timestamp:
Aug 19, 2015, 1:45:31 PM (9 years ago)
Author:
slebonnois
Message:

SL: bug corrections

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.TITAN/Tools/tmc.F90

    r1454 r1468  
    281281  do ilev=1,altlength
    282282! Need to be consistent with GCM computations
    283 !    if (za(ilon,ilat,ilev,itim).ne.miss_val) then
     283    if (vitu(ilon,ilat,ilev,itim).ne.miss_val) then
    284284     rayon(ilon,ilat,ilev,itim) = a0
    285285!     rayon(ilon,ilat,ilev,itim) = za(ilon,ilat,ilev,itim) + a0
    286286      grav(ilon,ilat,ilev,itim) = g0*a0*a0 &
    287287                 /(rayon(ilon,ilat,ilev,itim)*rayon(ilon,ilat,ilev,itim))
    288 !    else
    289 !     rayon(ilon,ilat,ilev,itim) = miss_val
    290 !      grav(ilon,ilat,ilev,itim) = miss_val
    291 !    endif
     288    else
     289     rayon(ilon,ilat,ilev,itim) = miss_val
     290      grav(ilon,ilat,ilev,itim) = miss_val
     291    endif
    292292  enddo
    293293 enddo
Note: See TracChangeset for help on using the changeset viewer.