Changeset 1427 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys

Timestamp:

May 18, 2015, 3:32:29 PM (10 years ago)

Author:

aslmd

Message:

LES VENUS. correction finale bug heating rates and parallelism

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (1 diff)

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

@@ -1499 +1499 @@
 !!!!!!!!!!!!!!!!!!!
 IF (firstcall .EQV. .true.) THEN
-  ALLOCATE(dtrad(nlayer))
-  dtrad(1:nlayer) = M_Q2(1,1:nlayer,1)
-  print *, "HEATING RATE", dtrad(kps:kpe)
+  ALLOCATE(dtrad(nlayer+1))
+  DO k=kps,kpe+1
+    !! the only solution to avoid 0 points in M_Q2
+    !! -- in each domain decomposition case
+    dtrad(k) = MAXVAL(M_Q2(:,k,:)) + MINVAL(M_Q2(:,k,:))
+  ENDDO
+  print*,"HEATING RATE",dtrad(kps:kpe+1)
 ENDIF
 DO i= 1,ngrid
-  pdt(i,kps:kpe) = dtrad(1:nlayer)
+  pdt(i,kps:kpe) = dtrad(kps:kpe)
 ENDDO
 !!!!!!!!!!!!!!!!!!!