Changeset 964 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Jul 26, 2012 9:14:24 AM (12 years ago)

Author:

raasch

Message:

old profil-parameters (cross_xtext, cross_normalized_x, etc. ) and respective code removed
(check_open, check_parameters, close_file, data_output_profiles, data_output_spectra, header, modules, parin)

reformatting (netcdf)

append feature removed from unit 14 (check_open)

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (4 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! check of old profil-parameters removed
 !
 ! Former revisions:
@@ -2518 +2518 @@
        END SELECT
 
-!
-!--    Check to which of the predefined coordinate systems the profile belongs
-       DO  k = 1, crmax
-          IF ( INDEX( cross_profiles(k), ' '//TRIM( data_output_pr(i) )//' ' ) &
-               /=0 ) &
-          THEN
-             dopr_crossindex(i) = k
-             EXIT
-          ENDIF
-       ENDDO
-!
-!--    Generate the text for the labels of the PROFIL output file. "-characters
-!--    must be substituted, otherwise PROFIL would interpret them as TeX
-!--    control characters
-       dopr_label(i) = data_output_pr(i)
-       position = INDEX( dopr_label(i) , '"' )
-       DO WHILE ( position /= 0 )
-          dopr_label(i)(position:position) = ''''
-          position = INDEX( dopr_label(i) , '"' )
-       ENDDO
-
     ENDDO
-
-!
-!-- y-value range of the coordinate system (PROFIL).
-!-- x-value range determined in plot_1d.
-    IF ( .NOT. ocean )  THEN
-       cross_uymin = 0.0
-       IF ( z_max_do1d == -1.0 )  THEN
-          cross_uymax = zu(nzt+1)
-       ELSEIF ( z_max_do1d < zu(nzb+1)  .OR.  z_max_do1d > zu(nzt+1) )  THEN
-          WRITE( message_string, * )  'z_max_do1d = ', z_max_do1d, ' must ', &
-                 'be >= ', zu(nzb+1), ' or <= ', zu(nzt+1)
-          CALL message( 'check_parameters', 'PA0099', 1, 2, 0, 6, 0 )
-       ELSE
-          cross_uymax = z_max_do1d
-       ENDIF
-    ENDIF
-
-!
-!-- Check whether the chosen normalizing factor for the coordinate systems is 
-!-- permissible
-    DO  i = 1, crmax
-       SELECT CASE ( TRIM( cross_normalized_x(i) ) )  ! TRIM required on IBM
-
-          CASE ( '', 'wpt0', 'ws2', 'tsw2', 'ws3', 'ws2tsw', 'wstsw2' )
-             j = 0
-
-          CASE DEFAULT
-             message_string = 'unknown normalization method cross_normali' // &
-                              'zed_x = "' // TRIM( cross_normalized_x(i) ) // &
-                              '"'
-             CALL message( 'check_parameters', 'PA0100', 1, 2, 0, 6, 0 )
-
-       END SELECT
-       SELECT CASE ( TRIM( cross_normalized_y(i) ) )  ! TRIM required on IBM
-
-          CASE ( '', 'z_i' )
-             j = 0
-
-          CASE DEFAULT
-             message_string = 'unknown normalization method cross_normali' // &
-                              'zed_y = "' // TRIM( cross_normalized_y(i) ) // &
-                              '"'
-             CALL message( 'check_parameters', 'PA0101', 1, 2, 0, 6, 0 )
-
-       END SELECT
-    ENDDO
-!
-!-- Check normalized y-value range of the coordinate system (PROFIL)
-    IF ( z_max_do1d_normalized /= -1.0  .AND.  z_max_do1d_normalized <= 0.0 ) &
-    THEN
-       WRITE( message_string, * )  'z_max_do1d_normalized = ', &
-                                   z_max_do1d_normalized, ' must be >= 0.0'
-       CALL message( 'check_parameters', 'PA0101', 1, 2, 0, 6, 0 )
-    ENDIF
 
 
@@ -2864 +2789 @@
 
 !
-!-- Upper plot limit (grid point value) for 1D profiles
-    IF ( z_max_do1d == -1.0 )  THEN
-
-       nz_do1d = nzt+1
-
-    ELSE
-       DO  k = nzb+1, nzt+1
-          nz_do1d = k
-          IF ( zw(k) > z_max_do1d )  EXIT
-       ENDDO
-    ENDIF
-
-!
 !-- Upper plot limit for 2D vertical sections
     IF ( z_max_do2d == -1.0 )  z_max_do2d = zu(nzt)
@@ -2958 +2870 @@
              CASE ( 'iso2d' )
                 iso2d_output  = .TRUE.
-             CASE ( 'profil' )
-                profil_output = .TRUE.
              CASE ( 'avs' )
                 avs_output    = .TRUE.