Changeset 3906


Ignore:
Timestamp:
Aug 28, 2025, 11:51:20 PM (6 months ago)
Author:
llange
Message:

Mars PCM
Fixing index bug in soilwater (index kept constant in a loop).
LL

Location:
trunk/LMDZ.MARS
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/changelog.txt

    r3905 r3906  
    49644964== 28/08/2025 == LL
    49654965Cleaning of the 'soilwater' routine.
     4966
     4967== 28/08/2025 == LL
     4968Fixing index bug in soilwater (index kept constant in a loop).
     4969

  • trunk/LMDZ.MARS/libf/phymars/soilwater.F90

    r3905 r3906  
    595595        ! Set loop flags
    596596        do ik = 1, nsoil
    597             condensation_flag = .false.
     597            condensation_flag(ik) = .false.
    598598        enddo
    599599       
    … …  
    704704                    znsoil(ig, ik) = zalpha(ig, ik) * znsoil(ig, ik + 1) + zbeta(ig, ik)
    705705                   
    706                     if (.not. over_mono_sat_flag(ig, nsoil)) then  ! modified 2024
     706                    if (.not. over_mono_sat_flag(ig, ik)) then  ! modified 2024
    707707                        adswater_temp(ig, ik) = (Ka(ig, ik) * ptimestep * znsoil(ig, ik) + adswprev(ig, ik)) &
    708708                                              / (1.D0 + ptimestep * Kd(ig, ik))
    … …  
    969969                        znsoil(ig, ik) = zalpha(ig, ik) * znsoil(ig, ik + 1) + zbeta(ig, ik)
    970970                       
    971                         if (.not. over_mono_sat_flag(ig, nsoil)) then  ! modified 2024
     971                        if (.not. over_mono_sat_flag(ig, ik)) then  ! modified 2024
    972972                            adswater_temp(ig, ik) = (Ka(ig, ik) * ptimestep * znsoil(ig, ik) + adswprev(ig, ik)) &
    973973                                                  / (1.D0 + ptimestep * Kd(ig, ik))
Note: See TracChangeset for help on using the changeset viewer.