Changeset 766 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

Aug 31, 2012, 3:22:47 PM (12 years ago)

Author:

acolaitis

Message:

LES. Added mars mode 3 (active dust with dust_mass and dust_number) and rectified pd_scalar default state.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC

Files:

• LES/modif_mars/Registry.EM (modified) (1 diff)
• WRFV2/phys/module_lmd_driver.F (modified) (5 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM

@@ -1400 +1400 @@
 rconfig   logical     pd_moist_dfi        namelist,dynamics     max_domains    .false.  rh    "pd_moist_dfi"    "positive-definite RK3 transport switch"      ""
 rconfig   logical     pd_chem             namelist,dynamics     max_domains    .false. rh    "pd_chem"          "positive-definite RK3 transport switch"      ""
-rconfig   logical     pd_scalar           namelist,dynamics     max_domains    .false. rh    "pd_scalar"        "positive-definite RK3 transport switch"      ""
+rconfig   logical     pd_scalar           namelist,dynamics     max_domains    .true. rh    "pd_scalar"        "positive-definite RK3 transport switch"      ""
 rconfig   logical     pd_tke              namelist,dynamics     max_domains    .false. rh    "pd_tke"           "positive-definite RK3 transport switch"      ""
 rconfig   logical top_radiation           namelist,dynamics     max_domains    .false. rh    "top_radiation"         ""      ""

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

@@ -69 +69 @@
    USE module_model_constants
    USE module_wrf_error
+   USE module_init_utilities
 ! add new modules here, if needed ...
 
@@ -188 +189 @@
    INTEGER :: lctindex,ziindex
    LOGICAL :: end_of_file
+   ! *** ----------------------- ***
+
+   ! *** for Mars Mode 3         ***
+   INTEGER nm3
+   PARAMETER(nm3=1000)
+   REAL, DIMENSION(nm3) :: qdust,qndust,dust_p_level
    ! *** ----------------------- ***
 
@@ -878 +885 @@
    ENDIF
   ENDIF
+  !!! DUST AND DUSTN
+  IF (MARS_MODE .EQ. 3 .and. firstcall .EQV. .true. .and. (.not. restart)) THEN
+      call read_dust(qdust,qndust,dust_p_level,nm3)
+      if (dust_p_level(1) .lt. p_prof(1)) then !input profile needs rescaling to avoid a plateau. This happens when you use different sources to initialize the pressure of the LES and the input_dust. Typicaly: you use the MCD for input_therm and different runs for input_dust, that dont have the same surface pressure ! This trick is ok because we dont want to initialize with a very precise profile of dust, just a realistic one.
+         dust_p_level = dust_p_level * p_prof(1)/dust_p_level(1)         
+      endif
+      DO k=1, nlayer
+        q_prof(k,1) = interp_0(qdust, dust_p_level, p_prof(k), nm3)
+        q_prof(k,2) = interp_0(qndust, dust_p_level, p_prof(k), nm3)
+      ENDDO
+      SCALAR(i,kps:kpe,j,2:3)=q_prof(:,:)  ! because we need to initialize the value of scalar! 
+  ENDIF
 #else
 SELECT CASE (MARS_MODE)
@@ -1057 +1076 @@
 #ifdef NEWPHYS
     CASE(3)
-    qsurf_val(:)=0.                !!! temporaire, a definir            
+    qsurf_val(1)=q_prof(1,1)                !!! temporaire, a definir       
+    qsurf_val(2)=q_prof(1,2)     
     CASE(11) 
     qsurf_val(1)=0.
@@ -1697 +1717 @@
 !!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 
+
+!!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine read_dust(qdust,qndust,dust_p_level,n)
+      implicit none
+      integer n
+      real qdust(n),qndust(n),dust_p_level(n)
+      logical end_of_file
+
+      integer k
+
+! first element is the surface
+
+      open(unit=11,file='input_dust',form='formatted',status='old')
+      rewind(11)
+      end_of_file = .false.
+      k = 0
+      do while (.not. end_of_file)
+
+        read(11,*,end=102) qdust(k+1), qndust(k+1), dust_p_level(k+1)
+        write(*,*) k, qdust(k+1), qndust(k+1), dust_p_level(k+1)
+        k = k+1
+        go to 113
+ 102    end_of_file = .true.
+ 113    continue
+      enddo
+
+      close(unit=11,status = 'keep')
+
+      end subroutine read_dust
+!!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
 END MODULE module_lmd_driver