- Timestamp:
- 2010-07-09T14:40:45+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/SBC/fldread.F90
r1970 r2004 838 838 WRITE(aname,'(a3,i2.2)') 'src',jn 839 839 data_tmp(:,:) = 0 840 CALL iom_get ( inum, jpdom_unknown, aname, data_tmp(1:nlci,1:nlcj), & 841 kstart=(/nimpp,njmpp/), kcount=(/nlci,nlcj/) ) 840 CALL iom_get ( inum, jpdom_data, aname, data_tmp(:,:) ) 842 841 data_src(:,:) = INT(data_tmp(:,:)) 843 842 ref_wgts(nxt_wgt)%data_jpj(:,:,jn) = 1 + (data_src(:,:)-1) / ref_wgts(nxt_wgt)%ddims(1) … … 848 847 aname = ' ' 849 848 WRITE(aname,'(a3,i2.2)') 'wgt',jn 850 ref_wgts(nxt_wgt)%data_wgt(1:nlci,1:nlcj,jn) = 0.0 851 CALL iom_get ( inum, jpdom_unknown, aname, ref_wgts(nxt_wgt)%data_wgt(1:nlci,1:nlcj,jn), & 852 kstart=(/nimpp,njmpp/), kcount=(/nlci,nlcj/) ) 849 ref_wgts(nxt_wgt)%data_wgt(:,:,jn) = 0.0 850 CALL iom_get ( inum, jpdom_data, aname, ref_wgts(nxt_wgt)%data_wgt(:,:,jn) ) 853 851 END DO 854 852 CALL iom_close (inum) 855 853 856 854 ! find min and max indices in grid 857 ref_wgts(nxt_wgt)%botleft(1) = MINVAL(ref_wgts(nxt_wgt)%data_jpi( 1:nlci,1:nlcj,:))858 ref_wgts(nxt_wgt)%botleft(2) = MINVAL(ref_wgts(nxt_wgt)%data_jpj( 1:nlci,1:nlcj,:))859 ref_wgts(nxt_wgt)%topright(1) = MAXVAL(ref_wgts(nxt_wgt)%data_jpi( 1:nlci,1:nlcj,:))860 ref_wgts(nxt_wgt)%topright(2) = MAXVAL(ref_wgts(nxt_wgt)%data_jpj( 1:nlci,1:nlcj,:))855 ref_wgts(nxt_wgt)%botleft(1) = MINVAL(ref_wgts(nxt_wgt)%data_jpi(:,:,:)) 856 ref_wgts(nxt_wgt)%botleft(2) = MINVAL(ref_wgts(nxt_wgt)%data_jpj(:,:,:)) 857 ref_wgts(nxt_wgt)%topright(1) = MAXVAL(ref_wgts(nxt_wgt)%data_jpi(:,:,:)) 858 ref_wgts(nxt_wgt)%topright(2) = MAXVAL(ref_wgts(nxt_wgt)%data_jpj(:,:,:)) 861 859 862 860 ! and therefore dimensions of the input box
Note: See TracChangeset
for help on using the changeset viewer.