Changeset 3438 for trunk/LMDZ.PLUTO/libf

Timestamp:

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

Author:

afalco

Message:

Pluto PCM: Minor fixes.
AF

Location:

trunk/LMDZ.PLUTO/libf

Files:

• dynphy_lonlat/phypluto/newstart.F (modified) (1 diff)
• phypluto/condense_n2.F90 (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F

@@ -1454 +1454 @@
      &              (lonfi(ig)*180./pi.ge.val4)  .and.
      &              (lonfi(ig)*180./pi.lt.val5)  .and. 
-     &              (qsurf(ig,igcm_n2).lt.val7))  then     
+     &              (qsurf(ig,igcm_n2).gt.val7))  then     
                     qsurf(ig,igcm_co_ice)=qsurf(ig,igcm_co_ice)*val3
                     qsurf(ig,igcm_co_ice)=qsurf(ig,igcm_co_ice)+val6

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)