Changeset 3004 for palm/trunk/SOURCE/uv_exposure_model_mod.f90

Timestamp:

Apr 27, 2018 12:33:25 PM (6 years ago)

Author:

Giersch

Message:

precipitation_rate removed, further allocation checks for data output of averaged quantities implemented, double CALL of flow_statistics at the beginning of time_integration removed, further minor bugfixes, comments added

File:

• palm/trunk/SOURCE/uv_exposure_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/uv_exposure_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Further allocation checks implemented (averaged data will be assigned to fill
+! values if no allocation happened so far) 
+! 
+! 2932 2018-03-26 09:39:22Z maronga
 ! renamed uvexposure_par to biometeorology_parameters
 ! 
@@ -331 +335 @@
     LOGICAL      ::  two_d !< flag parameter that indicates 2D variables (horizontal cross sections)
 
+    REAL(wp) ::  fill_value = -999.0_wp    !< value for the _FillValue attribute
+
     REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf !< 
 
@@ -353 +359 @@
 
        CASE ( 'uvem_vitd3dose*_xy' )        ! 2d-array
+          IF ( .NOT. ALLOCATED( vitd3_exposure_av ) ) THEN
+             ALLOCATE( vitd3_exposure_av(nysg:nyng,nxlg:nxrg) )
+             vitd3_exposure_av = REAL( fill_value, KIND = wp )
+          ENDIF
           IF ( av == 1 )  THEN
              DO  i = nxl, nxr
@@ -502 +512 @@
 
           CASE ( 'uvem_vitd3dose*' )
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   vitd3_exposure_av(j,i) = vitd3_exposure_av(j,i) + vitd3_exposure(j,i)
+             IF ( ALLOCATED( vitd3_exposure_av ) ) THEN
+                DO  i = nxlg, nxrg
+                   DO  j = nysg, nyng
+                      vitd3_exposure_av(j,i) = vitd3_exposure_av(j,i) + vitd3_exposure(j,i)
+                   ENDDO
                 ENDDO
-             ENDDO
+             ENDIF