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

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/WRFV2/phys/module_lmd_driver.F (modified) (2 diffs)

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