Changeset 324


Ignore:
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
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/README

    r321 r324  
    10931093  Also adapted 'inichim_readcallphys': removed some obsolescent tests on
    10941094  number of tracers for given combinations of options.
     1095
     1096== 21/10/2011 == AS
     1097- Added possibility for CH4 tracer in tracer.h and initracer.F

  • trunk/LMDZ.MARS/libf/phymars/initracer.F

    r283 r324  
    342342          count=count+1
    343343        endif
     344        if (noms(iq).eq."ch4") then
     345          igcm_ch4=iq
     346          mmol(igcm_ch4)=16.
     347          count=count+1
     348        endif
    344349        if (noms(iq).eq."ar") then
    345350          igcm_ar=iq
    … …  
    347352          count=count+1
    348353        endif
    349         if (noms(iq).eq."n") then
     354        if (noms(iq).eq."n") then
    350355          igcm_n=iq
    351356          mmol(igcm_n)=14.
    … …  
    402407          count=count+1
    403408        endif
    404         if (noms(iq).eq."n2plus") then
     409        if (noms(iq).eq."n2plus") then
    405410          igcm_n2plus=iq
    406411          mmol(igcm_n2plus)=28.
    … …  
    433438          count=count+1
    434439        endif
    435        
     440
    436441
    437442      enddo ! of do iq=1,nqmx

  • trunk/LMDZ.MARS/libf/phymars/tracer.h

    r171 r324  
    5050      integer :: igcm_no2
    5151      integer :: igcm_n2d
     52      integer :: igcm_ch4
    5253      ! Ions
    5354      integer :: igcm_co2plus
    … …  
    7475     & igcm_o2,igcm_o3,igcm_h,igcm_h2,igcm_oh,igcm_ho2,igcm_h2o2,       &
    7576     & igcm_n2,igcm_ar,igcm_n,igcm_no,igcm_no2,igcm_n2d,                &
     77     & igcm_ch4,                                                        &
    7678     & igcm_co2plus,igcm_oplus,igcm_o2plus,igcm_coplus,igcm_cplus,      &
    7779     & igcm_nplus,igcm_noplus,igcm_n2plus,igcm_hplus,igcm_elec,         &

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

    r279 r324  
    178178      3) tracers='2' ;;
    179179     11) tracers='4' ;;
     180     42) tracers='18' ;;
    180181      *) tracers='1' ;;
    181182  esac

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/compile_pgf

    r235 r324  
    88-L$NETCDF/lib -lnetcdf \
    99-I$NETCDF/include \
    10 -o readmeteo.exe
     10-o readmeteo.exe #-DPHOTOCHEM
    1111
    1212pgf90 create_readmeteo.F90 \
    1313-L$NETCDF/lib -lnetcdf \
    1414-I$NETCDF/include \
    15 -o create_readmeteo.exe
     15-o create_readmeteo.exe #-DPHOTOCHEM \
    1616
    1717\rm fix_no_info.inc 2> /dev/null

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/readmeteo.F90

    r73 r324  
    9898character*13, dimension(:), allocatable :: date_out
    9999character*19, dimension(:), allocatable :: date_out2
     100
     101#ifdef PHOTOCHEM
     102real, dimension(:,:,:,:,:), allocatable :: chemtrac
     103integer :: nchemtrac,i
     104CHARACTER*20,DIMENSION(:),ALLOCATABLE :: wtnom
     105#endif
     106
     107
    100108!***************************************************************************
    101109!***************************************************************************
    … …  
    230238allocate(interm(lonlen,latlen))
    231239allocate(gwparam(lonlen,latlen,5))
    232 allocate(ghtsfile(lonlen,latlen))       !! no scan axis
     240allocate(ghtsfile(lonlen,latlen))    !! no scan axis
    233241allocate(vide(lonlen,latlen))
    234242allocate(ones(lonlen,latlen))
    235243allocate(lat(latlen), lon(lonlen), alt(altlen), time(timelen))
    236244allocate(aps(altlen),bps(altlen),levels(altlen))
     245#ifdef PHOTOCHEM
     246nchemtrac = 14
     247allocate(wtnom(nchemtrac))
     248wtnom(1)  = "c_co2"
     249wtnom(2)  = "c_co"
     250wtnom(3)  = "c_o"
     251wtnom(4)  = "c_o1d"
     252wtnom(5)  = "c_o2"
     253wtnom(6)  = "c_o3"
     254wtnom(7)  = "c_h"
     255wtnom(8)  = "c_h2"
     256wtnom(9)  = "c_oh"
     257wtnom(10) = "c_ho2"
     258wtnom(11) = "c_h2o2"
     259wtnom(12) = "c_ch4"
     260wtnom(13) = "c_n2"
     261wtnom(14) = "c_ar"
     262allocate(chemtrac(lonlen,latlen,altlen,timelen,nchemtrac))
     263chemtrac(:,:,:,:,:)=0
     264#endif
    237265
    238266tfile(:,:,:,:)=0
    … …  
    523551   IF (ierr .NE. NF_NOERR) THEN
    524552      ierr = NF_INQ_VARID (nid,"t",nvarid)
    525         IF (ierr .NE. NF_NOERR) THEN
     553        IF (ierr .NE. NF_NOERR) THEN
    526554          PRINT *, "Error: Readmeteo <t> not found"
    527555          stop
    528         ENDIF
     556        ENDIF
    529557   ENDIF
    530558#ifdef NC_DOUBLE
    … …  
    627655!! special dust stuff
    628656!!------------------------
    629 
    630657
    631658!SELECT CASE(ident)
    … …  
    677704!!!!!!!!
    678705!!!!!!!! new physics
     706
     707
     708!!!!!!!!!!!!!!!!!!!!!!!!NEW PHYSICS + PHOTOCHEM
     709!!!!!!!!!!!!!!!!!!!!!!!!NEW PHYSICS + PHOTOCHEM
     710#ifdef PHOTOCHEM
     711    print *,'photochem'
     712    DO i=1,nchemtrac
     713     print *,wtnom(i)
     714     ierr=NF_INQ_VARID(nid,wtnom(i),nvarid)
     715     if (ierr.ne.NF_NOERR) then
     716       write(*,*) "...No ",wtnom(i), " - set to 0"
     717       chemtrac(:,:,:,:,i)=0.
     718     else
     719       ierr=NF_GET_VAR_REAL(nid,nvarid,chemtrac(:,:,:,:,i))
     720     endif
     721    ENDDO
     722#endif
     723!!!!!!!!!!!!!!!!!!!!!!!!NEW PHYSICS + PHOTOCHEM
     724!!!!!!!!!!!!!!!!!!!!!!!!NEW PHYSICS + PHOTOCHEM
     725
    679726
    680727
    … …  
    11921239!print *,'The field '//DESC//' was written to '//output
    11931240
     1241!------------------------!
     1242! >>> Write a variable   !
     1243!     PHOTOCHEMISTRY     !
     1244!------------------------!
     1245#ifdef PHOTOCHEM
     1246    DO i=1,nchemtrac
     1247       FIELD=wtnom(i)
     1248       UNITS='units'
     1249       DESC='desc'
     1250       XLVL=200100.
     1251       SLAB=chemtrac(:,:,1,time_out(l),i)
     1252       ! And now put everything in the destination file
     1253       ! ... Header
     1254       write(1) IFV
     1255       write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1256       write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1257       ! ... Data
     1258       write(1) SLAB
     1259    ENDDO
     1260#endif
    11941261
    11951262!!----------------------------------------------------
    … …  
    14551522END DO
    14561523!print *,'The field '//DESC//' was written to '//output
     1524
     1525
     1526!------------------------!
     1527! >>> Write a variable   !
     1528!     PHOTOCHEMISTRY     !
     1529!------------------------!
     1530#ifdef PHOTOCHEM
     1531    DO i=1,nchemtrac
     1532       FIELD=wtnom(i)
     1533       UNITS='units'
     1534       DESC='desc'
     1535       DO k = 1,altlen
     1536         XLVL=levels(k)
     1537         SLAB=chemtrac(:,:,k,time_out(l),i)
     1538         ! And now put everything in the destination file
     1539         ! ... Header
     1540         write(1) IFV
     1541         write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1542         write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1543         ! ... Data
     1544         write(1) SLAB
     1545       END DO
     1546    ENDDO
     1547#endif
    14571548
    14581549print *,'****done file '//output, int(100.*float(l)/float(FILES)), ' % '
    … …  
    14981589deallocate(aps,bps,levels)
    14991590
     1591#ifdef PHOTOCHEM
     1592deallocate(chemtrac)
     1593deallocate(wtnom)
     1594#endif
    15001595
    15011596print *, '------------------------'

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

    r316 r324  
    130130state  real  qdustn    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QDUSTN"        "Dust_number mixing ratio"   "kg kg-1"
    131131state  real  qco2      ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "QCO2"          "CO2 mixing ratio"           "kg kg-1"
     132state  real  chem_co   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_co"     ""   ""
     133state  real  chem_o    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_o"      ""   ""
     134state  real  chem_o1d  ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_o1d"    ""   ""
     135state  real  chem_o2   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_o2"     ""   ""
     136state  real  chem_o3   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_o3"     ""   ""
     137state  real  chem_h    ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_h"      ""   ""
     138state  real  chem_h2   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_h2"     ""   ""
     139state  real  chem_oh   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_oh"     ""   ""
     140state  real  chem_ho2  ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_ho2"    ""   ""
     141state  real  chem_h2o2 ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_h2o2"   ""   ""
     142state  real  chem_ch4  ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_ch4"    ""   ""
     143state  real  chem_n2   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_n2"     ""   ""
     144state  real  chem_ar   ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt)  "chem_ar"     ""   ""
    132145####
    133146####
    … …  
    13531366package   dust2eq      mars==3                      -              scalar:qdust,qdustn
    13541367package   newwater     mars==11                     -              scalar:qh2o,qh2o_ice,qdust,qdustn
     1368package   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
    13551369##### MARS OPTIONS
    13561370##### MARS OPTIONS

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

    r315 r324  
    650650
    651651SELECT CASE (MARS_MODE) !! ONLY ALLOW FOR MODES DEFINED IN Registry.EM
    652    CASE(4-10,12-19,22:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
     652   CASE(4-10,12-19,22-41,43:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
    653653   PRINT *, 'NOT SUPPORTED, to be done'
    654654   STOP
    … …  
    662662!package   radioac      mars==20                     -              scalar:qtrac1
    663663!package   radioac2     mars==21                     -              scalar:upward,downward
     664!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
    664665!!!!!!!!!!!!!!!!!!! FOR REFERENCE
    665666
    … …  
    672673      wtnom(nq) = 'co2'
    673674    CASE(1)
    674       wtnom(1)  = 'h2o_vap'     
    675       wtnom(2)  = 'h2o_ice'    
     675      wtnom(1)  = 'h2o_vap'
     676      wtnom(2)  = 'h2o_ice'
    676677    CASE(2)
    677678      wtnom(1)  = 'dust01'     
    678679    CASE(3)
    679680      wtnom(1)  = 'dust_mass'
    680       wtnom(2)  = 'dust_number'
     681      wtnom(2)  = 'dust_number' 
    681682    CASE(11)
    682683      wtnom(1)  = 'h2o_vap'
    683       wtnom(2)  = 'h2o_ice'
     684      wtnom(2)  = 'h2o_ice' 
    684685      wtnom(3)  = 'dust_mass'
    685686      wtnom(4)  = 'dust_number'
    … …  
    688689    CASE(21)
    689690      wtnom(1) = 'upward'
    690       wtnom(2) = 'downward'
     691      wtnom(2) = 'downward'
     692    CASE(42)
     693      wtnom(1)  = 'co2'
     694      wtnom(2)  = 'co'
     695      wtnom(3)  = 'o'
     696      wtnom(4)  = 'o1d'
     697      wtnom(5)  = 'o2'
     698      wtnom(6)  = 'o3'
     699      wtnom(7)  = 'h'
     700      wtnom(8)  = 'h2'
     701      wtnom(9)  = 'oh'
     702      wtnom(10)  = 'ho2'
     703      wtnom(11)  = 'h2o2'
     704      wtnom(12)  = 'ch4'
     705      wtnom(13)  = 'n2'
     706      wtnom(14)  = 'ar'
     707      wtnom(15)  = 'h2o_ice'
     708      wtnom(16)  = 'h2o_vap'
     709      wtnom(17)  = 'dust_mass'
     710      wtnom(18)  = 'dust_number'
    691711END SELECT
    692712#endif
    … …  
    729749q_prof(:,1:nq) = SCALAR(i,kps:kpe,j,2:nq+1)  !! the names were set above !! one dummy tracer in WRF
    730750  !!! CAS DU CO2
    731 DO iii=1,nq
    732  IF ( wtnom(iii) .eq. 'co2' ) q_prof(:,iii) = 0.95
    733 ENDDO
    734 
    735 IF ((MARS_MODE .EQ. 20) .OR. (MARS_MODE .EQ. 21)) THEN
     751  DO iii=1,nq
     752   IF ( wtnom(iii) .eq. 'co2' ) q_prof(:,iii) = 0.95
     753  ENDDO
     754  IF ((MARS_MODE .EQ. 20) .OR. (MARS_MODE .EQ. 21)) THEN
    736755   IF (firstcall .EQV. .true.) THEN
    737756      q_prof(:,:) = 0.95
    738757   ENDIF
    739 ENDIF
    740 
     758  ENDIF
    741759#else
    742760SELECT CASE (MARS_MODE)
    … …  
    841859!!!! ADDITIONAL SECURITY. THIS MIGHT HAPPEN WITH OLD INIT FILES.
    842860IF (z0_val == 0.) THEN
    843    PRINT *, 'WELL, z0 is 0, this is no good. Setting to old defaults value 0.01 m'
     861   IF ( (i == ips) .AND. (j == jps) ) PRINT *, 'WELL, z0 is 0, this is no good. Setting to old defaults value 0.01 m'
    844862   z0_val = 0.01
    845863ENDIF

  • trunk/MESOSCALE_DEV/PLOT/SPEC/GW/gravitwaveprof.pro

    r131 r324  
    5858nts=1  & nte=100  & nxs=50   & nxe=50
    5959nts=1  & nte=100  & nxs=60   & nxe=60
    60 nts=1  & nte=100  & nxs=20   & nxe=20
     60;nts=1  & nte=100  & nxs=20   & nxe=20
     61nts=29 & nte=29   & nxs=60   & nxe=60
    6162;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    6263;
    … …  
    150151     alt = [40.,80.]
    151152     alt = [0.,120.]
     153     alt = [10.,110.]
    152154     ;;;;;;;;;;;
    153155     ;;;;;;;;;;;
    … …  
    156158     minfield_init = 85.  ;100.
    157159     maxfield_init = 185. ;250.
     160     minfield_init = 90.
     161     maxfield_init = 230.
    158162     title_axis = ['Temperature (K)','Altitude (km)']
    159163     overplot = mcd_tpot
    … …  
    161165     SPAWN, "echo 'intervaly=10.' >> param_plot.idl"
    162166     SPAWN, "echo 'intervalx=5.' >> param_plot.idl"
     167     SPAWN, "echo 'intervalx=10.' >> param_plot.idl"
    163168     ;;;;;;;;;;;
    164169
    … …  
    216221        ;;oplot, yeye, vert/1000., linestyle=1
    217222
     223
     224loadct, 4
     225for i=16,37 do begin  ;10-37
     226for j=0,40 do begin
     227dec = j - 20 
     228overoverplot = reform(invar(nx+dec,middle+dec,*,i))
     229oplot, overoverplot, column, color=220  ;psym=3
     230endfor
     231endfor
     232
     233!p.thick = 4.0
     234oplot, overplot, overplot_column, linestyle=2
     235!p.thick = 2.0
     236
     237!p.thick = 4.0
     238oplot, reform(invar(nx,middle,*,28)), column
     239!p.thick = 2.0
     240
    218241@tempcond.inc
    219242overplot = reform(overplot(nx,middle,*,nt))
     243w = where(overplot le 0.) & overplot[w] = !VALUES.F_NAN
    220244oplot, overplot, overplot_column, linestyle=1
    221245
Note: See TracChangeset for help on using the changeset viewer.