Changeset 1425 for trunk/MESOSCALE/LMD_MM_MARS/SRC

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

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC

Files:

• LES/modif_mars/module_initialize_les.F.venus (modified) (1 diff)
• WRFV2/phys/module_lmd_driver.F (modified) (2 diffs)

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)
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

@@ -166 +166 @@
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)  :: PSFC,TSK
    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(OUT)  :: RTHBLTEN,RUBLTEN,RVBLTEN
-   REAL,DIMENSION(:),ALLOCATABLE,SAVE :: dtsw,dtlw
+   REAL,DIMENSION(:),ALLOCATABLE,SAVE :: dtrad
 #endif
 !-------------------------------------------
@@ -1499 +1499 @@
 !!!!!!!!!!!!!!!!!!!
 IF (firstcall .EQV. .true.) THEN
-  ALLOCATE(dtsw(nlayer))
-  ALLOCATE(dtlw(nlayer))
-  dtsw(1:nlayer) = M_Q2(1,1:nlayer,1)
-  dtlw(1:nlayer) = M_Q2(2,1:nlayer,2)
-ENDIF
-!print *, "SW", dtsw(kps:kpe) ! sw
-!print *, "LW", dtlw(kps:kpe) ! lw
+  ALLOCATE(dtrad(nlayer))
+  dtrad(1:nlayer) = M_Q2(1,1:nlayer,1)
+  print *, "HEATING RATE", dtrad(kps:kpe)
+ENDIF
 DO i= 1,ngrid
-  pdt(i,kps:kpe) = dtsw(1:nlayer) + dtlw(1:nlayer) 
+  pdt(i,kps:kpe) = dtrad(1:nlayer)
 ENDDO
 !!!!!!!!!!!!!!!!!!!