Ignore:
Timestamp:
Mar 27, 2019 3:56:16 PM (2 years ago)
Author:
pavelkrc
Message:

Code review of radiation_model_mod.f90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/chem_photolysis_mod.f90

    r3614 r3824  
    234234
    235235       USE radiation_model_mod,                                                   &
    236            ONLY:  calc_zenith, zenith
     236           ONLY:  calc_zenith, cos_zenith
    237237
    238238       IMPLICIT NONE
     
    247247       CALL calc_zenith
    248248
    249        IF ( zenith(0) > 0.0_wp ) THEN
    250           coszi = 1. / zenith(0)
     249       IF ( cos_zenith > 0.0_wp ) THEN
     250          coszi = 1. / cos_zenith
    251251
    252252          DO iphot = 1,nphot
     
    254254                IF ( TRIM( names_s(iav) ) == TRIM( phot_names(iphot) ) ) then
    255255                   phot_frequen(iphot)%freq(nzb+1:nzt,:,:) =    &
    256                                 par_l(iav) * zenith(0)**par_m(iav) *  EXP( -par_n(iav) * coszi )
     256                                par_l(iav) * cos_zenith**par_m(iav) *  EXP( -par_n(iav) * coszi )
    257257                ENDIF
    258258             ENDDO
Note: See TracChangeset for help on using the changeset viewer.