Changeset 72d3a5a in flexpart.git


Ignore:
Timestamp:
May 31, 2017, 10:11:57 AM (2 years ago)
Author:
Espen Sollum ATMOS <eso@…>
Branches:
dev, release-10, univie
Children:
e52967c
Parents:
6d420c3
Message:

Bugfix in verttransform (call to function ew)

Location:
src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • src/verttransform.f90

    r6481010 r72d3a5a  
    5050! added the variable cloud for use with scavenging - descr. in com_mod
    5151!*****************************************************************************
     52! Date: 2017-05-30 modification of a bug in ew. Don Morton (CTBTO project)   *
     53!*****************************************************************************
    5254!                                                                            *
    5355! Variables:                                                                 *
     
    155157
    156158
    157     tvold(ixm,jym)=tt2(ixm,jym,1,n)*(1.+0.378*ew*(td2(ixm,jym,1,n))/ &
     159    tvold(ixm,jym)=tt2(ixm,jym,1,n)*(1.+0.378*ew(td2(ixm,jym,1,n))/ &
    158160         ps(ixm,jym,1,n))
    159161    pold(ixm,jym)=ps(ixm,jym,1,n)
     
    215217  do jy=0,nymin1
    216218    do ix=0,nxmin1
    217       tvold(ix,jy)=tt2(ix,jy,1,n)*(1.+0.378*ew*(td2(ix,jy,1,n))/ &
     219      tvold(ix,jy)=tt2(ix,jy,1,n)*(1.+0.378*ew(td2(ix,jy,1,n))/ &
    218220           ps(ix,jy,1,n))
    219221    enddo
  • src/verttransform_nests.f90

    r341f4b7 r72d3a5a  
    5656!  the actual field dimensions
    5757!*****************************************************************************
     58! Date: 2017-05-30 modification of a bug in ew. Don Morton (CTBTO project)   *
     59!*****************************************************************************
    5860!                                                                            *
    5961! Variables:                                                                 *
     
    107109    do jy=0,nyn(l)-1
    108110      do ix=0,nxn(l)-1
    109         tvold(ix,jy)=tt2n(ix,jy,1,n,l)*(1.+0.378*ew*(td2n(ix,jy,1,n,l))/ &
     111        tvold(ix,jy)=tt2n(ix,jy,1,n,l)*(1.+0.378*ew(td2n(ix,jy,1,n,l))/ &
    110112             psn(ix,jy,1,n,l))
    111113      end do
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG