Changeset 983 for palm/trunk/SOURCE/netcdf.f90

Timestamp:

Aug 21, 2012 2:17:57 PM (12 years ago)

Author:

hoffmann

Message:

Bugfix in netcdf.f90 and improvements in palmplot.

File:

• palm/trunk/SOURCE/netcdf.f90 (modified) (6 diffs)

palm/trunk/SOURCE/netcdf.f90

@@ -7 +7 @@
 ! Current revisions:
 ! ------------------
-!
+! Bugfix in cross_profiles.
 !
 ! Former revisions:
@@ -3166 +3166 @@
           cross_profiles_adj  = ADJUSTL( cross_profiles )
           cross_profiles_numb = 999999
+          cross_profiles_char = ''
 
 !
@@ -3223 +3224 @@
           ENDIF
 
+          DO i = 1, crmax
+             IF ( cross_profiles_numb(i) == 999999 ) THEN
+                DO j = i + 1, crmax
+                   IF ( cross_profiles_numb(j) /= 999999 ) THEN
+                      cross_profiles_char(i) = cross_profiles_char(j)
+                      cross_profiles_numb(i) = cross_profiles_numb(j)
+                      cross_profiles_numb(j) = 999999
+                      EXIT
+                   ENDIF
+                ENDDO
+             ENDIF
+          ENDDO
+
+          DO i = 1, crmax-1
+             IF ( cross_profiles_numb(i + 1) == 999999 ) THEN
+                cross_profiles_count = i
+                EXIT
+             ENDIF
+          ENDDO
 !
 !--       Check if all profiles defined in data_output_pr are defined in 
@@ -3237 +3257 @@
                 ENDIF
 
-                IF ( j == cross_profiles_count )  THEN
+                IF (( j == cross_profiles_count ) .AND.                        &
+                    ( cross_profiles_count <= crmax - 1))  THEN
                    cross_profiles_count = cross_profiles_count + 1
                    cross_profiles_maxi  = cross_profiles_maxi  + 1
@@ -3260 +3281 @@
           IF ( cross_profiles_count >= crmax )  THEN
              message_string = 'It is not allowed to arrange more than '&
-                              // '100 profiles with cross_profiles.'
+                              // '100 profiles with cross_profiles. Apart '&
+                              // 'from that, all profiles are saved to the '&
+                              // 'netCDF file.'
              CALL message( 'define_netcdf_header', 'PA0354', 0, 0, 0, 6, 0 )
           ENDIF
@@ -3270 +3293 @@
           char_cross_profiles = ';'
           id_last = 1
+          cross_profiles_count = MIN( cross_profiles_count, crmax )
 
           DO  i = 1, cross_profiles_count