Changeset 1691 for palm/trunk/SOURCE/data_output_2d.f90

Timestamp:

Oct 26, 2015 4:17:44 PM (8 years ago)

Author:

maronga

Message:

various bugfixes and modifications of the atmosphere-land-surface-radiation interaction. Completely re-written routine to calculate surface fluxes (surface_layer_fluxes.f90) that replaces prandtl_fluxes. Minor formatting corrections and renamings

File:

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

palm/trunk/SOURCE/data_output_2d.f90

@@ -14 +14 @@
 ! PALM. If not, see <http://www.gnu.org/licenses/>.
 !
-! Copyright 1997-2014 Leibniz Universitaet Hannover
+! Copyright 1997-2015 Leibniz Universitaet Hannover
 !--------------------------------------------------------------------------------!
 !
 ! Current revisions:
 ! -----------------
-! 
+! Added output of Obukhov length (ol) and radiative heating rates  for RRTMG. 
+! Formatting corrections.
 ! 
 ! Former revisions:
@@ -112 +113 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  dzw, e, nr, p, pt, q, qc, ql, ql_c, ql_v, ql_vp, qr, qsws,      &
+        ONLY:  dzw, e, nr, ol, p, pt, q, qc, ql, ql_c, ql_v, ql_vp, qr, qsws,  &
                rho, sa, shf, tend, ts, u, us, v, vpt, w, z0, z0h, zu, zw
         
@@ -162 +163 @@
     USE radiation_model_mod,                                                   &
         ONLY:  rad_net, rad_net_av, rad_sw_in, rad_sw_in_av, rad_sw_out,       &
-               rad_sw_out_av, rad_lw_in, rad_lw_in_av, rad_lw_out,             &
-               rad_lw_out_av 
+               rad_sw_out_av, rad_sw_cs_hr, rad_sw_cs_hr_av, rad_sw_hr,        &
+               rad_sw_hr_av, rad_lw_in, rad_lw_in_av, rad_lw_out,              &
+               rad_lw_out_av, rad_lw_cs_hr, rad_lw_cs_hr_av, rad_lw_hr,        &
+               rad_lw_hr_av
 
     IMPLICIT NONE
@@ -543 +546 @@
                 IF ( mode == 'xy' )  level_z = zu
 
+             CASE ( 'ol*_xy' )        ! 2d-array
+                IF ( av == 0 ) THEN
+                   DO  i = nxlg, nxrg
+                      DO  j = nysg, nyng
+                         local_pf(i,j,nzb+1) = ol(j,i)
+                      ENDDO
+                   ENDDO
+                ELSE
+                   DO  i = nxlg, nxrg
+                      DO  j = nysg, nyng
+                         local_pf(i,j,nzb+1) = ol_av(j,i)
+                      ENDDO
+                   ENDDO
+                ENDIF
+                resorted = .TRUE.
+                two_d = .TRUE.
+                level_z(nzb+1) = zu(nzb+1)
+
              CASE ( 'p_xy', 'p_xz', 'p_yz' )
                 IF ( av == 0 )  THEN
@@ -937 +958 @@
                 ENDIF
 
+             CASE ( 'rad_lw_cs_hr_xy', 'rad_lw_cs_hr_xz', 'rad_lw_cs_hr_yz' )
+                IF ( av == 0 )  THEN
+                   to_be_resorted => rad_lw_cs_hr
+                ELSE
+                   to_be_resorted => rad_lw_cs_hr_av
+                ENDIF
+
+             CASE ( 'rad_lw_hr_xy', 'rad_lw_hr_xz', 'rad_lw_hr_yz' )
+                IF ( av == 0 )  THEN
+                   to_be_resorted => rad_lw_hr
+                ELSE
+                   to_be_resorted => rad_lw_hr_av
+                ENDIF
+
              CASE ( 'rad_sw_in_xy', 'rad_sw_in_xz', 'rad_sw_in_yz' )
                 IF ( av == 0 )  THEN
@@ -949 +984 @@
                 ELSE
                    to_be_resorted => rad_sw_out_av
+                ENDIF
+
+             CASE ( 'rad_sw_cs_hr_xy', 'rad_sw_cs_hr_xz', 'rad_sw_cs_hr_yz' )
+                IF ( av == 0 )  THEN
+                   to_be_resorted => rad_sw_cs_hr
+                ELSE
+                   to_be_resorted => rad_sw_cs_hr_av
+                ENDIF
+
+             CASE ( 'rad_sw_hr_xy', 'rad_sw_hr_xz', 'rad_sw_hr_yz' )
+                IF ( av == 0 )  THEN
+                   to_be_resorted => rad_sw_hr
+                ELSE
+                   to_be_resorted => rad_sw_hr_av
                 ENDIF
 
@@ -1227 +1276 @@
 !--                Exit the loop for layers beyond the data output domain
 !--                (used for soil model)
-                   IF ( layer_xy .GT. nzt_do )  THEN 
+                   IF ( layer_xy > nzt_do )  THEN 
                       EXIT loop1
                    ENDIF