Ignore:
Timestamp:
Dec 19, 2011, 6:54:32 PM (13 years ago)
Author:
aslmd
Message:

MESOSCALE: corrected a bug with dust outputs through MESOINI in new physics. added mars=12 for scavenging runs (see r76).

Location:
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM

    r341 r481  
    129129state  real  qdust     ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QDUST"         "Dust mixing ratio"          "kg kg-1"
    130130state  real  qdustn    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QDUSTN"        "Dust_number mixing ratio"   "kg kg-1"
     131state  real  qccn     ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QCCN"         "CCN mixing ratio"          "kg kg-1"
     132state  real  qccnn    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QCCNN"        "CCN number mixing ratio"   "kg kg-1"
    131133state  real  qco2      ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QCO2"          "CO2 mixing ratio"           "kg kg-1"
    132134state  real  chem_co   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_co"     ""   ""
     
    13661368package   dust2eq      mars==3                      -              scalar:qdust,qdustn
    13671369package   newwater     mars==11                     -              scalar:qh2o,qh2o_ice,qdust,qdustn
     1370package   newwatersca  mars==12                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn
    13681371package   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
    13691372##### MARS OPTIONS
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F

    r315 r481  
    12931293                                    its , ite , jts , jte , kts , kte )
    12941294endif       
    1295 if ( ( config_flags%mars == 1 ) .OR. ( config_flags%mars == 11 ) ) then
     1295if ( ( config_flags%mars == 1 ) .OR. ( config_flags%mars == 11 ) .OR. ( config_flags%mars == 12 ) ) then
    12961296       print *, '**** INTERPOLATE HV **** RANK 2 in SCALAR'
    12971297       !print *, size(scalar(0,0,0,:)), P_QH2O, P_QH2O_ICE
     
    13131313endif
    13141314#ifdef NEWPHYS
    1315 if ( config_flags%mars == 11 ) then
     1315if ( config_flags%mars == 11 .OR. config_flags%mars == 12 ) then
    13161316       print *, '**** INTERPOLATE DUSTQ **** RANK 4 in SCALAR'             
    13171317       CALL vert_interp_old ( grid%em_dustq_gc , grid%em_pd_gc , scalar(:,:,:,4) , grid%em_pb , &
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F

    r315 r481  
    27092709! .OR. ( (config_flags%mars .eq. 11) .AND. (is .EQ. 2) ) &    !! 11  pass only water vapor
    27102710 .OR. (  config_flags%mars .eq. 11                    ) &    !! 11  pass EVERYTHING
     2711 .OR. (  config_flags%mars .eq. 12                    ) &    !! 12  pass EVERYTHING
    27112712 .OR. config_flags%nested ) THEN                             !! *   pass all tracers if nested
    27122713         CALL relax_bdy_scalar ( scalar_tend(ims,kms,jms,is),            &
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/main/real_em.F

    r76 r481  
    712712IF ( (config_flags%mars .eq. 1) .OR. &
    713713     (config_flags%mars .eq. 3) .OR. &
    714      (config_flags%mars .eq. 11) ) THEN
     714     (config_flags%mars .eq. 11) .OR. &
     715     (config_flags%mars .eq. 12) ) THEN
    715716      CALL couple ( grid%em_mu_2 , grid%em_mub , q2bdy3dtemp1 , grid%scalar(:,:,:,3) , 't' , grid%msft , &
    716717                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, ips, ipe, jps, jpe, kps, kpe )
    717718ENDIF
    718 IF (config_flags%mars .eq. 11) THEN
     719IF ((config_flags%mars .eq. 11) .OR. (config_flags%mars .eq. 12)) THEN
    719720      CALL couple ( grid%em_mu_2 , grid%em_mub , q3bdy3dtemp1 , grid%scalar(:,:,:,4) , 't' , grid%msft , &
    720721                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, ips, ipe, jps, jpe, kps, kpe )
     
    786787IF ( (config_flags%mars .eq. 1) .OR. &
    787788     (config_flags%mars .eq. 3) .OR. &
    788      (config_flags%mars .eq. 11) ) THEN
     789     (config_flags%mars .eq. 11) .OR. &
     790     (config_flags%mars .eq. 12) ) THEN
    789791      CALL stuff_bdy     ( q2bdy3dtemp1 , grid%scalar_b(:,:,:,:,3)   , 'T' , ijds , ijde , spec_bdy_width      , &
    790792                                                                 ids , ide , jds , jde , kds , kde , &
     
    792794                                                                 ips , ipe , jps , jpe , kps , kpe )
    793795ENDIF
    794 IF (config_flags%mars .eq. 11) THEN
     796IF ((config_flags%mars .eq. 11) .OR. (config_flags%mars .eq. 12)) THEN
    795797      CALL stuff_bdy     ( q3bdy3dtemp1 , grid%scalar_b(:,:,:,:,4)   , 'T' , ijds , ijde , spec_bdy_width      , &
    796798                                                                 ids , ide , jds , jde , kds , kde , &
     
    863865IF ( (config_flags%mars .eq. 1) .OR. &
    864866     (config_flags%mars .eq. 3) .OR. &
    865      (config_flags%mars .eq. 11) ) THEN
     867     (config_flags%mars .eq. 11) .OR. &
     868     (config_flags%mars .eq. 12) ) THEN
    866869      CALL couple ( grid%em_mu_2 , grid%em_mub , q2bdy3dtemp2 , grid%scalar(:,:,:,3) , 't' , grid%msft , &
    867870                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, ips, ipe, jps, jpe, kps, kpe )
    868871ENDIF
    869 IF (config_flags%mars .eq. 11) THEN
     872IF ((config_flags%mars .eq. 11) .OR. (config_flags%mars .eq. 12)) THEN
    870873      CALL couple ( grid%em_mu_2 , grid%em_mub , q3bdy3dtemp2 , grid%scalar(:,:,:,4) , 't' , grid%msft , &
    871874                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, ips, ipe, jps, jpe, kps, kpe )
     
    941944IF ( (config_flags%mars .eq. 1) .OR. &
    942945     (config_flags%mars .eq. 3) .OR. &
    943      (config_flags%mars .eq. 11) ) THEN
     946     (config_flags%mars .eq. 11) .OR. &
     947     (config_flags%mars .eq. 12) ) THEN
    944948      CALL stuff_bdytend ( q2bdy3dtemp2 , q2bdy3dtemp1 , REAL(interval_seconds) ,grid%scalar_bt(:,:,:,:,3) , 'T' , &
    945949                                                            ijds , ijde , spec_bdy_width      , &
     
    948952                                                            ips , ipe , jps , jpe , kps , kpe )
    949953ENDIF
    950 IF (config_flags%mars .eq. 11) THEN
     954IF ((config_flags%mars .eq. 11) .OR. (config_flags%mars .eq. 12)) THEN
    951955      CALL stuff_bdytend ( q3bdy3dtemp2 , q3bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,4) , 'T' , &
    952956                                                            ijds , ijde , spec_bdy_width      , &
     
    11121116IF ( (config_flags%mars .eq. 1) .OR. &
    11131117     (config_flags%mars .eq. 3) .OR. &
    1114      (config_flags%mars .eq. 11) ) THEN
     1118     (config_flags%mars .eq. 11) .OR. &
     1119     (config_flags%mars .eq. 12) ) THEN
    11151120        CALL stuff_bdy     ( q2bdy3dtemp1 , grid%scalar_b(:,:,:,:,3)   , 'T', ijds , ijde , spec_bdy_width      , &
    11161121                                                                    ids , ide , jds , jde , kds , kde , &
     
    11181123                                                                    ips , ipe , jps , jpe , kps , kpe )
    11191124ENDIF
    1120 IF (config_flags%mars .eq. 11) THEN
     1125IF ((config_flags%mars .eq. 11) .OR. (config_flags%mars .eq. 12)) THEN
    11211126        CALL stuff_bdy     ( q3bdy3dtemp1 , grid%scalar_b(:,:,:,:,4)   , 'T', ijds , ijde , spec_bdy_width      , &
    11221127                                                                    ids , ide , jds , jde , kds , kde , &
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

    r341 r481  
    649649!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    650650
    651 SELECT CASE (MARS_MODE) !! ONLY ALLOW FOR MODES DEFINED IN Registry.EM
    652    CASE(4-10,12-19,22-41,43:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
    653    PRINT *, 'NOT SUPPORTED, to be done'
    654    STOP
    655 END SELECT
     651!SELECT CASE (MARS_MODE) !! ONLY ALLOW FOR MODES DEFINED IN Registry.EM
     652!   CASE(4-10,13-19,22-41,43:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
     653!   PRINT *, 'NOT SUPPORTED, to be done'
     654!   STOP
     655!END SELECT
    656656!!!!!!!!!!!!!!!!!!! FOR REFERENCE ; FROM REGISTRY.EM
    657657!package   default      mars==0                      -              -
     
    685685      wtnom(3)  = 'dust_mass'
    686686      wtnom(4)  = 'dust_number'
     687    CASE(12)
     688      wtnom(1)  = 'h2o_vap'
     689      wtnom(2)  = 'h2o_ice'
     690      wtnom(3)  = 'dust_mass'
     691      wtnom(4)  = 'dust_number'
     692      wtnom(5)  = 'ccn_mass'
     693      wtnom(6)  = 'ccn_number'
    687694    CASE(20)
    688695      wtnom(1) = 'qtrac1'
     
    920927    CASE(3)
    921928    qsurf_val(:)=0.                !!! temporaire, a definir           
    922     CASE(11)
     929    CASE(11-12)
    923930    qsurf_val(1)=0.
    924931    qsurf_val(2)=MARS_WICE(i,j)    !! logique avec wtnom(2) = 'h2o_ice' defini ci-dessus
Note: See TracChangeset for help on using the changeset viewer.