Changes in src/verttransform_ecmwf.f90 [79e0349:47f96e5] in flexpart.git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/verttransform_ecmwf.f90
r79e0349 r47f96e5 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/users/sec/kleinprojekte/hertlfit/' 111 character(len=60) :: fnameH,fnameI,fnameJ 110 ! character(len=60) :: zhgpath='/xnilu_wrk/flex_wrk/zhg/' 112 111 ! character(len=60) :: fnameA,fnameB,fnameC,fnameD,fnameE,fnameF,fnameG,fnameH 113 CHARACTER(LEN=3) :: aspec114 integer :: virr=0112 ! CHARACTER(LEN=3) :: aspec 113 ! integer :: virr=0 115 114 !real :: tot_cloud_h 116 115 !real :: dbg_height(nzmax) … … 733 732 !********* TEST ************'** 734 733 !teller(:)=0 735 virr=virr+1736 WRITE(aspec, '(i3.3)'), virr734 !virr=virr+1 735 !WRITE(aspec, '(i3.3)'), virr 737 736 738 737 !if (readclouds) then … … 771 770 !fnameF=trim(zhgpath)//trim(aspec)//'lsp.txt' 772 771 !fnameG=trim(zhgpath)//trim(aspec)//'convp.txt' 773 if (1.eq.2) then774 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: ',fnameH778 772 !if (readclouds) then 779 773 !OPEN(UNIT=111, FILE=fnameA,FORM='FORMATTED',STATUS = 'UNKNOWN') … … 782 776 !OPEN(UNIT=114, FILE=fnameD,FORM='FORMATTED',STATUS = 'UNKNOWN') 783 777 !else 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')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') 787 781 !endif 788 782 ! 789 do ix=0,nxmin1783 !do ix=0,nxmin1 790 784 !if (readclouds) then 791 785 !write(111,*) (icloud_stats(ix,jy,1,n),jy=0,nymin1) … … 794 788 !write(114,*) (icloud_stats(ix,jy,4,n),jy=0,nymin1) 795 789 !else 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) 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 799 793 !endif 800 end do794 !end do 801 795 ! 802 796 !if (readclouds) then … … 806 800 !CLOSE(114) 807 801 !else 808 CLOSE(115) 809 CLOSE(116) 810 CLOSE(117) 811 endif 802 !CLOSE(115) 803 !CLOSE(116) 804 !CLOSE(117) 812 805 !endif 813 806 !
Note: See TracChangeset
for help on using the changeset viewer.