Changeset 4034

Timestamp:

Nov 30, 2021, 7:36:27 PM (3 years ago)

Author:

lguez

Message:

Bug fix: define fields written to startphy

When compiled in debug mode, CE0L wrote NaNs? in startphy for variables
awake_dens, cv_gen, delta_tsurf, beta_aridity, ratqs_inter
and reading this startphy with gcm crashed in procedure
phyetat0_get, when computing min and max.

Location:

LMDZ6/trunk

Files:

• . (modified) (1 prop)
• libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90

@@ -93 +93 @@
   USE conf_phys_m, ONLY: conf_phys
   USE init_ssrf_m, ONLY: start_init_subsurf
+  USE phys_state_var_mod, ONLY: beta_aridity, delta_tsurf, awake_dens, cv_gen, &
+       ratqs_inter
   !use ioipsl_getincom
   IMPLICIT NONE
@@ -127 +129 @@
   LOGICAL :: filtre_oro=.false.
 
+  INCLUDE "compbl.h"
+  INCLUDE "thermcell.h"
+  
   deg2rad= pi/180.0
   iml=assert_eq(SIZE(masque,1),SIZE(phis,1),TRIM(modname)//" iml")
@@ -230 +235 @@
   wake_delta_pbl_TKE(:,:,:)=0
   wake_dens(:)=0
+  awake_dens = 0.
+  cv_gen = 0.
   ale_bl = 0.
   ale_bl_trig =0.
@@ -272 +279 @@
   CALL fonte_neige_init(run_off_lic_0)
   CALL pbl_surface_init( fder, snsrf, qsurf, tsoil )
+
+  IF (iflag_pbl>1 .AND. iflag_wake>=1  .AND. iflag_pbl_split >=1) then
+     delta_tsurf = 0.
+     beta_aridity = 0.
+  end IF
+
+  ratqs_inter = 0.002
   CALL phyredem( "startphy.nc" )