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

Timestamp:

Oct 21, 2011, 11:31:28 AM (14 years ago)

Author:

aslmd

Message:

MESOSCALE : Preparatory commit for the ultimate option mars=42 which

would allow mesoscale modeling with photochemistry.

[see r76 method to add 'mars' options]
Modified module_lmd_driver.F and Registry.EM and runmeso
Modified readmeteo.F90 and introduced an option -DPHOTOCHEM
...
Transparent to the casual user
Option mars=42 not yet finished though -- so do not use!

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Files:

• Registry/Registry.EM (modified) (2 diffs)
• phys/module_lmd_driver.F (modified) (6 diffs)

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

@@ -130 +130 @@
 state  real  qdustn    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QDUSTN"        "Dust_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  chem_co   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_co"     ""   ""
+state  real  chem_o    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_o"      ""   ""
+state  real  chem_o1d  ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_o1d"    ""   ""
+state  real  chem_o2   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_o2"     ""   ""
+state  real  chem_o3   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_o3"     ""   ""
+state  real  chem_h    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_h"      ""   ""
+state  real  chem_h2   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_h2"     ""   ""
+state  real  chem_oh   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_oh"     ""   ""
+state  real  chem_ho2  ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_ho2"    ""   ""
+state  real  chem_h2o2 ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_h2o2"   ""   ""
+state  real  chem_ch4  ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_ch4"    ""   ""
+state  real  chem_n2   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_n2"     ""   ""
+state  real  chem_ar   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_ar"     ""   ""
 ####
 ####
@@ -1353 +1366 @@
 package   dust2eq      mars==3                      -              scalar:qdust,qdustn
 package   newwater     mars==11                     -              scalar:qh2o,qh2o_ice,qdust,qdustn
+package   photochem    mars==42                     -              scalar:qco2,chem_co,chem_o,chem_o1d,chem_o2,chem_o3,chem_h,chem_h2,chem_oh,chem_ho2,chem_h2o2,chem_ch4,chem_n2,chem_ar,qh2o_ice,qh2o,qdust,qdustn
 ##### MARS OPTIONS
 ##### MARS OPTIONS

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

@@ -650 +650 @@
 
 SELECT CASE (MARS_MODE) !! ONLY ALLOW FOR MODES DEFINED IN Registry.EM
-   CASE(4-10,12-19,22:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
+   CASE(4-10,12-19,22-41,43:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
    PRINT *, 'NOT SUPPORTED, to be done'
    STOP
@@ -662 +662 @@
 !package   radioac      mars==20                     -              scalar:qtrac1
 !package   radioac2     mars==21                     -              scalar:upward,downward
+!package   photochem    mars==42                     -              scalar:qco2,chem_co,chem_o,chem_o1d,chem_o2,chem_o3,chem_h,chem_h2,chem_oh,chem_ho2,chem_h2o2,chem_ch4,chem_n2,chem_ar,qh2o_ice,qh2o,qdust,qdustn
 !!!!!!!!!!!!!!!!!!! FOR REFERENCE
 
@@ -672 +673 @@
       wtnom(nq) = 'co2'
     CASE(1)
-      wtnom(1)  = 'h2o_vap'     
-      wtnom(2)  = 'h2o_ice'     
+      wtnom(1)  = 'h2o_vap'
+      wtnom(2)  = 'h2o_ice' 
     CASE(2)
       wtnom(1)  = 'dust01'      
     CASE(3)
       wtnom(1)  = 'dust_mass'
-      wtnom(2)  = 'dust_number'
+      wtnom(2)  = 'dust_number' 
     CASE(11) 
       wtnom(1)  = 'h2o_vap'
-      wtnom(2)  = 'h2o_ice'
+      wtnom(2)  = 'h2o_ice' 
       wtnom(3)  = 'dust_mass'
       wtnom(4)  = 'dust_number'
@@ -688 +689 @@
     CASE(21)
       wtnom(1) = 'upward'
-      wtnom(2) = 'downward'
+      wtnom(2) = 'downward' 
+    CASE(42)
+      wtnom(1)  = 'co2'
+      wtnom(2)  = 'co'
+      wtnom(3)  = 'o'
+      wtnom(4)  = 'o1d'
+      wtnom(5)  = 'o2'
+      wtnom(6)  = 'o3'
+      wtnom(7)  = 'h'
+      wtnom(8)  = 'h2'
+      wtnom(9)  = 'oh'
+      wtnom(10)  = 'ho2'
+      wtnom(11)  = 'h2o2'
+      wtnom(12)  = 'ch4'
+      wtnom(13)  = 'n2'
+      wtnom(14)  = 'ar'
+      wtnom(15)  = 'h2o_ice'
+      wtnom(16)  = 'h2o_vap'
+      wtnom(17)  = 'dust_mass'
+      wtnom(18)  = 'dust_number'
 END SELECT
 #endif
@@ -729 +749 @@
 q_prof(:,1:nq) = SCALAR(i,kps:kpe,j,2:nq+1)  !! the names were set above !! one dummy tracer in WRF
   !!! CAS DU CO2
-DO iii=1,nq 
- IF ( wtnom(iii) .eq. 'co2' ) q_prof(:,iii) = 0.95
-ENDDO
-
-IF ((MARS_MODE .EQ. 20) .OR. (MARS_MODE .EQ. 21)) THEN
+  DO iii=1,nq 
+   IF ( wtnom(iii) .eq. 'co2' ) q_prof(:,iii) = 0.95
+  ENDDO
+  IF ((MARS_MODE .EQ. 20) .OR. (MARS_MODE .EQ. 21)) THEN
    IF (firstcall .EQV. .true.) THEN
       q_prof(:,:) = 0.95
    ENDIF
-ENDIF
-
+  ENDIF
 #else
 SELECT CASE (MARS_MODE)
@@ -841 +859 @@
 !!!! ADDITIONAL SECURITY. THIS MIGHT HAPPEN WITH OLD INIT FILES.
 IF (z0_val == 0.) THEN
-   PRINT *, 'WELL, z0 is 0, this is no good. Setting to old defaults value 0.01 m'
+   IF ( (i == ips) .AND. (j == jps) ) PRINT *, 'WELL, z0 is 0, this is no good. Setting to old defaults value 0.01 m'
    z0_val = 0.01
 ENDIF