Changeset 2753a5c in flexpart.git


Ignore:
Timestamp:
Dec 20, 2018, 2:55:22 PM (6 months ago)
Author:
Ignacio Pisso <ip@…>
Branches:
release-10
Children:
7123c70
Parents:
df96ea65 (diff), 4ad96c5 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.
Message:

resolved merging conflicts with GFS branch

Files:
35 added
179 deleted
28 edited
16 moved

Legend:

Unmodified
Added
Removed
  • options/COMMAND

    r6d73c4b r2753a5c  
    1616 ITSPLIT=        99999999, ! Interval of particle splitting (s)
    1717 LSYNCTIME=           900, ! All processes are synchronized to this time interval (s)
    18  CTL=          -5.0000000, ! Times step smaller than Lagr. time scale, LSYNCTIME if CTL<0, but can be shorter 
     18 CTL=          -5.0000000, ! CTL>1, ABL time step = (Lagrangian timescale (TL))/CTL, uses LSYNCTIME if CTL<0
    1919 IFINE=                 4, ! Reduction for time step in vertical transport, used only if CTL>1
    2020 IOUT=                  1, ! Output type: [1]mass 2]pptv 3]1&2 4]plume 5]1&4, +8 for NetCDF output     
    2121 IPOUT=                 0, ! Particle position output: 0]no 1]every output 2]only at end   
    2222 LSUBGRID=              0, ! Increase of ABL heights due to sub-grid scale orographic variations;[0]off 1]on
    23  LCONVECTION=           0, ! Switch for convection parameterization;[0]off 1]on   
     23 LCONVECTION=           1, ! Switch for convection parameterization;0]off [1]on   
    2424 LAGESPECTRA=           0, ! Switch for calculation of age spectra (needs AGECLASSES);[0]off 1]on 
    25  IPIN=                  0, ! Warm start from particle dump (needs previous partposit_end file); 0]no [1]yes 
     25 IPIN=                  0, ! Warm start from particle dump (needs previous partposit_end file); [0]no 1]yes 
    2626 IOUTPUTFOREACHRELEASE= 1, ! Separate output fields for each location in the RELEASE file; [0]no 1]yes
    2727 IFLUX=                 0, ! Output of mass fluxes through output grid box boundaries
    2828 MDOMAINFILL=           0, ! Switch for domain-filling, if limited-area particles generated at boundary
    29  IND_SOURCE=            1, ! Unit to be used at the ource   ;  [1]mass 2]mass mixing ratio
    30  IND_RECEPTOR=          1, ! Unit to be used at the receptor;  [1]mass 2]mass mixing ratio
     29 IND_SOURCE=            1, ! Unit to be used at the source   ;  [1]mass 2]mass mixing ratio
     30 IND_RECEPTOR=          1, ! Unit to be used at the receptor; [1]mass 2]mass mixing ratio 3]wet depo. 4]dry depo.
    3131 MQUASILAG=             0, ! Quasi-Lagrangian mode to track individual numbered particles
    3232 NESTED_OUTPUT=         0, ! Output also for a nested domain
    33  LINIT_COND=            0, ! Output sensitivity to initial conditions (backwards mode only)
     33 LINIT_COND=            0, ! Output sensitivity to initial conditions (bkw mode only) [0]off 1]conc 2]mmr
    3434 SURF_ONLY=             0, ! Output only for the lowest model layer, used w/ LINIT_COND=1 or 2
    3535 CBLFLAG=               0, ! Skewed, not Gaussian turbulence in the convective ABL, need large CTL and IFINE
  • options/RELEASES

    r6d73c4b r2753a5c  
    1 *************************************************************************
    2 *                                                                       *
    3 *                                                                       *
    4 *                                                                       *
    5 *   Input file for the Lagrangian particle dispersion model FLEXPART    *
    6 *                        Please select your options                     *
    7 *                                                                       *
    8 *                                                                       *
    9 *                                                                       *
    10 *************************************************************************
     1***************************************************************************************************************
     2*                                                                                                             *
     3*                                                                                                             *
     4*                                                                                                             *
     5*   Input file for the Lagrangian particle dispersion model FLEXPART                                          *
     6*                        Please select your options                                                           *
     7*                                                                                                             *
     8*                                                                                                             *
     9*                                                                                                             *
     10***************************************************************************************************************
    1111&RELEASES_CTRL
    12  NSPEC      =           1,
    13  SPECNUM_REL=          24,
     12 NSPEC      =           1, ! Total number of species
     13 SPECNUM_REL=          24, ! Species numbers in directory SPECIES
    1414 /
    15 &RELEASE
    16  IDATE1  =       20120101,
    17  ITIME1  =         090000,
    18  IDATE2  =       20120101,
    19  ITIME2  =         090000,
    20  LON1    =          0.000,
    21  LON2    =          0.000,
    22  LAT1    =         20.000,
    23  LAT2    =         20.000,
    24  Z1      =         50.000,
    25  Z2      =         50.000,
    26  ZKIND   =              1,
    27  MASS    =       1.0000E8,
    28  PARTS   =          10000,
    29  COMMENT =    "RELEASE 1",
     15&RELEASE                   ! For each release
     16 IDATE1  =       20120101, ! Release start date, YYYYMMDD: YYYY=year, MM=month, DD=day
     17 ITIME1  =         090000, ! Release start time in UTC HHMISS: HH hours, MI=minutes, SS=seconds
     18 IDATE2  =       20120101, ! Release end date, same as IDATE1
     19 ITIME2  =         090000, ! Release end time, same as ITIME1
     20 LON1    =          0.000, ! Left longitude of release box -180 < LON1 <180
     21 LON2    =          0.000, ! Right longitude of release box, same as LON1
     22 LAT1    =         20.000, ! Lower latitude of release box, -90 < LAT1 < 90
     23 LAT2    =         20.000, ! Upper latitude of release box same format as LAT1
     24 Z1      =         50.000, ! Lower height of release box meters/hPa above reference level
     25 Z2      =         50.000, ! Upper height of release box meters/hPa above reference level
     26 ZKIND   =              1, ! Reference level 1=above ground, 2=above sea level, 3 for pressure in hPa
     27 MASS    =       1.0000E8, ! Total mass emitted, only relevant for fwd simulations
     28 PARTS   =          10000, ! Total number of particles to be released
     29 COMMENT =    "RELEASE 1", ! Comment, written in the outputfile
    3030 /
  • options/SPECIES/SPECIES_008

    r8594fb8 r79996be9  
    1919 POHDCONST=-9.9,              ! OH Reaction rate - D [K]
    2020 POHNCONST=2.0,               !  OH Reaction rate -  N (no unit)
    21  PAREA_HOUR=0.578, 0.491, 0.428, 0.329, 0.384, 0.485, 0.763, 1.103, 1.084, 1.047, 1.096, 1.196, 1.298, 1.357, 1.447, 1.565, 1.636, 1.662, 1.401, 1.168, 1.031, 0.926, 0.816, 0.709
    22  PAREA_DOW=1.060, 1.060, 1.060, 1.060, 1.060, 0.900, 0.000
     21 PAREA_HOUR=0.578, 0.491, 0.428, 0.329, 0.384, 0.485, 0.763, 1.103, 1.084, 1.047, 1.096, 1.196, 1.298, 1.357, 1.447, 1.565, 1.636, 1.662, 1.401, 1.168, 1.031, 0.926, 0.816, 0.709,
     22 PPOINT_HOUR=0.845, 0.806, 0.786, 0.779, 0.793, 0.832, 0.895, 0.977, 1.031, 1.071, 1.105, 1.118, 1.131, 1.136, 1.143, 1.141, 1.133, 1.118, 1.097, 1.091, 1.079, 1.036, 0.966, 0.892,
     23 PAREA_DOW=1.060, 1.060, 1.060, 1.060, 1.060, 0.900, 0.000,
     24 PPOINT_DOW=1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000,
    2325 /
  • src/FLEXPART.f90

    r50958b8 r2753a5c  
    455455  call timemanager(metdata_format)
    456456
     457  if (verbosity.gt.0) then
    457458! NIK 16.02.2005
    458   do i=1,nspec
    459     write(*,*) '**********************************************'
    460     write(*,*) 'Scavenging statistics for species ', species(i), ':'
    461     write(*,*) 'Total number of occurences of below-cloud scavenging', &
    462          & tot_blc_count(i)
    463     write(*,*) 'Total number of occurences of in-cloud    scavenging', &
    464          & tot_inc_count(i)
    465     write(*,*) '**********************************************'
    466   end do
     459    do i=1,nspec
     460      if (tot_inc_count(i).gt.0) then
     461         write(*,*) '**********************************************'
     462         write(*,*) 'Scavenging statistics for species ', species(i), ':'
     463         write(*,*) 'Total number of occurences of below-cloud scavenging', &
     464           & tot_blc_count(i)
     465         write(*,*) 'Total number of occurences of in-cloud    scavenging', &
     466           & tot_inc_count(i)
     467         write(*,*) '**********************************************'
     468      endif
     469    end do
     470    write (*,*) 'timemanager> call wetdepo'
     471  endif
    467472 
    468473  write(*,*) 'CONGRATULATIONS: YOU HAVE SUCCESSFULLY COMPLETED A FLE&
  • src/com_mod.f90

    rc2bd55e re9e0f06  
    174174  real :: ri(5,numclass),rac(5,numclass),rcl(maxspec,5,numclass)
    175175  real :: rgs(maxspec,5,numclass),rlu(maxspec,5,numclass)
    176   real :: rm(maxspec),dryvel(maxspec),kao(maxspec)
     176  real :: rm(maxspec),dryvel(maxspec)
    177177  real :: ohcconst(maxspec),ohdconst(maxspec),ohnconst(maxspec)
    178   ! se  it is possible to associate a species with a second one to make transfer from gas to aerosol
    179   integer :: spec_ass(maxspec)
    180178
    181179  real :: area_hour(maxspec,24),point_hour(maxspec,24)
  • src/get_wetscav.f90

    r79e0349 re9e0f06  
    151151      if (height(il).gt.ztra1(jpart)) then
    152152        hz=il-1
    153 !        goto 26
    154153        exit
    155154      endif
    156155    end do
    157 !26  continue
    158 
    159156
    160157    if (ngrid.eq.0) then
     
    203200
    204201
    205 !ZHG oct 2014 : Calculated for 1) both 2) lsp 3) convp
     202!ZHG oct 2014 : Calculated for 1) both 2) lsp 3) convp - 2 and 3 not used removed by SE
    206203! Tentatively differentiate the grfraction for lsp and convp for treating differently the two forms
    207204! for now they are treated the same
    208205    grfraction(1)=max(0.05,cc*(lsp*lfr(i)+convp*cfr(j))/(lsp+convp))
    209     grfraction(2)=max(0.05,cc*(lfr(i)))
    210     grfraction(3)=max(0.05,cc*(cfr(j)))
    211 
    212206
    213207! 2) Computation of precipitation rate in sub-grid cell
    214208!******************************************************
    215209    prec(1)=(lsp+convp)/grfraction(1)
    216     prec(2)=(lsp)/grfraction(2)
    217     prec(3)=(convp)/grfraction(3)
    218 
    219210
    220211! 3) Computation of scavenging coefficients for all species
    221212!    Computation of wet deposition
    222213!**********************************************************
    223 
    224214
    225215      if (ngrid.gt.0) then
     
    237227!******************************************************************
    238228        if ((dquer(ks).le.0.).and.(weta_gas(ks).gt.0..or.wetb_gas(ks).gt.0.)) then
    239           !        if (weta(ks).gt.0. .or. wetb(ks).gt.0.) then
    240229          blc_count(ks)=blc_count(ks)+1
    241230          wetscav=weta_gas(ks)*prec(1)**wetb_gas(ks)
     
    271260          endif
    272261         
    273 !             write(*,*) 'bl-cloud, act_temp=',act_temp, ',prec=',prec(1),',wetscav=', wetscav, ', jpart=',jpart
    274 
    275262        endif ! gas or particle
    276263!      endif ! positive below-cloud scavenging parameters given in Species file
     
    285272        if ((ccn_aero(ks).gt.0. .or. in_aero(ks).gt.0.).or.(henry(ks).gt.0.and.dquer(ks).le.0)) then
    286273          inc_count(ks)=inc_count(ks)+1
    287 !          write(*,*) 'Incloud: ',inc_count
    288274! if negative coefficients (turned off) set to zero for use in equation
    289275          if (ccn_aero(ks).lt.0.) ccn_aero(ks)=0.
     
    300286! sec test
    301287!           cl=1E6*1E-7*prec(1)**0.3 !Sec GFS new
    302             cl=1E6*2E-7*prec(1)**0.36 !Sec ECMWF new
     288            cl=1E6*2E-7*prec(1)**0.36 !Sec ECMWF new, is also suitable for GFS
    303289!           cl=2E-7*prec(1)**0.36 !Andreas
    304290!           cl=1.6E-6*prec(1)**0.36 !Henrik
     
    322308          if (dquer(ks).gt.0.) then
    323309            S_i= frac_act/cl
    324 !           write(*,*) 'Si: ',S_i
    325 
    326310! GAS
    327311!****
    328312          else
    329 
    330313            cle=(1-cl)/(henry(ks)*(r_air/3500.)*act_temp)+cl
    331 !REPLACE to switch old/ new scheme
    332           ! S_i=frac_act/cle
    333314            S_i=1/cle
    334315          endif ! gas or particle
    335316
    336317! scavenging coefficient based on Hertel et al 1995 - using the S_i for either gas or aerosol
    337 !OLD
    338           if ((readclouds.and.ngrid.eq.0).or.(readclouds_this_nest.and.ngrid.gt.0)) then
     318!SEC wetscav fix, the cloud height is no longer needed, it gives wrong results
    339319            wetscav=incloud_ratio*S_i*(prec(1)/3.6E6)
    340           else
    341 !SEC wetscav fix
    342              wetscav=incloud_ratio*S_i*(prec(1)/3.6E6)
    343 !            wetscav=incloud_ratio*S_i*(prec(1)/3.6E6)/clouds_h
    344           endif
    345320        endif ! positive in-cloud scavenging parameters given in Species file
    346321      endif !incloud
  • src/netcdf_output_mod.f90

    rae43937 r4ad96c5  
    6767                       ccn_aero,in_aero, & ! wetc_in,wetd_in, &
    6868                       reldiff,henry,f0,density,dquer,dsigma,dryvel,&
    69 !                       weightmolar,ohreact,spec_ass,kao,vsetaver,&
    70                        weightmolar,ohcconst,ohdconst,spec_ass,kao,vsetaver,&
     69                       weightmolar,ohcconst,ohdconst,vsetaver,&
    7170                       ! for concoutput_netcdf and concoutput_nest_netcdf
    7271                       nxmin1,nymin1,nz,oro,oron,rho,rhon,&
     
    114113  logical, parameter :: write_area = .false.
    115114
     115  ! coordinate transformation from internal to world coord
     116  real :: xp1,yp1,xp2,yp2
    116117contains
    117118
     
    512513        call nf90_err(nf90_put_att(ncid, sID, 'ohcconst', ohcconst(i)))
    513514        call nf90_err(nf90_put_att(ncid, sID, 'ohdconst', ohdconst(i)))
    514         call nf90_err(nf90_put_att(ncid, sID, 'kao', kao(i)))
    515515        call nf90_err(nf90_put_att(ncid, sID, 'vsetaver', vsetaver(i)))
    516         call nf90_err(nf90_put_att(ncid, sID, 'spec_ass', spec_ass(i)))
    517516
    518517        if (lnest) then
     
    535534        call nf90_err(nf90_put_att(ncid, sID, 'ohcconst', ohcconst(i)))
    536535        call nf90_err(nf90_put_att(ncid, sID, 'ohdconst', ohdconst(i)))
    537         call nf90_err(nf90_put_att(ncid, sID, 'kao', kao(i)))
    538536        call nf90_err(nf90_put_att(ncid, sID, 'vsetaver', vsetaver(i)))
    539         call nf90_err(nf90_put_att(ncid, sID, 'spec_ass', spec_ass(i)))
    540537
    541538        if (lnest) then
     
    673670       call nf90_err(nf90_put_var(ncid, relendID, ireleaseend(i), (/i/)))
    674671       call nf90_err(nf90_put_var(ncid, relkindzID, kindz(i), (/i/)))
    675        call nf90_err(nf90_put_var(ncid, rellng1ID, xpoint1(i), (/i/)))
    676        call nf90_err(nf90_put_var(ncid, rellng2ID, xpoint2(i), (/i/)))
    677        call nf90_err(nf90_put_var(ncid, rellat1ID, ypoint1(i), (/i/)))
    678        call nf90_err(nf90_put_var(ncid, rellat2ID, ypoint2(i), (/i/)))
     672       xp1=xpoint1(i)*dx+xlon0
     673       yp1=ypoint1(i)*dy+ylat0
     674       xp2=xpoint2(i)*dx+xlon0
     675       yp2=ypoint2(i)*dy+ylat0
     676       call nf90_err(nf90_put_var(ncid, rellng1ID, xp1, (/i/)))
     677       call nf90_err(nf90_put_var(ncid, rellng2ID, xp2, (/i/)))
     678       call nf90_err(nf90_put_var(ncid, rellat1ID, yp1, (/i/)))
     679       call nf90_err(nf90_put_var(ncid, rellat2ID, yp2, (/i/)))
    679680       call nf90_err(nf90_put_var(ncid, relzz1ID, zpoint1(i), (/i/)))
    680681       call nf90_err(nf90_put_var(ncid, relzz2ID, zpoint2(i), (/i/)))
  • src/readspecies.f90

    raa8c34a rd6a245b  
    6767  character(len=16) :: pspecies
    6868  real :: pdecay, pweta_gas, pwetb_gas, preldiff, phenry, pf0, pdensity, pdquer
    69   real :: pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, pkao
     69  real :: pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst
    7070  real :: pcrain_aero, pcsnow_aero, pccn_aero, pin_aero
    71   integer :: readerror, pspec_ass
     71  real :: parea_dow(7), parea_hour(24), ppoint_dow(7), ppoint_hour(24)
     72  integer :: readerror
    7273
    7374! declare namelist
     
    7677       pcrain_aero, pcsnow_aero, pccn_aero, pin_aero, &
    7778       preldiff, phenry, pf0, pdensity, pdquer, &
    78        pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, pspec_ass, pkao
     79       pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, &
     80       parea_dow, parea_hour, ppoint_dow, ppoint_hour
    7981
    8082  pspecies="" ! read failure indicator value
     
    9698  pohdconst=-9.9E-09
    9799  pohnconst=2.0
    98   pspec_ass=-9
    99   pkao=-99.99
    100100  pweightmolar=-999.9
    101101
     102  do j=1,24           ! initialize everything to no variation
     103    area_hour(pos_spec,j)=1.
     104    point_hour(pos_spec,j)=1.
     105  end do
     106  do j=1,7
     107    area_dow(pos_spec,j)=1.
     108    point_dow(pos_spec,j)=1.
     109  end do
     110
     111  if (readerror.ne.0) then ! text format input
    102112! Open the SPECIES file and read species names and properties
    103113!************************************************************
     
    164174    read(unitspecies,'(f8.2)',end=22) ohnconst(pos_spec)
    165175!  write(*,*) ohnconst(pos_spec)
    166     read(unitspecies,'(i18)',end=22) spec_ass(pos_spec)
    167 !  write(*,*) spec_ass(pos_spec)
    168     read(unitspecies,'(f18.2)',end=22) kao(pos_spec)
    169 !       write(*,*) kao(pos_spec)
     176
     177! Read in daily and day-of-week variation of emissions, if available
     178!*******************************************************************
     179
     180    read(unitspecies,*,end=22)
     181    do j=1,24     ! 24 hours, starting with 0-1 local time
     182      read(unitspecies,*) ihour,area_hour(pos_spec,j),point_hour(pos_spec,j)
     183    end do
     184    read(unitspecies,*)
     185    do j=1,7      ! 7 days of the week, starting with Monday
     186      read(unitspecies,*) idow,area_dow(pos_spec,j),point_dow(pos_spec,j)
     187    end do
    170188
    171189    pspecies=species(pos_spec)
     
    188206    pohdconst=ohdconst(pos_spec)
    189207    pohnconst=ohnconst(pos_spec)
    190     pspec_ass=spec_ass(pos_spec)
    191     pkao=kao(pos_spec)
    192 
    193   else
     208
     209
     210    do j=1,24     ! 24 hours, starting with 0-1 local time
     211      parea_hour(j)=area_hour(pos_spec,j)
     212      ppoint_hour(j)=point_hour(pos_spec,j)
     213    end do
     214    do j=1,7      ! 7 days of the week, starting with Monday
     215      parea_dow(j)=area_dow(pos_spec,j)
     216      ppoint_dow(j)=point_dow(pos_spec,j)
     217    end do
     218
     219  else ! namelist available
    194220
    195221    species(pos_spec)=pspecies
     
    212238    ohdconst(pos_spec)=pohdconst
    213239    ohnconst(pos_spec)=pohnconst
    214     spec_ass(pos_spec)=pspec_ass
    215     kao(pos_spec)=pkao
     240
     241    do j=1,24     ! 24 hours, starting with 0-1 local time
     242      area_hour(pos_spec,j)=parea_hour(j)
     243      point_hour(pos_spec,j)=ppoint_hour(j)
     244    end do
     245    do j=1,7      ! 7 days of the week, starting with Monday
     246      area_dow(pos_spec,j)=parea_dow(j)
     247      point_dow(pos_spec,j)=ppoint_dow(j)
     248    end do
    216249
    217250  endif
     
    303336  end if
    304337
    305   if (spec_ass(pos_spec).gt.0) then
    306     spec_found=.FALSE.
    307     do j=1,pos_spec-1
    308       if (spec_ass(pos_spec).eq.specnum(j)) then
    309         spec_ass(pos_spec)=j
    310         spec_found=.TRUE.
    311         ASSSPEC=.TRUE.
    312       endif
    313     end do
    314     if (spec_found.eqv..false.) then
    315       goto 997
    316     endif
    317   endif
    318 
    319   if (dsigma(i).eq.1.) dsigma(i)=1.0001   ! avoid floating exception
    320338  if (dsigma(i).eq.0.) dsigma(i)=1.0001   ! avoid floating exception
    321339
     
    330348
    331349
    332 ! Read in daily and day-of-week variation of emissions, if available
    333 !*******************************************************************
    334 ! HSO: This is not yet implemented as namelist parameters
    335 
    336   do j=1,24           ! initialize everything to no variation
    337     area_hour(i,j)=1.
    338     point_hour(i,j)=1.
    339   end do
    340   do j=1,7
    341     area_dow(i,j)=1.
    342     point_dow(i,j)=1.
    343   end do
    344 
    345   if (readerror.ne.0) then ! text format input
    346 
    347     read(unitspecies,*,end=22)
    348     do j=1,24     ! 24 hours, starting with 0-1 local time
    349       read(unitspecies,*) ihour,area_hour(i,j),point_hour(i,j)
    350     end do
    351     read(unitspecies,*)
    352     do j=1,7      ! 7 days of the week, starting with Monday
    353       read(unitspecies,*) idow,area_dow(i,j),point_dow(i,j)
    354     end do
    355 
    356350  endif
    357351
  • src/timemanager.f90

    ra9cf4b1 rc7d1052  
    449449
    450450        !write(*,46) float(itime)/3600,itime,numpart
    451 45      format(i13,' SECONDS SIMULATED: ',i13, ' PARTICLES:    Uncertainty: ',3f7.3)
     45145      format(i13,' Seconds simulated: ',i13, ' Particles:    Uncertainty: ',3f7.3)
    45245246      format(' Simulated ',f7.1,' hours (',i13,' s), ',i13, ' particles')
    453453        if (ipout.ge.1) call partoutput(itime)    ! dump particle positions
  • .gitignore

    r027e844 rf976c5d  
    55.DS_Store
    66output
     7FLEXPART
     8FLEXPART_MPI
  • README.md

    r3d04845 r6616daf  
    2121```
    2222> cd src
     23
     24change the appropriate paths to libraries in the makefile
     25
    2326> make
    2427```
     
    2629* Deployment instructions
    2730
    28    FLEXPART is a standalone executable   
     31   FLEXPART is a standalone executable 
     32   The necessary ECMWF wind fields can be obtained testing flex_ecmwf
     33   The AVAILABLE file works with the default ERA 5 retrieved winds
     34   In the winds are available in flex_ecmwf/work it should suffice to execute
     35   ./src/FLEXPART in the main directory 
    2936
    3037### Contribution guidelines ###
  • create_tarball.sh

    r6d420c3 r6616daf  
    11#!/bin/bash
    22
    3 version=10.2
    4 tarball_tmp=flexpart$version
     3#define version number
     4version=10.3beta
    55
     6# define tarball name
     7tarball_tmp=flexpart_v$version
     8
     9# clean old package
     10rm -r $tarball_tmp
     11
     12# create basic dir structure
    613mkdir $tarball_tmp
    714mkdir $tarball_tmp/src
     15
     16# copy fortran source files
    817cp src/*.f90 $tarball_tmp/src
     18cp -r src/gributils $tarball_tmp/src
     19
     20# copy makefile
    921cp src/makefile $tarball_tmp/src
    10 cp src/makefile.gfs $tarball_tmp/src
    11 cp -r options $tarball_tmp
     22#cp src/makefile.gfs $tarball_tmp/src
     23
     24# copy default options (for the distribution they work with the defult flex_ecmwf test winds)
     25# cp -r options $tarball_tmp
     26cp -r options_flex_ecmwf_EA $tarball_tmp/options
     27
     28# copy default pathnames
     29cp pathnames_distribution $tarball_tmp/pathnames
     30
     31
     32# add ECMWF retrieve routines
     33mkdir $tarball_tmp/preprocess
     34#mkdir $tarball_tmp/preprocess/flex_ecmwf
     35mkdir $tarball_tmp/preprocess/flex_extract
     36#cp -r flex_ecmwf_src/* $tarball_tmp/preprocess/flex_ecmwf/
     37cp -r flex_ecmwf_src/* $tarball_tmp/preprocess/flex_extract/
     38
     39# copy default AVAILABLE
     40cp AVAILABLE_flex_ecmwf_EA $tarball_tmp/AVAILABLE
     41
     42
     43# directory for reading routines
     44mkdir $tarball_tmp/postprocess
     45
     46# add fortran reading routines
     47mkdir $tarball_tmp/postprocess/flex_read_fortran
     48cp postprocess/flex_read_fortran/*.f $tarball_tmp/postprocess/flex_read_fortran
     49cp postprocess/flex_read_fortran/*.f90 $tarball_tmp/postprocess/flex_read_fortran
     50cp postprocess/flex_read_fortran/makefile $tarball_tmp/postprocess/flex_read_fortran
     51
     52# add matlab reading routines
     53mkdir $tarball_tmp/postprocess/flex_read_matlab
     54cp postprocess/flex_read_matlab/*.m $tarball_tmp/postprocess/flex_read_matlab
     55
     56# examples
     57cp -r examples $tarball_tmp/
     58
    1259
    1360mkdir $tarball_tmp/tests
     61
    1462cp -r tests/NILU/test_1 $tarball_tmp/tests/
     63cp -r tests/flex_gen_cases $tarball_tmp/tests/
     64
    1565
    1666#return
    1767
    18 tar -cvf flexpart$version.tar $tarball_tmp/*
    19 
     68#tar -cvf flexpart$version.tar $tarball_tmp/*
     69echo now can run:  "tar -cvf $tarball_tmp.tar $tarball_tmp/*"
     70echo preliminary: scp flexpart_v10.3beta.tar njord:public_html/FLEXPART
  • options.reference/AGECLASSES

    • Property mode changed from 100755 to 100644
    rdc69f6a r0ff3b23  
    1212*                                              *
    1313************************************************
    14 1          Integer        Number of age classes
    15 1728000 
     141          Integer        Number of age classes 
     151728000
  • options.reference/AGECLASSES.oldformat

    • Property mode changed from 100755 to 100644
    rca350ba r0ff3b23  
    1212*                                              *
    1313************************************************
    14 1          Integer        Number of age classes
    15 1728000 
     141          Integer        Number of age classes 
     151728000
  • options.reference/COMMAND

    r2bfec12 r0ff3b23  
    11&COMMAND
    22 LDIRECT=          1,
    3  IBDATE=   20140911,
    4  IBTIME=     050000,
    5  IEDATE=   20140911,
    6  IETIME=     150000,
     3 IBDATE=   20170102,
     4 IBTIME=     060000,
     5 IEDATE=   20170102,
     6 IETIME=     120000,
    77 LOUTSTEP=       3600,
    88 LOUTAVER=       3600,
    99 LOUTSAMPLE=        900,
    10  ITSPLIT=  99997200,
     10 ITSPLIT=  99999999,
    1111 LSYNCTIME=        900,
    1212 CTL= -5.0000000    ,
    1313 IFINE=          4,
    14  IOUT=          1,
     14 IOUT=    1,
    1515 IPOUT=          0,
    16  LSUBGRID=          1,
    17  LCONVECTION=          1,
     16 LSUBGRID=          0,
     17 LCONVECTION=          0,
    1818 LAGESPECTRA=          0,
    1919 IPIN=          0,
     
    2424 IND_RECEPTOR=          1,
    2525 MQUASILAG=          0,
    26  NESTED_OUTPUT=          1,
     26 NESTED_OUTPUT=          0,
    2727 LINIT_COND=          0,
    2828 SURF_ONLY=          0,
  • options.reference/IGBP_int1.dat

    • Property mode changed from 100755 to 100644
  • options.reference/OH_7lev_agl.dat

    • Property mode changed from 100755 to 100644
  • options.reference/OUTGRID

    • Property mode changed from 100755 to 100644
    r4bf4a69 r0ff3b23  
    1515
    16163.  -----             4X,I5
    17       360               NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
     17     720               NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
    1818    NUMXGRID
    1919
    20204.  -----             4X,I5
    21       180             NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
     21     360             NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
    2222    NUMYGRID
    2323
    24245.  ------.---        4X,F10.3
    25          1.000        GRID DISTANCE IN X DIRECTION
     25         0.50        GRID DISTANCE IN X DIRECTION
    2626    DXOUTLON
    2727
    28286.  ------.---        4X,F10.3
    29          1.000        GRID DISTANCE IN Y DIRECTION
     29         0.50        GRID DISTANCE IN Y DIRECTION
    3030    DYOUTLAT
    3131
    323210. -----.-           4X, F7.1
    3333      100.0
     34    LEVEL 1           HEIGHT OF LEVEL (UPPER BOUNDARY)
     35
     3610. -----.-           4X, F7.1
     37      500.0
     38    LEVEL 2           HEIGHT OF LEVEL (UPPER BOUNDARY)
     39
     4010. -----.-           4X, F7.1
     41     1000.0
     42    LEVEL 3           HEIGHT OF LEVEL (UPPER BOUNDARY)
     43
     4410. -----.-           4X, F7.1
     45    10000.0
    3446    LEVEL 4           HEIGHT OF LEVEL (UPPER BOUNDARY)
    3547
    364810. -----.-           4X, F7.1
    37      3000.0
    38     LEVEL 4           HEIGHT OF LEVEL (UPPER BOUNDARY)
    39 
    40 10. -----.-           4X, F7.1
    41     50000.0
    42     LEVEL 4           HEIGHT OF LEVEL (UPPER BOUNDARY)
    43 
     49    40000.0
     50    LEVEL 5           HEIGHT OF LEVEL (UPPER BOUNDARY)
  • options.reference/OUTGRID_NEST

    rca350ba r0ff3b23  
    77
    881.  ------.----       4X,F11.4
    9          1.0000       GEOGRAFICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
     9        -2.0000       GEOGRAFICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
    1010    OUTLONLEFT        (left boundary of the first grid cell - not its centre)
    1111
    12122.  ------.----       4X,F11.4
    13         39.5000       GEOGRAFICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
     13         1.0000       GEOGRAFICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
    1414    OUTLATLOWER       (lower boundary of the first grid cell - not its centre)
    1515
    16163.  -----             4X,I5
    17       220             NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
     17      400             NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
    1818    NUMXGRID
    1919
    20204.  -----             4X,I5
    21       220             NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
     21      200             NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
    2222    NUMYGRID
    2323
    24245.  ------.-----      4X,F12.5
    25          0.02000      GRID DISTANCE IN X DIRECTION
     25         0.01000      GRID DISTANCE IN X DIRECTION
    2626    DXOUTLON
    2727
    28286.  ------.-----      4X,F12.5
    29          0.02000      GRID DISTANCE IN Y DIRECTION
     29         0.01000      GRID DISTANCE IN Y DIRECTION
    3030    DYOUTLAT
  • options.reference/RECEPTORS

    • Property mode changed from 100755 to 100644
    rca350ba r0ff3b23  
    1111
    12122.  ------.----       4X,F11.4
    13          6.1333       GEOGRAFICAL LONGITUDE
     13        21.1333       GEOGRAFICAL LONGITUDE
    1414    XRECEPTOR
    1515
    16163.  ------.----       4X,F11.4
    17         49.0833       GEOGRAFICAL LATITUDE
     17        38.0833       GEOGRAFICAL LATITUDE
    1818    YRECEPTOR 
    1919================================================================================
  • options.reference/RECEPTORS.oldformat

    • Property mode changed from 100755 to 100644
    rca350ba r0ff3b23  
    1111
    12122.  ------.----       4X,F11.4
    13          6.1333       GEOGRAFICAL LONGITUDE
     13        21.1333       GEOGRAFICAL LONGITUDE
    1414    XRECEPTOR
    1515
    16163.  ------.----       4X,F11.4
    17         49.0833       GEOGRAFICAL LATITUDE
     17        38.0833       GEOGRAFICAL LATITUDE
    1818    YRECEPTOR 
    1919================================================================================
  • options.reference/surfdata.t

    • Property mode changed from 100755 to 100644
  • options.reference/surfdepo.t

    • Property mode changed from 100755 to 100644
  • options/AGECLASSES

    rf9ce123 r547b188  
    1212*                                              *
    1313************************************************
    14 1          Integer        Number of age classes
    15 1728000
     14&AGECLASS
     15 NAGECLASS=          1,
     16 LAGE=    1728000,
     17 /
  • options/OUTGRID

    r5f42c27 r9397d0f  
    1 ********************************************************************************
    2 *                                                                              *
    3 *      Input file for the Lagrangian particle dispersion model FLEXPART        *
    4 *                       Please specify your output grid                        *
    5 *                                                                              *
    6 ********************************************************************************
    7 
    8 1.  ------.----       4X,F11.4
    9       -179.0000       GEOGRAFICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
    10     OUTLONLEFT        (left boundary of the first grid cell - not its centre)
    11 
    12 2.  ------.----       4X,F11.4
    13        -90.0000       GEOGRAFICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
    14     OUTLATLOWER       (lower boundary of the first grid cell - not its centre)
    15 
    16 3.  -----             4X,I5
    17      720               NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
    18     NUMXGRID
    19 
    20 4.  -----             4X,I5
    21      360             NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
    22     NUMYGRID
    23 
    24 5.  ------.---        4X,F10.3
    25          0.50        GRID DISTANCE IN X DIRECTION
    26     DXOUTLON
    27 
    28 6.  ------.---        4X,F10.3
    29          0.50        GRID DISTANCE IN Y DIRECTION
    30     DYOUTLAT
    31 
    32 10. -----.-           4X, F7.1
    33      1000.0
    34     LEVEL 1           HEIGHT OF LEVEL (UPPER BOUNDARY)
    35 
    36 10. -----.-           4X, F7.1
    37      2000.0
    38     LEVEL 2           HEIGHT OF LEVEL (UPPER BOUNDARY)
    39 
    40 10. -----.-           4X, F7.1
    41      3000.0
    42     LEVEL 3           HEIGHT OF LEVEL (UPPER BOUNDARY)
    43 
    44 10. -----.-           4X, F7.1
    45      4000.0
    46     LEVEL 4           HEIGHT OF LEVEL (UPPER BOUNDARY)
    47 
    48 10. -----.-           4X, F7.1
    49      5000.0
    50     LEVEL 5           HEIGHT OF LEVEL (UPPER BOUNDARY)
    51 
    52 10. -----.-           4X, F7.1
    53      6000.0
    54     LEVEL 6           HEIGHT OF LEVEL (UPPER BOUNDARY)
    55 
    56 10. -----.-           4X, F7.1
    57      7000.0
    58     LEVEL 7           HEIGHT OF LEVEL (UPPER BOUNDARY)
    59 
    60 10. -----.-           4X, F7.1
    61      8000.0
    62     LEVEL 8           HEIGHT OF LEVEL (UPPER BOUNDARY)
    63 
    64 10. -----.-           4X, F7.1
    65      9000.0
    66     LEVEL 9           HEIGHT OF LEVEL (UPPER BOUNDARY)
    67 
    68 10. -----.-           4X, F7.1
    69     10000.0
    70     LEVEL10           HEIGHT OF LEVEL (UPPER BOUNDARY)
    71 
    72 10. -----.-           4X, F7.1
    73     11000.0
    74     LEVEL11           HEIGHT OF LEVEL (UPPER BOUNDARY)
    75 
    76 10. -----.-           4X, F7.1
    77     12000.0
    78     LEVEL12           HEIGHT OF LEVEL (UPPER BOUNDARY)
    79 
    80 10. -----.-           4X, F7.1
    81     13000.0
    82     LEVEL13           HEIGHT OF LEVEL (UPPER BOUNDARY)
    83 
    84 10. -----.-           4X, F7.1
    85     14000.0
    86     LEVEL14           HEIGHT OF LEVEL (UPPER BOUNDARY)
    87 
    88 10. -----.-           4X, F7.1
    89     15000.0
    90     LEVEL15           HEIGHT OF LEVEL (UPPER BOUNDARY)
    91 
    92 10. -----.-           4X, F7.1
    93     16000.0
    94     LEVEL16           HEIGHT OF LEVEL (UPPER BOUNDARY)
    95 
    96 10. -----.-           4X, F7.1
    97     17000.0
    98     LEVEL17           HEIGHT OF LEVEL (UPPER BOUNDARY)
    99 
    100 10. -----.-           4X, F7.1
    101     18000.0
    102     LEVEL18           HEIGHT OF LEVEL (UPPER BOUNDARY)
    103 
    104 10. -----.-           4X, F7.1
    105     19000.0
    106     LEVEL19           HEIGHT OF LEVEL (UPPER BOUNDARY)
    107 
    108 10. -----.-           4X, F7.1
    109     20000.0
    110     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    111 
    112 10. -----.-           4X, F7.1
    113     21000.0
    114     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    115 
    116 10. -----.-           4X, F7.1
    117     22000.0
    118     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    119 
    120 10. -----.-           4X, F7.1
    121     23000.0
    122     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    123 
    124 10. -----.-           4X, F7.1
    125     24000.0
    126     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    127 
    128 10. -----.-           4X, F7.1
    129     25000.0
    130     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    131 
    132 10. -----.-           4X, F7.1
    133     26000.0
    134     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    135 
    136 10. -----.-           4X, F7.1
    137     27000.0
    138     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    139 
    140 10. -----.-           4X, F7.1
    141     28000.0
    142     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    143 
    144 10. -----.-           4X, F7.1
    145     29000.0
    146     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    147 
    148 10. -----.-           4X, F7.1
    149     30000.0
    150     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    151 
    152 10. -----.-           4X, F7.1
    153     40000.0
    154     LEVEL20           HEIGHT OF LEVEL (UPPER BOUNDARY)
    155 
    156 
     1!*******************************************************************************
     2!                                                                              *
     3!      Input file for the Lagrangian particle dispersion model FLEXPART        *
     4!                       Please specify your output grid                        *
     5!                                                                              *
     6! OUTLON0    = GEOGRAPHYICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID     *
     7! OUTLAT0    = GEOGRAPHYICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID      *
     8! NUMXGRID   = NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)       *
     9! NUMYGRID   = NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)       *
     10! DXOUT      = GRID DISTANCE IN X DIRECTION                                    *
     11! DYOUN      = GRID DISTANCE IN Y DIRECTION                                    *
     12! OUTHEIGHTS = HEIGHT OF LEVELS (UPPER BOUNDARY)                               *
     13!*******************************************************************************
     14&OUTGRID
     15 OUTLON0=    -25.00,
     16 OUTLAT0=     10.00,
     17 NUMXGRID=       85,
     18 NUMYGRID=       65,
     19 DXOUT=        1.00,
     20 DYOUT=        1.00,
     21 OUTHEIGHTS=  100.0, 500.0, 1000.0 , 50000.0,
     22 /
  • options/OUTGRID_NEST

    r2bfec12 r9397d0f  
    1 ********************************************************************************
    2 *                                                                              *
    3 *      Input file for the Lagrangian particle dispersion model FLEXPART        *
    4 *                       Please specify your output grid                        *
    5 *                                                                              *
    6 ********************************************************************************
    7 
    8 1.  ------.----       4X,F11.4
    9         -2.0000       GEOGRAFICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
    10     OUTLONLEFT        (left boundary of the first grid cell - not its centre)
    11 
    12 2.  ------.----       4X,F11.4
    13         -1.0000       GEOGRAFICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
    14     OUTLATLOWER       (lower boundary of the first grid cell - not its centre)
    15 
    16 3.  -----             4X,I5
    17       400             NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
    18     NUMXGRID
    19 
    20 4.  -----             4X,I5
    21       200             NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
    22     NUMYGRID
    23 
    24 5.  ------.-----      4X,F12.5
    25          0.01000      GRID DISTANCE IN X DIRECTION
    26     DXOUTLON
    27 
    28 6.  ------.-----      4X,F12.5
    29          0.01000      GRID DISTANCE IN Y DIRECTION
    30     DYOUTLAT
     1!*******************************************************************************
     2!                                                                              *
     3!      Input file for the Lagrangian particle dispersion model FLEXPART        *
     4!                       Please specify your output grid                        *
     5!                                                                              *
     6! OUTLON0    = GEOGRAPHYICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID     *
     7! OUTLAT0    = GEOGRAPHYICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID      *
     8! NUMXGRID   = NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)       *
     9! NUMYGRID   = NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)       *
     10! DXOUT      = GRID DISTANCE IN X DIRECTION                                    *
     11! DYOUN      = GRID DISTANCE IN Y DIRECTION                                    *
     12!*******************************************************************************
     13&OUTGRIDN
     14 OUTLON0N=    -25.00,
     15 OUTLAT0N=     10.00,
     16 NUMXGRIDN=       85,
     17 NUMYGRIDN=       65,
     18 DXOUTN=        1.00,
     19 DYOUTN=        1.00,
     20 /
  • options/RECEPTORS

    r2bfec12 r4bc2654  
    11&RECEPTORS
    2  RECEPTOR="receptor 1      ",
    3  LON=   0.0000000    ,
    4  LAT=   0.0000000    ,
     2 RECEPTOR="receptor 1",
     3 LON     =   0.0000000,
     4 LAT     =   0.0000000,
    55 /
    66&RECEPTORS
    7  RECEPTOR="receptor 2      ",
    8  LON=  1.00000000    ,
    9  LAT=  0.00000000    ,
     7 RECEPTOR="receptor 2",
     8 LON     =  1.00000000,
     9 LAT     =  0.00000000,
    1010 /
  • options/SPECIES/SPECIES_002

    r8d70e43 rce1e95b  
    11&SPECIES_PARAMS
    2  PSPECIES="Ozone",            ! Tracer name
     2 PSPECIES="O3",               ! Tracer name
    33 PDECAY=-9.9,                 ! Species half life
    44 PWETA_GAS=-0.9E-9,           ! Below-cloud scavenging (gases) - A (weta_gas)
  • options/SPECIES/SPECIES_003

    r8d70e43 rce1e95b  
    1313 PDRYVEL=-9.99,               ! Alternative: dry deposition velocity
    1414 PRELDIFF=1.2,                ! Dry deposition (gases) - D
    15  PHENRY=3.0E-03,              ! Dry deposition (gases) - Henrys const.
     15 PHENRY=2.0E-03,              ! Dry deposition (gases) - Henrys const.
    1616 PF0=0.0,                     ! Dry deposition (gases) - f0 (reactivity)
    1717 PWEIGHTMOLAR=30.0,           ! molweight
  • options/SPECIES/SPECIES_005

    r8d70e43 rce1e95b  
    22 PSPECIES="HNO3",             ! Tracer name
    33 PDECAY=-9.9,                 ! Species half life
    4  PWETA_GAS=5.0E-05,           ! Below-cloud scavenging (gases) - A (weta_gas)
     4 PWETA_GAS=8.0E-04,           ! Below-cloud scavenging (gases) - A (weta_gas)
    55 PWETB_GAS=0.62,              ! Below-cloud scavenging (gases) - B (wetb_gas)
    66 PCRAIN_AERO=-9.9,            ! Below-cloud scavenging (particles) - Crain (crain_aero)
     
    1414 PRELDIFF=1.9,                ! Dry deposition (gases) - D
    1515 PHENRY=1E14,                 ! Dry deposition (gases) - Henrys const.
    16  PF0=-9,                      ! Dry deposition (gases) - f0 (reactivity)
     16 PF0=0,                       ! Dry deposition (gases) - f0 (reactivity)
    1717 PWEIGHTMOLAR=63.0,           ! molweight
    1818 POHCCONST=-0.9E+9,           ! OH Reaction rate - C [cm^3/molecule/sec]
  • options/SPECIES/SPECIES_011

    r8d70e43 rce1e95b  
    1414 PRELDIFF=1.1,                ! Dry deposition (gases) - D
    1515 PHENRY=2.0E14,               ! Dry deposition (gases) - Henrys const.
    16  PF0=-9,                      ! Dry deposition (gases) - f0 (reactivity)
     16 PF0=0.0,                     ! Dry deposition (gases) - f0 (reactivity)
    1717 PWEIGHTMOLAR=17.0,           ! molweight
    1818 POHCCONST=-0.9E+9,           ! OH Reaction rate - C [cm^3/molecule/sec]
  • options/SPECIES/SPECIES_012

    r8d70e43 rce1e95b  
    1515 PHENRY=-0.9E+9,              ! Dry deposition (gases) - Henrys const.
    1616 PF0=-9,                      ! Dry deposition (gases) - f0 (reactivity)
    17  PWEIGHTMOLAR=0.0,            ! molweight
     17 PWEIGHTMOLAR=-9.9,           ! molweight
    1818 POHCCONST=-0.9E+9,           ! OH Reaction rate - C [cm^3/molecule/sec]
    1919 POHDCONST=-9.9,              ! OH Reaction rate - D [K]
  • options/SPECIES/SPECIES_013

    r8d70e43 rce1e95b  
    88 PCCN_AERO=0.9,               ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
    99 PIN_AERO=0.1,                ! In-cloud scavenging (particles) - INeff (in_aero)
    10  PDENSITY=-0.9E+9,            ! Dry deposition (particles) - rho
     10 PDENSITY=2.0E+3,             ! Dry deposition (particles) - rho
    1111 PDQUER=4.0E-07,              ! Dry deposition (particles) - dquer
    1212 PDSIGMA=0.3,                 ! Dry deposition (particles) - dsig
     
    1515 PHENRY=-0.9E+9,              ! Dry deposition (gases) - Henrys const.
    1616 PF0=-9,                      ! Dry deposition (gases) - f0 (reactivity)
    17  PWEIGHTMOLAR=2000.0,         ! molweight
     17 PWEIGHTMOLAR=-9.9,           ! molweight
    1818 POHCCONST=-0.9E+9,           ! OH Reaction rate - C [cm^3/molecule/sec]
    1919 POHDCONST=-9.9,              ! OH Reaction rate - D [K]
  • options/SPECIES/SPECIES_017

    r8d70e43 rce1e95b  
    1414 PRELDIFF=-9.9,               ! Dry deposition (gases) - D
    1515 PHENRY=-0.9E-9,              ! Dry deposition (gases) - Henrys const.
    16  PF0=-9,                      ! Dry deposition (gases) - f0 (reactivity)
     16 PF0=-9.0,                      ! Dry deposition (gases) - f0 (reactivity)
    1717 PWEIGHTMOLAR=350.0,          ! molweight
    1818 POHCCONST=-0.9E-9,           ! OH Reaction rate - C [cm^3/molecule/sec]
  • options/SPECIES/SPECIES_019

    r8d70e43 rce1e95b  
    11&SPECIES_PARAMS
    2  PSPECIES="KR-85",            ! Tracer name
     2 PSPECIES="Kr-85",            ! Tracer name
    33 PDECAY=-9.9,                 ! Species half life
    44 PWETA_GAS=-0.9E-9,           ! Below-cloud scavenging (gases) - A (weta_gas)
     
    1515 PHENRY=-0.9E-9,              ! Dry deposition (gases) - Henrys const.
    1616 PF0=-9,                      ! Dry deposition (gases) - f0 (reactivity)
    17  PWEIGHTMOLAR=-9.9,           ! molweight
     17 PWEIGHTMOLAR=83.8,           ! molweight
    1818 POHCCONST=-0.9E-9,           ! OH Reaction rate - C [cm^3/molecule/sec]
    1919 POHDCONST=-9.9,              ! OH Reaction rate - D [K]
  • options/SPECIES/SPECIES_025

    r8d70e43 rce1e95b  
    1515 PHENRY=-0.9E-9,              ! Dry deposition (gases) - Henrys const.
    1616 PF0=-9,                      ! Dry deposition (gases) - f0 (reactivity)
    17  PWEIGHTMOLAR=29.0,           ! molweight
     17 PWEIGHTMOLAR=-9.9,           ! molweight
    1818 POHCCONST=-0.9E-9,           ! OH Reaction rate - C [cm^3/molecule/sec]
    1919 POHDCONST=-9.9,              ! OH Reaction rate - D [K]
  • options/SPECIES/SPECIES_040

    r8d70e43 rce1e95b  
    11&SPECIES_PARAMS
    22 PSPECIES="BC",               ! Tracer name
    3  PDECAY=-999.99,              ! Species half life
     3 PDECAY=-9.9,                 ! Species half life
    44 PWETA_GAS=-0.9E-9,           ! Below-cloud scavenging (gases) - A (weta_gas)
    55 PWETB_GAS=-9.9,              ! Below-cloud scavenging (gases) - B (wetb_gas)
  • pathnames_ECMWF_OPER

    rca350ba r898f494  
    44/xnilu_wrk/flex_wrk/WIND_FIELDS/AVAILABLE_ECMWF_OPER_fields_global
    55============================================
    6 
     6                                           
  • src/par_mod.f90

    r79e0349 rdf96ea65  
    148148  ! ECMWF
    149149! integer,parameter :: nxmax=361,nymax=181,nuvzmax=92,nwzmax=92,nzmax=92,nxshift=359 ! 1.0 degree 92 level
     150!  integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0 ! 1.0 degree 138 level
    150151!   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359 ! 1.0 degree 138 level
    151152! integer,parameter :: nxmax=721,nymax=361,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359  ! 0.5 degree 138 level
     
    153154
    154155! GFS
    155    integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0
     156   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138
     157   integer :: nxshift=0
    156158
    157159
  • src/pathnames

    rf28aa0a ra712161  
    11../options/
    2 ../output/
     2./output/
    33/
    44/xnilu_wrk/flex_wrk/WIND_FIELDS/AVAILABLE_ECMWF_OPER_fields_global
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG