Changeset 126 for trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2

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

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

@@ -155 +155 @@
    INTEGER ::    i,j,k,its,ite,jts,jte,ij
    INTEGER ::    subs,iii
+   REAL ::    tau_decay
 
    ! *** for LMD physics
@@ -620 +621 @@
 
 SELECT CASE (MARS_MODE) !! ONLY ALLOW FOR MODES DEFINED IN Registry.EM
-   CASE(4-10,12:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
+   CASE(4-10,12-19,21:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
    PRINT *, 'NOT SUPPORTED, to be done'
    STOP
@@ -630 +631 @@
 !package   dust2eq      mars==3                      -              scalar:qdust,qdustn
 !package   newwater     mars==11                     -              scalar:qh2o,qh2o_ice,qdust,qdustn
+!package   radioac      mars==20                     -              scalar:qtrac1
 !!!!!!!!!!!!!!!!!!! FOR REFERENCE
 
@@ -652 +654 @@
       wtnom(3)  = 'dust_mass'
       wtnom(4)  = 'dust_number'
+    CASE(20) 
+      wtnom(1) = 'qtrac1'
 END SELECT
 #endif
@@ -696 +700 @@
    IF ( wtnom(iii) .eq. 'co2' ) q_prof(:,iii) = 0.95
   ENDDO
+
+    !! Mars mode 20 : add a passive tracer with radioactive-like decay
+    !! INIT HERE
+    IF (firstcall .EQV. .true.) THEN
+    IF (MARS_MODE .EQ. 20) THEN
+      q_prof(:,:) = 0.95
+    ENDIF
+    ENDIF
+
 #else
 SELECT CASE (MARS_MODE)
@@ -840 +853 @@
     qsurf_val(2)=MARS_WICE(i,j)    !! logique avec wtnom(2) = 'h2o_ice' defini ci-dessus
     qsurf_val(3)=0.                !! not coupled with lifting for the moment [non remobilise]
+    CASE(20)
+    qsurf_val(:)=0.
 #else
     CASE(3:)
@@ -1195 +1210 @@
 #ifdef NEWPHYS
 SCALAR(i,kps:kpe,j,1)=0.
-SCALAR(i,kps:kpe,j,2:nq+1)=SCALAR(i,kps:kpe,j,2:nq+1)+pdq(subs,kps:kpe,1:nq)*dt  !!! here dt is needed
+!! Mars mode 20 : add a passive tracer with radioactive-like decay
+IF (MARS_MODE .EQ. 20) THEN
+   print *, 'RADIOACTIVE-LIKE TRACER WITH SOURCE AT SURFACE LAYER.'
+   tau_decay=60.*10. !! why not make it a namelist argument?
+   SCALAR(i,kps:kpe,j,2) = SCALAR(i,kps:kpe,j,2)*exp(-dt/tau_decay)
+   SCALAR(i,1,j,2) = SCALAR(i,1,j,2) + 1. !! this tracer is emitted in the surface layer
+ELSE
+   SCALAR(i,kps:kpe,j,2:nq+1)=SCALAR(i,kps:kpe,j,2:nq+1)+pdq(subs,kps:kpe,1:nq)*dt  !!! here dt is needed
+ENDIF
 #else
 SELECT CASE (MARS_MODE)