Changeset 3712 for trunk/LMDZ.MARS/libf

Timestamp:

Apr 3, 2025, 10:45:25 AM (3 months ago)

Author:

emillour

Message:

Mars PCM:
Add extra checking in nirco2abs to ensure having nircorr==1 when o and
co2 tracers are available.
EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• conf_phys.F (modified) (1 diff)
• nirco2abs.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/conf_phys.F ¶

@@ -224 +224 @@
      $              "1-> new correction",
      $              "(matters only if callnirco2=T)"
-#ifdef MESOSCALE
+
          nircorr=0      !default value. this is OK below 60 km.
-#else
-         nircorr=0      !default value
-#endif
+         ! nircorr=1 should be prefered, but requires that co2
+         ! and o tracers be available
          call getin_p("nircorr",nircorr)
          write(*,*) " nircorr = ",nircorr

trunk/LMDZ.MARS/libf/phymars/nirco2abs.F ¶

@@ -126 +126 @@
             call abort_physic("nirco2abs","need an O tracer",1)
           endif
-        endif
+        else ! when nircorr==0
+          ! Make a sanity check; if both o and co2 tracers are available
+          ! then there is no reason to have nircorr==0
+          ico2=igcm_co2
+          io=igcm_o
+          if ((ico2/=0).and.(io/=0)) then
+            write(*,*) "nirco2abs error: O and CO2 tracers available"
+            write(*,*) " but nircorr=0; we recommended using nircorr=1"
+            call abort_physic("nirco2abs","need nircorr=1",1)
+          endif
+        endif ! of if (nircorr.eq.1)
         firstcall=.false.
       endif