Changeset 2112

Timestamp:

Mar 6, 2019, 4:43:27 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:

• correction of the Van Leer scheme in vlz_fi.F, used by the sedimentation routine (particular case evoked previously)

MV

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/vlz_fi.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -2703 +2703 @@
 == 14/02/2019 == MV
 - handling of a particular case in Van Leer advection in rocketduststorm_mod.F90: when the vertical velocity is very large at a particular layer and null below or above, mass fluxes must not entrain more tracers than available in the layer
+
+== 06/03/2019 == MV
+- correction of the Van Leer scheme in vlz_fi.F, used by the sedimentation routine (particular case evoked previously)

trunk/LMDZ.MARS/libf/phymars/vlz_fi.F

@@ -122 +122 @@
                 wq(ij,l) = Mqtot 
             end if
+
+c            it cannot entrain more than available mass !
+             if ( masse(ij,l)*q(ij,l) .lt. -(wq(ij,l+1)-wq(ij,l)) ) then
+               wq(ij,l) = wq(ij,l+1)+masse(ij,l)*q(ij,l)
+             end if
+
           end if
          end if
@@ -169 +175 @@
                 stop
              end if
+
+c            it cannot entrain more than available mass !
+             if ( masse(ij,l)*q(ij,l) .lt. -(wq(ij,l+1)-wq(ij,l)) ) then
+               wq(ij,l+1) = wq(ij,l)-masse(ij,l)*q(ij,l)
+             end if
+
           endif
          endif