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

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added checks for the urban surface model,
+! increased counter in DO WHILE loop over data_output (for urban surface output)
 ! 
 ! Former revisions:
@@ -423 +424 @@
     USE profil_parameter
     USE radiation_model_mod,                                                   &
-        ONLY: radiation, radiation_check_data_output,                          &
-              radiation_check_data_output_pr, radiation_check_parameters
+        ONLY:  radiation, radiation_check_data_output,                         &
+               radiation_check_data_output_pr, radiation_check_parameters
     USE spectra_mod,                                                           &
         ONLY:  calculate_spectra, spectra_check_parameters
@@ -431 +432 @@
     USE statistics
     USE transpose_indices
+    USE urban_surface_mod,                                                     &
+        ONLY:  urban_surface, usm_check_data_output, usm_check_parameters
     USE wind_turbine_model_mod,                                                &
-        ONLY: wtm_check_parameters, wind_turbine
+        ONLY:  wtm_check_parameters, wind_turbine
 
 
@@ -1100 +1103 @@
 !-- When land surface model is used, perform additional checks
     IF ( land_surface )  CALL lsm_check_parameters
+
+!
+!-- When urban surface model is used, perform additional checks
+    IF ( urban_surface )  CALL usm_check_parameters
 
 !
@@ -1616 +1623 @@
     ENDIF
 
+!
+!   This IF clause needs revision, got too complex!!
     IF ( surface_heatflux == 9999999.9_wp  )  THEN
        constant_heatflux = .FALSE.
-       IF ( large_scale_forcing  .OR.  land_surface )  THEN
+       IF ( large_scale_forcing  .OR.  land_surface  .OR.  urban_surface )  THEN
           IF ( ibc_pt_b == 0 )  THEN
              constant_heatflux = .FALSE.
@@ -2834 +2843 @@
 !-- Check and set steering parameters for 2d/3d data output and averaging
     i   = 1
-    DO  WHILE ( data_output(i) /= ' '  .AND.  i <= 100 )
+    DO  WHILE ( data_output(i) /= ' '  .AND.  i <= 500 )
 !
 !--    Check for data averaging
@@ -3053 +3062 @@
              IF ( unit == 'illegal' )  THEN
                 CALL radiation_check_data_output( var, unit, i, ilen, k )
+             ENDIF
+
+!
+!--          Block of urban surface model outputs
+             IF ( unit == 'illegal' .AND. urban_surface .AND. var(1:3) == 'us_' ) THEN
+                 CALL usm_check_data_output( var, unit )
              ENDIF
 
@@ -3718 +3733 @@
     ENDIF
 
-    IF ( large_scale_forcing  .AND.  topography /= 'flat' )  THEN
+    IF ( large_scale_forcing  .AND.  topography /= 'flat'                      &
+                              .AND.  .NOT.  urban_surface )  THEN
        message_string = 'The usage of large scale forcing from external &'//   & 
                         'file LSF_DATA is not implemented for non-flat topography'