Changeset 1764 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Aug 7, 2017, 4:38:48 PM (7 years ago)

Author:

jvatant

Message:

Fixed unitialized radius and qext in initracer

-> These arrays were printed and even tested in callsedim without being initialized !
-> Please stop handling uninitialized values !!!

Also updated README for modif from r1725

JVO

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callsedim.F (modified) (1 diff)
• initracer.F (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/callsedim.F

@@ -106 +106 @@
  
       do iq=1,nq
-       if((radius(iq).gt.1.e-9).and.(iq.ne.igcm_co2_ice)) then   
+       if((radius(iq).gt.1.e-9).and.(iq.ne.igcm_co2_ice)) then ! JVO 08/2017 : be careful radius was tested uninitialized (fixed) ... 
+       
 !         (no sedim for gases, and co2_ice sedim is done in condense_co2)      
 

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

@@ -61 +61 @@
        alpha_lift(:)=0.
        alpha_devil(:)=0.
+       
+       ! Added by JVO 2017 : these arrays are handled later
+       ! -> initialization is the least we can do, please !!!
+       radius(:)=0.
+       qext(:)=0.
 
 ! Initialization: get tracer names from the dynamics and check if we are