convert data to chunk stucture in cam4

[wangxuejiao14@lzu_edu_cn]

Dear all,      I have read some variables(eg:ua1(144,96,26,1))  from a netcdf file,these variables should be converted to the "chunk" structure,so how can I  convert  these variables(eg:u(144,96,26,1)) to the chunk structure(u(i,k)) ?I have read this poster http://bb.cgd.ucar.edu/node/1002286, but I can't understand better,can you enxplain me in detail ? Thank you ! 

 

[wangxuejiao14@lzu_edu_cn]

or should I use this in physpkg.F90do c = begchunk,endchunk        ncol = phys_state(c)%ncol        call get_lon_all_p(c, ncol, lons)        call get_lat_all_p(c, ncol, lats)        do k=1,pver          do i=1,ncol             ic = lons(i)             jc = lats(i)             ua1(i,k) = u_A1(ic,jc,k,1)             va1(i,k) = v_A1(ic,jc,k,1)             ta1(i,k) = t_A1(ic,jc,k,1)             qa1(i,k) = q_A1(ic,jc,k,1)          enddo        enddo     enddoThank you!

 

[wangxuejiao14@lzu_edu_cn]

or should I use this in physpkg.F90do c = begchunk,endchunk        ncol = phys_state(c)%ncol        call get_lon_all_p(c, ncol, lons)        call get_lat_all_p(c, ncol, lats)        do k=1,pver          do i=1,ncol             ic = lons(i)             jc = lats(i)             ua1(i,k) = u_A1(ic,jc,k,1)             va1(i,k) = v_A1(ic,jc,k,1)             ta1(i,k) = t_A1(ic,jc,k,1)             qa1(i,k) = q_A1(ic,jc,k,1)          enddo        enddo     enddoThank you!

 

[wangxuejiao14@lzu_edu_cn]

or should I use this in physpkg.F90do c = begchunk,endchunk        ncol = phys_state(c)%ncol        call get_lon_all_p(c, ncol, lons)        call get_lat_all_p(c, ncol, lats)        do k=1,pver          do i=1,ncol             ic = lons(i)             jc = lats(i)             ua1(i,k) = u_A1(ic,jc,k,1)             va1(i,k) = v_A1(ic,jc,k,1)             ta1(i,k) = t_A1(ic,jc,k,1)             qa1(i,k) = q_A1(ic,jc,k,1)          enddo        enddo     enddoThank you!