Changeset 3312 for LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Apr 11, 2018, 10:27:28 AM (6 years ago)

Author:

Laurent Fairhead

Message:

Continuing phasing of DYNAMICO and LMDZ physics

File:

• LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phys_state_var_mod.F90 (modified) (7 diffs)

LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phys_state_var_mod.F90

@@ -279 +279 @@
       REAL,ALLOCATABLE,SAVE :: total_rain(:), nday_rain(:)  
 !$OMP THREADPRIVATE(total_rain,nday_rain)
+      REAL,ALLOCATABLE,SAVE :: paire_ter(:)
+!$OMP THREADPRIVATE(paire_ter)
 ! albsol1: albedo du sol total pour SW visible
 ! albsol2: albedo du sol total pour SW proche IR
@@ -410 +412 @@
       ! tendencies on wind due to gravity waves
 
+      LOGICAL,SAVE :: is_initialized=.FALSE.
+!$OMP THREADPRIVATE(is_initialized)    
+
 CONTAINS
 
@@ -430 +435 @@
 include "clesphys.h"
 
+      IF (is_initialized) RETURN
+      is_initialized=.TRUE.
       ALLOCATE(pctsrf(klon,nbsrf))
       ALLOCATE(ftsol(klon,nbsrf))
@@ -540 +547 @@
       ALLOCATE(pfrac_1nucl(klon,klev))
       ALLOCATE(total_rain(klon), nday_rain(klon))
+      ALLOCATE(paire_ter(klon))
       ALLOCATE(albsol1(klon), albsol2(klon))
 !albedo SB >>>
@@ -595 +603 @@
       ALLOCATE(ale_bl_trig(klon))
 !!! fin nrlmd le 10/04/2012
-      if (ok_gwd_rando) allocate(du_gwd_rando(klon, klev))
+      IF (ok_gwd_rando) THEN
+        allocate(du_gwd_rando(klon, klev))
+        du_gwd_rando(:,:)=0.
+      ENDIF
       if (.not. ok_hines .and. ok_gwd_rando) allocate(du_gwd_front(klon, klev))
 
@@ -680 +691 @@
       deallocate(pfrac_1nucl)
       deallocate(total_rain, nday_rain)
+      deallocate(paire_ter)
       deallocate(albsol1, albsol2)
 !albedo SB >>>
@@ -727 +739 @@
       deallocate(ale_bl_trig)
 !!! fin nrlmd le 10/04/2012
-
+      is_initialized=.FALSE.
 END SUBROUTINE phys_state_var_end