Ignore:
Timestamp:
Oct 5, 2012, 2:52:43 AM (12 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/dyn_em
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F

    r680 r802  
    13201320endif
    13211321#ifdef NEWPHYS
     1322if ( (config_flags%mars == 10) ) then
     1323       print *, '**** INTERPOLATE CO2 **** RANK 2 in SCALAR'
     1324       CALL vert_interp_old ( grid%em_co2_gc , grid%em_pd_gc , scalar(:,:,:,2) , grid%em_pb , &
     1325                                    num_metgrid_levels , 'Q' , &
     1326                                    interp_type , lagrange_order , lowest_lev_from_sfc , &
     1327                                    zap_close_levels , force_sfc_in_vinterp , &
     1328                                    ids , ide , jds , jde , kds , kde , &
     1329                                    ims , ime , jms , jme , kms , kme , &
     1330                                    its , ite , jts , jte , kts , kte )
     1331endif
    13221332if ( (config_flags%mars == 11) .OR. (config_flags%mars == 12) ) then
    13231333       print *, '**** INTERPOLATE DUSTQ **** RANK 4 in SCALAR'             
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F

    r678 r802  
    27112711 .OR. (  config_flags%mars .eq.  2                    ) &    !! 2   pass the dust tracer
    27122712 .OR. (  config_flags%mars .eq.  3                    ) &    !! 3   pass dust Q and N
     2713 .OR. (  config_flags%mars .eq.  10                   ) &    !! 10  pass the co2 tracer
    27132714! .OR. ( (config_flags%mars .eq. 11) .AND. (is .EQ. 2) ) &    !! 11  pass only water vapor
    27142715 .OR. (  config_flags%mars .eq. 11                    ) &    !! 11  pass EVERYTHING
Note: See TracChangeset for help on using the changeset viewer.