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

Timestamp:

Jul 23, 2017, 9:52:26 PM (7 years ago)

Author:

aslmd

Message:

MESOSCALE VENUS. commented cp=f(T) to easily find where changes are. removed obsolete commented tests or prints.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em

Files:

• module_big_step_utilities_em.F (modified) (9 diffs)
• solve_em.F (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_big_step_utilities_em.F

@@ -18 +18 @@
    USE module_configure
    USE module_wrf_error
+
 CONTAINS
 
@@ -568 +569 @@
       ENDDO
       ENDDO
-      !print*,'muu(i,j)',muu(1,:)
+
    ELSE IF (name .EQ. 'v')THEN
 
@@ -587 +588 @@
       ENDDO
       ENDDO
-      !print*,'muv(i,j)',muv(1,:)
+
    ELSE IF (name .EQ. 'w')THEN
       itf=MIN(ite,ide-1)
@@ -620 +621 @@
       ENDDO
       ENDDO
-      !print*,'field(i,k,j)',field(1,:,1)
-      !print*,'(mu(i,j)+mub(i,j))',(mu(1,:)+mub(1,:))   
+
    ENDIF
 
@@ -1052 +1052 @@
 
   IF (non_hydrostatic) THEN
-      print*,'n_moist',n_moist
-      print*,'PARAM_FIRST_SCALAR',PARAM_FIRST_SCALAR
+
       IF (n_moist >= PARAM_FIRST_SCALAR ) THEN  
-        print*,'p(i,k,j)_before',p(1,:,1)
+
         DO j=jts,jtf
         DO k=kts,ktf
@@ -1076 +1075 @@
         ENDDO
         ENDDO
-        print*,'p(i,k,j)_after',p(1,:,1)
+
       ELSE
 
@@ -4076 +4075 @@
    INTEGER :: i_start, i_end, j_start, j_end, k_start, k_end, i_startu, j_startv
    INTEGER :: i, j, k
-   REAL    :: w1, w2, z0, z1, z2,temp
+   REAL    :: w1, w2, z0, z1, z2
 
 !-----------------------------------------------------------------------
@@ -4100 +4099 @@
 
 !  compute thermodynamics and velocities at pressure points
+
     do j = j_start,j_end
     do k = k_start, k_end
@@ -4107 +4107 @@
       p_phy(i,k,j) = p(i,k,j) + pb(i,k,j)
       pi_phy(i,k,j) = (p_phy(i,k,j)/p1000mb)**rcp
+      !! TAKE INTO ACCOUNT cp=f(T) on Venus
       IF (planet .eq. "venus" ) THEN
         t_phy(i,k,j)= (th_phy(i,k,j)**nu - nu*(TT00**nu)*log((p1000mb/p_phy(i,k,j))**rcp))**(1/nu)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F

@@ -160 +160 @@
   ijds = min(ids, jds)
   ijde = max(ide, jde)
-  !num_moist = 3
+
   num_3d_m        = num_moist
   num_3d_c        = num_chem