Changeset 3479 for palm/trunk

Timestamp:

Nov 1, 2018 4:00:30 PM (6 years ago)

Author:

gronemeier

Message:

reworked check for output quantities; assign new palm-error numbers, set unit according to data standard

File:

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

palm/trunk/SOURCE/biometeorology_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! - reworked check for output quantities
+! - assign new palm-error numbers
+! - set unit according to data standard.
+! 
+! 3475 2018-10-30 21:16:31Z kanani
 ! Add option for using MRT from RTM instead of simplified global value
 ! 
@@ -482 +487 @@
  SUBROUTINE biom_check_data_output( var, unit )
 
-       USE control_parameters,                                                 &
-           ONLY: data_output, message_string
-
-       IMPLICIT NONE
-
-       CHARACTER (LEN=*) ::  unit     !< The unit for the variable var
-       CHARACTER (LEN=*) ::  var      !< The variable in question
-
-
-                unit = '°C'
-                IF ( .NOT. biometeorology ) THEN
-                  message_string = 'output of "' // TRIM( var ) // '" req' //  &
-                        'uires biometeorology = .TRUE. in initialisati' &
-                        // 'on_parameters'
-                  CALL message( 'check_parameters', 'PA0114', 1, 2, 0, 6, 0 )
-                  unit = 'illegal'
-                ENDIF
-                IF ( .NOT.  radiation ) THEN
-                   message_string = 'output of "' // TRIM( var ) // '" require'&
-                                    // 's radiation = .TRUE.'
-                   CALL message( 'check_parameters', 'PA0509', 1, 2, 0, 6, 0 )
-                   unit = 'illegal'
-                ENDIF
-                IF ( mrt_nlevels == 0 ) THEN
-                   message_string = 'output of "' // TRIM( var ) // '" require'&
-                                    // 's mrt_nlevels > 0'
-                   CALL message( 'check_parameters', 'PA0510', 1, 2, 0, 6, 0 )
-                   unit = 'illegal'
-                ENDIF
+    USE control_parameters,                                                    &
+        ONLY: data_output, message_string
+
+    IMPLICIT NONE
+
+    CHARACTER (LEN=*) ::  unit     !< The unit for the variable var
+    CHARACTER (LEN=*) ::  var      !< The variable in question
+
+
+    SELECT CASE ( TRIM( var ) )
+     
+       CASE( 'biom_tmrt', 'biom_mrt', 'biom_pet', 'biom_pt', 'biom_utci' )
+          unit = 'degree_C'
+
+       CASE DEFAULT
+          unit = 'illegal'
+
+    END SELECT
+
+    IF ( unit /= 'illegal' )  THEN
+!
+!--    Futher checks if output belongs to biometeorology
+       IF ( .NOT. biometeorology .OR. .NOT. radiation ) THEN
+         message_string = 'output of "' // TRIM( var ) // '" req' //  &
+               'uires biometeorology = .TRUE. and radiation = .TRUE. ' &
+               // 'in initialisation_parameters'
+         CALL message( 'biom_check_data_output', 'PA0561', 1, 2, 0, 6, 0 )
+       ENDIF
+       IF ( mrt_nlevels == 0 ) THEN
+          message_string = 'output of "' // TRIM( var ) // '" require'&
+                           // 's mrt_nlevels > 0'
+          CALL message( 'biom_check_data_output', 'PA0562', 1, 2, 0, 6, 0 )
+       ENDIF
+
+    ENDIF
 
  END SUBROUTINE biom_check_data_output