Changeset 55 for trunk/mars/libf

Timestamp:

Feb 4, 2011, 1:46:39 AM (14 years ago)

Author:

aslmd

Message:

LMD_MM_MARS: element manquant pour runs traceurs avec nouvelle physique
--> avec mars=0 et non actif, le modele fonctionne toujours bien [avec starts JBM et conditions de ses runs]
--> mettre les traceurs avec mars=11 [actifs ou non] semble faire crasher le modele: pourquoi ?

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq1.inc
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq2.inc
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq3.inc
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq4.inc
M 54 mars/libf/phymars/meso_physiq.F
le tableau tnom est defini dans module_lmd_driver en fonction de config_flag%mars [MARS_MODE]

et passe dans meso_physiq [pour ensuite servir dans initracer]

--> le COMMON advtrac.h est alors necessaire

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
config_flag%mars defini une serie de traceurs ordonnee dans le tableau scalar

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/share/module_model_constants.F
on passe au cp et R du GCM martien pour une complete correspondance

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq.inc
ce fichier est en fait inutile et pourrait etre supprime

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/main/real_em.F
modifications pour prise en compte des traceurs avec la nouvelle physique
--> avantage des traceurs generiques dans la physique
--> pour l'instant dans SCALAR l'indice 1 est dummy, 2/3 water vapor/ice et dernier CO2

M 54 mesoscale/LMD_MM_MARS/SRC/PREP_MARS/readmeteo_newphys.F90
M 54 mesoscale/LMD_MM_MARS/SRC/WPS/wps_mars/metgrid/METGRID.TBL.ARW_MarsBase_newphys
modifications pour possibilites de passer les tableaux necessaires pour

initialiser et guider les traceurs dans la nouvelle physique

M mars/libf/phymars/dimradmars.h
M mars/libf/phymars/callradite.F
version traceurs non actifs pour tests basiques

Location:

trunk/mars/libf/phymars

Files:

• callradite.F (modified) (1 diff)
• dimradmars.h (modified) (1 diff)
• meso_physiq.F (modified) (5 diffs)

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
-          name_iaer(2) = "h2o_ice"
+c          name_iaer(2) = "h2o_ice"
 c        ----------------------------------------------------------
 

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=2)
+      parameter (naerkind=1)
+!      parameter (naerkind=2)
 
 ! Reference wavelengths used to compute reference optical depth (m)

trunk/mars/libf/phymars/meso_physiq.F

@@ -4 +4 @@
      $            pday,ptime,ptimestep,
      $            pplev,pplay,pphi,
-     $            pu,pv,pt,pq,
-     $            pw,
+     $            pu,pv,pt,pq,pw,
+     $            wtnom,
      $            pdu,pdv,pdt,pdq,pdpsrf,tracerdyn,
      $            wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,
@@ -156 +156 @@
 #include "wrf_output_3d.h"
 
+#include "advtrac.h"   !!! this is necessary for tracers (in dyn3d)
 
 c Arguments :
@@ -193 +194 @@
       real alpha,lay1 ! coefficients for building layers
       integer iloop
+      INTEGER tracerset    !!! this corresponds to config%mars
 c ****WRF
       REAL pday
       REAL ptime 
       logical tracerdyn
+      CHARACTER (len=20) :: wtnom(nqmx) ! tracer name
 
 c   outputs:
@@ -392 +395 @@
             enddo
             !!!!!!!!!!!!!!!!! DONE in soil_setting.F
+      tnom(:)=wtnom(:)   !! est rempli dans advtrac.h
+      PRINT*,'check: tracernames ', tnom
      !!!new physics
      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -422 +427 @@
               !!donc ecritphys est passe en argument ...
       PRINT*,'Write LMD physics each:',ecritphy,' seconds'
+c
+      DO iq=1, nq 
+        PRINT*, tnom(iq), pq(:,:,iq)
+      ENDDO
 c
 c ****WRF