Changeset 1372 for LMDZ4


Ignore:
Timestamp:
Apr 27, 2010, 12:56:22 PM (14 years ago)
Author:
musat
Message:

Bug fix: overflow for pplay field
IM

Location:
LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcellV0_main.F90

    r1371 r1372  
    593593      pcon(ig)=pplay(ig,1)*(zo(ig,1)/zqsat(ig,1))**CHI
    594594      enddo
    595       do k=1,nlay
     595!IM   do k=1,nlay
     596      do k=1,nlay-1
    596597         do ig=1,ngrid
    597598         if ((pcon(ig).le.pplay(ig,k))  &
    … …  
    600601         endif
    601602         enddo
     603      enddo
     604!IM
     605      do ig=1,ngrid
     606        if (pcon(ig).le.pplay(ig,nlay)) then
     607           zcon2(ig)=zlay(ig,nlay)-(pcon(ig)-pplay(ig,nlay))/(RG*rho(ig,nlay))/100.
     608           abort_message = 'thermcellV0_main: les thermiques vont trop haut '
     609           CALL abort_gcm (modname,abort_message,1)
     610        endif
    602611      enddo
    603612      if (prt_level.ge.1) print*,'14b OK convect8'

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_main.F90

    r1338 r1372  
    610610      pcon(ig)=pplay(ig,1)*(zo(ig,1)/zqsat(ig,1))**CHI
    611611      enddo
    612       do k=1,nlay
     612!IM   do k=1,nlay
     613      do k=1,nlay-1
    613614         do ig=1,ngrid
    614615         if ((pcon(ig).le.pplay(ig,k))  &
    … …  
    617618         endif
    618619         enddo
     620      enddo
     621!IM
     622      do ig=1,ngrid
     623        if (pcon(ig).le.pplay(ig,nlay)) then
     624           zcon2(ig)=zlay(ig,nlay)-(pcon(ig)-pplay(ig,nlay))/(RG*rho(ig,nlay))/100.
     625           abort_message = 'thermcellV0_main: les thermiques vont trop haut '
     626           CALL abort_gcm (modname,abort_message,1)
     627        endif
    619628      enddo
    620629      if (prt_level.ge.1) print*,'14b OK convect8'
Note: See TracChangeset for help on using the changeset viewer.