Changeset 1425 for trunk/MESOSCALE/LMD_MM_MARS/SRC/LES

Timestamp:

May 11, 2015, 7:54:05 AM (10 years ago)

Author:

aslmd

Message:

nophys LES. corrected a bug in parallel in heating rates prescription

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F.venus (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F.venus

@@ -877 +877 @@
     DO k=1,kte-1
        p_level = grid%znu(k)*(pd_surf - grid%p_top) + grid%p_top
-       !!! we use Q2 as a vehicle for heating rates! sick!
-       grid%m_q2(1,k,1) = interp_0( profdustq, pd_in, p_level, nl_in )
-       grid%m_q2(2,k,2) = interp_0( profdustn, pd_in, p_level, nl_in )
+       DO j = jts, jte
+       DO i = its, ite
+         !!! we use Q2 as a vehicle for heating rates! sick!
+         grid%m_q2(i,k,j) = interp_0( profdustq, pd_in, p_level, nl_in ) &
+                          + interp_0( profdustn, pd_in, p_level, nl_in )
+       ENDDO
+       ENDDO 
     ENDDO
-    print *, "HEATING RATE SW", grid%m_q2(1,:,1)
-    print *, "HEATING RATE LW", grid%m_q2(2,:,2)
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!