Changeset f28aa0a in flexpart.git


Ignore:
Timestamp:
Jun 22, 2016, 1:22:20 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:
462f74b, 0f7835d, 6b3dad4
Parents:
79abee9
Message:

netCDF output; added (hardcoded) option for writing grid cell volume/area

Files:
3 edited

Legend:

Unmodified
Added
Removed
  • options/COMMAND

    r79abee9 rf28aa0a  
    88 LOUTAVER=       3600,
    99 LOUTSAMPLE=        900,
    10  ITSPLIT=  99997200,
     10 ITSPLIT=  99999999,
    1111 LSYNCTIME=        900,
    1212 CTL= -5.0000000    ,
  • src/netcdf_output_mod.f90

    r341f4b7 rf28aa0a  
    4040  !  - Deposition fields can be calculated in double precision, see variable
    4141  !    'dep_prec' in par_mod
     42  !  - Hardcoded options 'write_vol' and 'write_area' for grid cell
     43  !    volume and area
    4244  !*****************************************************************************
    4345
     
    106108  ! switch output of release point information on/off
    107109  logical, parameter :: write_releases = .true.
     110
     111  ! switch output of grid cell volume and area on/off
     112  logical, parameter :: write_vol = .false.
     113  logical, parameter :: write_area = .false.
    108114
    109115contains
     
    390396
    391397  ! volume
    392   call nf90_err(nf90_def_var(ncid, 'volume', nf90_float, (/ lonDimID, latDimID, levDimID /), &
    393        volID))
    394 
     398  if (write_vol) call nf90_err(nf90_def_var(ncid, 'volume', nf90_float, &
     399       &(/ lonDimID, latDimID, levDimID /), volID))
    395400  ! area
    396   call nf90_err(nf90_def_var(ncid, 'area', nf90_float, (/ lonDimID, latDimID /), &
    397        areaID))
     401  if (write_area) call nf90_err(nf90_def_var(ncid, 'area', nf90_float, &
     402       &(/ lonDimID, latDimID /), areaID))
    398403
    399404
     
    633638
    634639  ! volume
    635   if (lnest) then
    636     call nf90_err(nf90_put_var(ncid, volID, volumen(:,:,:)))
    637   else
    638     call nf90_err(nf90_put_var(ncid, volID, volume(:,:,:)))
     640  if (write_vol) then
     641    if (lnest) then
     642      call nf90_err(nf90_put_var(ncid, volID, volumen(:,:,:)))
     643    else
     644      call nf90_err(nf90_put_var(ncid, volID, volume(:,:,:)))
     645    end if
    639646  end if
    640647
    641648  ! area
    642   if (lnest) then
    643     call nf90_err(nf90_put_var(ncid, areaID, arean(:,:)))
    644   else
    645     call nf90_err(nf90_put_var(ncid, areaID, area(:,:)))
     649  if (write_area) then
     650    if (lnest) then
     651      call nf90_err(nf90_put_var(ncid, areaID, arean(:,:)))
     652    else
     653      call nf90_err(nf90_put_var(ncid, areaID, area(:,:)))
     654    end if
    646655  end if
    647656
  • src/pathnames

    rd6a0977 rf28aa0a  
    44/xnilu_wrk/flex_wrk/WIND_FIELDS/AVAILABLE_ECMWF_OPER_fields_global
    55============================================
    6 
    7 /xnilu_wrk/flex_wrk/zhg/ECMWF_CLWC/Availables
    8 
    9 
    10  
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG