Changeset 3620

Timestamp:

Dec 11, 2018 12:29:43 PM (6 years ago)

Author:

moh.hefny

Message:

fix date and time and update 3d rad_lw_out

Location:

palm/trunk/SOURCE

Files:

• land_surface_model_mod.f90 (modified) (3 diffs)
• radiation_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! update the 3d rad_lw_out array
+!
+! 3597 2018-12-04 08:40:18Z maronga 
 ! Added pt_2m / theta_2m. Removed unncessary _eb strings.
 ! 
@@ -1755 +1758 @@
 
     USE pegrid
+    USE radiation_model_mod,  ONLY:  rad_lw_out
 
     IMPLICIT NONE
@@ -2191 +2195 @@
 
        surf%shf(m) = - f_shf * ( surf%pt1(m) - surf%pt_surface(m) ) / c_p
+       ! update the 3d field of rad_lw_ou array
+       IF ( horizontal ) rad_lw_out(k+k_off,j+j_off,i+i_off) = surf%rad_lw_out(m)
 
        IF ( humidity )  THEN

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! - fix initialization of date and time after calling zenith
+! - fix a bug in radiation_solar_pos
+!
+! 3616 2018-12-10 09:44:36Z Salim
 ! fix manipulation of time variables in radiation_presimulate_solar_pos
 !
@@ -2441 +2445 @@
 !--       whether the sun is up
           CALL calc_zenith     
+          ! readjust date and time to its initial value
+          CALL init_date_and_time
 !
 !--       Calculate initial surface albedo for different surfaces
@@ -2634 +2640 @@
 
        END SELECT
+
+! readjust date and time to its initial value
+       CALL init_date_and_time
 
        RETURN
@@ -8012 +8021 @@
       DO  it = 0, CEILING(( end_time - spinup_time ) / dt_radiation)
          time_since_reference_point = REAL(it, wp) * dt_radiation
-         simulated_time = simulated_time + dt_spinup
+         simulated_time = simulated_time + dt_radiation
          CALL simulate_pos
       ENDDO