Changeset 3537 for palm/trunk/UTIL/inifor/src/inifor_io.f90

Timestamp:

Nov 20, 2018 10:53:14 AM (5 years ago)

Author:

eckhard

Message:

inifor: COSMO-D2 support

File:

• palm/trunk/UTIL/inifor/src/inifor_io.f90 (modified) (3 diffs)

palm/trunk/UTIL/inifor/src/inifor_io.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! New routine get_netcdf_dim_vector()
+! 
+! 
+! 3534 2018-11-19 15:35:16Z raasch
 ! bugfix: INTENT attribute changed
 ! 
@@ -173 +177 @@
 
 
+    SUBROUTINE get_netcdf_dim_vector(filename, varname, array)
+
+       CHARACTER(LEN=*), INTENT(IN)         ::  filename
+       CHARACTER(LEN=*), INTENT(IN)         ::  varname
+       REAL(dp), ALLOCATABLE, INTENT(INOUT) ::  array(:)
+
+       INTEGER ::  ncid, varid, dimlen
+       INTEGER ::  dimids(NF90_MAX_VAR_DIMS)
+
+       IF ( nf90_open( TRIM(filename), NF90_NOWRITE, ncid ) .EQ. NF90_NOERR .AND. &
+            nf90_inq_varid( ncid, varname, varid ) .EQ. NF90_NOERR )  THEN
+
+          CALL check(nf90_inquire_variable( ncid, varid, dimids = dimids ))
+          CALL check(nf90_inquire_dimension( ncid, dimids(1), len = dimlen ))
+
+          ALLOCATE(array(dimlen))
+          CALL check(nf90_get_var( ncid, varid, array ))
+
+       ELSE
+
+          message = "Failed to read '" // TRIM(varname) // &
+             "' from file '" // TRIM(filename) // "'."
+          CALL abort('get_netcdf_dim_vector', message)
+
+       END IF
+
+    END SUBROUTINE get_netcdf_dim_vector
+
+
     SUBROUTINE get_input_dimensions(in_var, ncid)
 
@@ -209 +242 @@
 
           message = "Failed reading NetCDF variable " //                       &
-             TRIM(in_var % name) // " with " // TRIM(str(in_var%ndim)) //      &
+             TRIM(in_var % name) // " with " // TRIM(str(in_var % ndim)) //    &
              " dimensions because only two- and and three-dimensional" //      &
              " variables are supported."