Changeset 519 for palm/trunk/SOURCE/netcdf.f90
- Timestamp:
- Mar 19, 2010 5:30:02 AM (14 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/netcdf.f90 (modified) (38 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/netcdf.f90
r494 r519 7 7 ! Current revisions: 8 8 ! ------------------ 9 ! 9 ! Particle number defined as unlimited dimension in case of NetCDF4 output, 10 ! special characters like * and " are now allowed for NetCDF variable names, 11 ! replacement of these characters removed, routine clean_netcdf_varname 12 ! removed 10 13 ! 11 14 ! Former revisions: … … 95 98 CHARACTER (LEN=4) :: grid_x, grid_y, grid_z 96 99 CHARACTER (LEN=6) :: mode 97 CHARACTER (LEN=10) :: netcdf_var_name, netcdf_var_name_base, & 98 precision, var 100 CHARACTER (LEN=10) :: netcdf_var_name, precision, var 99 101 CHARACTER (LEN=80) :: time_average_text 100 102 CHARACTER (LEN=2000) :: var_list, var_list_old … … 1556 1558 !-- surface variable. Define it with id_dim_zu1_xy. 1557 1559 IF ( INDEX( do2d(av,i), '*' ) /= 0 ) THEN 1558 ! 1559 !-- First, remove those characters not allowed by NetCDF 1560 netcdf_var_name = do2d(av,i) 1561 CALL clean_netcdf_varname( netcdf_var_name ) 1562 1563 nc_stat = NF90_DEF_VAR( id_set_xy(av), netcdf_var_name, & 1560 1561 nc_stat = NF90_DEF_VAR( id_set_xy(av), do2d(av,i), & 1564 1562 nc_precision(1), & 1565 1563 (/ id_dim_x_xy(av), id_dim_y_xy(av),& … … 1568 1566 id_var_do2d(av,i) ) 1569 1567 1570 var_list = TRIM(var_list) // TRIM( netcdf_var_name) // ';'1568 var_list = TRIM(var_list) // TRIM( do2d(av,i) ) // ';' 1571 1569 1572 1570 ELSE … … 1827 1825 DO WHILE ( do2d(av,i)(1:1) /= ' ' ) 1828 1826 IF ( INDEX( do2d(av,i), 'xy' ) /= 0 ) THEN 1829 netcdf_var_name = do2d(av,i) 1830 CALL clean_netcdf_varname( netcdf_var_name ) 1831 var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';' 1827 var_list = TRIM( var_list ) // TRIM( do2d(av,i) ) // ';' 1832 1828 ENDIF 1833 1829 i = i + 1 … … 1967 1963 DO WHILE ( do2d(av,i)(1:1) /= ' ' ) 1968 1964 IF ( INDEX( do2d(av,i), 'xy' ) /= 0 ) THEN 1969 netcdf_var_name = do2d(av,i) 1970 CALL clean_netcdf_varname( netcdf_var_name ) 1971 nc_stat = NF90_INQ_VARID( id_set_xy(av), netcdf_var_name, & 1965 nc_stat = NF90_INQ_VARID( id_set_xy(av), do2d(av,i), & 1972 1966 id_var_do2d(av,i) ) 1973 1967 CALL handle_netcdf_error( 'netcdf', 138 ) … … 2401 2395 DO WHILE ( do2d(av,i)(1:1) /= ' ' ) 2402 2396 IF ( INDEX( do2d(av,i), 'xz' ) /= 0 ) THEN 2403 netcdf_var_name = do2d(av,i) 2404 CALL clean_netcdf_varname( netcdf_var_name ) 2405 var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';' 2397 var_list = TRIM( var_list ) // TRIM( do2d(av,i) ) // ';' 2406 2398 ENDIF 2407 2399 i = i + 1 … … 2541 2533 DO WHILE ( do2d(av,i)(1:1) /= ' ' ) 2542 2534 IF ( INDEX( do2d(av,i), 'xz' ) /= 0 ) THEN 2543 netcdf_var_name = do2d(av,i) 2544 CALL clean_netcdf_varname( netcdf_var_name ) 2545 nc_stat = NF90_INQ_VARID( id_set_xz(av), netcdf_var_name, & 2535 nc_stat = NF90_INQ_VARID( id_set_xz(av), do2d(av,i), & 2546 2536 id_var_do2d(av,i) ) 2547 2537 CALL handle_netcdf_error( 'netcdf', 177 ) … … 2989 2979 DO WHILE ( do2d(av,i)(1:1) /= ' ' ) 2990 2980 IF ( INDEX( do2d(av,i), 'yz' ) /= 0 ) THEN 2991 netcdf_var_name = do2d(av,i) 2992 CALL clean_netcdf_varname( netcdf_var_name ) 2993 var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';' 2981 var_list = TRIM( var_list ) // TRIM( do2d(av,i) ) // ';' 2994 2982 ENDIF 2995 2983 i = i + 1 … … 3129 3117 DO WHILE ( do2d(av,i)(1:1) /= ' ' ) 3130 3118 IF ( INDEX( do2d(av,i), 'yz' ) /= 0 ) THEN 3131 netcdf_var_name = do2d(av,i) 3132 CALL clean_netcdf_varname( netcdf_var_name ) 3133 nc_stat = NF90_INQ_VARID( id_set_yz(av), netcdf_var_name, & 3119 nc_stat = NF90_INQ_VARID( id_set_yz(av), do2d(av,i), & 3134 3120 id_var_do2d(av,i) ) 3135 3121 CALL handle_netcdf_error( 'netcdf', 216 ) … … 3229 3215 var_list = ';' 3230 3216 DO i = 1, dopr_n 3231 !3232 !-- First, remove those characters not allowed by NetCDF3233 netcdf_var_name = data_output_pr(i)3234 CALL clean_netcdf_varname( netcdf_var_name )3235 3217 3236 3218 IF ( statistic_regions == 0 ) THEN … … 3238 3220 ! 3239 3221 !-- Define the z-axes (each variable gets its own z-axis) 3240 nc_stat = NF90_DEF_DIM( id_set_pr, 'z'//TRIM(netcdf_var_name), & 3222 nc_stat = NF90_DEF_DIM( id_set_pr, & 3223 'z' // TRIM( data_output_pr(i) ), & 3241 3224 nzt+2-nzb, id_dim_z_pr(i,0) ) 3242 3225 CALL handle_netcdf_error( 'netcdf', 223 ) 3243 3226 3244 nc_stat = NF90_DEF_VAR( id_set_pr, 'z'//TRIM(netcdf_var_name), & 3245 NF90_DOUBLE, id_dim_z_pr(i,0), & 3227 nc_stat = NF90_DEF_VAR( id_set_pr, & 3228 'z' // TRIM( data_output_pr(i) ), & 3229 NF90_DOUBLE, id_dim_z_pr(i,0), & 3246 3230 id_var_z_pr(i,0) ) 3247 3231 CALL handle_netcdf_error( 'netcdf', 224 ) … … 3252 3236 ! 3253 3237 !-- Define the variable 3254 nc_stat = NF90_DEF_VAR( id_set_pr, netcdf_var_name,&3238 nc_stat = NF90_DEF_VAR( id_set_pr, data_output_pr(i), & 3255 3239 nc_precision(5), (/ id_dim_z_pr(i,0), & 3256 3240 id_dim_time_pr /), id_var_dopr(i,0) ) … … 3265 3249 CALL handle_netcdf_error( 'netcdf', 228 ) 3266 3250 3267 var_list = TRIM( var_list) // TRIM(netcdf_var_name) //';'3251 var_list = TRIM( var_list ) // TRIM( data_output_pr(i) ) // ';' 3268 3252 3269 3253 ELSE … … 3277 3261 !-- Define the z-axes (each variable gets it own z-axis) 3278 3262 nc_stat = NF90_DEF_DIM( id_set_pr, & 3279 'z'//TRIM(netcdf_var_name)//suffix, &3263 'z'//TRIM(data_output_pr(i))//suffix, & 3280 3264 nzt+2-nzb, id_dim_z_pr(i,j) ) 3281 3265 CALL handle_netcdf_error( 'netcdf', 229 ) 3282 3266 3283 3267 nc_stat = NF90_DEF_VAR( id_set_pr, & 3284 'z'//TRIM(netcdf_var_name)//suffix, &3268 'z'//TRIM(data_output_pr(i))//suffix, & 3285 3269 nc_precision(5), id_dim_z_pr(i,j), & 3286 3270 id_var_z_pr(i,j) ) … … 3294 3278 !-- Define the variable 3295 3279 nc_stat = NF90_DEF_VAR( id_set_pr, & 3296 TRIM( netcdf_var_name) // suffix, &3280 TRIM(data_output_pr(i)) // suffix, & 3297 3281 nc_precision(5), & 3298 3282 (/ id_dim_z_pr(i,j), & … … 3310 3294 CALL handle_netcdf_error( 'netcdf', 234 ) 3311 3295 3312 var_list = TRIM( var_list) // TRIM(netcdf_var_name) // &3296 var_list = TRIM( var_list ) // TRIM( data_output_pr(i) ) // & 3313 3297 suffix // ';' 3314 3298 … … 3376 3360 DO i = 1, dopr_n 3377 3361 3378 netcdf_var_name = data_output_pr(i)3379 CALL clean_netcdf_varname( netcdf_var_name )3380 3381 3362 IF ( statistic_regions == 0 ) THEN 3382 var_list = TRIM( var_list) // TRIM(netcdf_var_name) // ';'3363 var_list = TRIM( var_list ) // TRIM( data_output_pr(i) ) // ';' 3383 3364 ELSE 3384 3365 DO j = 0, statistic_regions 3385 3366 WRITE ( suffix, '(''_'',I1)' ) j 3386 var_list = TRIM( var_list) // TRIM(netcdf_var_name) // &3367 var_list = TRIM( var_list ) // TRIM( data_output_pr(i) ) // & 3387 3368 suffix // ';' 3388 3369 ENDDO … … 3445 3426 DO i = 1, dopr_n 3446 3427 3447 netcdf_var_name_base = data_output_pr(i)3448 CALL clean_netcdf_varname( netcdf_var_name_base )3449 3450 3428 IF ( statistic_regions == 0 ) THEN 3451 nc_stat = NF90_INQ_VARID( id_set_pr, netcdf_var_name_base, &3429 nc_stat = NF90_INQ_VARID( id_set_pr, data_output_pr(i), & 3452 3430 id_var_dopr(i,0) ) 3453 3431 CALL handle_netcdf_error( 'netcdf', 245 ) … … 3455 3433 DO j = 0, statistic_regions 3456 3434 WRITE ( suffix, '(''_'',I1)' ) j 3457 netcdf_var_name = TRIM( netcdf_var_name_base) // suffix3435 netcdf_var_name = TRIM( data_output_pr(i) ) // suffix 3458 3436 nc_stat = NF90_INQ_VARID( id_set_pr, netcdf_var_name, & 3459 3437 id_var_dopr(i,j) ) … … 3525 3503 var_list = ';' 3526 3504 DO i = 1, dots_num 3527 !3528 !-- First, remove those characters not allowed by NetCDF3529 netcdf_var_name = dots_label(i)3530 CALL clean_netcdf_varname( netcdf_var_name )3531 3505 3532 3506 IF ( statistic_regions == 0 ) THEN 3533 3507 3534 nc_stat = NF90_DEF_VAR( id_set_ts, netcdf_var_name,&3508 nc_stat = NF90_DEF_VAR( id_set_ts, dots_label(i), & 3535 3509 nc_precision(6), (/ id_dim_time_ts /), & 3536 3510 id_var_dots(i,0) ) … … 3545 3519 CALL handle_netcdf_error( 'netcdf', 255 ) 3546 3520 3547 var_list = TRIM( var_list) // TRIM(netcdf_var_name) // ';'3521 var_list = TRIM( var_list ) // TRIM( dots_label(i) ) // ';' 3548 3522 3549 3523 ELSE … … 3555 3529 3556 3530 nc_stat = NF90_DEF_VAR( id_set_ts, & 3557 TRIM( netcdf_var_name ) // suffix,&3531 TRIM( dots_label(i) ) // suffix, & 3558 3532 nc_precision(6), & 3559 3533 (/ id_dim_time_ts /), & … … 3571 3545 CALL handle_netcdf_error( 'netcdf', 347 ) 3572 3546 3573 var_list = TRIM( var_list) // TRIM(netcdf_var_name) // &3547 var_list = TRIM( var_list ) // TRIM( dots_label(i) ) // & 3574 3548 suffix // ';' 3575 3549 … … 3607 3581 DO i = 1, dots_num 3608 3582 3609 netcdf_var_name = dots_label(i)3610 CALL clean_netcdf_varname( netcdf_var_name )3611 3612 3583 IF ( statistic_regions == 0 ) THEN 3613 var_list = TRIM( var_list) // TRIM(netcdf_var_name) // ';'3584 var_list = TRIM( var_list ) // TRIM( dots_label(i) ) // ';' 3614 3585 ELSE 3615 3586 DO j = 0, statistic_regions 3616 3587 WRITE ( suffix, '(''_'',I1)' ) j 3617 var_list = TRIM( var_list) // TRIM(netcdf_var_name) // &3588 var_list = TRIM( var_list ) // TRIM( dots_label(i) ) // & 3618 3589 suffix // ';' 3619 3590 ENDDO … … 3676 3647 DO i = 1, dots_num 3677 3648 3678 netcdf_var_name_base = dots_label(i)3679 CALL clean_netcdf_varname( netcdf_var_name_base )3680 3681 3649 IF ( statistic_regions == 0 ) THEN 3682 nc_stat = NF90_INQ_VARID( id_set_ts, netcdf_var_name_base, &3650 nc_stat = NF90_INQ_VARID( id_set_ts, dots_label(i), & 3683 3651 id_var_dots(i,0) ) 3684 3652 CALL handle_netcdf_error( 'netcdf', 265 ) … … 3686 3654 DO j = 0, statistic_regions 3687 3655 WRITE ( suffix, '(''_'',I1)' ) j 3688 netcdf_var_name = TRIM( netcdf_var_name_base) // suffix3656 netcdf_var_name = TRIM( dots_label(i) ) // suffix 3689 3657 nc_stat = NF90_INQ_VARID( id_set_ts, netcdf_var_name, & 3690 3658 id_var_dots(i,j) ) … … 3842 3810 !-- Define the variable 3843 3811 netcdf_var_name = TRIM( data_output_sp(i) ) // '_x' 3844 CALL clean_netcdf_varname( netcdf_var_name )3845 3812 IF ( TRIM( grid_z ) == 'zw' ) THEN 3846 3813 nc_stat = NF90_DEF_VAR( id_set_sp, netcdf_var_name, & … … 3873 3840 !-- Define the variable 3874 3841 netcdf_var_name = TRIM( data_output_sp(i) ) // '_y' 3875 CALL clean_netcdf_varname( netcdf_var_name )3876 3842 IF ( TRIM( grid_z ) == 'zw' ) THEN 3877 3843 nc_stat = NF90_DEF_VAR( id_set_sp, netcdf_var_name, & … … 4174 4140 ! 4175 4141 !-- Define particle coordinate (maximum particle number) 4176 nc_stat = NF90_DEF_DIM( id_set_prt, 'particle_number', & 4177 maximum_number_of_particles, id_dim_prtnum ) 4142 IF ( netcdf_data_format < 3 ) THEN 4143 nc_stat = NF90_DEF_DIM( id_set_prt, 'particle_number', & 4144 maximum_number_of_particles, id_dim_prtnum) 4145 ELSE 4146 ! 4147 !-- NetCDF4 allows more than one unlimited dimension 4148 nc_stat = NF90_DEF_DIM( id_set_prt, 'particle_number', & 4149 NF90_UNLIMITED, id_dim_prtnum) 4150 ENDIF 4178 4151 CALL handle_netcdf_error( 'netcdf', 314 ) 4179 4152 … … 4311 4284 DO i = 1, dopts_num 4312 4285 4313 !4314 !-- First, remove those characters not allowed by NetCDF4315 netcdf_var_name = dopts_label(i)4316 CALL clean_netcdf_varname( netcdf_var_name )4317 4318 4286 DO j = 0, number_of_particle_groups 4319 4287 … … 4325 4293 4326 4294 nc_stat = NF90_DEF_VAR( id_set_pts, & 4327 TRIM( netcdf_var_name ) // suffix1,&4295 TRIM( dopts_label(i) ) // suffix1, & 4328 4296 nc_precision(6), & 4329 4297 (/ id_dim_time_pts /), & … … 4347 4315 CALL handle_netcdf_error( 'netcdf', 402 ) 4348 4316 4349 var_list = TRIM( var_list) // TRIM(netcdf_var_name) // &4317 var_list = TRIM( var_list ) // TRIM( dopts_label(i) ) // & 4350 4318 suffix1 // ';' 4351 4319 … … 4385 4353 DO i = 1, dopts_num 4386 4354 4387 netcdf_var_name = dopts_label(i)4388 CALL clean_netcdf_varname( netcdf_var_name )4389 4390 4355 DO j = 0, number_of_particle_groups 4391 4356 … … 4396 4361 ENDIF 4397 4362 4398 var_list = TRIM( var_list) // TRIM(netcdf_var_name) // &4363 var_list = TRIM( var_list ) // TRIM( dopts_label(i) ) // & 4399 4364 suffix1 // ';' 4400 4365 … … 4459 4424 DO i = 1, dopts_num 4460 4425 4461 netcdf_var_name_base = dopts_label(i)4462 CALL clean_netcdf_varname( netcdf_var_name_base )4463 4464 4426 DO j = 0, number_of_particle_groups 4465 4427 … … 4470 4432 ENDIF 4471 4433 4472 netcdf_var_name = TRIM( netcdf_var_name_base) // suffix14434 netcdf_var_name = TRIM( dopts_label(i) ) // suffix1 4473 4435 4474 4436 nc_stat = NF90_INQ_VARID( id_set_pts, netcdf_var_name, & … … 4546 4508 #endif 4547 4509 END SUBROUTINE handle_netcdf_error 4548 4549 4550 4551 SUBROUTINE clean_netcdf_varname( string )4552 #if defined( __netcdf )4553 4554 !------------------------------------------------------------------------------!4555 !4556 ! Description:4557 ! ------------4558 ! Replace those characters in string which are not allowed by NetCDF.4559 !------------------------------------------------------------------------------!4560 4561 USE netcdf_control4562 4563 IMPLICIT NONE4564 4565 CHARACTER (LEN=10), INTENT(INOUT) :: string4566 4567 INTEGER :: i, ic4568 4569 DO i = 1, replace_num4570 DO4571 ic = INDEX( string, replace_char(i) )4572 IF ( ic /= 0 ) THEN4573 string(ic:ic) = replace_by(i)4574 ELSE4575 EXIT4576 ENDIF4577 ENDDO4578 ENDDO4579 4580 #endif4581 END SUBROUTINE clean_netcdf_varname
Note: See TracChangeset
for help on using the changeset viewer.
