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

Timestamp:

Sep 5, 2012, 5:12:50 PM (12 years ago)

Author:

acolaitis

Message:

LES. Reorganized mode 3. dust_profile is now read by ideal.exe and stored in mars_qdust and mars_qndust. Added -mcmodel for large executable in the debug options of makegcm.

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (7 diffs)

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

@@ -41 +41 @@
         MARS_Z0, &
         CST_Z0, &
+        MARS_QDUST, &
+        MARS_QNDUST, &
 #endif
         MARS_GW, &
@@ -69 +71 @@
    USE module_model_constants
    USE module_wrf_error
-   USE module_init_utilities
 ! add new modules here, if needed ...
 
@@ -145 +146 @@
 REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT   )  :: &
      MARS_FLUXRAD,MARS_WSTAR
+REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN   )  :: &
+     MARS_QDUST, MARS_QNDUST
 #endif
 REAL, DIMENSION( ims:ime, 5, jms:jme ), INTENT(IN   )  :: &
@@ -189 +192 @@
    INTEGER :: lctindex,ziindex
    LOGICAL :: end_of_file
-   ! *** ----------------------- ***
-
-   ! *** for Mars Mode 3         ***
-   INTEGER nm3
-   PARAMETER(nm3=1000)
-   REAL, DIMENSION(nm3) :: qdust,qndust,dust_p_level
    ! *** ----------------------- ***
 
@@ -871 +868 @@
 !--------------------------------!
 #ifdef NEWPHYS
+!!! Initialize SCALAR for MODE 3.
+IF ((MARS_MODE .EQ. 3) .and. (firstcall .EQV. .true.) .and. (.not. restart)) THEN
+     SCALAR(i,kps:kpe,j,2)=MARS_QDUST(i,kps:kpe,j)
+     SCALAR(i,kps:kpe,j,3)=MARS_QNDUST(i,kps:kpe,j)
+ENDIF
 IF (MARS_MODE .EQ. 0) THEN
     q_prof(:,1)=0.95
@@ -884 +886 @@
       q_prof(:,:) = 0.95
    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
@@ -1717 +1707 @@
 !!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