Ignore:
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:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Ocean_skin/libf/phylmd/phys_output_ctrlout_mod.F90

    r3744 r3747  
    19481948
    19491949   type(ctrl_out), save:: o_delta_sst &
    1950         = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'delta_sst', &
     1950        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'delta_SST', &
    19511951        "ocean-air interface temperature minus bulk SST", "K", '')
    19521952
  • LMDZ6/branches/Ocean_skin/libf/phylmd/phys_output_write_mod.F90

    r3744 r3747  
    4040         o_q2m, o_ustar, o_u10m, o_v10m, &
    4141         o_wind10m, o_wind10max, o_gusts, o_sicf, &
    42          o_psol, o_mass, o_qsurf, o_qsol, o_tsoil_deep_land, &
     42         o_psol, o_mass, o_qsurf, o_qsol, &
    4343         o_precip, o_rain_fall, o_rain_con, o_ndayrain, o_plul, o_pluc, o_plun, &
    4444         o_snow, o_msnow, o_fsnow, o_evap, o_ep,o_epmax_diag, & ! epmax_cape
     
    365365    USE ocean_slab_mod, ONLY: nslay, tslab, slab_bilg, tice, seaice, &
    366366        slab_ekman,slab_hdiff,slab_gm,dt_ekman, dt_hdiff, dt_gm, dt_qflux
    367     USE pbl_surface_mod, ONLY: snow, ftsoil
     367    USE pbl_surface_mod, ONLY: snow
    368368    USE indice_sol_mod, ONLY: nbsrf
    369369    USE infotrac_phy, ONLY: nqtot, nqo, type_trac, tname, niadv
     
    396396    INCLUDE "compbl.h"
    397397    INCLUDE "YOMCST.h"
    398     INCLUDE "dimsoil.h"
    399398
    400399    ! Input
     
    716715       IF (.NOT. ok_veget) THEN
    717716          CALL histwrite_phy(o_qsol, qsol)
    718           CALL histwrite_phy(o_tsoil_deep_land, ftsoil(:, nsoilmx, is_ter))
    719717       ENDIF
    720718
Note: See TracChangeset for help on using the changeset viewer.