Changeset f13406c in flexpart.git for src/readreleases.f90
- Timestamp:
- Dec 23, 2013, 6:23:38 PM (10 years ago)
- Branches:
- master, 10.4.1_pesei, FPv9.3.1, FPv9.3.1b_testing, FPv9.3.2, GFS_025, NetCDF, bugfixes+enhancements, deposition, dev, fp9.3.1-20161214-nc4, grib2nc4_repair, inputlist, laptop, release-10, release-10.4.1, scaling-bug, svn-petra, svn-trunk, univie
- Children:
- 4fbe7a5
- Parents:
- e200b7a
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/readreleases.f90
re200b7a rf13406c 57 57 ! weta, wetb parameters to determine the wet scavenging coefficient * 58 58 ! zpoint1,zpoint2 height range, over which release takes place * 59 ! num_min_discrete if less, release cannot be randomized and happens at * 60 ! time mid-point of release interval * 59 61 ! * 60 62 !***************************************************************************** … … 68 70 69 71 integer :: numpartmax,i,j,id1,it1,id2,it2,specnum_rel,idum,stat 72 integer,parameter :: num_min_discrete=100 70 73 real :: vsh(ni),fracth(ni),schmih(ni),releaserate,xdum,cun 71 real(kind=dp) :: jul1,jul2,jul date74 real(kind=dp) :: jul1,jul2,julm,juldate 72 75 character(len=50) :: line 73 76 logical :: old … … 376 379 jul1=juldate(id1,it1) 377 380 jul2=juldate(id2,it2) 381 julm=(jul1+jul2)/2. 378 382 if (jul1.gt.jul2) then 379 383 write(*,*) 'FLEXPART MODEL ERROR' … … 391 395 stop 392 396 endif 393 ireleasestart(numpoint)=int((jul1-bdate)*86400.) 394 ireleaseend(numpoint)=int((jul2-bdate)*86400.) 397 if (npart(numpoint).gt.num_min_discrete) then 398 ireleasestart(numpoint)=int((jul1-bdate)*86400.) 399 ireleaseend(numpoint)=int((jul2-bdate)*86400.) 400 else 401 ireleasestart(numpoint)=int((julm-bdate)*86400.) 402 ireleaseend(numpoint)=int((julm-bdate)*86400.) 403 endif 395 404 else if (ldirect.eq.-1) then 396 405 if ((jul1.lt.edate).or.(jul2.gt.bdate)) then … … 401 410 stop 402 411 endif 403 ireleasestart(numpoint)=int((jul1-bdate)*86400.) 404 ireleaseend(numpoint)=int((jul2-bdate)*86400.) 412 if (npart(numpoint).gt.num_min_discrete) then 413 ireleasestart(numpoint)=int((jul1-bdate)*86400.) 414 ireleaseend(numpoint)=int((jul2-bdate)*86400.) 415 else 416 ireleasestart(numpoint)=int((julm-bdate)*86400.) 417 ireleaseend(numpoint)=int((julm-bdate)*86400.) 418 endif 405 419 endif 406 420 endif
Note: See TracChangeset
for help on using the changeset viewer.