Hi Sudheer,
There is no problem with dynamically allocating an array in fortran90 and
then passing that array to the netcdf routines that read a variable. The
array does *not* have to be dimensioned with parameters. I've done this
many, many times with several different fortran compilers, including
portland group, g95, gfortran, and intel fortran.
It's hard to tell what the problem you are running into is, since you say
"it will not accept as it is not defined as a parameter." *What* will not
accept it? The fortran compiler? Which compiler are you using? It's
completely standard, run of the mill fortran90 to be able to use
dynamically allocated arrays everywhere where you can use an array defined
using parameters. Without seeing the problematic code and the error
message it gives, my best suggestion is that you simply have an ordinary
fortran bug somewhere, and it doesn't have anything directly to do with
netcdf.
Regards,
--Dave
Sudheer Joseph wrote:
> Thank you Dave,
> The subroutine you mentioned is very
> useful. But my question was slightly different. suppose if you get the
> dimension size before reading the 3 dimensional variable. we cannot
> use that information to read the netcdf 3rd dimension at a stretch.
> I mean
> if I use the
>
> NTIM= nf_get_dim(time) ( just for illustration sake )
> and use this for allocating the dimension for 3D variable it will not
> accept as it is not defined as a parameter.
>
> allocate(pres_in(NLONS, NLATS, NTIM))
> count = (/ NLONS, NLATS,NTIM/)
> start = (/ 1, 1, 1/)
>
> One possibility came to my mind yesterday night was to read the data
> as xy slabs and use a do loop to read the time dimension. It appears
> there is no way to read the 3D variable at one call if we are not
> defining the dimension as parameter.
> with best regards,
> Sudheer
>
> On Fri, Jun 11, 2010 at 3:42 AM, Dave Allured <dave.allured@xxxxxxxx>
> wrote:
>> Sudheer,
>>
>> If the names of the dimensions are always the same, then use
>> nf90_inq_dimid followed by nf90_inquire_dimension to get the length
>> of each dimension. See the example on this page:
>>
>> http://www.unidata.ucar.edu/software/netcdf/docs/netcdf-f90/NF90_005fINQUIRE_005fDIMENSION.html
>>
>> If all the dimensions have coordinate variables, then the following
>> subroutine will get several dimension lengths in one call:
>>
>> http://www.esrl.noaa.gov/psd/people/dave.allured/data/f90/lib/getnc_dims.f90
>>
>> If the dimension names change between Netcdf files, then the problem
>> is more complicated. We can address that if neither of the above
>> works for you.
>>
>> --Dave A.
>> NOAA/PSD/CIRES
>>
>> On 6/10/2010 7:55 AM, Sudheer Joseph wrote:
>>> Dear Users,
>>> I have been coding a program to read multiple
>>> netcdf files having different time dimensions. I want to know if there
>>> is a way to read netcdf varibale with out defining the dimension as
>>> parameter. with in fortran 90 program.
>>> normally we need to define like below
>>> ! timesteps of data.
>>> integer, parameter :: NDIMS = 3, NRECS =1
>>> integer, parameter :: NTIM=1464, NLATS = 94, NLONS = 192
>>> allocate(pres_in(NLONS, NLATS, NTIM))count = (/ NLONS, NLATS,NTIM/)
>>> start = (/ 1, 1, 1/)
>>>
>>> ! Read the surface pressure and temperature data from the file,
>>> one
>>> ! record at a time.
>>> call check( nf90_get_var(ncid, pres_varid, pres_in, start =
>>> start,count = count) )
>>>
>>> In my case ntim may be different for different files. I wanted to use
>>> the reading part as a subroutine where ntim is provided as a variable
>>> and passed to subroutine
>>> but if I declare ntim as varible then netcdf gives an error.
>>> If any one have a solution please help. I am not looking for solutions
>>> like replacing the parameter statement using shell script and re
>>> compiling it etc.
>>>
>>> with best regards,
>>> Sudheer
>>>
>>> _______________________________________________
>>> netcdfgroup mailing list
>>> netcdfgroup@xxxxxxxxxxxxxxxx
>>> For list information or to unsubscribe, visit:
>>> http://www.unidata.ucar.edu/mailing_lists/
>>
>> _______________________________________________
>> netcdfgroup mailing list
>> netcdfgroup@xxxxxxxxxxxxxxxx
>> For list information or to unsubscribe, visit:
>> http://www.unidata.ucar.edu/mailing_lists/
>>
>
>
>
> --
> with best regards
>
> Sudheer
>
> **********************************************************************************
> Sudheer Joseph
> Scientist
> Indian National Centre for Ocean Information Services (INCOIS)
> Ocean Valley, Post Box No# 21,
> IDA Jeedimetla P.O.
> Hyderabad, Ranga Reddy District - 500 055
> Andhra Pradesh, India.
> TEl:+91-40-23044600(R),Tel:+91-9440832534(Mobile)
> Tel:+91-40-23886047(O),Fax:+91-40-23892910(O)
> E-mail:sjo@xxxxxxxxxxxxxx; sudheer.joseph@xxxxxxxxx; sjo@xxxxxxxxxxxxx.
> Web- http://oppamthadathil.tripod.com
> --------------* ---------------
> "The ultimate measure of a man is
> not where he stands in moments of
> comfort and convenience, but where
> he stands at times of challenge and
> controversy."
> Martin Luther King, Jr.
> "What we have done for ourselves alone dies with us.
> What we have done for others and the world remains and is immortal."
> - Albert Pines
>
> _______________________________________________
> netcdfgroup mailing list
> netcdfgroup@xxxxxxxxxxxxxxxx
> For list information or to unsubscribe, visit:
> http://www.unidata.ucar.edu/mailing_lists/
>
>
-------------------------------------------------------------------
David W. Pierce
Division of Climate, Atmospheric Science, and Physical Oceanography
Scripps Institution of Oceanography
(858) 534-8276 (voice) / (858) 534-8561 (fax) dpierce@xxxxxxxx
-------------------------------------------------------------------
