Changeset 3435 for LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90
- Timestamp:
- Jan 22, 2019, 4:21:59 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90
r3208 r3435 18 18 INTEGER, SAVE :: radpas ! radiation is called every "radpas" step 19 19 INTEGER, SAVE :: cvpas ! convection is called every "cvpas" step 20 INTEGER, SAVE :: cvpas_0 ! reference value for cvpas20 INTEGER, SAVE :: cvpas_0 = 1 ! reference value for cvpas 21 21 INTEGER, SAVE :: wkpas ! wake scheme is called every "wkpas" step 22 22 REAL, PARAMETER :: missing_val_nf90=nf90_fill_real … … 25 25 !$OMP THREADPRIVATE(cvpas_0) 26 26 !$OMP THREADPRIVATE(wkpas) 27 REAL, SAVE :: dtime, solaire_etat028 !$OMP THREADPRIVATE( dtime, solaire_etat0)27 REAL, SAVE :: phys_tstep=0, solaire_etat0 28 !$OMP THREADPRIVATE(phys_tstep, solaire_etat0) 29 29 30 30 REAL, ALLOCATABLE, SAVE :: pctsrf(:,:) … … 286 286 REAL,ALLOCATABLE,SAVE :: total_rain(:), nday_rain(:) 287 287 !$OMP THREADPRIVATE(total_rain,nday_rain) 288 REAL,ALLOCATABLE,SAVE :: paire_ter(:) 289 !$OMP THREADPRIVATE(paire_ter) 288 290 ! albsol1: albedo du sol total pour SW visible 289 291 ! albsol2: albedo du sol total pour SW proche IR … … 417 419 ! tendencies on wind due to gravity waves 418 420 421 LOGICAL,SAVE :: is_initialized=.FALSE. 422 !$OMP THREADPRIVATE(is_initialized) 423 419 424 CONTAINS 420 425 … … 437 442 include "clesphys.h" 438 443 444 IF (is_initialized) RETURN 445 is_initialized=.TRUE. 439 446 ALLOCATE(pctsrf(klon,nbsrf)) 440 447 ALLOCATE(ftsol(klon,nbsrf)) … … 452 459 ALLOCATE(snow_fall(klon)) 453 460 ALLOCATE(solsw(klon), sollw(klon)) 461 sollw=0.0 454 462 ALLOCATE(radsol(klon)) 455 463 ALLOCATE(swradcorr(klon)) … … 541 549 ALLOCATE(wake_deltat(klon,klev), wake_deltaq(klon,klev)) 542 550 ALLOCATE(wake_s(klon), awake_dens(klon), wake_dens(klon)) 551 awake_dens = 0. 543 552 ALLOCATE(wake_Cstar(klon)) 544 553 ALLOCATE(wake_pe(klon), wake_fip(klon)) … … 549 558 ALLOCATE(pfrac_1nucl(klon,klev)) 550 559 ALLOCATE(total_rain(klon), nday_rain(klon)) 560 ALLOCATE(paire_ter(klon)) 551 561 ALLOCATE(albsol1(klon), albsol2(klon)) 552 562 !albedo SB >>> … … 566 576 ALLOCATE(topsw(klon), toplw(klon)) 567 577 ALLOCATE(sollwdown(klon), sollwdownclr(klon)) 578 sollwdown = 0. 568 579 ALLOCATE(toplwdown(klon), toplwdownclr(klon)) 569 580 ALLOCATE(topsw0(klon),toplw0(klon),solsw0(klon),sollw0(klon)) 581 sollw0 = 0. 570 582 ALLOCATE(albpla(klon)) 571 583 !IM ajout variables CFMIP2/CMIP5 … … 604 616 ALLOCATE(ale_bl_trig(klon)) 605 617 !!! fin nrlmd le 10/04/2012 606 if (ok_gwd_rando) allocate(du_gwd_rando(klon, klev)) 607 if (.not. ok_hines .and. ok_gwd_rando) allocate(du_gwd_front(klon, klev)) 608 618 IF (ok_gwd_rando) THEN 619 allocate(du_gwd_rando(klon, klev)) 620 du_gwd_rando(:,:)=0. 621 ENDIF 622 IF (.not. ok_hines .and. ok_gwd_rando) THEN 623 ALLOCATE(du_gwd_front(klon, klev)) 624 du_gwd_front(:,:) = 0 !ym missing init 625 ENDIF 609 626 END SUBROUTINE phys_state_var_init 610 627 … … 691 708 deallocate(pfrac_1nucl) 692 709 deallocate(total_rain, nday_rain) 710 deallocate(paire_ter) 693 711 deallocate(albsol1, albsol2) 694 712 !albedo SB >>> … … 738 756 deallocate(ale_bl_trig) 739 757 !!! fin nrlmd le 10/04/2012 740 758 is_initialized=.FALSE. 741 759 END SUBROUTINE phys_state_var_end 742 760
Note: See TracChangeset
for help on using the changeset viewer.