Changeset 354 for palm/trunk/SOURCE/read_3d_binary.f90

Timestamp:

Jul 13, 2009 1:14:25 PM (15 years ago)

Author:

maronga

Message:

xy cross section output of sensible and latent heatflux is now available

File:

• palm/trunk/SOURCE/read_3d_binary.f90 (modified) (2 diffs)

palm/trunk/SOURCE/read_3d_binary.f90

@@ -521 +521 @@
                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
 
+                CASE ( 'qsws_av' )
+                   IF ( .NOT. ALLOCATED( qsws_av ) )  THEN
+                      ALLOCATE( qsws_av(nys-1:nyn+1,nxl-1:nxr+1) )
+                   ENDIF  
+                   IF ( k == 1 )  READ ( 13 )  tmp_2d
+                   qsws_av(nysc-1:nync+1,nxlc-1:nxrc+1) = &
+                                          tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
+
                 CASE ( 'qswst' )
                    IF ( k == 1 )  READ ( 13 )  tmp_2d
@@ -610 +618 @@
                    shf_m(nysc-1:nync+1,nxlc-1:nxrc+1) = &
                                           tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
-
+                CASE ( 'shf_av' )
+                   IF ( .NOT. ALLOCATED( shf_av ) )  THEN
+                      ALLOCATE( shf_av(nys-1:nyn+1,nxl-1:nxr+1) )
+                   ENDIF
+                   IF ( k == 1 )  READ ( 13 )  tmp_2d
+                   shf_av(nysc-1:nync+1,nxlc-1:nxrc+1) = &
+                                          tmp_2d(nysf-1:nynf+1,nxlf-1:nxrf+1)
                 CASE ( 'spectrum_x' )
                    IF ( k == 1 )  THEN