Changeset 553d0a7 in flexpart.git for src/readwind_ecmwf_mpi.f90
- Timestamp:
- Nov 7, 2022, 10:31:29 AM (18 months ago)
- Branches:
- scaling-bug
- Parents:
- 49e63b9
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/readwind_ecmwf_mpi.f90
r92fab65 r553d0a7 211 211 elseif ((parCat.eq.0).and.(parNum.eq.6).and.(typSurf.eq.103)) then ! 2D 212 212 isec1(6)=168 ! indicatorOfParameter 213 elseif ((parCat.eq.1).and.(parNum.eq.11).and.(typSurf.eq.1)) then ! SD 213 elseif ((parCat.eq.1).and.(parNum.eq.11).and.(typSurf.eq.1)) then ! SDE 214 214 isec1(6)=141 ! indicatorOfParameter 215 conversion_factor=1000. 215 conversion_factor=0.1 ! convert snow depth to water equivalent 216 elseif ((parCat.eq.1).and.(parNum.eq.254).and.(typSurf.eq.1)) then ! SD 217 isec1(6)=141 ! indicatorOfParameter 216 218 elseif ((parCat.eq.6).and.(parNum.eq.1) .or. parId .eq. 164) then ! CC 217 219 isec1(6)=164 ! indicatorOfParameter 218 220 elseif ((parCat.eq.1).and.(parNum.eq.9) .or. parId .eq. 142) then ! LSP 219 221 isec1(6)=142 ! indicatorOfParameter 220 elseif ((parCat.eq.1).and.(parNum.eq.10)) then ! CP222 elseif ((parCat.eq.1).and.(parNum.eq.10)) then ! ACPCP 221 223 isec1(6)=143 ! indicatorOfParameter 222 conversion_factor=1000.223 224 elseif ((parCat.eq.0).and.(parNum.eq.11).and.(typSurf.eq.1)) then ! SHF 224 225 isec1(6)=146 ! indicatorOfParameter … … 315 316 zsec4(nxfield*(ny-j-1)+i+1) 316 317 if(isec1(6).eq.141) sd(i,j,1,n)= &!! SNOW DEPTH 317 zsec4(nxfield*(ny-j-1)+i+1) /conversion_factor318 zsec4(nxfield*(ny-j-1)+i+1)*conversion_factor 318 319 if(isec1(6).eq.151) msl(i,j,1,n)= &!! SEA LEVEL PRESS. 319 320 zsec4(nxfield*(ny-j-1)+i+1) … … 333 334 endif 334 335 if(isec1(6).eq.143) then !! CONVECTIVE PREC. 335 convprec(i,j,1,n)=zsec4(nxfield*(ny-j-1)+i+1) /conversion_factor336 convprec(i,j,1,n)=zsec4(nxfield*(ny-j-1)+i+1) 336 337 if (convprec(i,j,1,n).lt.0.) convprec(i,j,1,n)=0. 337 338 endif
Note: See TracChangeset
for help on using the changeset viewer.