Changeset 964 for palm/trunk/SOURCE/close_file.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/close_file.f90 (modified) (3 diffs)

palm/trunk/SOURCE/close_file.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! old profil-units (40:49) and respective code removed
 !
 ! Former revisions:
@@ -67 +67 @@
     CHARACTER (LEN=2)   ::  suffix
     CHARACTER (LEN=10)  ::  datform = 'lit_endian'
-    CHARACTER (LEN=80)  ::  rtext, title, utext = '', xtext = '', ytext = ''
-
-    INTEGER ::  av, anzzeile, cranz, cross_count, cross_numbers, dimx, dimy, &
-                fid, file_id, j, k, legpos = 1, planz, timodex = 1
-    INTEGER, DIMENSION(100) ::  klist, lstyle, cucol
-
-    LOGICAL ::  checkuf = .TRUE., datleg = .TRUE., dp = .FALSE., &
-                grid = .TRUE., rand = .TRUE., swap, twoxa = .TRUE., &
-                twoya = .TRUE.
-
-    REAL    ::  ansx = -999.999, ansy = -999.999, gwid = 0.1, rlegfak, &
-                sizex, sizey, texfac, utmove = 50.0, uxmax, uxmin, uymax, &
-                uymin, yright
-    REAL, DIMENSION(100) ::  lwid, normx, normy
-
-    NAMELIST /CROSS/   ansx, ansy, cucol, grid, gwid, klist, legpos, lstyle, &
-                       lwid, normx, normy, rand, rlegfak, sizex, sizey, &
-                       texfac, timodex, twoxa, twoya, utext, utmove, uxmax, &
-                       uxmin, uymax, uymin, xtext, ytext
-    NAMELIST /GLOBAL/  checkuf, datform, dimx, dimy, dp, planz, sizex, sizey, &
-                       title, yright
-    NAMELIST /RAHMEN/  anzzeile, cranz, datleg, rtext, swap
+    CHARACTER (LEN=80)  ::  rtext, title
+
+    INTEGER ::  av, anzzeile, dimx, dimy, &
+                fid, file_id, j, k, planz
+
+    LOGICAL ::  checkuf = .TRUE., datleg = .TRUE., dp = .FALSE., swap
+
+    REAL ::  sizex, sizey, yright
+
+    NAMELIST /GLOBAL/  checkuf, datform, dimx, dimy, dp, planz, &
+                       title
+    NAMELIST /RAHMEN/  anzzeile, datleg, rtext, swap
 
 !
@@ -190 +180 @@
                    WRITE ( 32, 3200)  ' ', TRIM( run_description_header ), &
                                       nx+2, ny+2, nz_do3d+1, do3d_avs_n
-                ENDIF
-
-             CASE ( 40:49 )
-!
-!--             Write PROFIL namelist parameters for 1D profiles.
-!--             First determine, how many crosses are to be drawn.
-                IF ( myid == 0 )  THEN
-                   cross_numbers = 0
-                   DO  j = 1, crmax
-                      IF ( cross_profile_number_count(j) /= 0 )  THEN
-                         cross_numbers = cross_numbers + 1
-                      ENDIF
-                   ENDDO
-
-                   IF ( cross_numbers /= 0 )  THEN
-!
-!--                   Determine remaining RAHMEN parameters
-                      swap = .FALSE.
-                      rtext = '\0.5 ' // TRIM( run_description_header ) // &
-                              '    ' // TRIM( region( fid - 40 ) )
-!
-!--                   Write RAHMEN parameters
-                      IF ( statistic_regions == 0  .AND.  fid == 40 )  THEN
-                         suffix = ''
-                      ELSE
-                         WRITE ( suffix, '(''_'',I1)' )  fid - 40
-                      ENDIF
-                      OPEN ( 90, FILE='PLOT1D_PAR' // TRIM( suffix ), &
-                                 FORM='FORMATTED', DELIM='APOSTROPHE' )
-!
-!--                   Subtitle for crosses with time averaging
-                      IF ( averaging_interval_pr /= 0.0 )  THEN
-                         WRITE ( utext, 4000 )  averaging_interval_pr
-                      ENDIF
-!
-!--                   Determine and write CROSS parameters for each individual 
-!--                   cross
-                      cross_count = 0
-                      DO  j = 1, crmax
-                         k = cross_profile_number_count(j)
-                         IF ( k /= 0 )  THEN
-                            cross_count = cross_count + 1
-!
-!--                         Write RAHMEN parameters
-                            IF ( MOD( cross_count-1, &
-                                      profile_rows*profile_columns ) == 0 ) &
-                            THEN
-!
-!--                            Determine number of crosses still to be drawn
-                               cranz = MIN( cross_numbers - cross_count + 1, &
-                                            profile_rows * profile_columns )
-!
-!--                            If the first line cannot be filled with crosses
-!--                            completely, the default number of crosses per 
-!--                            line has to be reduced.
-                               IF ( cranz < profile_columns )  THEN
-                                  anzzeile = cranz
-                               ELSE
-                                  anzzeile = profile_columns
-                               ENDIF
-
-                               WRITE ( 90, RAHMEN )
-
-                            ENDIF
-!
-!--                         Store graph numbers
-                            klist(1:k) = cross_profile_numbers(1:k,j)
-                            klist(k+1:100) = 999999
-!
-!--                         Store graph attributes
-                            cucol  = cross_linecolors(:,j)
-                            lstyle = cross_linestyles(:,j)
-                            lwid = 0.6
-!
-!--                         Sizes, text etc.
-                            sizex = 100.0; sizey = 120.0
-                            rlegfak = 0.7; texfac = 1.0
-!
-!--                         Determine range of x-axis values
-                            IF ( cross_normalized_x(j) == ' ' )  THEN
-!
-!--                            Non-normalized profiles
-                               IF ( cross_uxmin(j) == 0.0  .AND. &
-                                    cross_uxmax(j) == 0.0 )  THEN
-                                  uxmin = cross_uxmin_computed(j)
-                                  uxmax = cross_uxmax_computed(j)
-                                  IF ( uxmin == uxmax )  uxmax = uxmin + 1.0
-                               ELSE
-!
-!--                               Values set in check_parameters are used here
-                                  uxmin = cross_uxmin(j); uxmax = cross_uxmax(j)
-                               ENDIF
-                            ELSE
-!
-!--                            Normalized profiles
-                               IF ( cross_uxmin_normalized(j) == 0.0  .AND. &
-                                    cross_uxmax_normalized(j) == 0.0 )  THEN
-                                  uxmin = cross_uxmin_normalized_computed(j)
-                                  uxmax = cross_uxmax_normalized_computed(j)
-                                  IF ( uxmin == uxmax )  uxmax = uxmin + 1.0
-                               ELSE
-!
-!--                               Values set in check_parameters are used here
-                                  uxmin = cross_uxmin_normalized(j)
-                                  uxmax = cross_uxmax_normalized(j)
-                               ENDIF
-                            ENDIF
-!
-!--                         Range of y-axis values
-!--                         may be re-adjusted during normalization if required
-                            uymin = cross_uymin(j); uymax = cross_uymax(j)
-                            ytext = 'height in m'
-!
-!--                         Normalization of the axes
-                            normx = cross_normx_factor(:,j)
-                            normy = cross_normy_factor(:,j)
-!
-!--                         Labelling of the axes
-                            IF ( cross_normalized_x(j) == ' ' )  THEN
-                               xtext = cross_xtext(j)
-                            ELSE
-                               xtext = TRIM( cross_xtext(j) ) // ' / ' // &
-                                       cross_normalized_x(j)
-                            ENDIF
-                            IF ( cross_normalized_y(j) == ' ' )  THEN
-                               ytext = 'height in m'
-                            ELSE
-                               ytext = 'height in m' // ' / ' // &
-                                       cross_normalized_y(j)
-!
-!--                            Determine upper limit of value range
-                               IF ( z_max_do1d_normalized /= -1.0 )  THEN
-                                  uymax = z_max_do1d_normalized
-                               ENDIF
-                            ENDIF
-
-                            WRITE ( 90, CROSS )
-
-                         ENDIF
-                      ENDDO
-
-                      CLOSE ( 90 )
-                   ENDIF
-                ENDIF
-
-             CASE ( 50:59 )
-!
-!--             Write PROFIL namelist parameters for time series
-!--             first determine number of crosses to be drawn
-                IF ( myid == 0 )  THEN
-                   cranz = 0
-                   DO  j = 1, 12
-                      IF ( cross_ts_number_count(j) /= 0 )  cranz = cranz+1
-                   ENDDO
-
-                   IF ( cranz /= 0 )  THEN
-!
-!--                   Determine RAHMEN parameters
-                      anzzeile = 1
-                      swap = .TRUE.
-                      rtext = '\1.0 ' // TRIM( run_description_header ) // &
-                              '    ' // TRIM( region( fid - 50 ) )
-!
-!--                   Write RAHMEN parameters
-                      IF ( statistic_regions == 0  .AND.  fid == 50 )  THEN
-                         suffix = ''
-                      ELSE
-                         WRITE ( suffix, '(''_'',I1)' )  fid - 50
-                      ENDIF
-                      OPEN ( 90, FILE='PLOTTS_PAR' // TRIM( suffix ), &
-                                 FORM='FORMATTED', DELIM='APOSTROPHE' )
-                      WRITE ( 90, RAHMEN )
-!
-!--                   Determine and write CROSS parameters for each individual 
-!--                   cross
-                      DO  j = 1, 12
-                         k = cross_ts_number_count(j)
-                         IF ( k /= 0 )  THEN
-!
-!--                         Store graph numbers
-                            klist(1:k) = cross_ts_numbers(1:k,j)
-                            klist(k+1:100) = 999999
-!
-!--                         Store graph attributes
-                            cucol(1:k)  = linecolors(1:k)
-                            lstyle(1:k) = linestyles(1:k)
-                            lwid = 0.4
-!
-!--                         Sizes, text etc.
-                            sizex = 250.0; sizey = 40.0
-                            rlegfak = 1.5; texfac = 1.5
-                            xtext = 'time in s'
-                            ytext = ''
-                            utext = ''
-!
-!--                         Determine range of y-axis values
-                            IF ( cross_ts_uymin(j) == 999.999 )  THEN
-                               uymin = cross_ts_uymin_computed(j)
-                            ELSE
-                               uymin = cross_ts_uymin(j)
-                            ENDIF
-                            IF ( cross_ts_uymax(j) == 999.999 )  THEN
-                               uymax = cross_ts_uymax_computed(j)
-                            ELSE
-                               uymax = cross_ts_uymax(j)
-                            ENDIF
-                            IF ( uymin == uymax )  uymax = uymin + 1.0
-!
-!--                         Range of x-axis values
-                            uxmin = 0.0; uxmax = simulated_time
-!
-!--                         Normalizations
-                            normx = 1.0; normy = 1.0
-
-                            WRITE ( 90, CROSS )
-
-                         ENDIF
-                      ENDDO
-
-                      CLOSE ( 90 )
-                   ENDIF
                 ENDIF