Changeset 1757 for palm/trunk/SOURCE/land_surface_model.f90

Timestamp:

Feb 22, 2016 3:49:32 PM (8 years ago)

Author:

maronga

Message:

some changes in land surface model, radiation model, nudging and some minor updates

File:

• palm/trunk/SOURCE/land_surface_model.f90 (modified) (8 diffs)

palm/trunk/SOURCE/land_surface_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: set tm_soil_m to zero after allocation. Added parameter 
+! unscheduled_radiation_calls to control calls of the radiation model based on
+! the skin temperature change during one time step (preliminary version). Set 
+! qsws_soil_eb to zero at model start (previously set to qsws_eb). Removed MAX
+! function as it cannot be vectorized.
 ! 
 ! Former revisions:
@@ -132 +136 @@
     USE radiation_model_mod,                                                   &
         ONLY:  force_radiation_call, radiation_scheme, rad_net, rad_sw_in,     &
-               rad_lw_out, rad_lw_out_change_0, sigma_sb
+               rad_lw_out, rad_lw_out_change_0, sigma_sb,                      &
+               unscheduled_radiation_calls
         
 #if defined ( __rrtmg )
@@ -320 +325 @@
 
     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::                                 &
-              tt_soil_m, & !< t_soil storage array 
-              tm_soil_m, & !< m_soil storage array 
+              tt_soil_m, & !< t_soil storage array
+              tm_soil_m, & !< m_soil storage array
               root_fr      !< root fraction (sum=1) 
 
@@ -668 +673 @@
        tt_surface_m = 0.0_wp
        tt_soil_m    = 0.0_wp
+       tm_soil_m    = 0.0_wp
        tm_liq_eb_m  = 0.0_wp
        c_liq        = 0.0_wp
@@ -681 +687 @@
 
        qsws_liq_eb  = 0.0_wp
-       qsws_soil_eb = qsws_eb
+       qsws_soil_eb = 0.0_wp
        qsws_veg_eb  = 0.0_wp
 
@@ -1235 +1241 @@
 
 !
-!--          In case of fast changes in the skin temperature, it is required to
-!--          update the radiative fluxes in order to keep the solution stable
-             IF ( ABS( t_surface_p(j,i) - t_surface(j,i) ) > 0.2_wp )  THEN
+!--          In case of fast changes in the skin temperature, it is possible to
+!--          update the radiative fluxes independently from the prescribed
+!--          radiation call frequency. This effectively prevents oscillations,
+!--          especially when setting skip_time_do_radiation /= 0. The threshold
+!--          value of 0.2 used here is just a first guess. This method should be
+!--          revised in the future as tests have shown that the threshold is 
+!--          often reached, when no oscillations would occur (causes immense
+!--          computing time for the radiation code).
+             IF ( ABS( t_surface_p(j,i) - t_surface(j,i) ) > 0.2_wp .AND.      &
+                  unscheduled_radiation_calls )  THEN
                 force_radiation_call_l = .TRUE.
              ENDIF
@@ -1388 +1401 @@
 !
 !--    Make a logical OR for all processes. Force radiation call if at
-!--    least one processor
-       IF ( intermediate_timestep_count == intermediate_timestep_count_max-1 ) &
-       THEN
+!--    least one processor reached the threshold change in skin temperature
+       IF ( unscheduled_radiation_calls .AND. intermediate_timestep_count      &
+            == intermediate_timestep_count_max-1 )  THEN
 #if defined( __parallel )
           IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
@@ -1628 +1641 @@
 !
 !--             Account for dry soils (find a better solution here!)
-                m_soil_p(:,j,i) = MAX(m_soil_p(:,j,i),0.0_wp)
+                DO  k = nzb_soil, nzt_soil
+                   IF ( m_soil_p(k,j,i) < 0.0_wp )  m_soil_p(k,j,i) = 0.0_wp
+                ENDDO
+
+
+
 
 !