Changeset 027e844 in flexpart.git for src/calcpar.f90
- Timestamp:
- Oct 13, 2017, 11:59:48 AM (6 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 02095e3
- Parents:
- fe32dca
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/calcpar.f90
r6ecb30a r027e844 76 76 real :: ttlev(nuvzmax),qvlev(nuvzmax),obukhov,scalev,ol,hmixplus 77 77 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 79 79 real :: uuh(0:nxmax-1,0:nymax-1,nuvzmax) 80 80 real :: vvh(0:nxmax-1,0:nymax-1,nuvzmax) … … 136 136 137 137 ! 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) 140 140 else 141 141 llev=0 142 142 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) 144 144 end if 145 145
Note: See TracChangeset
for help on using the changeset viewer.