Changeset 3859 for palm/trunk


Ignore:
Timestamp:
Apr 3, 2019 8:30:31 PM (2 years ago)
Author:
maronga
Message:

comments in radiation model updated, minor bugfix in palm_csd

Location:
palm/trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/palm_csd

    r3773 r3859  
    2525# -----------------
    2626# $Id$
     27# Bugfix: wrong variable naming for 'y'
     28#
     29# 3773 2019-03-01 08:56:57Z maronga
    2730# Unspecificed changes
    2831#
     
    566569   nc_write_attribute(filename[i], 'x', 'standard_name','projection_x_coordinate')   
    567570   nc_write_attribute(filename[i], 'x', 'units', 'm')
    568    nc_write_attribute(filename[i], 'y', 'long_name', 'x')
     571   nc_write_attribute(filename[i], 'y', 'long_name', 'y')
    569572   nc_write_attribute(filename[i], 'y', 'standard_name', 'projection_y_coordinate')   
    570573   nc_write_attribute(filename[i], 'y', 'units', 'm')
  • palm/trunk/SOURCE/radiation_model_mod.f90

    r3847 r3859  
    2828! -----------------
    2929! $Id$
     30! Added some descriptions
     31!
     32! 3847 2019-04-01 14:51:44Z suehring
    3033! Implement check for dt_radiation (must be > 0)
    3134!
     
    831834
    832835!
    833 !-- Flag parameters for RRTMGS (should not be changed)
     836!-- Flag parameters to be passed to RRTMG (should not be changed until ice phase in clouds is allowed)
    834837    INTEGER(iwp), PARAMETER :: rrtm_idrv     = 1, & !< flag for longwave upward flux calculation option (0,1)
    835838                               rrtm_inflglw  = 2, & !< flag for lw cloud optical properties (0,1,2)
     
    34923495
    34933496!
    3494 !--       Avoid temperature/humidity jumps at the top of the LES domain by
    3495 !--       linear interpolation from nzt+2 to nzt+7
     3497!--       Avoid temperature/humidity jumps at the top of the PALM domain by
     3498!--       linear interpolation from nzt+2 to nzt+7. Jumps are induced by
     3499!--       discrepancies between the values in the  domain and those above that
     3500!--       are prescribed in RRTMG
    34963501          DO k = nzt+2, nzt+7
    34973502             rrtm_tlay(0,k) = rrtm_tlay(0,nzt+1)                            &
     
    35073512          ENDDO
    35083513
    3509 !--       Linear interpolate to zw grid
     3514!--       Linear interpolate to zw grid. Loop reaches one level further up
     3515!--       due to the staggered grid in RRTMG
    35103516          DO k = nzb+2, nzt+8
    35113517             rrtm_tlev(0,k)   = rrtm_tlay(0,k-1) + (rrtm_tlay(0,k) -        &
     
    46144620! Description:
    46154621! ------------
    4616 !> Read trace gas data from file
     4622!> Read trace gas data from file and convert into trace gas paths / volume
     4623!> mixing ratios. If a user-defined input file is provided it needs to follow
     4624!> the convections used in RRTMG (see respective netCDF files shipped with
     4625!> RRTMG)
    46174626!------------------------------------------------------------------------------!
    46184627    SUBROUTINE read_trace_gas_data
     
    46384647                       id_var    !< NetCDf id ot the absorber
    46394648
    4640        REAL(wp) :: p_mls_l, p_mls_u, p_wgt_l, p_wgt_u, p_mls_m
     4649       REAL(wp) :: p_mls_l, &    !< pressure lower limit for interpolation
     4650                   p_mls_u, &    !< pressure upper limit for interpolation
     4651                   p_wgt_l, &    !< pressure weight lower limit for interpolation
     4652                   p_wgt_u, &    !< pressure weight upper limit for interpolation
     4653                   p_mls_m       !< mean pressure between upper and lower limits
    46414654
    46424655
Note: See TracChangeset for help on using the changeset viewer.