Write RIC file
Expected behavior
Output the trajectory in the RIC format.
Initial state
Variation: "GMT"
Input variables
| Variable name | Description |
|---|---|
| FLOAT.floatname | The floatname that is defined in the rfb file. |
| FLOAT.currentcycle | The first digit that is defined in the FLOAT.cycle variable in the rfb file. |
| time | rafos julian day |
| lat | latitude array |
| lon | longitude array |
| temperature | temperature array |
| pressure | pressure array |
Output parts
All characters are converted to lowercase!
| Shortname | Description | Value/Format | Default value |
|---|---|---|---|
| [PREFIX] | Header/Prefix of every line | FLOAT.floatname - cycle FLOAT.currentcycle | |
| [TIME] | 10 character char/number | YYMMDDHHmm | |
| [LAT] | Latitude | 7.3f | -99.0 |
| [LON] | Longitude | 8.3f | -999.0 |
| [TEMP] | Temperature | 7.2f | -999.0 |
| [PRES] | Pressure | 7.2f | -99.0 |
| [TIMEINCLUDINGOFFSET] | Input time | time + 2440000 |
Output of one line
[PREFIX] [TIME][LAT][LON] 0 -999.000 -999.000 [TEMP] [PRES] [TIMEINCLUDINGOFFSET]
Matlab code for conversion of time of one line (extracted from previous version of ARTOA3):
fdate=rd2dmy(time(ix));%fdate=[day,month,year]
if(fdate(3) >= 2000)
fdate(3)= fdate(3)-2000;
else
fdate(3)= fdate(3)-1900;
end
year=fdate(3);
mon=fdate(2);
day=fdate(1);
hour=24*(time(ix)-floor(time(ix)))+0.0000000001;
min=round(hour-floor(hour));
sec=round(min-floor(min));
hour=floor(hour);
min=round(min);
Edited by Olaf Boebel