Changeset 553d0a7 in flexpart.git for src/readwind_ecmwf.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.f90
r92fab65 r553d0a7 201 201 elseif ((parCat.eq.0).and.(parNum.eq.6).and.(typSurf.eq.103)) then ! 2D 202 202 isec1(6)=168 ! indicatorOfParameter 203 elseif ((parCat.eq.1).and.(parNum.eq.11).and.(typSurf.eq.1)) then ! SD 203 elseif ((parCat.eq.1).and.(parNum.eq.11).and.(typSurf.eq.1)) then ! SDE 204 204 isec1(6)=141 ! indicatorOfParameter 205 conversion_factor=1000. 205 conversion_factor=0.1 ! convert snow depth to water equivalent 206 elseif ((parCat.eq.1).and.(parNum.eq.254).and.(typSurf.eq.1)) then ! SD 207 isec1(6)=141 ! indicatorOfParameter 206 208 elseif ((parCat.eq.6).and.(parNum.eq.1) .or. parId .eq. 164) then ! CC 207 209 isec1(6)=164 ! indicatorOfParameter 208 210 elseif ((parCat.eq.1).and.(parNum.eq.9) .or. parId .eq. 142) then ! LSP 209 211 isec1(6)=142 ! indicatorOfParameter 210 elseif ((parCat.eq.1).and.(parNum.eq.10)) then ! CP212 elseif ((parCat.eq.1).and.(parNum.eq.10)) then ! ACPCP 211 213 isec1(6)=143 ! indicatorOfParameter 212 conversion_factor=1000.213 214 elseif ((parCat.eq.0).and.(parNum.eq.11).and.(typSurf.eq.1)) then ! SHF 214 215 isec1(6)=146 ! indicatorOfParameter … … 305 306 zsec4(nxfield*(ny-j-1)+i+1) 306 307 if(isec1(6).eq.141) sd(i,j,1,n)= &!! SNOW DEPTH 307 zsec4(nxfield*(ny-j-1)+i+1) /conversion_factor308 zsec4(nxfield*(ny-j-1)+i+1)*conversion_factor 308 309 if(isec1(6).eq.151) msl(i,j,1,n)= &!! SEA LEVEL PRESS. 309 310 zsec4(nxfield*(ny-j-1)+i+1) … … 323 324 endif 324 325 if(isec1(6).eq.143) then !! CONVECTIVE PREC. 325 convprec(i,j,1,n)=zsec4(nxfield*(ny-j-1)+i+1) /conversion_factor326 convprec(i,j,1,n)=zsec4(nxfield*(ny-j-1)+i+1) 326 327 if (convprec(i,j,1,n).lt.0.) convprec(i,j,1,n)=0. 327 328 endif
Note: See TracChangeset
for help on using the changeset viewer.