Changeset 210 for palm/trunk/UTIL

Timestamp:

Nov 6, 2008 8:54:02 AM (16 years ago)

Author:

raasch

Message:

updates in dvr routines for new dvr version

File:

• palm/trunk/UTIL/combine_plot_fields.f90 (modified) (6 diffs)

palm/trunk/UTIL/combine_plot_fields.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Size of pf3d adjusted to the required output size (1 gridpoint less, along
+! all three dimensions), because output of a subset of the data
+! (pf3d(nxa:nxe...) in the NF90_PUT_VAR statement caused segmentation fault
+! with the INTEL compiler.
 ! Subdomain data are read into temporary arrays pf_tmp/pf3d_tmp in order to
 ! avoid INTEL compiler warnings about (automatic) creation of temporary arrays
@@ -93 +97 @@
 
     INTEGER ::  av, danz, i, id,             &
-                j, model, models, nc_stat, nxa, nxag, nxe, nxeg, nya,   &
+                j, k, model, models, nc_stat, nxa, nxag, nxe, nxeg, nya,   &
                 nyag, nye, nyeg, nza, nzag, nze, nzeg, pos, time_step, xa, xe, &
-                ya, ye, za, ze
+                xxa, xxe, ya, ye, yya, yye, za, ze, zza, zze
 
     INTEGER, DIMENSION(0:1) ::  current_level, current_var, fanz, id_set, &
@@ -603 +607 @@
 !
 !--             Get variable ID from name
-!                print*, '*** find id for "',TRIM( var_name ),'" begin'
                 nc_stat = NF90_INQ_VARID( id_set(av), TRIM( var_name ), &
                      id_var(av,i) )
                 IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 16 )
-!                print*, '*** find id for "',TRIM( var_name ),'" end'
 
              ENDDO
@@ -636 +638 @@
                 READ ( id+110 )  nxag, nxeg, nyag, nyeg, nzag, nzeg
                 READ ( id+110 )  nxa, nxe, nya, nye, nza, nze
-                ALLOCATE ( pf3d(nxag:nxeg,nyag:nyeg,nzag:nzeg) )
+                ALLOCATE ( pf3d(nxa:nxe,nya:nye,nza:nze) )
                 IF ( avs_output )  THEN
                    OPEN ( 2, FILE='PLOT3D_DATA'//TRIM( model_string ), &
@@ -661 +663 @@
              ALLOCATE( pf3d_tmp(xa:xe,ya:ye,za:ze) )
              READ ( id+110 )  pf3d_tmp
-             pf3d(xa:xe,ya:ye,za:ze) = pf3d_tmp
+
+             xxa = MAX( nxa, xa )
+             xxe = MIN( nxe, xe )
+             yya = MAX( nya, ya )
+             yye = MIN( nye, ye )
+             zza = MAX( nza, za )
+             zze = MIN( nze, ze )
+             DO  k = zza, zze
+                DO  j = yya, yye
+                   DO  i = xxa, xxe
+                      pf3d(i,j,k) = pf3d_tmp(i,j,k)
+                   ENDDO
+                ENDDO
+             ENDDO
+
              DEALLOCATE( pf3d_tmp )
              IF ( id == 0 )  fanz(av) = fanz(av) + 1
@@ -688 +704 @@
 !--          Now write the data
              nc_stat = NF90_PUT_VAR( id_set(av), id_var(av,current_var(av)), &
-                                  pf3d(nxa:nxe,nya:nye,nza:nze),      &
-                                  start = (/ 1, 1, 1, time_step /),   &
+                                     pf3d, start = (/ 1, 1, 1, time_step /), &
                               count = (/ nxe-nxa+1, nye-nya+1, nze-nza+1, 1 /) )
              IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 18 )