Ticket #583: fldread.patch
File fldread.patch, 2.3 KB (added by sga, 14 years ago) |
---|
-
fldread.F90
814 814 aname = ' ' 815 815 WRITE(aname,'(a3,i2.2)') 'src',jn 816 816 data_tmp(:,:) = 0 817 CALL iom_get ( inum, jpdom_unknown, aname, data_tmp(1:nlci,1:nlcj), & 818 kstart=(/nimpp,njmpp/), kcount=(/nlci,nlcj/) ) 817 CALL iom_get ( inum, jpdom_autoglo, aname, data_tmp(1:jpi,1:jpj) ) 819 818 data_src(:,:) = INT(data_tmp(:,:)) 820 819 ref_wgts(nxt_wgt)%data_jpj(:,:,jn) = 1 + (data_src(:,:)-1) / ref_wgts(nxt_wgt)%ddims(1) 821 820 ref_wgts(nxt_wgt)%data_jpi(:,:,jn) = data_src(:,:) - ref_wgts(nxt_wgt)%ddims(1)*(ref_wgts(nxt_wgt)%data_jpj(:,:,jn)-1) … … 824 823 DO jn = 1, ref_wgts(nxt_wgt)%numwgt 825 824 aname = ' ' 826 825 WRITE(aname,'(a3,i2.2)') 'wgt',jn 827 ref_wgts(nxt_wgt)%data_wgt(1:nlci,1:nlcj,jn) = 0.0 828 CALL iom_get ( inum, jpdom_unknown, aname, ref_wgts(nxt_wgt)%data_wgt(1:nlci,1:nlcj,jn), & 829 kstart=(/nimpp,njmpp/), kcount=(/nlci,nlcj/) ) 826 ref_wgts(nxt_wgt)%data_wgt(1:jpi,1:jpj,jn) = 0.0 827 CALL iom_get ( inum, jpdom_autoglo, aname, ref_wgts(nxt_wgt)%data_wgt(1:jpi,1:jpj,jn) ) 830 828 END DO 831 829 CALL iom_close (inum) 832 830 833 831 ! find min and max indices in grid 834 ref_wgts(nxt_wgt)%botleft(1) = MINVAL(ref_wgts(nxt_wgt)%data_jpi(1: nlci,1:nlcj,:))835 ref_wgts(nxt_wgt)%botleft(2) = MINVAL(ref_wgts(nxt_wgt)%data_jpj(1: nlci,1:nlcj,:))836 ref_wgts(nxt_wgt)%topright(1) = MAXVAL(ref_wgts(nxt_wgt)%data_jpi(1: nlci,1:nlcj,:))837 ref_wgts(nxt_wgt)%topright(2) = MAXVAL(ref_wgts(nxt_wgt)%data_jpj(1: nlci,1:nlcj,:))832 ref_wgts(nxt_wgt)%botleft(1) = MINVAL(ref_wgts(nxt_wgt)%data_jpi(1:jpi,1:jpj,:)) 833 ref_wgts(nxt_wgt)%botleft(2) = MINVAL(ref_wgts(nxt_wgt)%data_jpj(1:jpi,1:jpj,:)) 834 ref_wgts(nxt_wgt)%topright(1) = MAXVAL(ref_wgts(nxt_wgt)%data_jpi(1:jpi,1:jpj,:)) 835 ref_wgts(nxt_wgt)%topright(2) = MAXVAL(ref_wgts(nxt_wgt)%data_jpj(1:jpi,1:jpj,:)) 838 836 839 837 ! and therefore dimensions of the input box 840 838 ref_wgts(nxt_wgt)%jpiwgt = ref_wgts(nxt_wgt)%topright(1) - ref_wgts(nxt_wgt)%botleft(1) + 1