Changeset c2bd55e in flexpart.git for src/netcdf_output_mod.f90

Timestamp:

Nov 24, 2017, 1:52:19 PM (6 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

0ecc1fe

Parents:

0c00f1f

Message:

Fixed an issue with nested input fields, and nested output fields that have more than 1000 points in x or y (netcdf output)

File:

• src/netcdf_output_mod.f90 (modified) (2 diffs)

src/netcdf_output_mod.f90

@@ -272 +272 @@
   character(len=10)           :: fprefix
   character(len=3)            :: anspec
-  character                   :: adate*8,atime*6,timeunit*32
-  real, dimension(1000)       :: coord
+  CHARACTER                   :: adate*8,atime*6,timeunit*32
+  ! ESO DBG: WHY IS THIS HARDCODED TO 1000?
+  !REAL, DIMENSION(1000)       :: coord
+  real, allocatable, dimension(:) :: coord
 
   integer                     :: cache_size
@@ -613 +615 @@
   ! longitudes (grid cell centers)
   if (lnest) then
+    if (.not.allocated(coord)) allocate(coord(numxgridn))
      do i = 1,numxgridn
         coord(i) = outlon0n + (i-0.5)*dxoutn
      enddo
      call nf90_err(nf90_put_var(ncid, lonID, coord(1:numxgridn)))
+     deallocate(coord)
   else
+    if (.not.allocated(coord)) allocate(coord(numxgrid))
      do i = 1,numxgrid
         coord(i) = outlon0 + (i-0.5)*dxout
      enddo
      call nf90_err(nf90_put_var(ncid, lonID, coord(1:numxgrid)))
+     deallocate(coord)
   endif
   ! latitudes (grid cell centers)
   if (lnest) then
+    if (.not.allocated(coord)) allocate(coord(numygridn))
      do i = 1,numygridn
         coord(i) = outlat0n + (i-0.5)*dyoutn
      enddo
      call nf90_err(nf90_put_var(ncid, latID, coord(1:numygridn)))
+     deallocate(coord)
   else
+    if (.not.allocated(coord)) allocate(coord(numygrid))
      do i = 1,numygrid
         coord(i) = outlat0 + (i-0.5)*dyout
      enddo
      call nf90_err(nf90_put_var(ncid, latID, coord(1:numygrid)))
+     deallocate(coord)
   endif
   ! levels