...
 
Commits (2)
......@@ -329,6 +329,10 @@ subroutine gridcheck_ecmwf
!***********************************************************
xauxa=abs(xaux2+dx-360.-xaux1)
!cyclicx
print*, 'xauxa=abs(xaux2+dx-360.-xaux1)'
print*, xauxa,xaux2,dx,xaux1
!stop
if (xauxa.lt.0.001) then
nx=nxfield+1 ! field is cyclic
xglobal=.true.
......
......@@ -156,10 +156,10 @@ module par_mod
! integer,parameter :: nxmax=181,nymax=91,nuvzmax=92,nwzmax=92,nzmax=92,nxshift=0 ! CERA 2.0 degree 92 level
! GFS
integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138
! integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138
integer :: nxshift=0 ! shift not fixed for the executable
integer,parameter :: nxmax=721,nymax=361,nuvzmax=138,nwzmax=138,nzmax=138
!*********************************************
! Maximum dimensions of the nested input grids
!*********************************************
......
......@@ -116,9 +116,21 @@ subroutine readoutgrid
yr=outlat0+real(numygrid)*dyout
xr1=xlon0+real(nxmin1)*dx
yr1=ylat0+real(nymin1)*dy
if ((outlon0+eps.lt.xlon0).or.(outlat0+eps.lt.ylat0) &
.or.(xr.gt.xr1+eps).or.(yr.gt.yr1+eps)) then
write(*,*) outlon0,outlat0
if ((outlat0+eps.lt.ylat0).or.(yr.gt.yr1+eps)) then
write(*,*) 'outlat0 =',outlat0
write(*,*) 'yr1,ylat0,yr,dyout =', yr1,ylat0,yr,dyout
write(*,*) ' #### FLEXPART MODEL ERROR! LAT OF OUTPUT ####'
write(*,*) ' #### GRID IS OUTSIDE MODEL DOMAIN. CHANGE ####'
write(*,*) ' #### FILE OUTGRID IN DIRECTORY ####'
write(*,'(a)') path(1)(1:length(1))
stop
endif
!cyclicx
print*, xglobal
! If cyclic input field allow output grid to cross border (do not stop)
if ((.not.xglobal).and.((outlon0+eps.lt.xlon0).or.(xr.gt.xr1+eps))) then
write(*,*) outlon0
write(*,*) xr1,yr1,xlon0,ylat0,xr,yr,dxout,dyout
write(*,*) ' #### FLEXPART MODEL ERROR! PART OF OUTPUT ####'
write(*,*) ' #### GRID IS OUTSIDE MODEL DOMAIN. CHANGE ####'
......