1. summary: Version 2 of gx3 topography Hudson's Bay was opened up and Borneo was separated from Sumatra. 2. created by: Matthew Hecht 3. date and time: November 12, 2001 4. created on: dataproc 5. dataset format: ieeei4 6. endian: big 7. archival copy of this file is located at /CCSM/inputdata/ocn/pop/gx3v2/grid/topography_20011112.ieeei4 8. creation summary: Ran fortran code modkmt_20011112.f90. Output of that code, which records all changes, including 4 1-point holes which were filled, was: upon input, max of kmt = 25 upon input, min of kmt = 0 looping with i fastest, we see changes in these places: at 41 51 kmt changed from 0 to 5 at 41 52 kmt changed from 0 to 5 at 41 53 kmt changed from 0 to 4 at 42 53 kmt changed from 0 to 4 at 41 54 kmt changed from 0 to 3 at 42 54 kmt changed from 0 to 3 at 41 55 kmt changed from 0 to 3 at 42 55 kmt changed from 0 to 3 at 41 56 kmt changed from 0 to 4 at 42 56 kmt changed from 0 to 4 at 41 57 kmt changed from 0 to 4 at 42 57 kmt changed from 0 to 4 at 41 58 kmt changed from 0 to 4 at 42 58 kmt changed from 0 to 4 at 87 100 kmt changed from 0 to 3 at 90 106 kmt changed from 0 to 7 at 90 107 kmt changed from 0 to 7 at 91 107 kmt changed from 0 to 9 looping on j fastest, we see those same changes (but maybe in a more convenient form?): at 41 51 kmt changed from 0 to 5 at 41 52 kmt changed from 0 to 5 at 41 53 kmt changed from 0 to 4 at 41 54 kmt changed from 0 to 3 at 41 55 kmt changed from 0 to 3 at 41 56 kmt changed from 0 to 4 at 41 57 kmt changed from 0 to 4 at 41 58 kmt changed from 0 to 4 at 42 53 kmt changed from 0 to 4 at 42 54 kmt changed from 0 to 3 at 42 55 kmt changed from 0 to 3 at 42 56 kmt changed from 0 to 4 at 42 57 kmt changed from 0 to 4 at 42 58 kmt changed from 0 to 4 at 87 100 kmt changed from 0 to 3 at 90 106 kmt changed from 0 to 7 at 90 107 kmt changed from 0 to 7 at 91 107 kmt changed from 0 to 9 1 grid pt hole corrected at 97 13 1 grid pt hole corrected at 53 47 1 grid pt hole corrected at 15 102 1 grid pt hole corrected at 85 112 upon output, max of kmt = 25 upon output, min of kmt = 0 9. derived from: /NJN01/POP/INPUTDATA/100x116x25/GREENLAND/topography.ieee 10. code, scripts, namelists, etc, used to create this file: program modkmt c--------------------------------------------------------- c--------------------------------------------------------- implicit none integer :: i, j, k, im1, ip1, max_k4 integer, parameter :: imt=100, jmt=116 integer*4, dimension(imt,jmt) :: kmt, kmt_old integer*4, dimension(3) :: kmt_size open(unit=10, file='topography.ieee', status='old', $ form='unformatted', access='direct', recl=imt*jmt*4) read(10, rec=1) kmt kmt_old = kmt write(*,*) 'upon input, max of kmt = ', maxval(kmt) write(*,*) 'upon input, min of kmt = ', minval(kmt) c---- James Bay kmt(87,100) = kmt(88,100) c---- Hudson Strait kmt(90,106) = 7 kmt(90,107) = 7 kmt(91,107) = 9 c---- Borneo--Sumatra kmt(41:42,56:58) = 4 kmt(41:42,54:55) = 3 kmt(41:42,53) = 4 kmt(41:42,51:52) = 5 c---- check print*, ' ' print*, 'looping with i fastest, we see changes in these places:' do j=1,jmt do i=1,imt if (kmt_old(i,j) .ne. kmt(i,j)) $ write(*,1002) i, j, kmt_old(i,j), kmt(i,j) enddo enddo 1002 format('at ',i3,1x,i3,' kmt changed from ',i2,' to ',i2) print*, ' ' print*, 'looping on j fastest, we see those same changes' print*, '(but maybe in a more convenient form?):' do i=1,imt do j=1,jmt if (kmt_old(i,j) .ne. kmt(i,j)) $ write(*,1002) i, j, kmt_old(i,j), kmt(i,j) enddo enddo print*, ' ' c---- set min to 3 where(kmt .eq. 1) kmt = 3 where(kmt .eq. 2) kmt = 3 c---- fill 1-pt holes do j=2, jmt-1 do i=1,imt if (i .ne. 1) then im1 = i - 1 else im1 = imt endif if (i .ne. imt) then ip1 = i + 1 else ip1 = 1 endif max_k4 = max( $ kmt(im1, j ), $ kmt(ip1, j ), $ kmt(i, j-1), $ kmt(i, j+1) ) if (kmt(i,j) .gt. max_k4) then write(*, 1001) '1 grid pt hole corrected at ', i, j kmt(i,j) = max_k4 endif enddo enddo 1001 format(a, 2(i4, 1x)) write(*,*) 'upon output, max of kmt = ', maxval(kmt) write(*,*) 'upon output, min of kmt = ', minval(kmt) c---- write in POP's format open(unit=12, file='kmt_gx3v2_011112.da', $ status='new', form='unformatted', access='direct', $ recl=imt*jmt*4) write(12, rec=1) kmt stop end