Changeset 488 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

Dec 21, 2011, 3:15:08 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: finally adopted other strategy for mars=12. ccn tracers are initialized to 0 even at boundaries and we adopt same treatment for those than other tracers (flux from boundary conditions at boundaries).

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Files:

• Registry/Registry.EM (modified) (1 diff)
• dyn_em/module_initialize_real.F (modified) (2 diffs)
• dyn_em/solve_em.F (modified) (2 diffs)

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

@@ -129 +129 @@
 state  real  qdust     ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QDUST"         "Dust mixing ratio"          "kg kg-1"
 state  real  qdustn    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QDUSTN"        "Dust_number mixing ratio"   "kg kg-1"
-state  real  qccn     ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QCCN"         "CCN mixing ratio"          "kg kg-1"
-state  real  qccnn    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QCCNN"        "CCN number mixing ratio"   "kg kg-1"
+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  chem_co   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_co"     ""   ""

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

@@ -1313 +1313 @@
 endif
 #ifdef NEWPHYS
-if ( config_flags%mars == 11 .OR. config_flags%mars == 12 ) then
+if ( (config_flags%mars == 11) .OR. (config_flags%mars == 12) ) then
        print *, '**** INTERPOLATE DUSTQ **** RANK 4 in SCALAR'             
        CALL vert_interp_old ( grid%em_dustq_gc , grid%em_pd_gc , scalar(:,:,:,4) , grid%em_pb , &
@@ -1330 +1330 @@
                                     ims , ime , jms , jme , kms , kme , &
                                     its , ite , jts , jte , kts , kte )
+endif
+if ( config_flags%mars == 12 ) then
+       print *, "**** mars=12... CCN TRACERS SET TO 0"
+       DO k = 1, kte
+       DO j = jts, MIN(jte,jde-1)
+       DO i = its, MIN(ite,ide-1)
+         scalar(i,k,j,6) = 1.E-30
+         scalar(i,k,j,7) = 1.E-30
+       ENDDO
+       ENDDO
+       ENDDO
 endif
 #endif

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

@@ -2709 +2709 @@
 ! .OR. ( (config_flags%mars .eq. 11) .AND. (is .EQ. 2) ) &    !! 11  pass only water vapor
  .OR. (  config_flags%mars .eq. 11                    ) &    !! 11  pass EVERYTHING
- .OR. (  (config_flags%mars .eq. 12) .AND. (is .LE. 5) ) &    !! 12  pass EVERYTHING except 2 ccn tracers
+ .OR. (  config_flags%mars .eq. 12                    ) &    !! 12  pass EVERYTHING 
  .OR. config_flags%nested ) THEN                             !! *   pass all tracers if nested
          CALL relax_bdy_scalar ( scalar_tend(ims,kms,jms,is),            &
@@ -2770 +2770 @@
 !IF (      ((config_flags%mars .eq. 1 ) .and. (is .ne. 2))  &   
 !     .OR. ((config_flags%mars .eq. 11) .and. (is .ne. 2))  ) THEN
-IF  ( ((config_flags%mars .eq. 12 ) .and. (is .ge. 6)) &
-      .OR. ( config_flags%mars .gt. 50 ) ) THEN
+IF ( config_flags%mars .gt. 50 ) THEN
 
      !!! YOU HAVE TO ADD A CONDITION HERE IF YOU ADD A mars OPTION in REGISTRY  !!!