Deltime.c 3.77 KB
Newer Older
1
2
3
4
/*
  This file is part of CDO. CDO is a collection of Operators to
  manipulate and analyse Climate model Data.

5
  Copyright (C) 2003-2016 Uwe Schulzweida, <uwe.schulzweida AT mpimet.mpg.de>
6
7
8
9
10
11
12
13
14
15
16
17
  See COPYING file for copying and redistribution conditions.

  This program is free software; you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation; version 2 of the License.

  This program is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.
*/

Uwe Schulzweida's avatar
Uwe Schulzweida committed
18
#include <ctype.h>
19

Ralf Mueller's avatar
Ralf Mueller committed
20
#include <cdi.h>
21
#include "cdo.h"
Uwe Schulzweida's avatar
Uwe Schulzweida committed
22
#include "cdo_int.h"
23
24
25
#include "pstream.h"


26
void *Deltime(void *argument)
27
{
28
  int nrecs;
29
  int recID, varID, levelID;
30
  int vdate /*, vtime */;
31
32
33
34
  int copytimestep;
  int gridsize;
  int nmiss;
  int year, month, day;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
35
  int dday, dmon;
36
  double *array = NULL;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
37
  const char *cmons[]={"", "jan","feb","mar","apr","may","jun","jul","aug","sep","oct","nov","dec"};
38
39
40

  cdoInitialize(argument);

41
42
43
44
  bool lcopy = UNCHANGED_RECORD;

  int DELDAY   = cdoOperatorAdd("delday",   0, 0, NULL);
  int DEL29FEB = cdoOperatorAdd("del29feb", 0, 0, NULL);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
45

46
47
  UNUSED(DELDAY);

48
  int operatorID = cdoOperatorID();
Uwe Schulzweida's avatar
Uwe Schulzweida committed
49
50
51
52
53
54
55
56
57
58
59
60

  if ( operatorID == DEL29FEB )
    {
      dday = 29;
      dmon = 2;
    }
  else
    {
      int im;
      int nsel;
      char *sarg;
      nsel = operatorArgc();
Uwe Schulzweida's avatar
Uwe Schulzweida committed
61
      if ( nsel < 1 ) cdoAbort("Too few arguments!");
Uwe Schulzweida's avatar
Uwe Schulzweida committed
62
63
64
65
      if ( nsel > 1 ) cdoAbort("Too many arguments!");
      sarg = operatorArgv()[0];
      dday = atoi(sarg);
      dmon = 0;
66
      while ( isdigit(*sarg) ) sarg++;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
67
68
69
      if ( isalpha(*sarg) )
	{
	  char smon[32];
70
71
	  strncpy(smon, sarg, sizeof(smon)-1);
	  smon[sizeof(smon)-1] = 0;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
72
73
74
75
76
77
78
79
80
81
	  strtolower(smon);
	  for ( im = 0; im < 12; ++im )
	    if ( memcmp(smon, cmons[im+1], 3) == 0 ) break;

	  if ( im < 12 ) dmon = im + 1;
	}
    }

  if ( cdoVerbose ) cdoPrint("delete day %d%s", dday, cmons[dmon]);

82
  int streamID1 = streamOpenRead(cdoStreamName(0));
83

84
85
  int vlistID1 = streamInqVlist(streamID1);
  int vlistID2 = vlistDuplicate(vlistID1);
86

87
88
  int taxisID1 = vlistInqTaxis(vlistID1);
  int taxisID2 = taxisDuplicate(taxisID1);
89
90
91
  taxisDefCalendar(taxisID2, CALENDAR_365DAYS);
  vlistDefTaxis(vlistID2, taxisID2);

92
  int streamID2 = streamOpenWrite(cdoStreamName(1), cdoFiletype());
93
94
95
96
97
98

  streamDefVlist(streamID2, vlistID2);

  if ( ! lcopy )
    {
      gridsize = vlistGridsizeMax(vlistID1);
99
      array = (double*) Malloc(gridsize*sizeof(double));
100
101
    }
      
102
103
104
  int nfound = 0;
  int tsID  = 0;
  int tsID2 = 0;
105
106
107
  while ( (nrecs = streamInqTimestep(streamID1, tsID)) )
    {
      vdate = taxisInqVdate(taxisID1);
108
      // vtime = taxisInqVtime(taxisID1);
109

Uwe Schulzweida's avatar
Uwe Schulzweida committed
110
      cdiDecodeDate(vdate, &year, &month, &day);
111

Uwe Schulzweida's avatar
Uwe Schulzweida committed
112
      if ( day == dday && (month == dmon || dmon == 0) )
113
	{
Uwe Schulzweida's avatar
Uwe Schulzweida committed
114
	  nfound++;
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
	  copytimestep = FALSE;
	  if ( cdoVerbose )
	    cdoPrint("Delete %4.4d-%2.2d-%2.2d at timestep %d", year, month, day, tsID+1);
	}
      else
	copytimestep = TRUE;

      if ( copytimestep )
	{
	  taxisCopyTimestep(taxisID2, taxisID1);

	  streamDefTimestep(streamID2, tsID2++);

	  for ( recID = 0; recID < nrecs; recID++ )
	    {
	      streamInqRecord(streamID1, &varID, &levelID);
	      streamDefRecord(streamID2, varID, levelID);
	      if ( lcopy )
		{
		  streamCopyRecord(streamID2, streamID1);
		}
	      else
		{
		  streamReadRecord(streamID1, array, &nmiss);
		  streamWriteRecord(streamID2, array, nmiss);
		}
	    }
	}
       
      tsID++;
    }

  streamClose(streamID2);
  streamClose(streamID1);

Uwe Schulzweida's avatar
Uwe Schulzweida committed
150
151
152
  if ( nfound == 0 )
    cdoWarning("Day %d%s not found!", dday, cmons[dmon]);

153
  if ( ! lcopy )
154
    if ( array ) Free(array);
155
156
157
158
159

  vlistDestroy(vlistID2);

  cdoFinish();

160
  return 0;
161
}