Changeset 126 for trunk/mesoscale/LMD_LES_MARS/modif_mars

Timestamp:

May 23, 2011, 8:40:24 PM (14 years ago)

Author:

acolaitis

Message:

LMD_LES_MARS et LMD_MM_MARS

M 125 000-USERS
Entree d'Arnaud Colaitis chez les committers.

M 125 mars/libf/phymars/dimradmars.h
M 125 mars/libf/phymars/callradite.F
Cas simple conrath sans doubleq OK avec LES par defaut

M 125 mesoscale/LMD_LES_MARS/modif_mars/module_initialize_les.F
Initialisation du LES a partir de profils de temperature potentielle
calcules avec r_cp de module_initialize_les
[en lieu et place du calcul effectue dans wrf_sounding.sci... risque de mismatch !]

M 125 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
M 125 mesoscale/LMD_LES_MARS/modif_mars/Registry.EM
Mars mode 20 : add a passive tracer with radioactive-like decay

A 0 mesoscale/LMD_LES_MARS/modif_mars/makeles_gnome
A + - mesoscale/TESTS/newphys_tracers/callphys.def.mars0.ideal
Modifications mineures, ajout de fichiers pour cas specifiques

Location:

trunk/mesoscale/LMD_LES_MARS/modif_mars

Files:

• Registry.EM (modified) (2 diffs)
• makeles_gnome (added)
• module_initialize_les.F (modified) (3 diffs)

trunk/mesoscale/LMD_LES_MARS/modif_mars/Registry.EM

@@ -118 +118 @@
 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  qtrac1    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "qtrac1"     "Decaying tracer 1"          "kg kg-1"
 ####
 ####
@@ -1488 +1489 @@
 package   water        mars==1                      -              moist:qv;scalar:qh2o,qh2o_ice
 package   dust         mars==2                      -              moist:qv;scalar:qdust
+package   radioac      mars==20                     -              scalar:qtrac1
 ##### MARS OPTIONS
 ##### MARS OPTIONS

trunk/mesoscale/LMD_LES_MARS/modif_mars/module_initialize_les.F

@@ -798 +798 @@
       real h_input(n), th_input(n), qv_input(n), u_input(n), v_input(n)
 
+! input therm data (element 0 is the ground so it's n+1 but n is 1000 anyway so...)
+
+     real r_therm(n),cp_therm(n),p_therm(n),rho_therm(n),t_therm(n)
+
 ! diagnostics
 
@@ -816 +820 @@
       call read_sounding( p_surf, th_surf, qv_surf, &
                           h_input, th_input, qv_input, u_input, v_input,n, nl, debug )
+
+! and the therm :
+
+      call read_therm(r_therm,cp_therm,p_therm,rho_therm,t_therm,n)
+
+! To use r/cp as defined above, one has to recompute teta from T (default MCD computes
+! teta for a variable r/cp)
+
+      do k=1,nl
+        th_input(k) = t_therm(k)*(p1000mb/p_therm(k))**(r/cp)
+      enddo
+      th_surf = t_therm(1)*(p1000mb/p_therm(1))**(r/cp)
+! -----
 
       if(dry) then
@@ -957 +974 @@
       end subroutine read_sounding
 
+      subroutine read_therm(r,cp,p,rho,t,n)
+      implicit none
+      integer n
+      real r(n),cp(n),p(n),rho(n),t(n)
+      logical end_of_file
+
+      integer k
+
+! first element is the surface
+
+      open(unit=11,file='input_therm',form='formatted',status='old')
+      rewind(11)
+      end_of_file = .false.
+      k = 0
+      do while (.not. end_of_file)
+
+        read(11,*,end=101) r(k+1), cp(k+1), p(k+1), rho(k+1), t(k+1)
+        write(*,*) k, r(k+1), cp(k+1), p(k+1), rho(k+1), t(k+1)
+        k = k+1
+        go to 112
+ 101    end_of_file = .true.
+ 112    continue
+      enddo
+
+      close(unit=11,status = 'keep')
+
+      end subroutine read_therm
+
 END MODULE module_initialize_ideal