Changeset 3439 for trunk/LMDZ.PLUTO

Timestamp:

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

Author:

afalco

Message:

Pluto PCM: added some incompatibility checks.
AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• inifis_mod.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (2 diffs)

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

@@ -1438 +1438 @@
 
 ! Test of incompatibility:
-     if (condcosurf.and.no_n2frost) then
-       call abort_physic(rname, "CO surface condensation and no_n2frost are both active which may not be relevant", 1)
-     end if
+
+     if ((.not.tracer).and.(haze)) then
+       call abort_physic(rname, 'if haze are on, tracers must be on!', 1)
+
+     endif
+     if (callgasvis.and..not.callsoil) then
+         call abort_physic(rname, 'if callgasvis is set, callsoil must be true', 1)
+
+     endif        
+     if (paleo.and..not.fast) then
+         call abort_physic(rname, 'if paleo is set, fast must be true', 1)
+
+     endif        
+     if ((haze_proffix.or.haze_radproffix).and..not.aerohaze) then
+         call abort_physic(rname, 'for now, haze/rad proffix only works w aerohaze=T', 1)
+
+     endif    
+      if (condcosurf.and.no_n2frost) then
+        call abort_physic(rname, "CO surface condensation and no_n2frost are both active which may not be relevant", 1)
+      end if
 
      if ((cpp_mugaz_mode >= 1).and.(is_master).and.(ngrid>1)) then

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

@@ -376 +376 @@
       REAL zdqflow(ngrid,nq)
 
-      REAL zdteuv(ngrid,nlayer)    ! (K/s)
-      REAL zdtconduc(ngrid,nlayer) ! (K/s)
-      REAL zdumolvis(ngrid,nlayer)
-      REAL zdvmolvis(ngrid,nlayer)
-      real zdqmoldiff(ngrid,nlayer,nq)
-
       ! Haze relatated tendancies
       REAL zdqhaze(ngrid,nlayer,nq)
@@ -2444 +2438 @@
                call writediagfi(ngrid,"zdtch4cloud","ch4 cloud","T s-1",&
                            3,zdtch4cloud)
-             endif
+               call writediagfi(ngrid,"zdqch4cloud","ch4 cloud","T s-1",&
+                           3,zdqch4cloud(1,1,igcm_ch4_gas))
+            endif
 
          endif