Changeset 79e0349 in flexpart.git


Ignore:
Timestamp:
Dec 5, 2018, 11:50:52 AM (8 months ago)
Author:
Sabine <sabine.eckhardt@…>
Branches:
release-10
Children:
db91eb7
Parents:
2bfec12
Message:

Doing the fit for scavenging when using metfields with no clouds

Location:
src
Files:
3 edited

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

    r1f55de1 r79e0349  
    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=359 ! 1.0 degree 138 level
     150!   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359 ! 1.0 degree 138 level
    151151! integer,parameter :: nxmax=721,nymax=361,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359  ! 0.5 degree 138 level
    152152!  integer,parameter :: nxmax=181,nymax=91,nuvzmax=92,nwzmax=92,nzmax=92,nxshift=0  ! CERA 2.0 degree 92 level
    153153
    154154! GFS
    155 !  integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0
     155   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0
    156156
    157157
  • 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!
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG