Changeset db712a8 in flexpart.git for src/gridcheck_nests.f90
- Timestamp:
- Mar 3, 2016, 12:34:56 PM (8 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 38b7917
- Parents:
- b0434e1
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/gridcheck_nests.f90
rb0434e1 rdb712a8 172 172 elseif ((parCat.eq.1).and.(parNum.eq.83).and.(typSurf.eq.105)) then ! clwc 173 173 isec1(6)=246 ! indicatorOfParameter 174 ! readclouds=.true.175 ! sumclouds=.false.176 174 elseif ((parCat.eq.1).and.(parNum.eq.84).and.(typSurf.eq.105)) then ! ciwc 177 175 isec1(6)=247 ! indicatorOfParameter … … 180 178 elseif ((parCat.eq.201).and.(parNum.eq.31).and.(typSurf.eq.105)) then ! qc 181 179 isec1(6)=201031 ! indicatorOfParameter 182 ! readclouds=.true.183 ! sumclouds=.true.184 180 elseif ((parCat.eq.3).and.(parNum.eq.0).and.(typSurf.eq.1)) then !SP 185 181 isec1(6)=134 ! indicatorOfParameter 186 182 elseif ((parCat.eq.2).and.(parNum.eq.32)) then ! W, actually eta dot 187 183 isec1(6)=135 ! indicatorOfParameter 188 elseif ((parCat.eq.128).and.(parNum.eq.77)) then ! W, actually eta dot !added bymc to make it consistent with new gridche k.f90189 isec1(6)=135 ! indicatorOfParameter ! ! " " " " " " " " " " " " " " " " " " " " " " " " " " "184 elseif ((parCat.eq.128).and.(parNum.eq.77)) then ! W, actually eta dot !added bymc to make it consistent with new gridcheck.f90 185 isec1(6)=135 ! indicatorOfParameter ! 190 186 elseif ((parCat.eq.3).and.(parNum.eq.0).and.(typSurf.eq.101)) then !SLP 191 187 isec1(6)=151 ! indicatorOfParameter … … 258 254 endif ! ifield 259 255 256 if (nxn(l).gt.nxmaxn) then 257 write(*,*) 'FLEXPART error: Too many grid points in x direction.' 258 write(*,*) 'Reduce resolution of wind fields (file GRIDSPEC)' 259 write(*,*) 'for nesting level ',l 260 write(*,*) 'Or change parameter settings in file par_mod.' 261 write(*,*) nxn(l),nxmaxn 262 stop 263 endif 264 265 if (nyn(l).gt.nymaxn) then 266 write(*,*) 'FLEXPART error: Too many grid points in y direction.' 267 write(*,*) 'Reduce resolution of wind fields (file GRIDSPEC)' 268 write(*,*) 'for nesting level ',l 269 write(*,*) 'Or change parameter settings in file par_mod.' 270 write(*,*) nyn(l),nymaxn 271 stop 272 endif 273 260 274 !HSO get the second part of the grid dimensions only from GRiB1 messages 261 275 if (isec1(6) .eq. 167 .and. (gotGrib.eq.0)) then !added by mc to make it consistent with new gridchek.f90 note that gotGrid must be changed in gotGrib!! … … 333 347 if(nuvzn.eq.nlev_ec) nwzn=nlev_ecn+1 334 348 335 if (nxn(l).gt.nxmaxn) then336 write(*,*) 'FLEXPART error: Too many grid points in x direction.'337 write(*,*) 'Reduce resolution of wind fields (file GRIDSPEC)'338 write(*,*) 'for nesting level ',l339 write(*,*) 'Or change parameter settings in file par_mod.'340 write(*,*) nxn(l),nxmaxn341 stop342 endif343 344 if (nyn(l).gt.nymaxn) then345 write(*,*) 'FLEXPART error: Too many grid points in y direction.'346 write(*,*) 'Reduce resolution of wind fields (file GRIDSPEC)'347 write(*,*) 'for nesting level ',l348 write(*,*) 'Or change parameter settings in file par_mod.'349 write(*,*) nyn(l),nymaxn350 stop351 endif352 353 349 if ((nuvzn.gt.nuvzmax).or.(nwzn.gt.nwzmax)) then 354 350 write(*,*) 'FLEXPART error: Nested wind fields have too many'// &
Note: See TracChangeset
for help on using the changeset viewer.