Changeset 2706 for palm/trunk/SOURCE/surface_mod.f90
- Timestamp:
- Dec 18, 2017 6:33:49 PM (6 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/surface_mod.f90 (modified) (16 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_mod.f90
r2698 r2706 21 21 ! Current revisions: 22 22 ! ------------------ 23 ! 23 ! In case of restarts read and write pt_surface 24 24 ! 25 25 ! Former revisions: … … 2235 2235 IF ( ALLOCATED( surf_def_h(l)%rib ) ) & 2236 2236 surf_h(l)%rib(mm(l)) = surf_def_h(l)%rib(m) 2237 IF ( ALLOCATED( surf_def_h(l)%pt_surface ) ) & 2238 surf_h(l)%pt_surface(mm(l)) = surf_def_h(l)%pt_surface(m) 2237 2239 IF ( ALLOCATED( surf_def_h(l)%usws ) ) & 2238 2240 surf_h(l)%usws(mm(l)) = surf_def_h(l)%usws(m) … … 2290 2292 IF ( ALLOCATED( surf_lsm_h%rib ) ) & 2291 2293 surf_h(0)%rib(mm(0)) = surf_lsm_h%rib(m) 2294 IF ( ALLOCATED( surf_lsm_h%pt_surface ) ) & 2295 surf_h(l)%pt_surface(mm(l)) = surf_lsm_h%pt_surface(m) 2292 2296 IF ( ALLOCATED( surf_lsm_h%usws ) ) & 2293 2297 surf_h(0)%usws(mm(0)) = surf_lsm_h%usws(m) … … 2345 2349 IF ( ALLOCATED( surf_usm_h%rib ) ) & 2346 2350 surf_h(0)%rib(mm(0)) = surf_usm_h%rib(m) 2351 IF ( ALLOCATED( surf_usm_h%pt_surface ) ) & 2352 surf_h(l)%pt_surface(mm(l)) = surf_usm_h%pt_surface(m) 2347 2353 IF ( ALLOCATED( surf_usm_h%usws ) ) & 2348 2354 surf_h(0)%usws(mm(0)) = surf_usm_h%usws(m) … … 2436 2442 IF ( ALLOCATED( surf_def_v(l)%rib ) ) & 2437 2443 surf_v(l)%rib(mm(l)) = surf_def_v(l)%rib(m) 2444 IF ( ALLOCATED( surf_def_v(l)%pt_surface ) ) & 2445 surf_v(l)%pt_surface(mm(l)) = surf_def_v(l)%pt_surface(m) 2438 2446 IF ( ALLOCATED( surf_def_v(l)%shf ) ) & 2439 2447 surf_v(l)%shf(mm(l)) = surf_def_v(l)%shf(m) … … 2492 2500 IF ( ALLOCATED( surf_lsm_v(l)%rib ) ) & 2493 2501 surf_v(l)%rib(mm(l)) = surf_lsm_v(l)%rib(m) 2502 IF ( ALLOCATED( surf_lsm_v(l)%pt_surface ) ) & 2503 surf_v(l)%pt_surface(mm(l)) = surf_lsm_v(l)%pt_surface(m) 2494 2504 IF ( ALLOCATED( surf_lsm_v(l)%usws ) ) & 2495 2505 surf_v(l)%usws(mm(l)) = surf_lsm_v(l)%usws(m) … … 2552 2562 IF ( ALLOCATED( surf_usm_v(l)%rib ) ) & 2553 2563 surf_v(l)%rib(mm(l)) = surf_usm_v(l)%rib(m) 2564 IF ( ALLOCATED( surf_usm_v(l)%pt_surface ) ) & 2565 surf_v(l)%pt_surface(mm(l)) = surf_usm_v(l)%pt_surface(m) 2554 2566 IF ( ALLOCATED( surf_usm_v(l)%usws ) ) & 2555 2567 surf_v(l)%usws(mm(l)) = surf_usm_v(l)%usws(m) … … 2657 2669 WRITE ( 14 ) surf_h(l)%rib 2658 2670 ENDIF 2671 WRITE ( 14 ) 'surf_h(' // dum // ')%pt_surface ' 2672 IF ( ALLOCATED ( surf_h(l)%pt_surface ) ) THEN 2673 WRITE ( 14 ) surf_h(l)%pt_surface 2674 ENDIF 2659 2675 WRITE ( 14 ) 'surf_h(' // dum // ')%usws ' 2660 2676 IF ( ALLOCATED ( surf_h(l)%usws ) ) THEN … … 2757 2773 IF ( ALLOCATED ( surf_v(l)%rib ) ) THEN 2758 2774 WRITE ( 14 ) surf_v(l)%rib 2775 ENDIF 2776 WRITE ( 14 ) 'surf_v(' // dum // ')%pt_surface ' 2777 IF ( ALLOCATED ( surf_v(l)%pt_surface ) ) THEN 2778 WRITE ( 14 ) surf_v(l)%pt_surface 2759 2779 ENDIF 2760 2780 WRITE ( 14 ) 'surf_v(' // dum // ')%shf ' … … 2997 3017 IF ( ALLOCATED( surf_h(0)%rib ) .AND. kk == 1 ) & 2998 3018 READ ( 13 ) surf_h(0)%rib 3019 CASE ( 'surf_h(0)%pt_surface' ) 3020 IF ( ALLOCATED( surf_h(0)%pt_surface ) .AND. kk == 1 ) & 3021 READ ( 13 ) surf_h(0)%pt_surface 2999 3022 CASE ( 'surf_h(0)%usws' ) 3000 3023 IF ( ALLOCATED( surf_h(0)%usws ) .AND. kk == 1 ) & … … 3073 3096 IF ( ALLOCATED( surf_h(1)%rib ) .AND. kk == 1 ) & 3074 3097 READ ( 13 ) surf_h(1)%rib 3098 CASE ( 'surf_h(1)%pt_surface' ) 3099 IF ( ALLOCATED( surf_h(1)%pt_surface ) .AND. kk == 1 ) & 3100 READ ( 13 ) surf_h(1)%pt_surface 3075 3101 CASE ( 'surf_h(1)%usws' ) 3076 3102 IF ( ALLOCATED( surf_h(1)%usws ) .AND. kk == 1 ) & … … 3227 3253 IF ( ALLOCATED( surf_v(0)%rib ) .AND. kk == 1 ) & 3228 3254 READ ( 13 ) surf_v(0)%rib 3255 CASE ( 'surf_v(0)%pt_surface' ) 3256 IF ( ALLOCATED( surf_v(0)%pt_surface ) .AND. kk == 1 ) & 3257 READ ( 13 ) surf_v(0)%pt_surface 3229 3258 CASE ( 'surf_v(0)%shf' ) 3230 3259 IF ( ALLOCATED( surf_v(0)%shf ) .AND. kk == 1 ) & … … 3306 3335 IF ( ALLOCATED( surf_v(1)%rib ) .AND. kk == 1 ) & 3307 3336 READ ( 13 ) surf_v(1)%rib 3337 CASE ( 'surf_v(1)%pt_surface' ) 3338 IF ( ALLOCATED( surf_v(1)%pt_surface ) .AND. kk == 1 ) & 3339 READ ( 13 ) surf_v(1)%pt_surface 3308 3340 CASE ( 'surf_v(1)%shf' ) 3309 3341 IF ( ALLOCATED( surf_v(1)%shf ) .AND. kk == 1 ) & … … 3385 3417 IF ( ALLOCATED( surf_v(2)%rib ) .AND. kk == 1 ) & 3386 3418 READ ( 13 ) surf_v(2)%rib 3419 CASE ( 'surf_v(2)%pt_surface' ) 3420 IF ( ALLOCATED( surf_v(2)%pt_surface ) .AND. kk == 1 ) & 3421 READ ( 13 ) surf_v(2)%pt_surface 3387 3422 CASE ( 'surf_v(2)%shf' ) 3388 3423 IF ( ALLOCATED( surf_v(2)%shf ) .AND. kk == 1 ) & … … 3464 3499 IF ( ALLOCATED( surf_v(3)%rib ) .AND. kk == 1 ) & 3465 3500 READ ( 13 ) surf_v(3)%rib 3501 CASE ( 'surf_v(3)%pt_surface' ) 3502 IF ( ALLOCATED( surf_v(3)%pt_surface ) .AND. kk == 1 ) & 3503 READ ( 13 ) surf_v(3)%pt_surface 3466 3504 CASE ( 'surf_v(3)%shf' ) 3467 3505 IF ( ALLOCATED( surf_v(3)%shf ) .AND. kk == 1 ) & … … 3646 3684 ENDIF 3647 3685 3686 IF ( INDEX( TRIM( field_chr ), '%pt_surface' ) /= 0 ) THEN 3687 IF ( ALLOCATED( surf_target%pt_surface ) .AND. & 3688 ALLOCATED( surf_file%pt_surface ) ) & 3689 surf_target%pt_surface(m_target) = surf_file%pt_surface(m_file) 3690 ENDIF 3691 3648 3692 IF ( INDEX( TRIM( field_chr ), '%usws' ) /= 0 ) THEN 3649 3693 IF ( ALLOCATED( surf_target%usws ) .AND. &
Note: See TracChangeset
for help on using the changeset viewer.
