Changeset 101 in lmdz_wrf


Ignore:
Timestamp:
Jul 25, 2014, 3:47:31 PM (10 years ago)
Author:
lfita
Message:

Fixing zw2 < 0. Issue giving NaN due to SQRT(zw2)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_height.F90

    r100 r101  
    7575      CALL check_var3D(lfname, lvarname, zw2, ngrid, nlay+1, largest, .FALSE.)
    7676
     77! L. Fita, LMD July 2014. Fixing the If when zw2 == 0
    7778      do l=1,nlay
    7879         do ig=1,ngrid
    7980            if (l.le.lmax(ig)) then
    80                 if (zw2(ig,l).lt.0.)then
    81                   print*,'pb2 zw2<0'
    82                 endif
     81! Previous version
     82!                if (zw2(ig,l).lt.0.)then
     83!                  print*,'pb2 zw2<0'
     84!                endif
     85!                zw2(ig,l)=sqrt(zw2(ig,l))
     86!                wmax(ig)=max(wmax(ig),zw2(ig,l))
     87!            else
     88!                 zw2(ig,l)=0.
     89!            endif
     90! New version
     91              if (zw2(ig,l).lt.0.)then
     92                print*,'pb2 zw2<0'
     93                zw2(ig,l)=0.
     94                wmax(ig)=max(wmax(ig),zw2(ig,l))
     95              else
    8396                zw2(ig,l)=sqrt(zw2(ig,l))
    8497                wmax(ig)=max(wmax(ig),zw2(ig,l))
     98              endif 
    8599            else
    86100                 zw2(ig,l)=0.
    87101            endif
     102
    88103          enddo
    89104      enddo
Note: See TracChangeset for help on using the changeset viewer.