Changeset 126

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

Files:

• 000-USERS (modified) (1 diff)
• mars/libf/phymars/callradite.F (modified) (1 diff)
• mars/libf/phymars/dimradmars.h (modified) (2 diffs)
• mesoscale/LMD_LES_MARS/modif_mars/Registry.EM (modified) (2 diffs)
• mesoscale/LMD_LES_MARS/modif_mars/makeles_gnome (added)
• mesoscale/LMD_LES_MARS/modif_mars/module_initialize_les.F (modified) (3 diffs)
• mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (7 diffs)
• mesoscale/TESTS/newphys_tracers/callphys.def.mars0.ideal (copied) (copied from trunk/mesoscale/TESTS/newphys_tracers/callphys.def.mars0) (4 diffs)

trunk/000-USERS

@@ -6 +6 @@
 Ehouarn Millour                 LMD [CNRS]                      start                           All GCMs
 Sebastien Lebonnois             LMD [CNRS]                      start                           Venus and Titan GCMs
-Aymeric Spiga                   LMD [UPMC]                      17 - 11 - 2011                  Mars GCM, MMM and LES
+Aymeric Spiga                   LMD [UPMC]                      17 - 11 - 2010                  Mars GCM, MMM and LES
+Arnaud Colaitis                 LMD [CNRS]                      23 - 05 - 2011                  Mars GCM, MMM and LES
 
 ------------------------------------------------------------------------------------------------------------------------

trunk/mars/libf/phymars/callradite.F

@@ -263 +263 @@
 c        PLEASE MAKE SURE that you set up the right number of
 c          scatterers in dimradmars.h (naerkind);
-c         name_iaer(1) = "dust_conrath"   !! poussiere classique
-          name_iaer(1) = "dust_doubleq"
+         name_iaer(1) = "dust_conrath"   !! poussiere classique
+c          name_iaer(1) = "dust_doubleq"
 cc        name_iaer(2) = "dust_submicron" !! JB: experimental
 c          name_iaer(2) = "h2o_ice"

trunk/mars/libf/phymars/dimradmars.h

@@ -27 +27 @@
 ! (ex: naerkind=2 if you use one dust mode and active ice ...)
       integer naerkind
-      parameter (naerkind=1)
+       parameter (naerkind=1)
 !      parameter (naerkind=2)
+!!!!!! CHANGE CHANGE nsizemax BELOW BELOW
 
 ! Reference wavelengths used to compute reference optical depth (m)
@@ -76 +77 @@
 !   the radiative calculations.
 
-      INTEGER, PARAMETER :: nsizemax = 60
-!      INTEGER, PARAMETER :: nsizemax = 1
+!      INTEGER, PARAMETER :: nsizemax = 60
+      INTEGER, PARAMETER :: nsizemax = 1
 
 ! Various initialisation for LW radiative code

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

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)

trunk/mesoscale/TESTS/newphys_tracers/callphys.def.mars0.ideal

@@ -1 +1 @@
 ##General options
 ##~~~~~~~~~~~~~~~
+#Directory where external input files are
+datadir = /san/work/colaitis/datafile/ 
+
 #Run with or without tracer transport ?
 tracer=.false.
@@ -26 +29 @@
 #  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
 #  =26 Mars Year 26 from TES assimilation
-iaervar = 26 
+iaervar = 1 
 # Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
-tauvis=0.2
+tauvis = 0.5
 # Dust vertical distribution: 
 # (=0: old distrib. (Pollack90), =1: top set by "topdustref",
@@ -43 +46 @@
 callnlte = .false.
 # call CO2 NIR absorption ?   matters only if callrad=T
-callnirco2 = .true.
+callnirco2 = .false.
 # call turbulent vertical diffusion ?
 calldifv  = .true.
 # call convective adjustment ?
-calladj  = .true.
+calladj  = .false.
 # call CO2 condensation ?
 callcond =.true.
@@ -55 +58 @@
 calllott  = .false.
 # Impose polar cap surface albedos as observed by TES?
-TESicealbedo = .true.
+TESicealbedo = .false.
 
 ## Radiative transfer options :