Changeset 2693

Timestamp:

Jun 21, 2022, 11:05:32 AM (3 years ago)

Author:

aslmd

Message:

Added radius and rho option in modern traceur

Location:

trunk/LMDZ.GENERIC

Files:

• deftank/traceur.def.modern (modified) (1 diff)
• libf/phystd/initracer.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/deftank/traceur.def.modern

@@ -52 +52 @@
 #               is_recomb_qotf ! 1 if tracer recombination is done on-the-fly, else 0
 #                                (if 1, must have is_recomb_qset=0)
+#               is_generic     ! 1 if tracer is generic, else 0
+#               radius          ! radius value in meter
+#               rho             ! density in kg/m3
 #
 # Chemistry:    SF_mode        ! 1 if surface set up value, else 2 sedimentation velocity

trunk/LMDZ.GENERIC/libf/phystd/initracer.F90

@@ -617 +617 @@
               ' is_generic=', is_generic(iq)       
           endif
+          !option radius
+          if (index(tracline,'radius=') .ne. 0) then 
+            read(tracline(index(tracline,'radius=') &
+              +len('radius='):),*) radius(iq)
+            write(*,*)'Parameter value (traceur.def) :'// &
+              "radius=",radius(iq), "m"
+          else
+            write(*,*) ' Parameter value (default)     :'// &
+            ' radius=', radius(iq)," m"     
+        endif
+        !option rho
+        if (index(tracline,'rho=') .ne. 0) then 
+          read(tracline(index(tracline,'rho=') &
+          +len('rho='):),*) rho_q(iq)
+          write(*,*)'Parameter value (traceur.def) :'// &
+            "rho=",rho_q(iq)
+        else
+          write(*,*) ' Parameter value (default)     :'// &
+          ' rho=', rho_q(iq)       
+      endif
       end subroutine get_tracdat