Changeset 2100 for trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90

Timestamp:

Feb 14, 2019, 3:49:40 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:

• 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

MV

File:

• trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90

@@ -607 +607 @@
                 wq(l+1) = -MQtot 
              end if
+
+             if ( masse(l)*q(l) .lt. -(wq(l+1)-wq(l)) ) then  ! particular case when the vertical velocity is very strong in the layer and null below (wq(l)=0)
+               wq(l+1) = wq(l)-masse(l)*q(l)
+             end if
+
           endif
          endif  ! w<0 (up)
@@ -657 +662 @@
                 wq(l) = MQtot 
             end if
+
+             if ( masse(l)*q(l) .lt. -(wq(l+1)-wq(l)) ) then  ! particular case when the vertical velocity is very strong in the layer and null above (wq(l+1)=0)
+               wq(l) = wq(l+1)+masse(l)*q(l)
+             end if
+
           end if
          end if ! w>0 (down)