[netcdfgroup] data read and transfer

• To: netcdfgroup@xxxxxxxxxxxxxxxx
• Subject: [netcdfgroup] data read and transfer
• From: "H.Dang" <danghy@xxxxxxxxx>
• Date: Wed, 24 Mar 2010 14:14:32 -0500

Hello,

I was just wondering if you could find obvious mistakes in below code that I
read an NC file. Code:

     DO nc_S = 1, N_GC_SPC
        call nc_chk1_rd1(bg_ncid,GC_SPC(nc_s), chk1_CGRID, nc_s,
nstep_nreps)
     ENDDO

The subroutine is:

     28    subroutine nc_chk1_rd1(chk_ncid, nc_data_name, chk1_data, N_ic,
nstep_nreps)
     29
     30     use netcdf
     31     use nc_adj_module
     32     USE GRID_CONF
     33     USE CGRID_SPCS
     34     USE KPP_Integrator
     35
     36     implicit none
     37
     38       integer :: chk_ncid
     39       character*16 nc_data_name
     40       real, intent(out) :: chk1_data(:,:,:,:)
     41       integer :: N_ic,nstep_nreps
     42       integer :: start1(4), count1(4)
     43       integer :: nc_data_varid
     44
     45       start1 = (/ 1, 1, 1, nstep_nreps /)
     46       count1 = (/ NCOLS, NROWS, NLAYS, 1 /)
     47
     48        call nc_check( nf90_inq_varid(chk_ncid, nc_data_name,
nc_data_varid))
     49        call nc_check( nf90_get_var(chk_ncid, nc_data_varid &
     50                  , chk1_data(:,:,:,N_ic), start1, count1
))
     51
     52    endsubroutine nc_chk1_rd1

The program is running now, but I'm just afraid I can miss something in the
code. I'm especially worring about the data transfer of the variable in red
color. Please let me know if you think is okay. Any comments welcome!


-- 
Cordially,
Hongyan
鸿雁