Ignore:
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
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_big_step_utilities_em.F

    r1724 r1729  
    1414
    1515   USE module_domain, ONLY : domain
    16    USE module_configure
    1716   USE module_model_constants
    1817   USE module_state_description
     
    44764475      p_phy(i,k,j) = p(i,k,j) + pb(i,k,j)
    44774476      pi_phy(i,k,j) = (p_phy(i,k,j)/p1000mb)**rcp
     4477      !! TAKE INTO ACCOUNT cp=f(T) on Venus
    44784478      IF (planet.eq. "venus" ) THEN
    4479       !call tpot2t(1,th_phy(i,k,j), t_phy(i,k,j), pi_phy(i,k,j))
     4479        !call tpot2t(1,th_phy(i,k,j), t_phy(i,k,j), pi_phy(i,k,j))
    44804480        t_phy(i,k,j)= (th_phy(i,k,j)**nu - nu*(TT00**nu)*log((p1000mb/p_phy(i,k,j))**rcp))**(1/nu)
    44814481      ELSE
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_big_step_utilities_em.F

    r1724 r1729  
    1818   USE module_configure
    1919   USE module_wrf_error
     20
    2021CONTAINS
    2122
     
    568569      ENDDO
    569570      ENDDO
    570       !print*,'muu(i,j)',muu(1,:)
     571
    571572   ELSE IF (name .EQ. 'v')THEN
    572573
     
    587588      ENDDO
    588589      ENDDO
    589       !print*,'muv(i,j)',muv(1,:)
     590
    590591   ELSE IF (name .EQ. 'w')THEN
    591592      itf=MIN(ite,ide-1)
     
    620621      ENDDO
    621622      ENDDO
    622       !print*,'field(i,k,j)',field(1,:,1)
    623       !print*,'(mu(i,j)+mub(i,j))',(mu(1,:)+mub(1,:))   
     623
    624624   ENDIF
    625625
     
    10521052
    10531053  IF (non_hydrostatic) THEN
    1054       print*,'n_moist',n_moist
    1055       print*,'PARAM_FIRST_SCALAR',PARAM_FIRST_SCALAR
     1054
    10561055      IF (n_moist >= PARAM_FIRST_SCALAR ) THEN 
    1057         print*,'p(i,k,j)_before',p(1,:,1)
     1056
    10581057        DO j=jts,jtf
    10591058        DO k=kts,ktf
     
    10761075        ENDDO
    10771076        ENDDO
    1078         print*,'p(i,k,j)_after',p(1,:,1)
     1077
    10791078      ELSE
    10801079
     
    40764075   INTEGER :: i_start, i_end, j_start, j_end, k_start, k_end, i_startu, j_startv
    40774076   INTEGER :: i, j, k
    4078    REAL    :: w1, w2, z0, z1, z2,temp
     4077   REAL    :: w1, w2, z0, z1, z2
    40794078
    40804079!-----------------------------------------------------------------------
     
    41004099
    41014100!  compute thermodynamics and velocities at pressure points
     4101
    41024102    do j = j_start,j_end
    41034103    do k = k_start, k_end
     
    41074107      p_phy(i,k,j) = p(i,k,j) + pb(i,k,j)
    41084108      pi_phy(i,k,j) = (p_phy(i,k,j)/p1000mb)**rcp
     4109      !! TAKE INTO ACCOUNT cp=f(T) on Venus
    41094110      IF (planet .eq. "venus" ) THEN
    41104111        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

    r1724 r1729  
    160160  ijds = min(ids, jds)
    161161  ijde = max(ide, jde)
    162   !num_moist = 3
     162
    163163  num_3d_m        = num_moist
    164164  num_3d_c        = num_chem
Note: See TracChangeset for help on using the changeset viewer.