Changeset 2007 for palm/trunk/SOURCE/netcdf_interface_mod.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/netcdf_interface_mod.f90 (modified) (12 diffs)

palm/trunk/SOURCE/netcdf_interface_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Added support for new urban surface model (temporary modifications of 
+! SELECT CASE ( ) necessary, see variable trimvar),
+! increased DIMENSION of do2d_unit, do3d_unit, id_var_do2d, id_var_do3d,
+! increased LEN of char_cross_profiles, var_list, var_list_old
 ! 
 ! Former revisions:
@@ -240 +243 @@
     CHARACTER (LEN=9), DIMENSION(300) ::  dopr_unit = 'unknown'
 
-    CHARACTER (LEN=7), DIMENSION(0:1,100) ::  do2d_unit, do3d_unit
+    CHARACTER (LEN=7), DIMENSION(0:1,500) ::  do2d_unit, do3d_unit
 
     CHARACTER (LEN=16), DIMENSION(25) ::  prt_var_names = &
@@ -319 +322 @@
 
     INTEGER(iwp), DIMENSION(dopts_num,0:10) ::  id_var_dopts
-    INTEGER(iwp), DIMENSION(0:1,100)        ::  id_var_do2d, id_var_do3d
+    INTEGER(iwp), DIMENSION(0:1,500)        ::  id_var_do2d, id_var_do3d
     INTEGER(iwp), DIMENSION(100,0:9)        ::  id_dim_z_pr, id_var_dopr, &
                                                 id_var_z_pr
@@ -441 +444 @@
         ONLY:  hom, statistic_regions
 
+    USE urban_surface_mod,                                                     &
+        ONLY:  urban_surface, usm_define_netcdf_grid
+
 
     IMPLICIT NONE
@@ -454 +460 @@
     CHARACTER (LEN=10)             ::  precision             !<
     CHARACTER (LEN=10)             ::  var                   !<
+    CHARACTER (LEN=20)             ::  trimvar               !< TRIM of output-variable string
     CHARACTER (LEN=80)             ::  time_average_text     !<
-    CHARACTER (LEN=2000)           ::  char_cross_profiles   !<
-    CHARACTER (LEN=2000)           ::  var_list              !<
-    CHARACTER (LEN=2000)           ::  var_list_old          !<
+    CHARACTER (LEN=4000)           ::  char_cross_profiles   !<
+    CHARACTER (LEN=4000)           ::  var_list              !<
+    CHARACTER (LEN=4000)           ::  var_list_old          !<
 
     CHARACTER (LEN=100), DIMENSION(1:crmax) ::  cross_profiles_adj   !<
@@ -506 +513 @@
     REAL(wp), DIMENSION(:), ALLOCATABLE   ::  netcdf_data    !<
     REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  netcdf_data_2d !<
+
 
 !
@@ -746 +754 @@
 
           DO WHILE ( domask(mid,av,i)(1:1) /= ' ' )
-
+!
+!--          Temporary solution to account for data output within the new urban 
+!--          surface model (urban_surface_mod.f90), see also SELECT CASE ( trimvar )
+             trimvar = TRIM( domask(mid,av,i) )
+             IF ( urban_surface  .AND.  trimvar(1:3) == 'us_' )  THEN
+                trimvar = 'usm_output'
+             ENDIF
 !
 !--          Check for the grid
              found = .FALSE.
-             SELECT CASE ( domask(mid,av,i) )
+             SELECT CASE ( trimvar )
 !
 !--             Most variables are defined on the scalar grid
@@ -782 +796 @@
                    grid_z = 'zw'
 
+!             
+!--       Block of urban surface model outputs   
+                CASE ( 'usm_output' )
+
+                   CALL usm_define_netcdf_grid( domask(mid,av,i), found, &
+                                                        grid_x, grid_y, grid_z )
 
                 CASE DEFAULT
@@ -806 +826 @@
                                                     grid_x, grid_y, grid_z )
                    ENDIF
-
+                                                 
                    IF ( .NOT. found )  THEN
                       WRITE ( message_string, * ) 'no grid defined for',       &
@@ -1256 +1276 @@
 
           DO WHILE ( do3d(av,i)(1:1) /= ' ' )
-
+!
+!--          Temporary solution to account for data output within the new urban 
+!--          surface model (urban_surface_mod.f90), see also SELECT CASE ( trimvar )
+             trimvar = TRIM( do3d(av,i) )
+             IF ( urban_surface  .AND.  trimvar(1:3) == 'us_' )  THEN
+                trimvar = 'usm_output'
+             ENDIF
 !
 !--          Check for the grid
              found = .FALSE.
-             SELECT CASE ( do3d(av,i) )
+             SELECT CASE ( trimvar )
 !
 !--             Most variables are defined on the scalar grid
@@ -1292 +1318 @@
                    grid_z = 'zw'
 
+!             
+!--             Block of urban surface model outputs   
+                CASE ( 'usm_output' )
+                   CALL usm_define_netcdf_grid( do3d(av,i), found, &
+                                                   grid_x, grid_y, grid_z )
 
                 CASE DEFAULT
@@ -1315 +1346 @@
                                                     grid_y, grid_z )
                    ENDIF
-
+                                                 
                    IF ( .NOT. found )  THEN
                       WRITE ( message_string, * ) 'no grid defined for varia', &