Changeset 3956 for LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Jul 6, 2021, 9:16:14 AM (3 years ago)

Author:

jyg

Message:

Bug fixes concerning various variables ill-initialized, ill-used, ill-printed, or ill-placed.
+ cv_gen moved from phys_local_var_mod.F90 to phys_state_var_mod.F90; ==> changes in physiq_mod.F90
and phys_output_write.F90
+ awake_dens added in phys_state_var_mod.F90
+ cv_gen and awake_dens now initialized in phyetat0.F90 and written in phyredem.F90
+ cv_gen, awake_dens, and solswfdiff now initialized in old_lmdz1d.F90 and scm.F90
+ useless variables suppressed in pbl_surface_mod.F90.

File:

• LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -249 +249 @@
 ! awake_dens  : number of active wakes per unit area
 ! wake_dens   : number of wakes per unit area
+! cv_gen      : birth rate of cumulonimbus per unit area.
 ! wake_occ    : occurence of wakes (= 1 if wakes occur, =0 otherwise)
 ! wake_Cstar  : vitesse d'etalement de la poche
@@ -261 +262 @@
       REAL,ALLOCATABLE,SAVE :: awake_dens(:), wake_dens(:)
 !$OMP THREADPRIVATE(awake_dens, wake_dens)
+      REAL,ALLOCATABLE,SAVE :: cv_gen(:)
+!$OMP THREADPRIVATE(cv_gen)
       REAL,ALLOCATABLE,SAVE :: wake_Cstar(:)
 !$OMP THREADPRIVATE(wake_Cstar)
@@ -490 +493 @@
       print*, 'allocate falb'
       ALLOCATE(falb_dir(klon,nsw,nbsrf),falb_dif(klon,nsw,nbsrf))
-      print*, 'allocate falb good', falb_dir(1,1,1)
+!!      print*, 'allocate falb good', falb_dir(1,1,1)
       ALLOCATE(chl_con(klon))
 !albedo SB <<<
@@ -588 +591 @@
       ALLOCATE(wake_deltat(klon,klev), wake_deltaq(klon,klev))
       ALLOCATE(wake_s(klon), awake_dens(klon), wake_dens(klon))
-      awake_dens = 0.
+!!      awake_dens = 0.  ! initialized in phyetat0
+      ALLOCATE(cv_gen(klon))
       ALLOCATE(wake_Cstar(klon))
       ALLOCATE(wake_pe(klon), wake_fip(klon))
@@ -742 +746 @@
       DEALLOCATE(wake_deltat, wake_deltaq)
       DEALLOCATE(wake_s, awake_dens, wake_dens)
+      DEALLOCATE(cv_gen)
       DEALLOCATE(wake_Cstar, wake_pe, wake_fip)
 !jyg<