anchor.diff:              call check( nf90_def_var(ncid, "w", NF90_REAL, dimids, temp_id), trim(err)//"w" )
               call check( nf90_put_att(ncid,temp_id,"standard_name","upward_air_velocity"))
               call check( nf90_put_att(ncid,temp_id,"long_name","Vertical wind"))
  -            call check( nf90_put_att(ncid,temp_id,"WARNING","Grid relative (i.e. add u*dz/dx) and scaled by dx/dz"))
  +!             call check( nf90_put_att(ncid,temp_id,"WARNING","Grid relative (i.e. add u*dz/dx) and scaled by dx/dz"))
               call check( nf90_put_att(ncid,temp_id,"units","m s-1"))
               varid(11)=temp_id
               
  
positive.code:             varid(30)=temp_id
              call check( nf90_inq_varid(ncid, "canwat", temp_id), trim(err)//"canwat" )
              varid(31)=temp_id
          endif
          call check( nf90_inq_varid(ncid, "ta2m", temp_id), trim(err)//"ta2m" )
          varid(34)=temp_id
          call check( nf90_inq_varid(ncid, "hus2m", temp_id), trim(err)//"hus2m" )
          varid(35)=temp_id
          call check( nf90_inq_varid(ncid, "u10m", temp_id), trim(err)//"u10m" )
          varid(36)=temp_id
          call check( nf90_inq_varid(ncid, "v10m", temp_id), trim(err)//"v10m" )
  
negative.code:             call check( nf90_put_var(ncid, varid(29), domain%lwup,         start_two_D),  trim(filename)//":lwup" )
              call check( nf90_put_var(ncid, varid(30), domain%snow_swe,     start_two_D),  trim(filename)//":snow_swe" )
              call check( nf90_put_var(ncid, varid(31), domain%canopy_water, start_two_D),  trim(filename)//":canopy_water" )
              call check( nf90_put_var(ncid, varid(34), domain%T2m, start_two_D),  trim(filename)//":t2m" )
              call check( nf90_put_var(ncid, varid(35), domain%Q2m, start_two_D),  trim(filename)//":hus2m" )
          endif
          call check( nf90_put_var(ncid, varid(36), domain%u10, start_two_D),  trim(filename)//":u10m" )
          call check( nf90_put_var(ncid, varid(37), domain%v10, start_two_D),  trim(filename)//":v10m" )
          
  