Changeset 2703

Timestamp:

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

Author:

aslmd

Message:

fixed the initialisation of rho_q + some indentations

File:

• trunk/LMDZ.GENERIC/libf/phystd/initracer.F90 (modified) (3 diffs)

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

@@ -134 +134 @@
        qext(:)=0.
        is_generic(:)= 0 !LT
+       rho_q(:) = 0. !need to be init here if we want to read it from modern traceur with get_tracdat
 
 
@@ -439 +440 @@
 !     Initialisation tracers ....
 !------------------------------------------------------------
-      rho_q(1:nq)=0
+      ! rho_q(1:nq)=0
 
       rho_dust=2500.  ! Mars dust density (kg.m-3)
@@ -626 +627 @@
             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
+          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