program bathy use netcdf use mod_bathy_field implicit none include 'netcdf.inc' character*40,parameter :: sdir="/space3/ocean/seabottom/data/" character*40,parameter :: cdir="/space3/bathymetry/data/" character*40,parameter :: odir="/space3/levitas/newfortran/bathy_img/" character*12, parameter ::ofile='pippo.bat' character (len = *), parameter :: & fcu =trim(sdir)//"20150101_d-INGV--RFVL-MFSeas2-MEDATL.nc" character (len = 30), dimension(6) :: fg character(len=11) :: att_name integer :: ncid,ndims,nvars,ngatts,unlimdimid integer :: retval integer :: dimid,nrecs,xtype,natts,attid,tlen,vrlen integer :: dimids(NF_MAX_VAR_DIMS) integer ic,i,j type(Sbathy)::medb real(r8)::del fg(1) = trim(cdir)//"B3.mnt" fg(2) = trim(cdir)//"B4.mnt" fg(3) = trim(cdir)//"C3.mnt" fg(4) = trim(cdir)//"C4.mnt" fg(5) = trim(cdir)//"D3.mnt" fg(6) = trim(cdir)//"D4.mnt" medb=Sbathy(fcu,fg(this_image()),this_image()) del=(medb%gmed%box%p1%plon-medb%gmed%box%p0%plon)/medb%gmed%nlon if(this_image().eq.1) then write(*,'(i2,1x,a30,2(f8.4,1x,i5,f8.4),f8.3)'), this_image(),fg(this_image()),& medb%gmed%box%p0%plon,medb%gmed%nlon,medb%gmed%box%p1%plon,& medb%gmed%box%p0%plat,medb%gmed%nlat,medb%gmed%box%p1%plat,1./del end if del=(medb%gemo%box%p1%plon-medb%gemo%box%p0%plon)/medb%gemo%nlon print*,'-------------------------------------------' write(*,'(i2,1x,a30,2(f8.4,1x,i5,f8.4),f8.3)'), this_image(),fg(this_image()),& medb%gemo%box%p0%plon,medb%gemo%nlon,medb%gemo%box%p1%plon,& medb%gemo%box%p0%plat,medb%gemo%nlat,medb%gemo%box%p1%plat,1/del sync all call write_bat(medb,(trim(odir)//trim(ofile))) print*, this_image(), minval(medb%gmed%z), minloc(medb%gmed%z) print*,'-------------------------------------------' end program bathy