Changeset 3747 for LMDZ6/branches/Ocean_skin

Timestamp:

Jul 2, 2020, 4:58:11 PM (4 years ago)

Author:

lguez

Message:

Bug fix: revert revision [3601]

Bug fix: revert revision [3601]. Procedure phys_output_write is
called a first time by physiq before the call to phyetat0, so before
ftsoil is allocated. So we cannot reference `ftsoil(:, nsoilmx,
is_ter) in phys_output_write`.

Location:

LMDZ6/branches/Ocean_skin/libf/phylmd

Files:

• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (4 diffs)

LMDZ6/branches/Ocean_skin/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1948 +1948 @@
 
    type(ctrl_out), save:: o_delta_sst &
-        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'delta_sst', &
+        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'delta_SST', &
         "ocean-air interface temperature minus bulk SST", "K", '')
 

LMDZ6/branches/Ocean_skin/libf/phylmd/phys_output_write_mod.F90

@@ -40 +40 @@
          o_q2m, o_ustar, o_u10m, o_v10m, &
          o_wind10m, o_wind10max, o_gusts, o_sicf, &
-         o_psol, o_mass, o_qsurf, o_qsol, o_tsoil_deep_land, &
+         o_psol, o_mass, o_qsurf, o_qsol, &
          o_precip, o_rain_fall, o_rain_con, o_ndayrain, o_plul, o_pluc, o_plun, &
          o_snow, o_msnow, o_fsnow, o_evap, o_ep,o_epmax_diag, & ! epmax_cape
@@ -365 +365 @@
     USE ocean_slab_mod, ONLY: nslay, tslab, slab_bilg, tice, seaice, &
         slab_ekman,slab_hdiff,slab_gm,dt_ekman, dt_hdiff, dt_gm, dt_qflux
-    USE pbl_surface_mod, ONLY: snow, ftsoil
+    USE pbl_surface_mod, ONLY: snow
     USE indice_sol_mod, ONLY: nbsrf
     USE infotrac_phy, ONLY: nqtot, nqo, type_trac, tname, niadv
@@ -396 +396 @@
     INCLUDE "compbl.h"
     INCLUDE "YOMCST.h"
-    INCLUDE "dimsoil.h"
 
     ! Input
@@ -716 +715 @@
        IF (.NOT. ok_veget) THEN
           CALL histwrite_phy(o_qsol, qsol)
-          CALL histwrite_phy(o_tsoil_deep_land, ftsoil(:, nsoilmx, is_ter))
        ENDIF