source: lmdz_wrf/WRFV3/tools/data.c @ 1

Last change on this file since 1 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 4.3 KB
Line 
1#include <stdio.h>
2#include <stdlib.h>
3#include <string.h>
4#ifdef _WIN32
5#define rindex(X,Y) strrchr(X,Y)
6#define index(X,Y) strchr(X,Y)
7#define bzero(X,Y) memset(X,0,Y)
8#else
9#  include <strings.h>
10#endif
11
12#include "registry.h"
13#include "protos.h"
14#include "data.h"
15
16#define MAXTOKENS 30
17
18int
19init_dim_table()
20{
21  Dim = NULL ;
22  return(0) ;
23}
24
25node_t * 
26new_node ( int kind )
27{ node_t *p ; p = (node_t *)malloc(sizeof(node_t)) ; bzero(p,sizeof(node_t)); p->node_kind = kind ; return (p) ; }
28
29int
30add_node_to_end ( node_t * node , node_t ** list )
31{
32  node_t * p ;
33  if ( *list == NULL ) 
34    { *list = node ; }
35  else
36  {
37    for ( p = *list ; p->next != NULL ; p = p->next ) ;
38    p->next = node ;
39  }
40  return(0) ;
41}
42
43int
44add_node_to_end_4d ( node_t * node , node_t ** list )
45{
46  node_t * p ;
47  if ( *list == NULL ) 
48    { *list = node ; }
49  else
50  {
51    for ( p = *list ; p->next4d != NULL ; p = p->next4d ) ;
52    p->next4d = node ;
53  }
54  return(0) ;
55}
56
57#if 0
58int
59show_nodelist( node_t * p )
60{
61  show_nodelist1( p , 0 ) ;
62}
63
64show_nodelist1( node_t * p , int indent )
65{
66  if ( p == NULL ) return(0) ;
67  show_node1( p, indent) ;
68  show_nodelist1( p->next, indent ) ;
69}
70
71int
72show_node( node_t * p )
73{
74  return(show_node1(p,0)) ;
75}
76
77int
78show_node1( node_t * p, int indent )
79{
80  char spaces[] = "                           " ;
81  char tmp[25] , t1[25] , t2[25] ;
82  char * x, *ca, *ld, *ss, *se, *sg ;
83  char *nodekind ;
84  node_t * q ;
85  int nl ;
86  int i ;
87
88  if ( p == NULL ) return(1) ;
89  strcpy(tmp, spaces) ;
90  if ( indent >= 0 && indent < 20 ) tmp[indent] = '\0' ;
91
92// this doesn't make much sense any more, ever since node_kind was
93// changed to a bit mask
94  if      ( p->node_kind & RCONFIG ) nodekind = "RCONFIG" ;
95  else if ( p->node_kind & I1      ) nodekind = "I1" ;
96  else if ( p->node_kind & FIELD   ) nodekind = "FIELD" ;
97  else if ( p->node_kind & FOURD   ) nodekind = "FOURD" ;
98  else if ( p->node_kind & MEMBER  ) nodekind = "MEMBER" ;
99  else if ( p->node_kind & RCONFIG ) nodekind = "RCONFIG" ;
100
101  if ( !p->scalar_array_member )
102  {
103  switch ( p->node_kind )
104  {
105  case RCONFIG :
106  case I1      :
107  case FIELD   :
108  case FOURD   :
109  case MEMBER  :
110    fprintf(stderr,"%s%s : %10s ndims %1d\n",tmp,nodekind,p->name, p->ndims) ;
111    for ( i = 0 ; i < p->ndims ; i++ )
112    {
113      sg = "" ;
114      switch ( p->dims[i]->coord_axis ) {
115        case COORD_X : ca = "X" ; if ( p->stag_x ) sg = "*" ; break ;
116        case COORD_Y : ca = "Y" ; if ( p->stag_y ) sg = "*" ; break ;
117        case COORD_Z : ca = "Z" ; if ( p->stag_z ) sg = "*" ; break ;
118        case COORD_C : ca = "C" ; break ;
119      }
120      switch ( p->dims[i]->len_defined_how ) {
121        case DOMAIN_STANDARD : ld = "STANDARD" ;  ss = "" ; se = "" ; break ;
122        case NAMELIST        : ld = "NAMELIST" ;  ss = p->dims[i]->associated_namelist_variable ; se="" ; break ;
123        case CONSTANT        : ld = "CONSTANT" ;  sprintf(t1,"%d",p->dims[i]->coord_start) ; ss = t1 ;
124                                                  sprintf(t2,"%d",p->dims[i]->coord_end  ) ; se = t2 ;
125                                                  break ;
126      }
127      fprintf(stderr,"      dim %1d: %c %2s%s %10s %10s %10s\n",i,p->dims[i]->dim_name,ca,sg,ld,ss,se) ;
128    }
129    nl = 0 ;
130    if ( strlen( p->use     ) > 0 ) {
131       nl = 1 ; fprintf(stderr,"      use: %s",p->use) ;
132       if ( p->scalar_array_member ) fprintf(stderr,"(4D)") ;
133    }
134    if ( strlen( p->dname   ) > 0 ) { nl = 1 ; fprintf(stderr,"    dname: %s",p->dname) ;    }
135    if ( strlen( p->descrip ) > 0 ) { nl = 1 ; fprintf(stderr,"  descrip: %s",p->descrip) ;    }
136    if ( nl == 1 ) fprintf(stderr,"\n") ;
137    show_node1( p->type, indent+1 ) ;
138    break ;
139  case TYPE  :
140    x = "derived" ;
141    if ( p->type_type == SIMPLE ) x = "simple" ;
142    fprintf(stderr,"%sTYPE : %10s %s ndims %1d\n",tmp,p->name,x, p->ndims) ;
143    show_nodelist1( p->fields, indent+1 ) ;
144    break ;
145  case DIM   :
146    break ;
147  default :
148    break ;
149  }
150  }
151  show_nodelist1( p->members , indent+2 ) ;
152  return(0) ;
153}
154#endif
155
156int
157set_mark ( int val , node_t * lst )
158{
159  node_t * p ;
160  if ( lst == NULL ) return(0) ;
161  for ( p = lst ; p != NULL ; p = p->next )
162  {
163    p->mark = val ;
164    set_mark( val , p->fields ) ;
165    set_mark( val , p->members ) ;
166  }
167  return(0) ;
168}
169
170int
171set_mark_4d ( int val , node_t * lst )
172{
173  node_t * p ;
174  if ( lst == NULL ) return(0) ;
175  for ( p = lst ; p != NULL ; p = p->next4d )
176  {
177    p->mark = val ;
178    set_mark( val , p->fields ) ;
179    set_mark( val , p->members ) ;
180  }
181  return(0) ;
182}
183
Note: See TracBrowser for help on using the repository browser.