Custom Query (964 matches)

Several Matlab scripts were added:

• 4dvar/obs_extract.m: Extracts data from input data assimilation observations NetCDF file at the requested time interval (days) and creates new observations NetCDF files at such intervals. This script may be used to generate input observation NetCDF files for the Ensemble Kalman Filter (EnKF) from 4D-Var observations files. That is, they can be used for 3-dimensional First Guess at Appropriate Time (3D-FGAT) data assimilation with EnKF.

• utility/daynum.m: Converts requested date (year, month, day, ...) into a serial date number. It uses the Proleptic Gregorian Calendar, which extends backward the date preceding 15 October 1582 with a year length of 365.2425 days. It is similar to Matlab function datenum but differs on the origin (day=0), which corresponds to Mar 1, 0000.

• utility/dayvec.m: Converts a given date number as computed by daynum to a date vector (year, month, day, hour, minutes, seconds). It is the inverse function to daynum. It is equivalent the native datevec function.

• utility/inferno.m: A 256 colormap by Nathaniel J. Smith and Stefan van der Walt.

• utility/magma.m: A 256 colormap by Nathaniel J. Smith and Stefan van der Walt.

• utility/plasma.m: A 256 colormap by Nathaniel J. Smith and Stefan van der Walt.

• utility/viridis.m: A 256 colormap by Eric Firing.

• utility/yearday.m: Given (year, month, day) this function computes the day of the year.

I also updated 4dvar/obs_read.m, landmask/uvp_masks.m, and netcdf/roms_metadata.m scripts.

Try the new colormap palettes. Maybe some of you have used these palettes before.

• Corrected couple of cut/paste typos in ROMS/Functionals/ana_cloud.h and ROMS/Utility/read_asspar.F.

• Updated routine datevec in dateclock.F to make sure that seconds are positive:

  !
  !  Convert fraction of a day.
  !
        seconds=DayFraction*86400.0_r8
        CT=3.0_r8*EPSILON(seconds)             ! comparison tolerance
        seconds=ROUND(seconds, CT)             ! tolerant round function
  !
        hour=INT(seconds/3600.0_r8)
        seconds=ABS(seconds-REAL(hour*3600,r8))
        minutes=INT(seconds/60.0_r8)
        seconds=ABS(seconds-REAL(minutes*60,r8))
  

• Updated Matlab script dayvec.m so it exactly as ROMS routine datevec in module dateclock.F. Mow, it returs more information fields to the output structure. For example:

  % dayvec(daynum(2000,2,28,23,58,45))
  
  ans = 
  
         daynum: 730483.999131944
      leap_year: 'true'
           yday: 59
           year: 2000
          month: 2
            day: 28
           hour: 23
        minutes: 58
        seconds: 45
         string: '2000-02-28 23:58:45.00'
  

• Added new function tround.m to round floating-point numbers using Fuzzy or Tolerant Floor function.