Ignore:
Timestamp:
Aug 6, 2021, 3:12:49 PM (3 years ago)
Author:
fhourdin
Message:

Correction de l'appel a "soil" dans la version isotopique
suite à une modification de cette routine dans phhylmd.

Location:
LMDZ6/trunk/libf/phylmdiso
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmdiso/ocean_forced_mod.F90

    r3940 r3975  
    262262!
    263263    USE dimphy
     264    USE geometry_mod, ONLY: longitude,latitude
     265
    264266    USE calcul_fluxs_mod
    265267    USE surface_data,     ONLY : calice, calsno
     
    379381    IF (soil_model) THEN
    380382! update tsoil and calculate soilcap and soilflux
    381        CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, tsoil,soilcap, soilflux)
     383       CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, qsol, &
     384        & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil,soilcap, soilflux)
     385
    382386       cal(1:knon) = RCPD / soilcap(1:knon)
    383387       radsol(1:knon) = radsol(1:knon)  + soilflux(1:knon)
  • LMDZ6/trunk/libf/phylmdiso/surf_land_bucket_mod.F90

    r3927 r3975  
    3030    USE cpl_mod
    3131    USE dimphy
    32     USE geometry_mod, ONLY: latitude
     32    USE geometry_mod, ONLY: longitude,latitude
    3333    USE mod_grid_phy_lmdz
    3434    USE mod_phys_lmdz_para
     
    176176    IF (soil_model) THEN
    177177!       write(*,*) 'surf_land_bucket 185'
    178        CALL soil(dtime, is_ter, knon, snow, tsurf, tsoil, soilcap, soilflux)
     178       CALL soil(dtime, is_ter, knon, snow, tsurf, qsol,  &
     179      & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux)
     180
    179181!       write(*,*) 'surf_land_bucket 187'
    180182       DO i=1, knon
  • LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90

    r3940 r3975  
    4444#endif
    4545#endif
     46    USE geometry_mod,     ONLY : longitude,latitude
     47
    4648!FC
    4749    USE ioipsl_getin_p_mod, ONLY : getin_p
     
    5355
    5456    USE indice_sol_mod
     57
    5558
    5659!    INCLUDE "indicesol.h"
     
    337340    ! use soil model and recalculate properly cal
    338341    IF (soil_model) THEN
    339        CALL soil(dtime, is_lic, knon, snow, tsurf, tsoil, soilcap, soilflux)
     342       CALL soil(dtime, is_lic, knon, snow, tsurf, qsol, &
     343    & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux)
     344
    340345       cal(1:knon) = RCPD / soilcap(1:knon)
    341346       radsol(1:knon)  = radsol(1:knon) + soilflux(1:knon)
Note: See TracChangeset for help on using the changeset viewer.