anchor.diff:          call check( nf90_put_att(ncid,temp_id,"standard_name","eastward_10m_wind_speed"))
           call check( nf90_put_att(ncid,temp_id,"long_name","Eastward wind speed at 2m"))
           call check( nf90_put_att(ncid,temp_id,"units","m s-1"))
  +        call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
           varid(36)=temp_id
   
           call check( nf90_def_var(ncid, "v10m", NF90_REAL, dimtwo_time, temp_id), trim(err)//"v10m" )
           call check( nf90_put_att(ncid,temp_id,"standard_name","northward_10m_wind_speed"))
           call check( nf90_put_att(ncid,temp_id,"long_name","Northward wind speed at 2m"))
           call check( nf90_put_att(ncid,temp_id,"units","m s-1"))
  +        call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
           varid(37)=temp_id
   
           call check( nf90_def_var(ncid, "hfss", NF90_REAL, dimtwo_time, temp_id), trim(err)//"hfss" )
  
positive.code:             call check( nf90_put_att(ncid,temp_id,"standard_name","soil_moisture_content"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Column Soil Moisture"))
              call check( nf90_put_att(ncid,temp_id,"units","kg m-2"))
              varid(26)=temp_id
  
              call check( nf90_def_var(ncid, "soil_t", NF90_REAL, dimids, temp_id), trim(err)//"soil_t" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","soil_temperature"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Column Soil Temperature"))
              call check( nf90_put_att(ncid,temp_id,"units","K"))
              varid(27)=temp_id
  
              call check( nf90_def_var(ncid, "snw", NF90_REAL, dimtwo_time, temp_id), trim(err)//"snw" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","surface_snow_amount"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Snow water equivalent"))
              call check( nf90_put_att(ncid,temp_id,"units","kg m-2"))
              varid(30)=temp_id
  
              call check( nf90_def_var(ncid, "canwat", NF90_REAL, dimtwo_time, temp_id), trim(err)//"canwat" )
              call check( nf90_put_att(ncid,temp_id,"standard_name","canopy_water_amount"))
              call check( nf90_put_att(ncid,temp_id,"long_name","Canopy water content"))
              call check( nf90_put_att(ncid,temp_id,"units","kg m-2"))
              varid(31)=temp_id
  
          endif
  
negative.code:                 call check( nf90_put_att(ncid,temp_id,"long_name","Froude Number"))
                  call check( nf90_put_att(ncid,temp_id,"description", "Tendency for topographic blocking to occur = U / (Z * N)"))
                  call check( nf90_put_att(ncid,temp_id,"units","(m s-1) (m s-1)-1"))
                  call check( nf90_put_att(ncid,temp_id,"coordinate","lat lon"))
                  varid(42)=temp_id
              endif
  
  