Changeset 72 for trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

Timestamp:

Feb 19, 2011, 5:36:31 PM (14 years ago)

Author:

aslmd

Message:

LMD_MM_MARS: corrections pour runs cycle de l'eau. inclusion des traceurs de
surface [manquement grave au modele precedent] et importation depuis les runs
GCM. attention water ice est q(:,nqmx-1) mais surface water ice est
qsurf(:,nqmx).
NB: il reste a impacter ces changements pour le cas NEWPHYS

M 71 SRC/WRFV2/Registry/Registry.EM
M 71 SRC/WRFV2/dyn_em/module_initialize_real.F
M 71 SRC/WRFV2/dyn_em/solve_em.F
M 71 SRC/WRFV2/dyn_em/module_initialize_quarter_ss.F
M 71 SRC/WRFV2/phys/module_lmd_driver.F
M 71 SRC/WRFV2/mars_lmd/libf/phymars/meso_physiq.F
Modifies pour initialiser, calculer et sortir qsurf dans le cas mars=2
[cycle de l'eau avec deux traceurs vapour et ice]

M 71 SRC/PREP_MARS/readmeteo.F90
M 71 SRC/WPS/wps_mars/metgrid/METGRID.TBL.ARW_MarsBase
Routines d'initialisation modifiees pour qsurf.

D 71 SRC/PREP_MARS/readmeteo_newphys.F90
D 71 SRC/PREP_MARS/readmeteo.F90_backup
M 71 SRC/PREP_MARS/readmeteo.F90
A + - SRC/PREP_MARS/obsolete/readmeteo_newphys.F90
A + - SRC/PREP_MARS/obsolete/readmeteo.F90_backup
M 71 SRC/WPS/wps_mars/metgrid/METGRID.TBL.ARW_MarsBase_newphys
M 71 SRC/WPS/wps_mars/switch_WPS_Mars
Desormais les programmes "readmeteo" et "metgrid" fonctionnent
indifferemment de la physique LMD utilisee (ancienne ou nouvelle)
--> une seule version de reference
--> retrocompatible pour l'ancienne physique

(les parametres inconnus sont regles a des valeurs "dummy" et ne sont pas utilisees)

A 0 SIMU/deftank/dust
Ajout mineur de fichiers parametres mais qui sera utile pour le cycle de poussiere

File:

• trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (8 diffs)

trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

@@ -30 +30 @@
         MARS_MODE, &
         MARS_ALB,MARS_TI,MARS_CICE,MARS_EMISS, &
+        MARS_WICE, &
         MARS_TSOIL, &
 #ifdef NEWPHYS
@@ -109 +110 @@
      XLAT,XLONG,HT,  &
      MARS_ALB,MARS_TI,MARS_EMISS,MARS_CICE, &
+     MARS_WICE, &
      SLPX,SLPY
 ! 3D arrays 
@@ -440 +442 @@
 !ENDIF   
 
-print *, 'check dynamics'
-  !!! in some cases, weird values are displayed
-  !!! despite the fact that outputs are OK...
-  !print *, 'u', MAXVAL(u), MINVAL(u)
-  !print *, 'v', MAXVAL(v), MINVAL(v)
-  !print *, 'w', MAXVAL(w), MINVAL(w)
-  !print *, 't', MAXVAL(t), MINVAL(t, MASK = t > 0)
-print *, 'u', u(10,1,10), u(10,15,10)
-print *, 'v', v(10,1,10), v(10,15,10)
-print *, 'w', w(10,1,10), w(10,15,10) 
-print *, 't', t(10,1,10), t(10,15,10)
+!print *, 'check dynamics'
+!  !!! in some cases, weird values are displayed
+!  !!! despite the fact that outputs are OK...
+!  !print *, 'u', MAXVAL(u), MINVAL(u)
+!  !print *, 'v', MAXVAL(v), MINVAL(v)
+!  !print *, 'w', MAXVAL(w), MINVAL(w)
+!  !print *, 't', MAXVAL(t), MINVAL(t, MASK = t > 0)
+!print *, 'u', u(10,1,10), u(10,15,10)
+!print *, 'v', v(10,1,10), v(10,15,10)
+!print *, 'w', w(10,1,10), w(10,15,10) 
+!print *, 't', t(10,1,10), t(10,15,10)
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -786 +788 @@
 
 !-------------------!
+! Tracer at surface !
+!-------------------!
+#ifdef NEWPHYS
+  !!! a faire !!!
+  PRINT *, 'WARNING WARNING no tracer at surface'
+  qsurf_val(:)=0.
+#else
+SELECT CASE (MARS_MODE)
+    CASE(0)  !! NO TRACERS (mars=0)
+    qsurf_val(:)=0.
+    CASE(1)  !! WATER CYCLE (mars=1)
+    qsurf_val(2)=MARS_WICE(i,j)    !! attention... H2O ice is tracer nqmx in qsurf in LMD physics
+    qsurf_val(1)=0.                
+    CASE(2)  !! DUST CYCLE (mars=2)
+    qsurf_val(:)=0.
+END SELECT
+#endif
+
+!-------------------!
 ! Slope inclination !
 !-------------------!
@@ -832 +853 @@
 PRINT *,'v_prof ',v_prof
 PRINT *,'tsoil ',tsoil_val 
+PRINT *,'qsurf ',qsurf_val 
 #ifdef NEWPHYS
 PRINT *,'isoil ',isoil_val
@@ -844 +866 @@
 !-------------------------!
 !-------------------------!
-q2_val(:)=0      !PBL wind variance
-qsurf_val(:)=0   !Tracer on surface
+q2_val(:)=0.      !PBL wind variance
 
 !-----------------!
@@ -1167 +1188 @@
 !!PRINT *, 'min',MINVAL(RUBLTEN, MASK=RUBLTEN/=0.),&
 !!        ' at',MINLOC(RUBLTEN, MASK=RUBLTEN/=0.)
-PRINT *, RUBLTEN(10,1,10), RUBLTEN(10,15,10)
+!PRINT *, RUBLTEN(10,1,10), RUBLTEN(10,15,10)
 PRINT *, 'v non-zero tendencies'
 !!PRINT *, 'max',MAXVAL(RVBLTEN, MASK=RVBLTEN/=0.),&
@@ -1173 +1194 @@
 !!PRINT *, 'min',MINVAL(RVBLTEN, MASK=RVBLTEN/=0.),&
 !!        ' at',MINLOC(RVBLTEN, MASK=RVBLTEN/=0.)
-PRINT *, RVBLTEN(10,1,10), RVBLTEN(10,15,10)
+!PRINT *, RVBLTEN(10,1,10), RVBLTEN(10,15,10)
 !!! STOP IF CRASH
 !IF (MAXVAL(RUBLTEN, MASK=RUBLTEN/=0.) == 0.) STOP