Changeset 1427 for trunk/MESOSCALE

Timestamp:

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

Author:

aslmd

Message:

LES VENUS. correction finale bug heating rates and parallelism

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC

Files:

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

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

@@ -172 +172 @@
    stretch_grid = .true.
 !  FOR LES, set stretch to false
-   stretch_grid = .false.
+!   stretch_grid = .false.
    delt = 3.
 !   z_scale = .50
    z_scale = .40
-   z_scale = .25
-   z_scale = .15
+!   z_scale = .25
+!   z_scale = .15
    pi = 2.*asin(1.0)
    write(6,*) ' pi is ',pi
@@ -875 +875 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     call read_dust(profdustq,profdustn,nl_in)
-    DO k=1,kte-1
+    DO k=1,kte!-1
        p_level = grid%znu(k)*(pd_surf - grid%p_top) + grid%p_top
        DO j = jts, jte
@@ -884 +884 @@
        ENDDO
        ENDDO 
+       print*,k,grid%m_q2(1,k,1)
     ENDDO
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

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