Changeset 8220 for branches/ORCHIDEE_2_2


Ignore:
Timestamp:
2023-10-10T15:50:31+02:00 (15 months ago)
Author:
anne.cozic
Message:

Add initialization of ptnlev1 used in chemistry (specific case of VOC coupling between atm and surf) before it was calculated by thermosoil_main

Location:
branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/sechiba.f90

    r7710 r8220  
    517517         mc_layh,  mcl_layh,   soilmoist,       njsc ,     & 
    518518         frac_snow_veg,frac_snow_nobio,totfrac_nobio,     & 
    519          snowdz, snowrho, snowtemp, lambda_snow, cgrnd_snow, dgrnd_snow, pb) 
     519         snowdz, snowrho, snowtemp, lambda_snow, cgrnd_snow, dgrnd_snow, pb, & 
     520         ptnlev1) 
    520521 
    521522 

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/thermosoil.f90

    r7710 r8220  
    235235                                    mc_layh,       mcl_layh,   tmc_layh,        njsc,     & 
    236236                                    frac_snow_veg,frac_snow_nobio,totfrac_nobio, & 
    237                                     snowdz, snowrho, snowtemp, lambda_snow, cgrnd_snow, dgrnd_snow, pb) 
     237                                    snowdz, snowrho, snowtemp, lambda_snow, cgrnd_snow, dgrnd_snow, pb, & 
     238                                    ptnlev1) 
    238239 
    239240    !! 0. Variable and parameter declaration 
    … …  
    265266    REAL(r_std),DIMENSION (kjpindex,ngrnd), INTENT (out)  :: ftempdiag        !! temperature profile on full depth for stream temperature (K) 
    266267    REAL(r_std),DIMENSION (kjpindex),INTENT(out)          :: gtemp            !! First soil layer temperature 
     268    REAL(r_std),DIMENSION (kjpindex), INTENT (out)        :: ptnlev1          !! 1st level soil temperature    
    267269 
    268270    !! 0.3 Modified variables 
    … …  
    521523    stempdiag(:,:) = ptn(:,1:nslm) 
    522524    ftempdiag(:,:) = ptn(:,1:ngrnd) 
    523      
     525 
     526    !! Initialize output arguments to be used in sechiba 
     527    ptnlev1(:) = ptn(:,1) 
    524528 
    525529    !! 2.3. Computes cgrnd, dgrnd, soilflx and soilcap coefficients only if they were not found in restart file. 
Note: See TracChangeset for help on using the changeset viewer.