Changeset 3712 for trunk/LMDZ.MARS

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

Files:

• changelog.txt (modified) (1 diff)
• deftank/callphys.def.GCM6 (modified) (1 diff)
• deftank/callphys.def.MCD6 (modified) (1 diff)
• libf/phymars/conf_phys.F (modified) (1 diff)
• libf/phymars/nirco2abs.F (modified) (1 diff)

trunk/LMDZ.MARS/changelog.txt

@@ -4789 +4789 @@
 Change default scheme for molecular diffusion in the thermosphere to the
 MPF scheme. 
+Add extra checking in nirco2abs to ensure having nircorr==1 when o and
+co2 tracers are available.

trunk/LMDZ.MARS/deftank/callphys.def.GCM6

@@ -67 +67 @@
 callnirco2 = .true.
 # NIR NLTE correction ? matters only if callnirco2=T
+# nircorr should be set to 1 if co2 and o tracers are available, 0 otherwise
 nircorr=0
 # call turbulent vertical diffusion ?

trunk/LMDZ.MARS/deftank/callphys.def.MCD6

@@ -67 +67 @@
 callnirco2 = .true.
 # NIR NLTE correction ? matters only if callnirco2=T
-nircorr=0
+# nircorr should be set to 1 if co2 and o tracers are available, 0 otherwise
+nircorr=1
 # call turbulent vertical diffusion ?
 calldifv  = .true.

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