Changeset 3914 for palm/trunk/SOURCE

Timestamp:

Apr 17, 2019 4:02:02 PM (6 years ago)

Author:

suehring

Message:

In order to obtain correct surface temperature during spinup set window fraction to zero (only during spinup) instead of just disabling time-integration of window-surface temperature.

File:

• palm/trunk/SOURCE/urban_surface_mod.f90 (modified) (7 diffs)

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! In order to obtain correct surface temperature during spinup set window 
+! fraction to zero (only during spinup) instead of just disabling
+! time-integration of window-surface temperature. 
+! 
+! 3901 2019-04-16 16:17:02Z suehring
 ! Workaround - set green fraction to zero ( green-heat model crashes ).
 ! 
@@ -3753 +3758 @@
             lsf_surf = .FALSE.
         ENDIF
-
 !
 !--     Flag surface elements belonging to the ground floor level. Therefore, 
@@ -7636 +7640 @@
         LOGICAL                               :: spinup             !true during spinup
         
+        REAL(wp)                              :: frac_win           !< window fraction, used to restore original values during spinup
+        REAL(wp)                              :: frac_green         !< green fraction, used to restore original values during spinup
+        REAL(wp)                              :: frac_wall          !< wall fraction, used to restore original values during spinup
         REAL(wp)                              :: stend_wall         !< surface tendency
         
@@ -7701 +7708 @@
         DO  m = 1, surf_usm_h%ns
 !
+!--       During spinup set green and window fraction to zero and restore
+!--       at the end of the loop.
+!--       Note, this is a temporary fix and need to be removed later.  
+           IF ( spinup )  THEN
+              frac_win   = surf_usm_h%frac(ind_wat_win,m)
+              frac_wall  = surf_usm_h%frac(ind_veg_wall,m)
+              frac_green = surf_usm_h%frac(ind_pav_green,m)
+              surf_usm_h%frac(ind_wat_win,m)   = 0.0_wp
+              surf_usm_h%frac(ind_veg_wall,m)  = 1.0_wp
+              surf_usm_h%frac(ind_pav_green,m) = 0.0_wp
+           ENDIF
+!
 !--        Get indices of respective grid point
            i = surf_usm_h%i(m)
@@ -8141 +8160 @@
               surf_usm_h%r_s(m) = 1.0E10_wp
            ENDIF
+!
+!--        During spinup green and window fraction are set to zero. Here, the original
+!--        values are restored.
+           IF ( spinup )  THEN
+              surf_usm_h%frac(ind_wat_win,m)   = frac_win
+              surf_usm_h%frac(ind_veg_wall,m)  = frac_wall
+              surf_usm_h%frac(ind_pav_green,m) = frac_green
+           ENDIF
  
        ENDDO
@@ -8148 +8175 @@
        DO  l = 0, 3
            DO  m = 1, surf_usm_v(l)%ns
+!
+!--           During spinup set green and window fraction to zero and restore
+!--           at the end of the loop.
+!--           Note, this is a temporary fix and need to be removed later. 
+              IF ( spinup )  THEN
+                 frac_win   = surf_usm_v(l)%frac(ind_wat_win,m)
+                 frac_wall  = surf_usm_v(l)%frac(ind_veg_wall,m)
+                 frac_green = surf_usm_v(l)%frac(ind_pav_green,m)
+                 surf_usm_v(l)%frac(ind_wat_win,m)   = 0.0_wp
+                 surf_usm_v(l)%frac(ind_veg_wall,m)  = 1.0_wp
+                 surf_usm_v(l)%frac(ind_pav_green,m) = 0.0_wp
+              ENDIF
 !
 !--          Get indices of respective grid point
@@ -8517 +8556 @@
               ELSE
                  surf_usm_v(l)%r_s(m) = 1.0E10_wp
+              ENDIF
+!
+!--           During spinup green and window fraction are set to zero. Here, the original
+!--           values are restored.
+              IF ( spinup )  THEN
+                 surf_usm_v(l)%frac(ind_wat_win,m)   = frac_win
+                 surf_usm_v(l)%frac(ind_veg_wall,m)  = frac_wall
+                 surf_usm_v(l)%frac(ind_pav_green,m) = frac_green
               ENDIF