Changeset 2823 for trunk/LMDZ.MARS/libf/phymars/dyn1d

Timestamp:

Nov 18, 2022, 10:43:43 AM (3 years ago)

Author:

emillour

Message:

Mars GCM:
Remove the "tracer" (logical) flag as we now always run with at least
one tracer.
EM

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -221 +221 @@
       endif
 
-! check if we are going to run with or without tracers
-      write(*,*) "Run with or without tracer transport ?"
-      tracer=.false. ! default value
-      call getin("tracer",tracer)
-      write(*,*) " tracer = ",tracer
-
-! while we're at it, check if there is a 'traceur.def' file
+
+! check if there is a 'traceur.def' file
 ! and process it.
-      if (tracer) then
       ! load tracer names from file 'traceur.def'
         open(90,file='traceur.def',status='old',form='formatted',
@@ -246 +240 @@
             write(*,*) "testphys1d: error reading number of tracers"
             write(*,*) "   (first line of traceur.def) "
+            stop
+          endif
+          if (nq<1) then
+            write(*,*) "testphys1d: error number of tracers"
+            write(*,*) "is nq=",nq," but must be >=1!"
             stop
           endif
@@ -318 +317 @@
         write(*,*) "testphys1d: initializing tracers"
         call read_profile(nq, nlayer, qsurf, q)
-      else
-      ! we still need to set (dummy) tracer number and names for physdem1
-        nq=1
-        nqtot=nq ! set value of nqtot (in infotrac module) as nq
-        ! allocate arrays:
-        allocate(tname(nq))
-        allocate(q(nlayer,nq))
-        allocate(qsurf(nq))
-        allocate(dq(nlayer,nq))
-        allocate(dqdyn(nlayer,nq))
-        allocate(mqtot(nq))
-        do iq=1,nq
-          write(str7,'(a1,i2.2)')'t',iq
-          tname(iq)=str7
-        enddo
-      ! and just to be clean, also initialize tracers to zero for physdem1
-        q(:,:)=0 
-        qsurf(:)=0      
-      endif ! of if (tracer)
-      
-      !write(*,*) "testphys1d q", q(1,:)
-      !write(*,*) "testphys1d qsurf", qsurf
 
 c  Date and local time at beginning of run