Changeset db712a8 in flexpart.git for src/gridcheck_nests.f90


Ignore:
Timestamp:
Mar 3, 2016, 12:34:56 PM (8 years ago)
Author:
Espen Sollum ATMOS <eso@…>
Branches:
master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
Children:
38b7917
Parents:
b0434e1
Message:

Completed handling of nested wind fields with cloud water (for wet deposition).

File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/gridcheck_nests.f90

    rb0434e1 rdb712a8  
    172172  elseif ((parCat.eq.1).and.(parNum.eq.83).and.(typSurf.eq.105)) then ! clwc
    173173    isec1(6)=246         ! indicatorOfParameter
    174     ! readclouds=.true.
    175     ! sumclouds=.false.
    176174  elseif ((parCat.eq.1).and.(parNum.eq.84).and.(typSurf.eq.105)) then ! ciwc
    177175    isec1(6)=247         ! indicatorOfParameter
     
    180178  elseif ((parCat.eq.201).and.(parNum.eq.31).and.(typSurf.eq.105)) then ! qc
    181179    isec1(6)=201031      ! indicatorOfParameter
    182     ! readclouds=.true.
    183     ! sumclouds=.true.
    184180  elseif ((parCat.eq.3).and.(parNum.eq.0).and.(typSurf.eq.1)) then !SP
    185181    isec1(6)=134         ! indicatorOfParameter
    186182  elseif ((parCat.eq.2).and.(parNum.eq.32)) then ! W, actually eta dot
    187183    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 gridchek.f90
    189     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    !
    190186  elseif ((parCat.eq.3).and.(parNum.eq.0).and.(typSurf.eq.101)) then !SLP
    191187    isec1(6)=151         ! indicatorOfParameter
     
    258254  endif ! ifield
    259255
     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
    260274  !HSO  get the second part of the grid dimensions only from GRiB1 messages
    261275  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!!
     
    333347  if(nuvzn.eq.nlev_ec) nwzn=nlev_ecn+1
    334348
    335   if (nxn(l).gt.nxmaxn) then
    336   write(*,*) 'FLEXPART error: Too many grid points in x direction.'
    337   write(*,*) 'Reduce resolution of wind fields (file GRIDSPEC)'
    338   write(*,*) 'for nesting level ',l
    339   write(*,*) 'Or change parameter settings in file par_mod.'
    340   write(*,*) nxn(l),nxmaxn
    341   stop
    342   endif
    343 
    344   if (nyn(l).gt.nymaxn) then
    345   write(*,*) 'FLEXPART error: Too many grid points in y direction.'
    346   write(*,*) 'Reduce resolution of wind fields (file GRIDSPEC)'
    347   write(*,*) 'for nesting level ',l
    348   write(*,*) 'Or change parameter settings in file par_mod.'
    349   write(*,*) nyn(l),nymaxn
    350   stop
    351   endif
    352 
    353349  if ((nuvzn.gt.nuvzmax).or.(nwzn.gt.nwzmax)) then
    354350  write(*,*) 'FLEXPART error: Nested wind fields have too many'// &
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG