Changeset 3642 for trunk

Timestamp:

Feb 24, 2025, 1:48:19 PM (4 months ago)

Author:

afalco

Message:

Pluto: minor fixes & better message when stumbling onto error.
AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• aeroptproperties.F90 (modified) (1 diff)
• condense_n2.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (2 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/aeroptproperties.F90 ¶

@@ -328 +328 @@
         IF ((MINVAL(reffrad).LT.minrad).OR.(MAXVAL(reffrad).GT.maxrad)) then
            WRITE(*,*) 'Warning: particle size in grid box #'
-           WRITE(*,*) ig,' is larger than optical properties. '
+           WRITE(*,*) i,' is larger than optical properties. '
            WRITE(*,*) 'reffrad=',MINVAL(reffrad),'-',MAXVAL(reffrad)
            WRITE(*,*) 'radiustab=',minrad,'-',maxrad

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

@@ -820 +820 @@
             w(l+1) = Mtot
             qm(l+1) = Mqtot / Mtot
-            write(*,*) 'top layer is disapearing !'
-            stop
+            print*,"m", m
+            print*,"Mtot", Mtot
+            print*,"Mqtot", Mqtot
+            call abort_physic("condense_n2","top layer is disapearing !",1)
          end if
       else      ! if(w(l+1).lt.0)

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 ¶

@@ -477 +477 @@
       ! local variables for skin depth check
       real :: therm_inertia(ngrid,nsoilmx)
-    !   real :: tidat_out(ngrid,nsoilmx)
       real :: inertia_min,inertia_max
       real :: diurnal_skin ! diurnal skin depth (m)
@@ -1161 +1160 @@
                     zdqdif,zdqsdif,qsat_ch4,qsat_ch4_l1) !,zq1temp_ch4,qsat_ch4)
 
-            zdtdif(1:ngrid,1:nlayer)=zdhdif(1:ngrid,1:nlayer)*zpopsk(1:ngrid,1:nlayer) 
+            zdtdif(1:ngrid,1:nlayer)=zdhdif(1:ngrid,1:nlayer)*zpopsk(1:ngrid,1:nlayer)
 
             bcond=1./tcond1p4Pa