Changeset d8107c2 in flexpart.git for src/netcdf_output_mod.f90


Ignore:
Timestamp:
Jan 13, 2016, 5:02:07 PM (8 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:
41d8574
Parents:
f75967d
Message:

Added variables volume/area to netcdf output. Fixed a bug where output of dry deposition where wrong for the MPI version.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/netcdf_output_mod.f90

    rcbad0f1 rd8107c2  
    244244  integer :: nspecDimID, npointDimID, nageclassDimID, ncharDimID, pointspecDimID
    245245  integer :: tID, lonID, latID, levID, poleID, lageID, oroID
     246  integer :: volID, areaID
    246247  integer :: rellng1ID, rellng2ID, rellat1ID, rellat2ID, relzz1ID, relzz2ID
    247248  integer :: relcomID, relkindzID, relstartID, relendID, relpartID, relxmassID
     
    267268  ! otherwise give an error which can look confusing).
    268269  ! *********************************************************************
    269   ! open(newunit=test_unit,file=trim(path(2)(1:length(2)))//'test_dir.txt',status='replace',&
    270   !      &err=100)
    271   !  close (test_unit, status='delete')
    272270  open(unit=unittmp,file=trim(path(2)(1:length(2)))//'test_dir.txt',status='replace',&
    273271       &err=100)
     
    341339  call nf90_err(nf90_def_dim(ncid, 'numpoint', numpoint, npointDimID))
    342340
     341
    343342  ! create variables
    344343  !*************************
     
    378377  call nf90_err(nf90_put_att(ncid, levID, 'standard_name', 'height'))
    379378  call nf90_err(nf90_put_att(ncid, levID, 'long_name', 'height above ground'))
     379
     380  ! volume
     381  call nf90_err(nf90_def_var(ncid, 'volume', nf90_float, (/ lonDimID, latDimID, levDimID /), &
     382       volID))
     383
     384  ! area
     385  call nf90_err(nf90_def_var(ncid, 'area', nf90_float, (/ lonDimID, latDimID /), &
     386       areaID))
     387
    380388
    381389  if (write_releases.eqv..true.) then
     
    612620  ! levels
    613621  call nf90_err(nf90_put_var(ncid, levID, outheight(1:numzgrid)))
    614  
     622
     623  ! volume
     624  call nf90_err(nf90_put_var(ncid, volID, volume(:,:,:)))
     625
     626  ! area
     627  call nf90_err(nf90_put_var(ncid, areaID, area(:,:)))
     628
    615629  if (write_releases.eqv..true.) then
    616630    ! release point information
     
    889903              if (mpi_mode.gt.0) then
    890904                do l=1,nclassunc
    891                   auxgrid(l)=drygridunc(ix,jy,ks,kp,l,nage)
     905                  auxgrid(l)=drygridunc0(ix,jy,ks,kp,l,nage)
    892906                end do
    893907              else
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG