Changeset 802 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Timestamp:

Oct 5, 2012, 2:52:43 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE : adaptation for polar runs with new physics. co2 mr is now taken from GCM and given as an input (previously was initialized at 0.95). added mars=10 which is a case with CO2 being the only tracer.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Files:

• Registry/Registry.EM (modified) (3 diffs)
• dyn_em/module_initialize_real.F (modified) (1 diff)
• dyn_em/solve_em.F (modified) (1 diff)
• main/real_em.F (modified) (1 diff)
• phys/module_lmd_driver.F (modified) (2 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM ¶

@@ -132 +132 @@
 state  real  qccn      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCCN"         "CCN mixing ratio"          "kg kg-1"
 state  real  qccnn     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCCNN"        "CCN number mixing ratio"   "kg kg-1"
-state  real  qco2      ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QCO2"          "CO2 mixing ratio"           "kg kg-1"
+state  real  qco2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO2"          "CO2 mixing ratio"           "kg kg-1"
 state  real  chem_co   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_co"     ""   ""
 state  real  chem_o    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_o"      ""   ""
@@ -221 +221 @@
 state  real   isoil_gc       igj   dyn_em  1  Z  i1  "ISOIL"            "soil thermal inertia"    "tiu"
 state  real   dsoil_gc       igj   dyn_em  1  Z  i1  "DSOIL"            "soil depths"             "m"
+state  real   co2_gc         igj   dyn_em  1  Z  i1  "CO2"              "co2 mixing ratio"        "kg/kg"
 state  real   dustq_gc       igj   dyn_em  1  Z  i1  "DUSTQ"            "dust mixing ratio"       "kg/kg"
 state  real   dustn_gc       igj   dyn_em  1  Z  i1  "DUSTN"            "dust number density"     "part/kg"
@@ -1374 +1375 @@
 package   dust1        mars==2                      -              scalar:qdust
 package   dust2eq      mars==3                      -              scalar:qdust,qdustn
+package   newbasis     mars==10                     -              scalar:qco2
 package   newwater     mars==11                     -              scalar:qh2o,qh2o_ice,qdust,qdustn
 package   newwatersca  mars==12                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F ¶

@@ -1320 +1320 @@
 endif
 #ifdef NEWPHYS
+if ( (config_flags%mars == 10) ) then
+       print *, '**** INTERPOLATE CO2 **** RANK 2 in SCALAR'
+       CALL vert_interp_old ( grid%em_co2_gc , grid%em_pd_gc , scalar(:,:,:,2) , grid%em_pb , &
+                                    num_metgrid_levels , 'Q' , &
+                                    interp_type , lagrange_order , lowest_lev_from_sfc , &
+                                    zap_close_levels , force_sfc_in_vinterp , &
+                                    ids , ide , jds , jde , kds , kde , &
+                                    ims , ime , jms , jme , kms , kme , &
+                                    its , ite , jts , jte , kts , kte )
+endif
 if ( (config_flags%mars == 11) .OR. (config_flags%mars == 12) ) then
        print *, '**** INTERPOLATE DUSTQ **** RANK 4 in SCALAR'             

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F ¶

@@ -2711 +2711 @@
  .OR. (  config_flags%mars .eq.  2                    ) &    !! 2   pass the dust tracer
  .OR. (  config_flags%mars .eq.  3                    ) &    !! 3   pass dust Q and N 
+ .OR. (  config_flags%mars .eq.  10                   ) &    !! 10  pass the co2 tracer
 ! .OR. ( (config_flags%mars .eq. 11) .AND. (is .EQ. 2) ) &    !! 11  pass only water vapor
  .OR. (  config_flags%mars .eq. 11                    ) &    !! 11  pass EVERYTHING

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/main/real_em.F ¶

@@ -717 +717 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 IF (config_flags%mars .gt. 0) THEN
+!! ceci suffit avec les modes 2 et 10 qui ne transportent qu'un seul traceur
       CALL couple ( grid%em_mu_2 , grid%em_mub , qbdy3dtemp1 , grid%scalar(:,:,:,2) , 't' , grid%msft , &
                     ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, ips, ipe, jps, jpe, kps, kpe )

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F ¶

@@ -780 +780 @@
 !!!                               SAME NAMING   AS IN THE LMD PHYSICS !!!!
 SELECT CASE (MARS_MODE)
-    CASE(0)  
+    CASE(0,10)  
       wtnom(nq) = 'co2'
     CASE(1)
@@ -869 +869 @@
     q_prof(:,1:nq) = SCALAR(i,kps:kpe,j,2:nq+1)  !! the names were set above !! one dummy tracer in WRF
 ENDIF
-  !!! CAS DU CO2
-  DO iii=1,nq 
-   IF ( wtnom(iii) .eq. 'co2' .and. (.not. restart)) q_prof(:,iii) = 0.95
-  ENDDO
+  !!!! CAS DU CO2
+  !DO iii=1,nq 
+  ! IF ( wtnom(iii) .eq. 'co2' .and. (.not. restart)) q_prof(:,iii) = 0.95
+  !ENDDO
   IF ((MARS_MODE .EQ. 20) .OR. (MARS_MODE .EQ. 21)) THEN
    IF (firstcall .EQV. .true. .and. (.not. restart)) THEN