Changeset 1c3c778 in flexpart.git
- Timestamp:
- Mar 14, 2016, 11:06:31 AM (8 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 7e52e2e
- Parents:
- c2162ce
- Location:
- src
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
src/makefile
rdb712a8 r1c3c778 11 11 # To use gfortran version 4.9, add "gcc=4.9" to the make command, e.g. 12 12 # 'make -j ecmwf gcc=4.9', 13 # also set environment variable LD_LIBRARY_ FLAGS=13 # also set environment variable LD_LIBRARY_PATH to point to compiler libraries 14 14 # 15 15 # USAGE 16 # Compile serial FLEXPART (ECMWF) 16 # Compile serial FLEXPART (ECMWF)a 17 17 # make [-j] ecmwf 18 18 # -
src/mpi_mod.f90
r38b7917 r1c3c778 61 61 62 62 use mpi 63 use par_mod, only: dp, 63 use par_mod, only: dp,sp 64 64 use com_mod, only: lroot 65 65 … … 1409 1409 integer :: d2s1 = nxmax*nymax 1410 1410 integer :: d2s2 = nxmax*nymax*maxspec 1411 integer :: d1_size1 = maxwf1411 !integer :: d1_size1 = maxwf 1412 1412 1413 1413 ! integer :: d3s1,d3s2,d2s1,d2s2 … … 2154 2154 2155 2155 pv(:,:,:,li:ui)=1.0 2156 clouds(:,:,:,li:ui)=0 .02156 clouds(:,:,:,li:ui)=0 2157 2157 2158 2158 clwc(:,:,:,li:ui)=0.0 … … 2161 2161 ! 2D fields 2162 2162 2163 cloudsh(:,:,li:ui)=0 .02163 cloudsh(:,:,li:ui)=0 2164 2164 vdep(:,:,:,li:ui)=0.0 2165 2165 ps(:,:,:,li:ui)=1.0e5 -
src/par_mod.f90
r38b7917 r1c3c778 83 83 real,parameter :: d_trop=50., d_strat=0.1 84 84 real,parameter :: rho_water=1000. !ZHG 2015 [kg/m3] 85 !ZHG MAR2016 86 real,parameter :: incloud_ratio=1.6 87 85 88 ! karman Karman's constant 86 89 ! href [m] Reference height for dry deposition … … 216 219 integer,parameter :: maxpart=40000000 217 220 integer,parameter :: maxspec=1 218 real,parameter :: minmass=0.0 !0.0001221 real,parameter :: minmass=0.0001 219 222 220 223 ! maxpart Maximum number of particles -
src/timemanager_mpi.f90
r7999df47 r1c3c778 825 825 ! eso :TODO: this not implemented yet (transfer particles to PID 0 or rewrite) 826 826 ! the tools to do this are already in mpi_mod.f90 827 ! :CUR: check this828 827 if (lroot) then 829 828 do j=1,numpart -
src/wetdepo.f90
rdb712a8 r1c3c778 143 143 ! Determine nested grid coordinates 144 144 !********************************** 145 readclouds_this_nest=.false. 145 146 146 147 if (ngrid.gt.0) then … … 149 150 ix=int(xtn) 150 151 jy=int(ytn) 151 if (readclouds_nest(ngrid)) then 152 readclouds_this_nest=.true. 153 else 154 readclouds_this_nest=.false. 155 end if 152 if (readclouds_nest(ngrid)) readclouds_this_nest=.true. 156 153 else 157 154 ix=int(xtra1(jpart)) … … 364 361 !OLD 365 362 if ((readclouds.and.ngrid.eq.0).or.(readclouds_this_nest.and.ngrid.gt.0)) then 366 wetscav= S_i*(prec(1)/3.6E6)363 wetscav=incloud_ratio*S_i*(prec(1)/3.6E6) 367 364 else 368 wetscav= S_i*(prec(1)/3.6E6)/clouds_h365 wetscav=incloud_ratio*S_i*(prec(1)/3.6E6)/clouds_h 369 366 endif 370 367 -
src/wetdepokernel_nest.f90
rdb712a8 r1c3c778 53 53 real :: x,y,deposit(maxspec),ddx,ddy,xl,yl,wx,wy,w 54 54 integer :: ix,jy,ixp,jyp,ks,kp,nunc,nage 55 56 real :: dbg_dx, dbg_dy, dbg_xoutshiftn, dbg_youtshiftn, dbg_dxoutn, dbg_dyoutn,dbg_t57 55 58 56 xl=(x*dx+xoutshiftn)/dxoutn
Note: See TracChangeset
for help on using the changeset viewer.