Changeset 3904 for trunk/LMDZ.MARS/libf/phymars/initracer.F90

Timestamp:

Aug 26, 2025, 2:17:21 PM (5 months ago)

Author:

jbclement

Message:

Mars PCM:
Some adjustments of sanity checks for the compatibility of options 'doubleq' and 'active'.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/initracer.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/initracer.F90

@@ -593 +593 @@
 
 if (doubleq) then
+    ! verify that we indeed have dust_mass and dust_number tracers
+    if (igcm_dust_mass == 0) then
+           write(*,*) "initracer: error!"
+            write(*,*) "  Cannot use doubleq option without a dust_mass tracer!"
+           call abort_physic("initracer","doubleq issue",1)
+    endif
+    if (igcm_dust_number == 0) then
+        write(*,*) "initracer: error!"
+        write(*,*) "  Cannot use doubleq option without a dust_number tracer!"
+        call abort_physic("initracer","doubleq issue",1)
+    endif
+
     ! "doubleq" technique
     ! -------------------
@@ -785 +797 @@
 ! some extra (possibly redundant) sanity checks for tracers:
 ! ---------------------------------------------------------
-if (doubleq) then
-    ! verify that we indeed have dust_mass and dust_number tracers
-    if (igcm_dust_mass == 0) then
-           write(*,*) "initracer: error!"
-            write(*,*) "  Cannot use doubleq option without a dust_mass tracer!"
-           call abort_physic("initracer","doubleq issue",1)
-    endif
-    if (igcm_dust_number == 0) then
-        write(*,*) "initracer: error!"
-        write(*,*) "  Cannot use doubleq option without a dust_number tracer!"
-        call abort_physic("initracer","doubleq issue",1)
-        endif
-endif
-
 if ((.not.doubleq).and.(dustbin>0)) then
     ! verify that we indeed have 'dustbin' dust tracers