id	summary	reporter	owner	description	type	status	priority	milestone	component	version	resolution	keywords	cc
946	SetFillVal argument missing when writing vertical coordinates	wilkin		"* The variables associated with depth are special in ROMS because of wetting/drying and other algorithms. We cannot add a **_!FillValue** in output NetCDF files. Although the output variables **z_rho**, **z_u**, **z_v**, and **z_w** are defined with **!SetFillVal = .FALSE.** in module **def_his.F**, the same optional variable is missing when writing in **wrt_his.F** and **wrt_quick**. For example, we need to have in **wrt_his.F** instead:
 {{{
!
!  Write time-varying depths of RHO-points.
!
      IF (Hout(idpthR,ng)) THEN
        scale=1.0_dp
        gtype=gfactor*r3dvar
        status=nf_fwrite3d(ng, model, HIS(ng)%ncid,                     &
     &                     HIS(ng)%Vid(idpthR),                         &
     &                     HIS(ng)%Rindex, gtype,                       &
     &                     LBi, UBi, LBj, UBj, 1, N(ng), scale,         &
# ifdef MASKING
     &                     GRID(ng) % rmask,                            &
# endif
     &                     GRID(ng) % z_r,                              &
     &                     SetFillVal = .FALSE.)
        IF (FoundError(status, nf90_noerr, __LINE__, MyFile)) THEN
          IF (Master) THEN
            WRITE (stdout,10) TRIM(Vname(1,idpthR)), HIS(ng)%Rindex
          END IF
          exit_flag=3
          ioerror=status
          RETURN
        END IF
      END IF
}}}
 To avoid those land values being read in MATLAB as 1e+37. Those variables were written correctly in **wrt_dai.F**, which is used for EnKF and 4D-Var.

* Corrected typo in **def_his.F** when defining implicit omega vertical velocity.

* Updated metadata **varinfo.yaml** to include **omega_implicit** and added scale factor **1/rho0** for wave dissipation variables."	bug	closed	minor	Release ROMS/TOMS 4.2	Nonlinear	4.1	Fixed