change netcdf read for shelf sea mask

hamocc/mo_read_shelfmask.F90

@@ -48,6 +48,8 @@ subroutine read_shelfmask(kpie,kpje,kbnd,pbath,omask)
     use mod_nctools,    only: ncfopn,ncread,ncfcls
     use mo_control_bgc, only: use_shelfsea_res_time,io_stdo_bgc
     use mo_param_bgc,   only: shelfbreak_depth
+    use netcdf,         only: nf90_open,nf90_close,nf90_nowrite
+    use mo_netcdf_bgcrw,only: read_netcdf_var
 
     !Arguments
     integer, intent(in) :: kpie                                     ! 1st dimension of model grid
@@ -58,7 +60,7 @@ subroutine read_shelfmask(kpie,kpje,kbnd,pbath,omask)
 
     ! Local variables
     logical :: file_exists=.false.
-    integer :: i,j,errstat,dummymask(2)
+    integer :: i,j,errstat,dummymask(2),ncid,ncstat
     real,allocatable  :: mask(:,:)
 
 
@@ -98,9 +100,14 @@ subroutine read_shelfmask(kpie,kpje,kbnd,pbath,omask)
         write(io_stdo_bgc,*) ''
         write(io_stdo_bgc,'(a)') 'read_shelfmask: read mask from ',trim(shelfsea_maskfile)
       endif
-      call ncfopn(trim(shelfsea_maskfile),'r',' ',1,iotype)
-      call ncread('shelfmask',mask,dummymask,0,0.)
-      call ncfcls
+
+      ncstat=nf90_open(trim(shelfsea_maskfile),nf90_nowrite,ncid)
+      call read_netcdf_var(ncid,'shelfmask',mask,1,1,0)
+      ncstat=nf90_close(ncid)
+
+      !call ncfopn(trim(shelfsea_maskfile),'r',' ',1,iotype)
+      !call ncread('shelfmask',mask,dummymask,0,0.)
+      !call ncfcls
     else
       ! reconstruct shelf sea mask from internal bathymetry
       !$OMP DO PARALLEL PRIVATE (i,j)