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

Timestamp:

Mar 4, 2015 5:44:27 PM (9 years ago)

Author:

maronga

Message:

LSM output of r_a and r_s added

File:

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

palm/trunk/SOURCE/data_output_2d.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added output of r_a and r_s
 ! 
 ! Former revisions:
@@ -139 +139 @@
                m_soil_av, nzb_soil, nzt_soil, qsws_eb, qsws_eb_av,             &
                qsws_liq_eb, qsws_liq_eb_av, qsws_soil_eb, qsws_soil_eb_av,     &
-               qsws_veg_eb, qsws_veg_eb_av, shf_eb, shf_eb_av, t_soil,         &
-               t_soil_av, zs
+               qsws_veg_eb, qsws_veg_eb_av, r_a, r_a_av, r_s, r_s_av, shf_eb,  &
+               shf_eb_av, t_soil, t_soil_av, zs
     
     USE netcdf_control
@@ -934 +934 @@
                    to_be_resorted => rho_av
                 ENDIF
+
+             CASE ( 'r_a*_xy' )        ! 2d-array
+                IF ( av == 0 )  THEN
+                   DO  i = nxlg, nxrg
+                      DO  j = nysg, nyng
+                         local_pf(i,j,nzb+1) = r_a(j,i)
+                      ENDDO
+                   ENDDO
+                ELSE
+                   DO  i = nxlg, nxrg
+                      DO  j = nysg, nyng
+                         local_pf(i,j,nzb+1) = r_a_av(j,i)
+                      ENDDO
+                   ENDDO
+                ENDIF
+                resorted = .TRUE.
+                two_d = .TRUE.
+                level_z(nzb+1) = zu(nzb+1)
+
+             CASE ( 'r_s*_xy' )        ! 2d-array
+                IF ( av == 0 )  THEN
+                   DO  i = nxlg, nxrg
+                      DO  j = nysg, nyng
+                         local_pf(i,j,nzb+1) = r_s(j,i)
+                      ENDDO
+                   ENDDO
+                ELSE
+                   DO  i = nxlg, nxrg
+                      DO  j = nysg, nyng
+                         local_pf(i,j,nzb+1) = r_s_av(j,i)
+                      ENDDO
+                   ENDDO
+                ENDIF
+                resorted = .TRUE.
+                two_d = .TRUE.
+                level_z(nzb+1) = zu(nzb+1)
 
              CASE ( 's_xy', 's_xz', 's_yz' )