Changeset 1662 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

Feb 8, 2017, 5:27:32 PM (8 years ago)

Author:

aslmd

Message:

added mars=32 in initialize_real

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F (modified) (4 diffs)

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

@@ -1334 +1334 @@
                                     its , ite , jts , jte , kts , kte )
 endif        
-if ( ( config_flags%mars == 1 ) .OR. ( config_flags%mars == 11 ) .OR. ( config_flags%mars == 12 ) ) then
+if ( ( config_flags%mars == 1 ) .OR. ( config_flags%mars == 11 ) .OR. ( config_flags%mars == 12 ) .OR. ( config_flags%mars == 32 ) ) then
        print *, '**** INTERPOLATE HV **** RANK 2 in SCALAR'
        !print *, size(scalar(0,0,0,:)), P_QH2O, P_QH2O_ICE
@@ -1364 +1364 @@
                                     its , ite , jts , jte , kts , kte )
 endif
-if ( (config_flags%mars == 11) .OR. (config_flags%mars == 12) ) then
+if ( (config_flags%mars == 11) .OR. (config_flags%mars == 12) .OR. (config_flags%mars == 32) ) 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 , &
@@ -1382 +1382 @@
                                     its , ite , jts , jte , kts , kte )
 endif
-if ( config_flags%mars == 12 ) then
+if ( (config_flags%mars == 12) .OR. (config_flags%mars == 32) ) then
        print *, '**** INTERPOLATE CCNQ **** RANK 6 in SCALAR'
        CALL vert_interp_old ( grid%em_ccnq_gc , grid%em_pd_gc , scalar(:,:,:,6) , grid%em_pb , &
@@ -1399 +1399 @@
                                     ims , ime , jms , jme , kms , kme , &
                                     its , ite , jts , jte , kts , kte )
+endif
+if ( (config_flags%mars == 32) ) then
+       print *, '**** INTERPOLATE CO2 **** RANK 8 in SCALAR'
+       CALL vert_interp_old ( grid%em_co2_gc , grid%em_pd_gc , scalar(:,:,:,8) , 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 )
+       print *, '**** set other CO2 tracers to 0'
+       scalar(:,:,:,9) = 0.
+       scalar(:,:,:,10) = 0.
+       scalar(:,:,:,11) = 0.
 endif
 #endif