Changeset 3465 for LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90

Timestamp:

Mar 14, 2019, 10:34:31 AM (5 years ago)

Author:

Laurent Fairhead

Message:

Further modifications for DYNAMICO/LMDZ convergence. These are based
on Yann's LMDZ6_V2 sources. Compiles on irene and converges with revision 3459
in a bucket configuration
YM/LF

File:

• LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90 (modified) (3 diffs)

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

@@ -40 +40 @@
   USE phys_state_var_mod, ONLY: zmea, zstd, zsig, zgam, zthe, zpic, zval, z0m, &
           solsw, radsol, t_ancien, wake_deltat, wake_s,  rain_fall, qsol, z0h, &
-          sollw, rugoro, q_ancien, wake_deltaq, wake_pe, snow_fall, ratqs,w01, &
+          sollw,sollwdown, rugoro, q_ancien, wake_deltaq, wake_pe, snow_fall, ratqs,w01, &
     sig1, ftsol, clwcon, fm_therm, wake_Cstar,  pctsrf,  entr_therm,radpas, f0,&
-    zmax0,fevap, rnebcon,falb_dir, wake_fip,    agesno,  detr_therm, pbl_tke,  &
+    zmax0,fevap, rnebcon,falb_dir, falb_dif, wake_fip,    agesno,  detr_therm, pbl_tke,  &
     phys_state_var_init, ql_ancien, qs_ancien, prlw_ancien, prsw_ancien, &
-    prw_ancien, sollwdown
+    prw_ancien, u10m,v10m, treedrg, u_ancien, v_ancien, wake_delta_pbl_TKE, wake_dens, &
+    ale_bl, ale_bl_trig, alp_bl
   USE comconst_mod, ONLY: pi, dtvr
 
@@ -194 +195 @@
   falb_dir(:, :, is_oce) = 0.5
   falb_dir(:, :, is_sic) = 0.6
+
+!ym warning missing init for falb_dif => set to 0
+  falb_dif(:,:,:)=0
+
+  u10m(:,:)=0  
+  v10m(:,:)=0  
+  treedrg(:,:,:)=0
+
   fevap(:,:) = 0.
   DO i=1,nbsrf; qsolsrf(:,i)=150.; END DO
@@ -211 +220 @@
   prw_ancien = 0.
   agesno     = 0.
-
+  
+  u_ancien = 0.
+  v_ancien = 0.
+  wake_delta_pbl_TKE(:,:,:)=0
+  wake_dens(:)=0
+  ale_bl = 0.
+  ale_bl_trig =0.
+  alp_bl=0.
+  
+  z0m(:,:)=0 ! ym missing 5th subsurface initialization
+ 
   z0m(:,is_oce) = rugmer(:)
   z0m(:,is_ter) = MAX(1.0e-05,zstd(:)*zsig(:)/2.0)