#include #include #include "protos.h" #include "registry.h" #include "data.h" #include #include #include "sym.h" int gen_namelist_defines ( char * dirname , int sw_dimension ) { FILE * fp ; char fname[NAMELEN] ; char fn[NAMELEN] ; node_t *p ; sprintf( fn, "namelist_defines%s.inc", sw_dimension?"":"2" ) ; if ( strlen(dirname) > 0 ) { sprintf(fname,"%s/%s",dirname,fn) ; } if ((fp = fopen( fname , "w" )) == NULL ) return(1) ; print_warning(fp,fname) ; fprintf(fp,"integer :: first_item_in_struct\n") ; for ( p = Domain.fields ; p != NULL ; p = p-> next ) { if ( p->node_kind & RCONFIG ) { if ( sw_dimension ) { if ( !strcmp( p->nentries, "1" ) ) fprintf(fp,"%s :: %s\n",p->type->name ,p->name) ; else if ( strcmp( p->nentries, "-" ) ) /* if not equal to "-" */ fprintf(fp,"%s , DIMENSION(%s) :: %s\n",p->type->name ,p->nentries,p->name) ; } else { fprintf(fp,"%s :: %s\n",p->type->name ,p->name) ; } } } fprintf(fp,"integer :: last_item_in_struct\n") ; close_the_file( fp ) ; return(0) ; } int gen_namelist_defaults ( char * dirname ) { FILE * fp ; char fname[NAMELEN] ; char *fn = "namelist_defaults.inc" ; node_t *p ; if ( strlen(dirname) > 0 ) { sprintf(fname,"%s/%s",dirname,fn) ; } if ((fp = fopen( fname , "w" )) == NULL ) return(1) ; print_warning(fp,fname) ; for ( p = Domain.fields ; p != NULL ; p = p-> next ) { if ( p->node_kind & RCONFIG && strcmp(p->dflt,"-") && strcmp(p->dflt,"")) { if ( !strncmp ( p->type->name , "character", 9 ) ) { fprintf(fp,"%s = \"%s\"\n",p->name ,p->dflt) ; } else { fprintf(fp,"%s = %s\n",p->name ,p->dflt) ; } } } close_the_file( fp ) ; return(0) ; } int gen_namelist_statements ( char * dirname ) { FILE * fp ; char fname[NAMELEN] ; char * fn = "namelist_statements.inc" ; char howset[NAMELEN] ; char *p1, *p2 ; node_t *p ; strcpy( fname, fn ) ; if ( strlen(dirname) > 0 ) { sprintf(fname,"%s/%s",dirname,fn) ; } if ((fp = fopen( fname , "w" )) == NULL ) return(1) ; print_warning(fp,fname) ; for ( p = Domain.fields ; p != NULL ; p = p-> next ) { if ( p->node_kind & RCONFIG ) { strcpy(howset,p->howset) ; if (( p1 = strtok(howset,",")) != NULL ) { p2 = strtok(NULL,",") ; if ( !strcmp(p1,"namelist") ) { if ( p2 == NULL ) { fprintf(stderr, "Warning: no namelist section specified for nl %s\n",p->name) ; continue ; } fprintf(fp,"NAMELIST /%s/ %s\n",p2,p->name) ; } } } } close_the_file( fp ) ; return(0) ; } int gen_get_nl_config ( char * dirname ) { FILE * fp ; char fname[NAMELEN] ; char * fn = "get_nl_config.inc" ; char * gs, * intnt ; char howset[NAMELEN] ; node_t *p ; int sw ; strcpy( fname, fn ) ; if ( strlen(dirname) > 0 ) { sprintf(fname,"%s/%s",dirname,fn) ; } if ((fp = fopen( fname , "w" )) == NULL ) return(1) ; print_warning(fp,fname) ; for ( sw = 0 ; sw < 2 ; sw++ ) { if ( sw == 0 ) { gs = "get" ; intnt = "OUT" ; } else { gs = "set" ; intnt = "IN" ; } for ( p = Domain.fields ; p != NULL ; p = p-> next ) { if ( p->node_kind & RCONFIG ) { strcpy(howset,p->howset) ; fprintf(fp,"SUBROUTINE nl_%s_%s ( id_id , %s )\n",gs,p->name, p->name) ; if ( sw_ifort_kludge ) { fprintf(fp," USE module_configure\n") ; } fprintf(fp," %s , INTENT(%s) :: %s\n",p->type->name,intnt,p->name) ; fprintf(fp," INTEGER id_id\n") ; fprintf(fp," CHARACTER*80 emess\n") ; if ( sw == 0 ) /* get */ { if ( !strcmp( p->nentries, "1" )) { if ( ! sw_ifort_kludge ) { fprintf(fp," IF ( id_id .NE. 1 ) THEN\n") ; fprintf(fp," call wrf_debug(1,&\n'WARNING in nl_%s_%s: %s applies to all domains. First arg ignored.')\n", gs,p->name, p->name ) ; fprintf(fp," ENDIF\n" ) ; } if ( !strncmp(p->type->name,"character",9)) { fprintf(fp," %s = trim(model_config_rec%%%s)\n",p->name,p->name) ; }else{ fprintf(fp," %s = model_config_rec%%%s\n",p->name,p->name) ; } } else { if ( ! sw_ifort_kludge ) { if ( !strcmp( p->nentries, "max_domains" )) { fprintf(fp," IF ( id_id .LT. 1 .OR. id_id .GT. model_config_rec%%max_dom ) THEN\n") ; fprintf(fp," WRITE(emess,*)'nl_%s_%s: Out of range domain number: ',id_id\n",gs,p->name) ; } else if ( !strcmp( p->nentries, "max_moves" )) { fprintf(fp," IF ( id_id .LT. 1 .OR. id_id .GT. model_config_rec%%num_moves ) THEN\n") ; fprintf(fp," WRITE(emess,*)'nl_%s_%s: Out of range move number: ',id_id\n",gs,p->name) ; } else if ( !strcmp( p->nentries, "max_eta" )) { fprintf(fp," IF ( id_id .LT. 1 .OR. id_id .GT. model_config_rec%%e_vert(1) ) THEN\n") ; fprintf(fp," WRITE(emess,*)'nl_%s_%s: Out of range eta_level number: ',id_id\n",gs,p->name) ; } else { fprintf(stderr,"Registry WARNING: multi element rconfig entry must be either max_domains, max_moves, or max_eta \n") ; } fprintf(fp," CALL wrf_error_fatal(emess)\n") ; fprintf(fp," ENDIF\n" ) ; } fprintf(fp," %s = model_config_rec%%%s(id_id)\n",p->name,p->name) ; } } else /* set */ { if ( !strcmp( p->nentries, "1" )) { if ( ! sw_ifort_kludge ) { fprintf(fp," IF ( id_id .NE. 1 ) THEN\n") ; fprintf(fp," call wrf_debug(1,&\n'WARNING in nl_%s_%s: %s applies to all domains. First arg ignored.')\n", gs,p->name, p->name ) ; fprintf(fp," ENDIF\n" ) ; } if ( !strncmp(p->type->name,"character",9)) { fprintf(fp," model_config_rec%%%s = trim(%s) \n",p->name,p->name) ; }else{ fprintf(fp," model_config_rec%%%s = %s \n",p->name,p->name) ; } } else { if ( ! sw_ifort_kludge ) { if ( !strcmp( p->nentries, "max_domains" )) { fprintf(fp," IF ( id_id .LT. 1 .OR. id_id .GT. model_config_rec%%max_dom ) THEN\n") ; fprintf(fp," WRITE(emess,*)'nl_%s_%s: Out of range domain number: ',id_id\n",gs,p->name) ; } else if ( !strcmp( p->nentries, "max_moves" )) { fprintf(fp," IF ( id_id .LT. 1 .OR. id_id .GT. model_config_rec%%num_moves ) THEN\n") ; fprintf(fp," WRITE(emess,*)'nl_%s_%s: Out of range move number: ',id_id\n",gs,p->name) ; } else if ( !strcmp( p->nentries, "max_eta" )) { fprintf(fp," IF ( id_id .LT. 1 .OR. id_id .GT. model_config_rec%%e_vert(1) ) THEN\n") ; fprintf(fp," WRITE(emess,*)'nl_%s_%s: Out of range eta_level number: ',id_id\n",gs,p->name) ; } else { fprintf(stderr,"Registry WARNING: multi element rconfig entry must be either max_domains, max_moves, or max_eta \n") ; } fprintf(fp," CALL wrf_error_fatal(emess)\n") ; fprintf(fp," ENDIF\n" ) ; } fprintf(fp," model_config_rec%%%s(id_id) = %s\n",p->name,p->name) ; } } fprintf(fp," RETURN\n") ; fprintf(fp,"END SUBROUTINE nl_%s_%s\n",gs,p->name ) ; } } } close_the_file( fp ) ; return(0) ; } int gen_config_assigns ( char * dirname ) { FILE * fp ; char fname[NAMELEN] ; char * fn = "config_assigns.inc" ; char tmp[NAMELEN] ; node_t *p ; strcpy( fname, fn ) ; if ( strlen(dirname) > 0 ) { sprintf(fname,"%s/%s",dirname,fn) ; } if ((fp = fopen( fname , "w" )) == NULL ) return(1) ; print_warning(fp,fname) ; fprintf(fp,"! Contains config assign statements for module_domain.F.\n") ; fprintf(fp,"#ifndef SOURCE_RECORD\n") ; fprintf(fp,"# define SOURCE_RECORD cfg%%\n") ; fprintf(fp,"#endif\n") ; fprintf(fp,"#ifndef SOURCE_REC_DEX\n") ; fprintf(fp,"# define SOURCE_REC_DEX\n") ; fprintf(fp,"#endif\n") ; fprintf(fp,"#ifndef DEST_RECORD\n") ; fprintf(fp,"# define DEST_RECORD new_grid%%\n") ; fprintf(fp,"#endif\n") ; for ( p = Domain.fields ; p != NULL ; p = p-> next ) { if ( p->node_kind & RCONFIG ) { if ( !strcmp( p->nentries, "1" )) strcpy( tmp, "" ) ; else strcpy( tmp, "SOURCE_REC_DEX" ) ; fprintf(fp," DEST_RECORD %-26s = SOURCE_RECORD %s %s\n",p->name,p->name,tmp) ; } } close_the_file( fp ) ; return(0) ; } int gen_config_reads ( char * dirname ) { FILE * fp ; char fname[NAMELEN] ; char * fn = "config_reads.inc" ; char howset[NAMELEN] ; char *p1, *p2 ; node_t *p ; strcpy( fname, fn ) ; if ( strlen(dirname) > 0 ) { sprintf(fname,"%s/%s",dirname,fn) ; } if ((fp = fopen( fname , "w" )) == NULL ) return(1) ; print_warning(fp,fname) ; fprintf(fp,"! Contains namelist statements for module_config.F.\n") ; fprintf(fp,"#ifndef NAMELIST_READ_UNIT\n") ; fprintf(fp,"# define NAMELIST_READ_UNIT nml_unit\n") ; fprintf(fp,"#endif\n") ; fprintf(fp,"#ifndef NAMELIST_READ_ERROR_LABEL\n") ; fprintf(fp,"# define NAMELIST_READ_ERROR_LABEL 9200\n") ; fprintf(fp,"#endif\n") ; fprintf(fp,"!\n") ; sym_forget() ; for ( p = Domain.fields ; p != NULL ; p = p-> next ) { if ( p->node_kind & RCONFIG ) { strcpy(howset,p->howset) ; p1 = strtok(howset,",") ; p2 = strtok(NULL,",") ; if ( !strcmp(p1,"namelist") ) { if ( p2 == NULL ) { fprintf(stderr, "Warning: no namelist section specified for nl %s\n",p->name) ; continue ; } if (sym_get( p2 ) == NULL) /* not in table yet */ { fprintf(fp," REWIND ( UNIT = NAMELIST_READ_UNIT )\n") ; fprintf(fp," READ ( UNIT = NAMELIST_READ_UNIT , NML = %s , ERR = NAMELIST_READ_ERROR_LABEL , END = NAMELIST_READ_ERROR_LABEL )\n",p2) ; fprintf(fp,"#ifndef NO_NAMELIST_PRINT\n") ; fprintf(fp," WRITE ( UNIT = * , NML = %s )\n",p2) ; fprintf(fp,"#endif\n") ; sym_add(p2) ; } } } } close_the_file( fp ) ; return(0) ; }