Changeset 3824 for palm/trunk/SOURCE/chem_photolysis_mod.f90

Timestamp:

Mar 27, 2019 3:56:16 PM (5 years ago)

Author:

pavelkrc

Message:

Code review of radiation_model_mod.f90

File:

• palm/trunk/SOURCE/chem_photolysis_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/chem_photolysis_mod.f90

@@ -234 +234 @@
 
        USE radiation_model_mod,                                                   &
-           ONLY:  calc_zenith, zenith
+           ONLY:  calc_zenith, cos_zenith
 
        IMPLICIT NONE
@@ -247 +247 @@
        CALL calc_zenith
 
-       IF ( zenith(0) > 0.0_wp ) THEN
-          coszi = 1. / zenith(0)
+       IF ( cos_zenith > 0.0_wp ) THEN
+          coszi = 1. / cos_zenith
 
           DO iphot = 1,nphot
@@ -254 +254 @@
                 IF ( TRIM( names_s(iav) ) == TRIM( phot_names(iphot) ) ) then
                    phot_frequen(iphot)%freq(nzb+1:nzt,:,:) =    &
-                                par_l(iav) * zenith(0)**par_m(iav) *  EXP( -par_n(iav) * coszi ) 
+                                par_l(iav) * cos_zenith**par_m(iav) *  EXP( -par_n(iav) * coszi ) 
                 ENDIF
              ENDDO