Changeset 79e0349 in flexpart.git
- Timestamp:
- Dec 5, 2018, 11:50:52 AM (5 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug
- Children:
- db91eb7
- Parents:
- 2bfec12
- Location:
- src
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
src/get_wetscav.f90
r6985a98 r79e0349 298 298 else !parameterize cloudwater m2/m3 299 299 !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 301 305 endif 302 306 … … 318 322 if (dquer(ks).gt.0.) then 319 323 S_i= frac_act/cl 324 ! write(*,*) 'Si: ',S_i 320 325 321 326 ! GAS … … 334 339 wetscav=incloud_ratio*S_i*(prec(1)/3.6E6) 335 340 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 337 344 endif 338 345 endif ! positive in-cloud scavenging parameters given in Species file -
src/par_mod.f90
r1f55de1 r79e0349 148 148 ! ECMWF 149 149 ! 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 level150 ! integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359 ! 1.0 degree 138 level 151 151 ! integer,parameter :: nxmax=721,nymax=361,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359 ! 0.5 degree 138 level 152 152 ! integer,parameter :: nxmax=181,nymax=91,nuvzmax=92,nwzmax=92,nzmax=92,nxshift=0 ! CERA 2.0 degree 92 level 153 153 154 154 ! GFS 155 !integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0155 integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0 156 156 157 157 -
src/verttransform_ecmwf.f90
r47f96e5 r79e0349 108 108 ! real :: rcw(0:nxmax-1,0:nymax-1) 109 109 ! 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 111 112 ! character(len=60) :: fnameA,fnameB,fnameC,fnameD,fnameE,fnameF,fnameG,fnameH 112 !CHARACTER(LEN=3) :: aspec113 !integer :: virr=0113 CHARACTER(LEN=3) :: aspec 114 integer :: virr=0 114 115 !real :: tot_cloud_h 115 116 !real :: dbg_height(nzmax) … … 732 733 !********* TEST ************'** 733 734 !teller(:)=0 734 !virr=virr+1735 !WRITE(aspec, '(i3.3)'), virr735 virr=virr+1 736 WRITE(aspec, '(i3.3)'), virr 736 737 737 738 !if (readclouds) then … … 770 771 !fnameF=trim(zhgpath)//trim(aspec)//'lsp.txt' 771 772 !fnameG=trim(zhgpath)//trim(aspec)//'convp.txt' 773 if (1.eq.2) then 774 fnameH=trim(zhgpath)//trim(aspec)//'tcwc.txt' 775 fnameI=trim(zhgpath)//trim(aspec)//'prec.txt' 776 fnameJ=trim(zhgpath)//trim(aspec)//'cloudsh.txt' 777 write(*,*) 'Writing data to file: ',fnameH 772 778 !if (readclouds) then 773 779 !OPEN(UNIT=111, FILE=fnameA,FORM='FORMATTED',STATUS = 'UNKNOWN') … … 776 782 !OPEN(UNIT=114, FILE=fnameD,FORM='FORMATTED',STATUS = 'UNKNOWN') 777 783 !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')784 OPEN(UNIT=115, FILE=fnameH,FORM='FORMATTED',STATUS = 'UNKNOWN') 785 OPEN(UNIT=116, FILE=fnameI,FORM='FORMATTED',STATUS = 'UNKNOWN') 786 OPEN(UNIT=117, FILE=fnameJ,FORM='FORMATTED',STATUS = 'UNKNOWN') 781 787 !endif 782 788 ! 783 !do ix=0,nxmin1789 do ix=0,nxmin1 784 790 !if (readclouds) then 785 791 !write(111,*) (icloud_stats(ix,jy,1,n),jy=0,nymin1) … … 788 794 !write(114,*) (icloud_stats(ix,jy,4,n),jy=0,nymin1) 789 795 !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-02792 !write(117,*) (convprec(ix,jy,1,n),jy=0,nymin1) !5.38330078E-02 796 write(115,*) (ctwc(ix,jy,n),jy=0,nymin1) 797 write(116,*) (lsprec(ix,jy,1,n)+convprec(ix,jy,1,n),jy=0,nymin1) 798 write(117,*) (cloudsh(ix,jy,n),jy=0,nymin1) 793 799 !endif 794 !end do800 end do 795 801 ! 796 802 !if (readclouds) then … … 800 806 !CLOSE(114) 801 807 !else 802 !CLOSE(115) 803 !CLOSE(116) 804 !CLOSE(117) 808 CLOSE(115) 809 CLOSE(116) 810 CLOSE(117) 811 endif 805 812 !endif 806 813 !
Note: See TracChangeset
for help on using the changeset viewer.