Changeset 027e844 in flexpart.git for src


Ignore:
Timestamp:
Oct 13, 2017, 11:59:48 AM (6 years ago)
Author:
Espen Sollum ATMOS <eso@…>
Branches:
master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
Children:
02095e3
Parents:
fe32dca
Message:

Bugfix for calcpar.

Location:
src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • src/calcpar.f90

    r6ecb30a r027e844  
    7676  real :: ttlev(nuvzmax),qvlev(nuvzmax),obukhov,scalev,ol,hmixplus
    7777  real :: ulev(nuvzmax),vlev(nuvzmax),ew,rh,vd(maxspec),subsceff,ylat
    78   real :: altmin,tvold,pold,zold,pint,tv,zlev(nuvzmax),hmixdummy
     78  real :: altmin,tvold,pold,zold,pint,tv,zlev(nuvzmax),hmixdummy,akzdummy
    7979  real :: uuh(0:nxmax-1,0:nymax-1,nuvzmax)
    8080  real :: vvh(0:nxmax-1,0:nymax-1,nuvzmax)
     
    136136
    137137        ! calculate inverse Obukhov length scale with tth(llev)
    138       ol=obukhov(ps(ix,jy,1,n),tt2(ix,jy,1,n),td2(ix,jy,1,n), &
    139             tth(ix,jy,llev,n),ustar(ix,jy,1,n),sshf(ix,jy,1,n),akm,bkm,akz(llev),metdata_format)
     138        ol=obukhov(ps(ix,jy,1,n),tt2(ix,jy,1,n),td2(ix,jy,1,n), &
     139             tth(ix,jy,llev,n),ustar(ix,jy,1,n),sshf(ix,jy,1,n),akm,bkm,akz(llev),metdata_format)
    140140      else
    141141        llev=0
    142142        ol=obukhov(ps(ix,jy,1,n),tt2(ix,jy,1,n),td2(ix,jy,1,n), &
    143             tth(ix,jy,2,n),ustar(ix,jy,1,n),sshf(ix,jy,1,n),akm,bkm,akz(llev),metdata_format)
     143            tth(ix,jy,2,n),ustar(ix,jy,1,n),sshf(ix,jy,1,n),akm,bkm,akzdummy,metdata_format)
    144144      end if
    145145
  • src/par_mod.f90

    rfe32dca r027e844  
    218218  !**************************************************
    219219
    220   integer,parameter :: maxpart=300000000
     220  integer,parameter :: maxpart=3000000
    221221  integer,parameter :: maxspec=1
    222222
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG