[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Characters arrays with f90

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.

  • Subject: Re: Characters arrays with f90
  • Date: Thu, 10 Oct 2002 16:12:03 -0600
>To: address@hidden
>From: Francois Vandenberghe <address@hidden>
>Subject: Characters arrays with f90
>Organization: UCAR/Unidata
>Keywords: 200210091704.g99H48109566


Hi,

> I have a problem understanding the way NETCDF f90 interface can load
> an array of characters. I tried to follow page 66 (7.7 Reading and writing
> character string values) of the documentation vserion 3 (February 200)
> but it doesn't work as I would like.
> 
> The problem is the following: I have a station names arrays "Staname"
> which contains "recNum" names (22 in that case) of "stanamLen" characters
> (6 in that case).
> 
> After reading both dimension recNum and stanamLen, I have the following
> statements:
> 
> CHARACTER, DIMENSION (:,:), ALLOCATABLE :: staName
> 
> ALLOCATE (staName (recNum,StanamLen))
> 
> iret = NF90_GET_VAR (nc_id, var_id (n), staName , map = (/ stanamLen, 1 /))
> 
> The output however gives:
> 
> 10:31:06 staName (  1) = C
> 10:31:06 staName (  2) = F
> 10:31:06 staName (  3) = L
> 10:31:06 staName (  4) = M
> 10:31:06 staName (  5) = O
> 10:31:06 staName (  6) = P
> 10:31:06 staName (  7) = R
> 10:31:06 staName (  8) = S
> 10:31:06 staName (  9) = C
> 10:31:06 staName ( 10) = G
> 10:31:06 staName ( 11) = G
> 10:31:06 staName ( 12) = L
> 10:31:06 staName ( 13) = L
> 10:31:06 staName ( 14) = N
> 10:31:06 staName ( 15) = O
> 10:31:06 staName ( 16) = P
> 10:31:06 staName ( 17) = P
> 10:31:06 staName ( 18) = S
> 10:31:06 staName ( 19) = W
> 10:31:06 staName ( 20) = S
> 10:31:06 staName ( 21) = V
> 10:31:06 staName ( 22) = A
> 
> while I am trying to read (from ncdump):
> 
>  staName =
>   "CCDNH ",
>   "FMEMD ",
>   "LAXCA ",
>   "MCGQU ",
>   "ONTCA ",
>   "PLTCO ",
>   "RUTNJ ",
>   "STWMA ",
>   "CCOCA ",
>   "GLACA ",
>   "GVYCA ",
>   "LHSCA ",
>   "LVRCA ",
>   "NPTOR ",
>   "OREMA ",
>   "PSPNY ",
>   "PYMMA ",
>   "SCHNY ",
>   "WFDCA ",
>   "SACCA ",
>   "VISCA ",
> 
> Any ideas what I am not doing correctly?

I'm not a Fortran90 expert, but I think you have reversed the
dimension orders in the ALLOCATE statement, which I would write:

 ALLOCATE (staName (StanamLen,recNum))

because the StanamLen dimension is the most rapidly varying, so it
should appear first.  The output of ncdump uses the opposite dimension
order from Fortran, so that may the cause of the confusion.

--Russ

_____________________________________________________________________

Russ Rew                                         UCAR Unidata Program
address@hidden                     http://www.unidata.ucar.edu