Changeset 3274 for palm/trunk/SOURCE/read_restart_data_mod.f90
- Timestamp:
- Sep 24, 2018 3:42:55 PM (6 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/read_restart_data_mod.f90 (modified) (18 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/read_restart_data_mod.f90
r3241 r3274 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Modularization of all bulk cloud physics code components 28 ! 29 ! 3241 2018-09-12 15:02:00Z raasch 27 30 ! unused variables removed 28 31 ! … … 111 114 q_init, ref_state, s_init, sa_init, u_init, ug, v_init, vg 112 115 116 USE bulk_cloud_model_mod, & 117 ONLY: bcm_rrd_global 118 113 119 USE date_and_time_mod, & 114 120 ONLY: day_of_year_init, time_utc_init … … 125 131 USE indices, & 126 132 ONLY: nz, nx, nx_on_file, ny, ny_on_file 127 128 USE microphysics_mod, &129 ONLY: c_sedimentation, collision_turbulence, &130 cloud_water_sedimentation, limiter_sedimentation, &131 nc_const, ventilation_effect132 133 133 134 USE model_1d_mod, & … … 332 333 CASE ( 'building_wall_south' ) 333 334 READ ( 13 ) building_wall_south 334 CASE ( 'c_sedimentation' )335 READ ( 13 ) c_sedimentation336 335 CASE ( 'call_psolver_at_all_substeps' ) 337 336 READ ( 13 ) call_psolver_at_all_substeps … … 350 349 CASE ( 'cloud_droplets' ) 351 350 READ ( 13 ) cloud_droplets 352 CASE ( 'cloud_physics' )353 READ ( 13 ) cloud_physics354 CASE ( 'cloud_scheme' )355 READ ( 13 ) cloud_scheme356 351 CASE ( 'cloud_top_radiation' ) 357 352 READ ( 13 ) cloud_top_radiation 358 CASE ( 'cloud_water_sedimentation' )359 READ ( 13 ) cloud_water_sedimentation360 353 CASE ( 'collective_wait' ) 361 354 READ ( 13 ) collective_wait 362 CASE ( 'collision_turbulence' )363 READ ( 13 ) collision_turbulence364 355 CASE ( 'conserve_volume_flow' ) 365 356 READ ( 13 ) conserve_volume_flow … … 461 452 CASE ( 'latitude' ) 462 453 READ ( 13 ) latitude 463 CASE ( 'limiter_sedimentation' )464 READ ( 13 ) limiter_sedimentation465 454 CASE ( 'longitude' ) 466 455 READ ( 13 ) longitude … … 484 473 CASE ( 'most_method' ) 485 474 READ ( 13 ) most_method 486 CASE ( 'nc_const' )487 READ ( 13 ) nc_const488 475 CASE ( 'netcdf_precision' ) 489 476 READ ( 13 ) netcdf_precision … … 520 507 CASE ( 'prandtl_number' ) 521 508 READ ( 13 ) prandtl_number 522 CASE ( 'precipitation' )523 READ ( 13 ) precipitation524 509 CASE ( 'psolver' ) 525 510 READ ( 13 ) psolver … … 732 717 CASE ( 'v_max_ijk' ) 733 718 READ ( 13 ) v_max_ijk 734 CASE ( 'ventilation_effect' )735 READ ( 13 ) ventilation_effect736 719 CASE ( 'vg' ) 737 720 READ ( 13 ) vg … … 786 769 CASE DEFAULT 787 770 771 IF ( .NOT. found ) CALL bcm_rrd_global( found ) 772 788 773 IF ( .NOT. found ) CALL wtm_rrd_global( found ) 789 774 … … 1061 1046 1062 1047 USE arrays_3d, & 1063 ONLY: e, kh, km, p, pt, q, ql, qc, nc, nr, prr, precipitation_amount,&1064 qr,s, sa, u, u_m_l, u_m_n, u_m_r, u_m_s, v, v_m_l, v_m_n, &1048 ONLY: e, kh, km, p, pt, q, ql, & 1049 s, sa, u, u_m_l, u_m_n, u_m_r, u_m_s, v, v_m_l, v_m_n, & 1065 1050 v_m_r, v_m_s, vpt, w, w_m_l, w_m_n, w_m_r, w_m_s 1066 1051 … … 1073 1058 ONLY: cpu_log, log_point_s 1074 1059 1075 USE gust_mod, &1060 USE gust_mod, & 1076 1061 ONLY : gust_rrd_local 1077 1062 … … 1084 1069 USE land_surface_model_mod, & 1085 1070 ONLY: lsm_rrd_local 1071 1072 USE bulk_cloud_model_mod, & 1073 ONLY : bcm_rrd_local 1086 1074 1087 1075 USE particle_attributes, & … … 1469 1457 tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1470 1458 1471 CASE ( 'nc' )1472 IF ( k == 1 ) READ ( 13 ) tmp_3d1473 nc(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1474 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1475 1476 CASE ( 'nc_av' )1477 IF ( .NOT. ALLOCATED( nc_av ) ) THEN1478 ALLOCATE( nc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )1479 ENDIF1480 IF ( k == 1 ) READ ( 13 ) tmp_3d1481 nc_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1482 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1483 1484 1485 CASE ( 'nr' )1486 IF ( k == 1 ) READ ( 13 ) tmp_3d1487 nr(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1488 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1489 1490 CASE ( 'nr_av' )1491 IF ( .NOT. ALLOCATED( nr_av ) ) THEN1492 ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )1493 ENDIF1494 IF ( k == 1 ) READ ( 13 ) tmp_3d1495 nr_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1496 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1497 1498 1459 CASE ( 'p' ) 1499 1460 IF ( k == 1 ) READ ( 13 ) tmp_3d … … 1525 1486 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1526 1487 1527 CASE ( 'prr' )1528 IF ( .NOT. ALLOCATED( prr ) ) THEN1529 ALLOCATE( prr(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )1530 ENDIF1531 IF ( k == 1 ) READ ( 13 ) tmp_3d1532 prr(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1533 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1534 1535 CASE ( 'prr_av' )1536 IF ( .NOT. ALLOCATED( prr_av ) ) THEN1537 ALLOCATE( prr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )1538 ENDIF1539 IF ( k == 1 ) READ ( 13 ) tmp_3d1540 prr_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1541 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1542 1543 CASE ( 'precipitation_amount' )1544 IF ( k == 1 ) READ ( 13 ) tmp_2d1545 precipitation_amount(nysc-nbgp:nync+nbgp, &1546 nxlc-nbgp:nxrc+nbgp) = &1547 tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1548 1549 1488 CASE ( 'pt' ) 1550 1489 IF ( k == 1 ) READ ( 13 ) tmp_3d … … 1573 1512 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1574 1513 1575 CASE ( 'qc' )1576 IF ( k == 1 ) READ ( 13 ) tmp_3d1577 qc(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1578 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1579 1580 CASE ( 'qc_av' )1581 IF ( .NOT. ALLOCATED( qc_av ) ) THEN1582 ALLOCATE( qc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )1583 ENDIF1584 IF ( k == 1 ) READ ( 13 ) tmp_3d1585 qc_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1586 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1587 1588 1514 CASE ( 'ql' ) 1589 1515 IF ( k == 1 ) READ ( 13 ) tmp_3d … … 1597 1523 IF ( k == 1 ) READ ( 13 ) tmp_3d 1598 1524 ql_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = & 1599 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1600 1601 CASE ( 'qr' )1602 IF ( k == 1 ) READ ( 13 ) tmp_3d1603 qr(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1604 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)1605 1606 CASE ( 'qr_av' )1607 IF ( .NOT. ALLOCATED( qr_av ) ) THEN1608 ALLOCATE( qr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )1609 ENDIF1610 IF ( k == 1 ) READ ( 13 ) tmp_3d1611 qr_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &1612 1525 tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp) 1613 1526 … … 1968 1881 1969 1882 ! 1883 !-- Read microphysics module restart data 1884 IF ( .NOT. found ) CALL bcm_rrd_local( i, k, nxlf, & 1885 nxlc, nxl_on_file, nxrf, nxrc, & 1886 nxr_on_file, nynf, nync, & 1887 nyn_on_file, nysf, nysc, & 1888 nys_on_file, tmp_2d, tmp_3d, found ) 1889 1890 ! 1970 1891 !-- Read surface related variables 1971 1892 IF ( .NOT. found ) CALL surface_rrd_local( i, k, nxlf, &
Note: See TracChangeset
for help on using the changeset viewer.
