def_rst.F: Corrected the output dimension for wetting and drying variables wetdry_mask_u and wetdry_mask_v. We need to have instead:
!
! Define wet/dry mask on U-points.
!
Vinfo( 1)=Vname(1,idUwet)
Vinfo( 2)=Vname(2,idUwet)
Vinfo( 3)=Vname(3,idUwet)
Vinfo( 9)='land'
Vinfo(10)='water'
Vinfo(14)=Vname(4,idUwet)
Vinfo(16)=Vname(1,idtime)
Vinfo(22)='coordinates'
Aval(5)=REAL(Iinfo(1,idUwet,ng),r8)
status=def_var(ng, iNLM, RST(ng)%ncid, RST(ng)%Vid(idUwet), &
& NF_FOUT, nvd3, su2dgrd, Aval, Vinfo, ncname, &
& SetFillVal = .FALSE.)
IF (exit_flag.ne.NoError) RETURN
!
! Define wet/dry mask on V-points.
!
Vinfo( 1)=Vname(1,idVwet)
Vinfo( 2)=Vname(2,idVwet)
Vinfo( 3)=Vname(3,idVwet)
Vinfo(14)=Vname(4,idVwet)
Vinfo(16)=Vname(1,idtime)
Vinfo( 9)='land'
Vinfo(10)='water'
Vinfo(22)='coordinates'
Aval(5)=REAL(Iinfo(1,idVwet,ng),r8)
status=def_var(ng, iNLM, RST(ng)%ncid, RST(ng)%Vid(idVwet), &
& NF_FOUT, nvd3, sv2dgrd, Aval, Vinfo, ncname, &
& SetFillVal = .FALSE.)
IF (exit_flag.ne.NoError) RETURN
Many thanks to Lyon Lanerolle for bringing this to my attention.