Changeset df96ea65 in flexpart.git


Ignore:
Timestamp:
Dec 20, 2018, 1:34:38 PM (3 months ago)
Author:
Ignacio Pisso <ip@…>
Branches:
release-10
Children:
2753a5c
Parents:
6d73c4b (diff), db91eb7 (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:

resolve merge conflicts with clouds_in_GFS

Files:
36 added
179 deleted
26 edited
16 moved

Legend:

Unmodified
Added
Removed
  • src/get_wetscav.f90

    r6985a98 r79e0349  
    298298          else                                  !parameterize cloudwater m2/m3
    299299!ZHG updated parameterization of cloud water to better reproduce the values coming from ECMWF
    300             cl=1.6E-6*prec(1)**0.36
     300! sec test
     301!           cl=1E6*1E-7*prec(1)**0.3 !Sec GFS new
     302            cl=1E6*2E-7*prec(1)**0.36 !Sec ECMWF new
     303!           cl=2E-7*prec(1)**0.36 !Andreas
     304!           cl=1.6E-6*prec(1)**0.36 !Henrik
    301305          endif
    302306
     
    318322          if (dquer(ks).gt.0.) then
    319323            S_i= frac_act/cl
     324!           write(*,*) 'Si: ',S_i
    320325
    321326! GAS
     
    334339            wetscav=incloud_ratio*S_i*(prec(1)/3.6E6)
    335340          else
    336             wetscav=incloud_ratio*S_i*(prec(1)/3.6E6)/clouds_h
     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
    337344          endif
    338345        endif ! positive in-cloud scavenging parameters given in Species file
  • src/par_mod.f90

    rd944d39 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
     150!  integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0 ! 1.0 degree 138 level
     151!   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
    152153!  integer,parameter :: nxmax=181,nymax=91,nuvzmax=92,nwzmax=92,nzmax=92,nxshift=0  ! CERA 2.0 degree 92 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/readwind_gfs.f90

    r6ecb30a rdb91eb7  
    1 !**********************************************************************
     1
    22! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010         *
    33! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa,             *
     
    108108  character(len=24) :: gribErrorMsg = 'Error reading grib file'
    109109  character(len=20) :: gribFunction = 'readwind_gfs'
     110  character(len=20) :: shortname
    110111
    111112
     
    120121
    121122  !HSO
    122 5   call grib_open_file(ifile,path(3)(1:length(3)) &
     123  call grib_open_file(ifile,path(3)(1:length(3)) &
    123124         //trim(wfname(indj)),'r',iret)
    124125  if (iret.ne.GRIB_SUCCESS) then
     
    162163
    163164  !read the grib2 identifiers
     165  call grib_get_string(igrib,'shortName',shortname,iret)
     166
    164167  call grib_get_int(igrib,'discipline',discipl,iret)
    165168!  call grib_check(iret,gribFunction,gribErrorMsg)
     
    173176       valSurf,iret)
    174177!  call grib_check(iret,gribFunction,gribErrorMsg)
    175 
     178 
     179!  write(*,*) 'Field: ',ifield,parCat,parNum,typSurf,shortname
    176180  !convert to grib1 identifiers
    177181  isec1(6)=-1
     
    214218    isec1(7)=105         ! indicatorOfTypeOfLevel
    215219    isec1(8)=10
     220  elseif ((parCat.eq.1).and.(parNum.eq.22).and.(typSurf.eq.100)) then ! CLWMR Cloud Mixing Ratio [kg/kg]:
     221    isec1(6)=153         ! indicatorOfParameter
     222    isec1(7)=100         ! indicatorOfTypeOfLevel
     223    isec1(8)=valSurf/100 ! level, convert to hPa
    216224  elseif ((parCat.eq.3).and.(parNum.eq.1).and.(typSurf.eq.101)) then ! SLP
    217225    isec1(6)=2           ! indicatorOfParameter
     
    548556        endif
    549557      endif
     558! SEC & IP 12/2018 read GFS clouds
     559      if(isec1(6).eq.153) then  !! CLWCR  Cloud liquid water content [kg/kg]
     560        clwch(i,j,nlev_ec-k+2,n)=zsec4(nxfield*(ny-j-1)+i+1)
     561        readclouds=.true.
     562        sumclouds=.true.
     563      endif
     564
    550565
    551566    end do
     
    675690    call shift_field(vvh,nxfield,ny,nuvzmax,nuvz,1,1)
    676691    call shift_field(wwh,nxfield,ny,nwzmax,nwz,1,1)
     692! IP & SEC adding GFS Clouds 20181205
     693    call shift_field(clwch,nxfield,ny,nuvzmax,nuvz,2,n)
    677694  endif
    678695
  • src/verttransform_ecmwf.f90

    r47f96e5 r79e0349  
    108108  ! real :: rcw(0:nxmax-1,0:nymax-1)
    109109  ! real :: rpc(0:nxmax-1,0:nymax-1)
    110   ! character(len=60) :: zhgpath='/xnilu_wrk/flex_wrk/zhg/'
     110  character(len=60) :: zhgpath='/xnilu_wrk/users/sec/kleinprojekte/hertlfit/'
     111  character(len=60) :: fnameH,fnameI,fnameJ
    111112  ! character(len=60) :: fnameA,fnameB,fnameC,fnameD,fnameE,fnameF,fnameG,fnameH
    112   ! CHARACTER(LEN=3)  :: aspec
    113   ! integer :: virr=0
     113  CHARACTER(LEN=3)  :: aspec
     114  integer :: virr=0
    114115  !real :: tot_cloud_h
    115116  !real :: dbg_height(nzmax)
     
    732733     !********* TEST ************'**
    733734!teller(:)=0
    734 !virr=virr+1
    735 !WRITE(aspec, '(i3.3)'), virr
     735virr=virr+1
     736WRITE(aspec, '(i3.3)'), virr
    736737
    737738!if (readclouds) then
     
    770771!fnameF=trim(zhgpath)//trim(aspec)//'lsp.txt'
    771772!fnameG=trim(zhgpath)//trim(aspec)//'convp.txt'
     773if (1.eq.2) then
     774fnameH=trim(zhgpath)//trim(aspec)//'tcwc.txt'
     775fnameI=trim(zhgpath)//trim(aspec)//'prec.txt'
     776fnameJ=trim(zhgpath)//trim(aspec)//'cloudsh.txt'
     777write(*,*) 'Writing data to file: ',fnameH
    772778!if (readclouds) then
    773779!OPEN(UNIT=111, FILE=fnameA,FORM='FORMATTED',STATUS = 'UNKNOWN')
     
    776782!OPEN(UNIT=114, FILE=fnameD,FORM='FORMATTED',STATUS = 'UNKNOWN')
    777783!else
    778 !OPEN(UNIT=115, FILE=fnameE,FORM='FORMATTED',STATUS = 'UNKNOWN')
    779 !OPEN(UNIT=116, FILE=fnameF,FORM='FORMATTED',STATUS = 'UNKNOWN')
    780 !OPEN(UNIT=117, FILE=fnameG,FORM='FORMATTED',STATUS = 'UNKNOWN')
     784OPEN(UNIT=115, FILE=fnameH,FORM='FORMATTED',STATUS = 'UNKNOWN')
     785OPEN(UNIT=116, FILE=fnameI,FORM='FORMATTED',STATUS = 'UNKNOWN')
     786OPEN(UNIT=117, FILE=fnameJ,FORM='FORMATTED',STATUS = 'UNKNOWN')
    781787!endif
    782788!
    783 !do ix=0,nxmin1
     789do ix=0,nxmin1
    784790!if (readclouds) then
    785791!write(111,*) (icloud_stats(ix,jy,1,n),jy=0,nymin1)
     
    788794!write(114,*) (icloud_stats(ix,jy,4,n),jy=0,nymin1)
    789795!else
    790 !write(115,*) (cloudsh(ix,jy,n),jy=0,nymin1)    !integer
    791 !write(116,*) (lsprec(ix,jy,1,n),jy=0,nymin1)   !7.83691406E-02
    792 !write(117,*) (convprec(ix,jy,1,n),jy=0,nymin1) !5.38330078E-02
     796write(115,*) (ctwc(ix,jy,n),jy=0,nymin1) 
     797write(116,*) (lsprec(ix,jy,1,n)+convprec(ix,jy,1,n),jy=0,nymin1)
     798write(117,*) (cloudsh(ix,jy,n),jy=0,nymin1)
    793799!endif
    794 !end do
     800end do
    795801!
    796802!if (readclouds) then
     
    800806!CLOSE(114)
    801807!else
    802 !CLOSE(115)
    803 !CLOSE(116)
    804 !CLOSE(117)
     808CLOSE(115)
     809CLOSE(116)
     810CLOSE(117)
     811endif
    805812!endif
    806813!
  • src/verttransform_gfs.f90

    r6ecb30a rdb91eb7  
    7575  integer :: rain_cloud_above,kz_inv
    7676  real :: f_qvsat,pressure
    77   real :: rh,lsp,convp
     77  real :: rh,lsp,cloudh_min,convp,prec
    7878  real :: rhoh(nuvzmax),pinmconv(nzmax)
    7979  real :: ew,pint,tv,tvold,pold,dz1,dz2,dz,ui,vi
     
    224224      tt(ix,jy,1,n)=tth(ix,jy,llev,n)
    225225      qv(ix,jy,1,n)=qvh(ix,jy,llev,n)
     226! IP & SEC, 201812 add clouds
     227      if (readclouds) then
     228         clwc(ix,jy,1,n)=clwch(ix,jy,llev,n)
     229      endif
    226230      pv(ix,jy,1,n)=pvh(ix,jy,llev)
    227231      rho(ix,jy,1,n)=rhoh(llev)
     
    231235      tt(ix,jy,nz,n)=tth(ix,jy,nuvz,n)
    232236      qv(ix,jy,nz,n)=qvh(ix,jy,nuvz,n)
     237! IP & SEC, 201812 add clouds
     238      if (readclouds) then
     239         clwc(ix,jy,nz,n)=clwch(ix,jy,nuvz,n)
     240      endif
    233241      pv(ix,jy,nz,n)=pvh(ix,jy,nuvz)
    234242      rho(ix,jy,nz,n)=rhoh(nuvz)
     
    242250            tt(ix,jy,iz,n)=tt(ix,jy,nz,n)
    243251            qv(ix,jy,iz,n)=qv(ix,jy,nz,n)
     252! IP & SEC, 201812 add clouds
     253            if (readclouds) then
     254               clwc(ix,jy,iz,n)=clwc(ix,jy,nz,n)
     255            endif
    244256            pv(ix,jy,iz,n)=pv(ix,jy,nz,n)
    245257            rho(ix,jy,iz,n)=rho(ix,jy,nz,n)
     
    258270            qv(ix,jy,iz,n)=(qvh(ix,jy,kz-1,n)*dz2 &
    259271            +qvh(ix,jy,kz,n)*dz1)/dz
     272! IP & SEC, 201812 add clouds
     273            if (readclouds) then
     274               clwc(ix,jy,iz,n)=(clwch(ix,jy,kz-1,n)*dz2 &
     275               +clwch(ix,jy,kz,n)*dz1)/dz
     276            endif
    260277            pv(ix,jy,iz,n)=(pvh(ix,jy,kz-1)*dz2+pvh(ix,jy,kz)*dz1)/dz
    261278            rho(ix,jy,iz,n)=(rhoh(kz-1)*dz2+rhoh(kz)*dz1)/dz
     
    496513
    497514
     515
     516!***********************************************************************************
     517! IP & SEC, 201812 GFS clouds read
     518  if (readclouds) then
     519! The method is loops all grids vertically and constructs the 3D matrix for clouds
     520! Cloud top and cloud bottom gid cells are assigned as well as the total column
     521! cloud water. For precipitating grids, the type and whether it is in or below
     522! cloud scavenging are assigned with numbers 2-5 (following the old metod).
     523! Distinction is done for lsp and convp though they are treated the same in regards
     524! to scavenging. Also clouds that are not precipitating are defined which may be
     525! to include future cloud processing by non-precipitating-clouds.
     526!***********************************************************************************
     527    write(*,*) 'Global NCEP fields: using cloud water'
     528    clw(:,:,:,n)=0.0
     529    ctwc(:,:,n)=0.0
     530    clouds(:,:,:,n)=0
     531! If water/ice are read separately into clwc and ciwc, store sum in clwc
     532    do jy=0,nymin1
     533      do ix=0,nxmin1
     534        lsp=lsprec(ix,jy,1,n)
     535        convp=convprec(ix,jy,1,n)
     536        prec=lsp+convp
     537! Find clouds in the vertical
     538        do kz=1, nz-1 !go from top to bottom
     539          if (clwc(ix,jy,kz,n).gt.0) then
     540! assuming rho is in kg/m3 and hz in m gives: kg/kg * kg/m3 *m3/kg /m = m2/m3
     541            clw(ix,jy,kz,n)=(clwc(ix,jy,kz,n)*rho(ix,jy,kz,n))*(height(kz+1)-height(kz))
     542            ctwc(ix,jy,n) = ctwc(ix,jy,n)+clw(ix,jy,kz,n)
     543            cloudh_min=min(height(kz+1),height(kz))
     544          endif
     545        end do
     546
     547! If Precipitation. Define removal type in the vertical
     548        if ((lsp.gt.0.01).or.(convp.gt.0.01)) then ! cloud and precipitation
     549
     550          do kz=nz,1,-1 !go Bottom up!
     551            if (clw(ix,jy,kz,n).gt. 0) then ! is in cloud
     552              cloudsh(ix,jy,n)=cloudsh(ix,jy,n)+height(kz)-height(kz-1)
     553              clouds(ix,jy,kz,n)=1                               ! is a cloud
     554              if (lsp.ge.convp) then
     555                clouds(ix,jy,kz,n)=3                            ! lsp in-cloud
     556              else
     557                clouds(ix,jy,kz,n)=2                             ! convp in-cloud
     558              endif                                              ! convective or large scale
     559            elseif((clw(ix,jy,kz,n).le.0) .and. (cloudh_min.ge.height(kz))) then ! is below cloud
     560              if (lsp.ge.convp) then
     561                clouds(ix,jy,kz,n)=5                             ! lsp dominated washout
     562              else
     563                clouds(ix,jy,kz,n)=4                             ! convp dominated washout
     564              endif                                              ! convective or large scale
     565            endif
     566
     567            if (height(kz).ge. 19000) then                        ! set a max height for removal
     568              clouds(ix,jy,kz,n)=0
     569            endif !clw>0
     570          end do !nz
     571        endif ! precipitation
     572      end do
     573    end do
     574  else
     575  write(*,*) 'Global NCEP fields: using cloud water from Parameterization'
    498576  !   write (*,*) 'initializing clouds, n:',n,nymin1,nxmin1,nz
    499577  !   create a cloud and rainout/washout field, clouds occur where rh>80%
     
    534612    end do
    535613  end do
     614  endif  ! IP & SEC 201812, GFS clouds read
    536615
    537616
  • .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/COMMAND

    r2bfec12 r6d73c4b  
     1***************************************************************************************************************
     2*                                                                                                             *
     3*      Input file for the Lagrangian particle dispersion model FLEXPART                                       *
     4*                           Please select your options                                                        *
     5*                                                                                                             *
     6***************************************************************************************************************
    17&COMMAND
    2  LDIRECT=               1,
    3  IBDATE=         20170102,
    4  IBTIME=           060000,
    5  IEDATE=         20170102,
    6  IETIME=           120000,
    7  LOUTSTEP=           3600,
    8  LOUTAVER=           3600,
    9  LOUTSAMPLE=          900,
    10  ITSPLIT=        99999999,
    11  LSYNCTIME=           900,
    12  CTL=          -5.0000000,
    13  IFINE=                 4,
    14  IOUT=                  1,
    15  IPOUT=                 0,
    16  LSUBGRID=              0,
    17  LCONVECTION=           0,
    18  LAGESPECTRA=           0,
    19  IPIN=                  0,
    20  IOUTPUTFOREACHRELEASE= 0,
    21  IFLUX=                 0,
    22  MDOMAINFILL=           0,
    23  IND_SOURCE=            1,
    24  IND_RECEPTOR=          1,
    25  MQUASILAG=             0,
    26  NESTED_OUTPUT=         0,
    27  LINIT_COND=            0,
    28  SURF_ONLY=             0,
    29  CBLFLAG=               0,
    30  OHFIELDS_PATH= "../../flexin/",
     8 LDIRECT=               1, ! Simulation direction in time   ; 1 (forward) or -1 (backward)
     9 IBDATE=         20120101, ! Start date of the simulation   ; YYYYMMDD: YYYY=year, MM=month, DD=day 
     10 IBTIME=           060000, ! Start time of the simulation   ; HHMISS: HH=hours, MI=min, SS=sec; UTC
     11 IEDATE=         20120101, ! End date of the simulation     ; same format as IBDATE
     12 IETIME=           120000, ! End  time of the simulation    ; same format as IBTIME
     13 LOUTSTEP=           3600, ! Interval of model output; average concentrations calculated every LOUTSTEP (s) 
     14 LOUTAVER=           3600, ! Interval of output averaging (s)
     15 LOUTSAMPLE=          900, ! Interval of output sampling  (s), higher stat. accuracy with shorter intervals
     16 ITSPLIT=        99999999, ! Interval of particle splitting (s)
     17 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 
     19 IFINE=                 4, ! Reduction for time step in vertical transport, used only if CTL>1
     20 IOUT=                  1, ! Output type: [1]mass 2]pptv 3]1&2 4]plume 5]1&4, +8 for NetCDF output     
     21 IPOUT=                 0, ! Particle position output: 0]no 1]every output 2]only at end   
     22 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   
     24 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 
     26 IOUTPUTFOREACHRELEASE= 1, ! Separate output fields for each location in the RELEASE file; [0]no 1]yes
     27 IFLUX=                 0, ! Output of mass fluxes through output grid box boundaries
     28 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
     31 MQUASILAG=             0, ! Quasi-Lagrangian mode to track individual numbered particles
     32 NESTED_OUTPUT=         0, ! Output also for a nested domain
     33 LINIT_COND=            0, ! Output sensitivity to initial conditions (backwards mode only) 
     34 SURF_ONLY=             0, ! Output only for the lowest model layer, used w/ LINIT_COND=1 or 2
     35 CBLFLAG=               0, ! Skewed, not Gaussian turbulence in the convective ABL, need large CTL and IFINE
     36 OHFIELDS_PATH= "../../flexin/", ! Default path for OH file
    3137 /
  • 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/RELEASES

    r5b7ec80 r6d73c4b  
     1*************************************************************************
     2*                                                                       *
     3*                                                                       *
     4*                                                                       *
     5*   Input file for the Lagrangian particle dispersion model FLEXPART    *
     6*                        Please select your options                     *
     7*                                                                       *
     8*                                                                       *
     9*                                                                       *
     10*************************************************************************
    111&RELEASES_CTRL
    2  NSPEC=          1,
    3  SPECNUM_REL=        24,
     12 NSPEC      =           1,
     13 SPECNUM_REL=          24,
    414 /
    515&RELEASE
    6  IDATE1=   20170102,
    7  ITIME1=     090000,
    8  IDATE2=   20170102,
    9  ITIME2=     090000,
    10  LON1=    0.000    ,
    11  LON2=    0.000    ,
    12  LAT1=   0.000    ,
    13  LAT2=   0.000    ,
    14  Z1=    50.000       ,
    15  Z2=    50.000       ,
    16  ZKIND=     1,
    17  MASS=  1.0000E8    ,
    18  PARTS=   10000
    19  COMMENT="TEST1",
     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",
    2030 /
  • 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/FLEXPART.f90

    ra9cf4b1 r50958b8  
    137137           
    138138  if (verbosity.gt.0) then
     139    print*, 'nxmax=',nxmax
     140    print*, 'nymax=',nymax
     141    print*, 'nzmax=',nzmax
     142    print*,'nxshift=',nxshift
    139143    write(*,*) 'call readpaths'
    140144  endif
  • 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