tripolegrid.F90

Go to the documentation of this file.

module tripolegrid

  use gengrid_kinds, only: dbl_kind,int_kind,cm
  use grdvars,       only: ni,nj,nv,mastertask,nverts,ncoord
  use grdvars,       only: lonct,latct,lonct_vert,latct_vert
  use grdvars,       only: loncu,latcu,loncu_vert,latcu_vert
  use grdvars,       only: loncv,latcv,loncv_vert,latcv_vert
  use grdvars,       only: lonbu,latbu,lonbu_vert,latbu_vert
  use grdvars,       only: wet4,areact,anglet,dp4
  use charstrings,   only: logmsg,history
  use vartypedefs,   only: maxvars, fixvars, fixvars_typedefine
  use netcdf

  implicit none
  private

  public write_tripolegrid

  contains
  
  subroutine write_tripolegrid(fname)

  character(len=*), intent(in) :: fname

  ! local variables
  integer :: ii,id,rc, ncid, dim2(2),dim3(3)
  integer :: idimid,jdimid,kdimid

!---------------------------------------------------------------------
! create the netcdf file
!---------------------------------------------------------------------

  ! define the output variables and file name
  call fixvars_typedefine

  ! create the file
  ! 64_bit offset reqd for 008 grid
  ! produces b4b results for smaller grids
  rc = nf90_create(trim(fname), nf90_64bit_offset, ncid)
  if(mastertask) then
    logmsg = '==> writing tripole grid to '//trim(fname)
    print '(a)', trim(logmsg)
    if(rc .ne. 0)print '(a)', 'nf90_create = '//trim(nf90_strerror(rc))
  end if

  rc = nf90_def_dim(ncid, 'ni', ni, idimid)
  rc = nf90_def_dim(ncid, 'nj', nj, jdimid)
  rc = nf90_def_dim(ncid, 'nv', nv, kdimid)
  
  !mask
  dim2(:) = (/idimid, jdimid/)
   rc = nf90_def_var(ncid, 'wet',     nf90_int,   dim2, id)
   rc = nf90_put_att(ncid, id,     'units',           'nd')
  !area
   rc = nf90_def_var(ncid, 'area', nf90_double,   dim2, id)
   rc = nf90_put_att(ncid, id,     'units',           'm2')
  !angleT
   rc = nf90_def_var(ncid, 'anglet', nf90_double, dim2, id)
   rc = nf90_put_att(ncid, id,     'units',      'radians')
  !bathymetry
   rc = nf90_def_var(ncid,  'depth', nf90_float,  dim2, id)
   rc = nf90_put_att(ncid, id,     'units',            'm')

  dim2(:) = (/idimid, jdimid/)
  do ii = 1,ncoord
   rc = nf90_def_var(ncid, trim(fixvars(ii)%var_name), nf90_double, dim2, id)
   rc = nf90_put_att(ncid, id,     'units', trim(fixvars(ii)%unit_name))
   rc = nf90_put_att(ncid, id, 'long_name', trim(fixvars(ii)%long_name))
   if(trim(fixvars(ii)%var_name(1:3)) .eq. "lon")then
    rc = nf90_put_att(ncid, id,  'lon_bounds', trim(fixvars(ii)%vertices))
   else
    rc = nf90_put_att(ncid, id,  'lat_bounds', trim(fixvars(ii)%vertices))
   endif
  enddo

  dim3(:) = (/idimid, jdimid, kdimid/)
  do ii = ncoord+1,ncoord+nverts
   rc = nf90_def_var(ncid, trim(fixvars(ii)%var_name), nf90_double, dim3, id)
   rc = nf90_put_att(ncid, id,     'units', trim(fixvars(ii)%unit_name))
   rc = nf90_put_att(ncid, id, 'long_name', trim(fixvars(ii)%long_name))
  enddo

  rc = nf90_put_att(ncid, nf90_global, 'history', trim(history))
  rc = nf90_enddef(ncid)

  rc = nf90_inq_varid(ncid,   'wet',        id)
  rc = nf90_put_var(ncid,        id, int(wet4))

  rc = nf90_inq_varid(ncid,  'area',      id)
  rc = nf90_put_var(ncid,        id,  areact)

  rc = nf90_inq_varid(ncid,'anglet',      id)
  rc = nf90_put_var(ncid,        id,  anglet)

  rc = nf90_inq_varid(ncid, 'depth',      id)
  rc = nf90_put_var(ncid,        id,     dp4)

  rc = nf90_inq_varid(ncid,  'lonCt',     id)
  rc = nf90_put_var(ncid,        id,   lonct)

  rc = nf90_inq_varid(ncid,  'latCt',     id)
  rc = nf90_put_var(ncid,        id,   latct)

  rc = nf90_inq_varid(ncid, 'lonCv',      id)
  rc = nf90_put_var(ncid,        id,   loncv)

  rc = nf90_inq_varid(ncid, 'latCv',      id)
  rc = nf90_put_var(ncid,        id,   latcv)
  
  rc = nf90_inq_varid(ncid, 'lonCu',      id)
  rc = nf90_put_var(ncid,        id,   loncu)

  rc = nf90_inq_varid(ncid, 'latCu',      id)
  rc = nf90_put_var(ncid,        id,   latcu)

  rc = nf90_inq_varid(ncid, 'lonBu',      id)
  rc = nf90_put_var(ncid,        id,   lonbu)

  rc = nf90_inq_varid(ncid, 'latBu',      id)
  rc = nf90_put_var(ncid,        id,   latbu)

  ! vertices
  rc = nf90_inq_varid(ncid,  'lonCt_vert',     id)
  rc = nf90_put_var(ncid,         id,  lonct_vert)

  rc = nf90_inq_varid(ncid,  'latCt_vert',     id)
  rc = nf90_put_var(ncid,         id,  latct_vert)

  rc = nf90_inq_varid(ncid, 'lonCv_vert',      id)
  rc = nf90_put_var(ncid,        id,   loncv_vert)

  rc = nf90_inq_varid(ncid, 'latCv_vert',      id)
  rc = nf90_put_var(ncid,        id,   latcv_vert)

  rc = nf90_inq_varid(ncid, 'lonCu_vert',      id)
  rc = nf90_put_var(ncid,        id,   loncu_vert)

  rc = nf90_inq_varid(ncid, 'latCu_vert',      id)
  rc = nf90_put_var(ncid,        id,   latcu_vert)

  rc = nf90_inq_varid(ncid, 'lonBu_vert',      id)
  rc = nf90_put_var(ncid,        id,   lonbu_vert)

  rc = nf90_inq_varid(ncid, 'latBu_vert',      id)
  rc = nf90_put_var(ncid,        id,   latbu_vert)

  rc = nf90_close(ncid)

  end subroutine write_tripolegrid
end module tripolegrid