Changeset 772 for trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars

Timestamp:

Sep 5, 2012, 5:12:50 PM (13 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.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars

Files:

• Registry.EM (modified) (3 diffs)
• module_first_rk_step_part1.F (modified) (1 diff)
• module_initialize_les.F (modified) (5 diffs)

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

@@ -114 +114 @@
 state  real  QH2O      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "QH2O"       "Water vapor mixing ratio"   "kg kg-1"
 state  real  QH2O_ICE  ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "QH2O_ICE"   "Water ice mixing ratio"     "kg kg-1"
-state  real  QDUST     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "QDUST"      "Dust mixing ratio"          "kg kg-1"
+state  real  qdust     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "qdust"      "Dust mixing ratio"          "kg kg-1"
+state  real  qndust    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "qndust"     "Dust number density"        " m-3 "
 state  real  qtrac1    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "qtrac1"     "Decaying tracer 1"          "kg kg-1"
 state  real  upward    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "upward"     "Decaying tracer surf"       "kg kg-1"
@@ -232 +233 @@
 state  real   mars_dsoil     ilj   misc  1  -  i012rd  "MARS_DSOIL"  "soil depths"             "m"
 state  real   mars_gw        ilj   misc  1  -  i012rd  "MARS_GW"     "gwparam"                 " "
+state  real   mars_qdust     ikj   misc  1  -  i012rd  "MARS_QDUST"   "dust mass mixing ratio" "kg/kg"
+state  real   mars_qndust    ikj   misc  1  -  i012rd  "MARS_QNDUST"   "dust number density"    " "
 ####
 ####
@@ -1496 +1499 @@
 package   water        mars==1                      -              moist:qv;scalar:qh2o,qh2o_ice
 package   dust         mars==2                      -              moist:qv;scalar:qdust
+package   dust2eq      mars==3                      -              scalar:qdust,qndust
 package   radioac      mars==20                     -              scalar:qtrac1
 package   radioac2     mars==21                     -              scalar:upward,downward

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_first_rk_step_part1.F

@@ -631 +631 @@
      &        ,MARS_Z0=grid%mars_z0             &
      &        ,CST_Z0=config_flags%init_Z0      &
+     &        ,MARS_QDUST=grid%mars_qdust       &
+     &        ,MARS_QNDUST=grid%mars_qndust     &
 #endif
      &        ,MARS_GW=grid%mars_gw                                       &

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

@@ -114 +114 @@
 !!MARS
  REAL :: lon_input, lat_input, alt_input, tsurf_input
+ ! for mode 3
+ LOGICAL :: MARS_MODE_INIT
+ REAL, DIMENSION(nl_max) :: qdust,qndust,dust_p_level
 !!MARS
 
@@ -346 +349 @@
     grid%mars_tsurf(i,j)=0. 
 !! <<
+!! >> Used for mode 3:
+    grid%mars_qdust(i,:,j)=0.
+    grid%mars_qndust(i,:,j)=0.
+!! <<
     grid%slpx(i,j) = 0.
     grid%slpy(i,j) = 0.
@@ -690 +697 @@
   print *, 'constant meridional wind ', config_flags%init_V
 ENDIF
+
 !!!MARS MARS
 
@@ -754 +762 @@
   ENDDO
   ENDDO
+
+!!!!! MARS INITIALIZE MODE 3 FOR LES
+
+!! RESTART=config_flags%restart -> this routine is not used in restart
+
+  MARS_MODE_INIT=config_flags%mars
+
+  IF (MARS_MODE_INIT .EQ. 3)  THEN
+  write (*,*) 'MARS MODE 3 INITIALIZATION, READING INPUT_DUST (module_initialize_les.F)'
+    ! load a profile of dust (same for all points)
+    call read_dust(qdust,qndust,dust_p_level,nl_in)
+    p_level = grid%znu(1)*(pd_surf - grid%p_top) + grid%p_top
+      IF (dust_p_level(1) .lt. p_level) 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_level/dust_p_level(1)
+     ENDIF
+     DO k=1,kte-1
+         p_level = grid%znu(k)*(pd_surf - grid%p_top) + grid%p_top
+         grid%mars_qdust(its:ite,k,jts:jte)=interp_0( qdust, dust_p_level, p_level, nl_in )
+         grid%mars_qndust(its:ite,k,jts:jte)=interp_0( qndust, dust_p_level, p_level, nl_in )
+     ENDDO
+  ENDIF
+
+!!!!! MARS 
+
 
  END SUBROUTINE init_domain_rk
@@ -1008 +1040 @@
       end subroutine read_therm
 
+!!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_initialize_ideal