Changeset 2007 for palm/trunk/SOURCE/data_output_3d.f90

Timestamp:

Aug 24, 2016 3:47:17 PM (8 years ago)

Author:

kanani

Message:

changes in the course of urban surface model implementation

File:

• palm/trunk/SOURCE/data_output_3d.f90 (modified) (6 diffs)

palm/trunk/SOURCE/data_output_3d.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Added support for new urban surface model (temporary modifications of 
+! SELECT CASE ( ) necessary, see variable trimvar)
 ! 
 ! Former revisions:
@@ -178 +179 @@
         ONLY:  radiation, radiation_data_output_3d
 
+    USE urban_surface_mod,                                                     &
+        ONLY:  nzub, nzut, urban_surface, usm_data_output_3d
+
 
     IMPLICIT NONE
@@ -200 +204 @@
 
     REAL(wp), DIMENSION(:,:,:), POINTER ::  to_be_resorted  !< 
+
+    CHARACTER (LEN=20) ::  trimvar  !< TRIM of output-variable string
 
 !
@@ -260 +266 @@
 
     DO  WHILE ( do3d(av,if)(1:1) /= ' ' )
-!
+
+!
+!--    Temporary solution to account for data output within the new urban 
+!--    surface model (urban_surface_mod.f90), see also SELECT CASE ( trimvar ).
 !--    Store the array chosen on the temporary array.
-       resorted = .FALSE.
-       nzb_do = nzb
-       nzt_do = nz_do3d
+       trimvar = TRIM( do3d(av,if) )
+       IF ( urban_surface  .AND.  trimvar(1:3) == 'us_' )  THEN
+          trimvar = 'usm_output'
+          resorted = .TRUE.
+          nzb_do   = nzub
+          nzt_do   = nzut
+       ELSE
+          resorted = .FALSE.
+          nzb_do   = nzb
+          nzt_do   = nz_do3d
+       ENDIF
 !
 !--    Set flag to steer output of radiation, land-surface, or user-defined
@@ -273 +290 @@
        ALLOCATE( local_pf(nxlg:nxrg,nysg:nyng,nzb_do:nzt_do) )
 
-       SELECT CASE ( TRIM( do3d(av,if) ) )
+       SELECT CASE ( trimvar )
 
           CASE ( 'e' )
@@ -553 +570 @@
                 to_be_resorted => w_av
              ENDIF
+!             
+!--       Block of urban surface model outputs   
+          CASE ( 'usm_output' )
+             CALL usm_data_output_3d( av, do3d(av,if), found, local_pf,     &
+                                         nzb_do, nzt_do )
 
           CASE DEFAULT