Changeset 242571d in flexpart.git
- Timestamp:
- Mar 1, 2015, 3:57:22 PM (9 years ago)
- Branches:
- master, 10.4.1_pesei, FPv9.3.1, FPv9.3.1b_testing, FPv9.3.2, GFS_025, bugfixes+enhancements, dev, fp9.3.1-20161214-nc4, grib2nc4_repair, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 18c1336
- Parents:
- 43c8684
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/readreleases.f90
rb4d29ce r242571d 126 126 ! prepare namelist output if requested 127 127 if (nmlout.eqv..true.) then 128 open(unitreleasesout,file=path(2)(1:length(2))//'RELEASES.namelist',access='append',status='new',err=1000) 128 !open(unitreleasesout,file=path(2)(1:length(2))//'RELEASES.namelist',access='append',status='new',err=1000) 129 open(unitreleasesout,file=path(2)(1:length(2))//'RELEASES.namelist',err=1000) 129 130 endif 130 131 … … 258 259 if (stat.ne.0) write(*,*)'ERROR: could not allocate xmasssave' 259 260 260 write (*,*) ' Releasepoints : ', numpoint261 write (*,*) 'readreleases> Releasepoints : ', numpoint 261 262 262 263 do i=1,numpoint … … 293 294 294 295 do i=1,nspec 296 if (verbosity.gt.0) then 297 print*, 'readreleases> call readspecies', i 298 endif 299 295 300 if (readerror.ne.0) then 296 301 read(unitreleases,*,err=998) specnum_rel(i) … … 494 499 endif ! if namelist format 495 500 501 502 if (verbosity.gt.1 .and. numpoint.eq.1) then ! verbosity 2 or larger 503 write(*,*) 'numpoint=', numpoint 504 print*, id1,it1 505 print*, id2,it2 506 print*, xpoint1(numpoint) 507 print*, ypoint1(numpoint) 508 print*, xpoint2(numpoint) 509 print*, ypoint2(numpoint) 510 print*, 'kindz=' , kindz(numpoint) 511 print*, zpoint1(numpoint) 512 print*, zpoint2(numpoint) 513 print*, npart(numpoint) 514 do i=1,nspec 515 !mass(i)= 516 print*, 'xmass=', xmass(numpoint,i) 517 end do 518 print*, compoint(numpoint) 519 endif 520 521 496 522 ! If a release point contains no particles, stop and issue error message 497 523 !*********************************************************************** … … 536 562 write(*,*) 'after simulation stops.' 537 563 write(*,*) 'Make files COMMAND and RELEASES consistent.' 564 write(*,*) jul1, ' < ' , bdate 565 write(*,*) ' .or. ' 566 write(*,*) jul2 , ' > ', edate 567 538 568 stop 539 569 endif … … 563 593 endif 564 594 595 if (verbosity.gt.1 .and. numpoint.eq.1) then ! verbosity 2 or larger 596 print*, 'ireleasestart(',numpoint,')', ireleasestart(numpoint) 597 print*, 'ireleaseend(',numpoint,')', ireleaseend(numpoint) 598 endif 599 565 600 ! Determine the release rate (particles per second) and total number 566 601 ! of particles released during the simulation
Note: See TracChangeset
for help on using the changeset viewer.