Changeset 2112 for trunk/LMDZ.MARS/libf/phymars

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

File:

• trunk/LMDZ.MARS/libf/phymars/vlz_fi.F (modified) (2 diffs)

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