Changeset 3438 for trunk/LMDZ.PLUTO/libf/phypluto

Timestamp:

Sep 25, 2024, 4:43:31 PM (4 months ago)

Author:

afalco

Message:

Pluto PCM: Minor fixes.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/condense_n2.F90 (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/condense_n2.F90

@@ -171 +171 @@
 
      ! calculate global mean surface pressure for the fast mode
+     IF (.not. ALLOCATED(kp)) ALLOCATE(kp(klon))
+     DO ig=1,klon
+        kp(ig) = exp(-phisfi(ig)/(r*38.))
+     ENDDO
      IF (fast) THEN
-        IF (.not. ALLOCATED(kp)) ALLOCATE(kp(klon))
-        DO ig=1,klon
-           kp(ig) = exp(-phisfi(ig)/(r*38.))
-        ENDDO
         p00=glob_average2d(kp) ! mean pres at ref level
      ENDIF
@@ -382 +382 @@
         pdpsrf(ig)   = -pdicen2(ig)*g
     !    pdpsrf(ig)   = 0. ! OPTION to check impact N2 sub/cond
-        IF (zplev(ig)+pdpsrf(ig)*subtimestep.le.0.0000001) then
+        IF (fast.and.(zplev(ig)+pdpsrf(ig)*subtimestep.le.0.0000001)) then
             pdpsrf(ig)=(0.0000001*kp(ig)/p00-zplev(ig))/subtimestep
             pdicen2(ig)=-pdpsrf(ig)/g
@@ -855 +855 @@
          Mtot = masse(m)
          MQtot = masse(m)*q(m)
-         !if (m.lt.klev) then ! because some compilers will have problems
-         !                      ! evaluating masse(klev+1)
+         if (m.lt.klev) then ! because some compilers will have problems
+                              ! evaluating masse(klev+1)
          do while ((m.lt.klev).and.(w(l+1).gt.(Mtot+masse(m+1))))
             m=m+1
             Mtot = Mtot + masse(m)
             MQtot = MQtot + masse(m)*q(m)
-         !   if (m.eq.klev) exit
+           if (m.eq.klev) exit
          end do
-         !endif
+         endif
          if (m.lt.klev) then
             sigw=(w(l+1)-Mtot)/masse(m+1)