Changeset 2100 for trunk/LMDZ.MARS

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

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/rocketduststorm_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -2700 +2700 @@
 == 06/02/2019 == MV
 - follow-up of the last change in rocketduststorm_mod.F90: corrections in vector definitions
+
+== 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

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)