debugprint.F90

Go to the documentation of this file.

module debugprint

  use esmf,    only : esmf_logfounderror, esmf_logerr_passthru
  use grdvars, only : ni,nj,ipole,angle,anglet
  use grdvars, only : htn,latct,lonct,latcv,loncv,latcu,loncu,latbu,lonbu
  use grdvars, only : xlatct,xlonct,xlatcu,xloncu
  use grdvars, only : latbu_vert,lonbu_vert,latcv_vert,loncv_vert
  use grdvars, only : latct_vert,lonct_vert,latcu_vert,loncu_vert

  implicit none
  private

  public :: checkseam, checkxlatlon, checkpoint

contains

  subroutine checkseam

    ! local variables
    integer :: j,i1,i2

    j = nj
    i1 = ipole(1); i2 = ipole(2)+1

    !htn must be the same along seam
    j = nj
    i1 = ipole(1); i2 = ipole(2)+1
    print *,'HTN across seam '
    print *,htn(i1-2,j),htn(i2+2,j)
    print *,htn(i1-1,j),htn(i2+1,j)
    print *,htn(i1,  j),htn(i2,  j)
    print *,htn(i1+1,j),htn(i2-1,j)
    print *,htn(i1+2,j),htn(i2-2,j)

    print *,'latCv across seam '
    print *,latcv(i1-2,j),latcv(i2+2,j)
    print *,latcv(i1-1,j),latcv(i2+1,j)
    print *,latcv(i1,  j),latcv(i2,  j)
    print *,latcv(i1+1,j),latcv(i2-1,j)
    print *,latcv(i1+2,j),latcv(i2-2,j)

    print *,'lonCv across seam '
    print *,loncv(i1-2,j),loncv(i2+2,j)
    print *,loncv(i1-1,j),loncv(i2+1,j)
    print *,loncv(i1,  j),loncv(i2,  j)
    print *,loncv(i1+1,j),loncv(i2-1,j)
    print *,loncv(i1+2,j),loncv(i2-2,j)

    print *,'angleT across seam '
    print *,anglet(i1-2,j),anglet(i2+2,j)
    print *,anglet(i1-1,j),anglet(i2+1,j)
    print *,anglet(i1,  j),anglet(i2,  j)
    print *,anglet(i1+1,j),anglet(i2-1,j)
    print *,anglet(i1+2,j),anglet(i2-2,j)

    print *,'latCu across seam '
    print *,latcu(i1-3,j),latcu(i2+2,j),latcu(i1-3,j)-latcu(i2+2,j)
    print *,latcu(i1-2,j),latcu(i2+1,j)
    print *,latcu(i1-1,j),latcu(i2+0,j)
    print *,latcu(i1,  j),latcu(i2-1,j)
    print *,latcu(i1+1,j),latcu(i2-2,j)
    print *,latcu(i1+2,j),latcu(i2-3,j)
    print *,latcu(i1+3,j),latcu(i2-4,j)

    print *,'lonCu across seam '
    print *,loncu(i1-3,j),loncu(i2+2,j),loncu(i1-3,j)+loncu(i2+2,j)
    print *,loncu(i1-2,j),loncu(i2+1,j)
    print *,loncu(i1-1,j),loncu(i2+0,j)
    print *,loncu(i1,  j),loncu(i2-1,j)
    print *,loncu(i1+1,j),loncu(i2-2,j)
    print *,loncu(i1+2,j),loncu(i2-3,j)
    print *,loncu(i1+3,j),loncu(i2-4,j)

    print *,'latCt across seam '
    print *,latct(i1-3,j),latct(i2+3,j),latct(i1-3,j)-latct(i2+3,j)
    print *,latct(i1-2,j),latct(i2+2,j)
    print *,latct(i1-1,j),latct(i2+1,j)
    print *,latct(i1,  j),latct(i2,  j)
    print *,latct(i1+1,j),latct(i2-1,j)
    print *,latct(i1+2,j),latct(i2-2,j)
    print *,latct(i1+3,j),latct(i2-3,j)

    print *,'lonCt across seam '
    print *,lonct(i1-3,j),lonct(i2+3,j),lonct(i1-3,j)+lonct(i2+3,j)
    print *,lonct(i1-2,j),lonct(i2+2,j)
    print *,lonct(i1-1,j),lonct(i2+1,j)
    print *,lonct(i1,  j),lonct(i2,  j)
    print *,lonct(i1+1,j),lonct(i2-1,j)
    print *,lonct(i1+2,j),lonct(i2-2,j)
    print *,lonct(i1+3,j),lonct(i2-3,j)
    print *
  end subroutine checkseam


  subroutine checkxlatlon

    ! local variables
    integer :: i

    print *,'============== Ct grid ==============='
    print *,'============== Left pole ============'
    do i = ipole(1)-3,ipole(1)+3
       print '(i5,6f12.5)',i,lonct(i,nj),xlonct(i),lonct(i,nj)+xlonct(i),latct(i,nj),xlatct(i),latct(i,nj)-xlatct(i)
    enddo
    print *

    print *,'============ Right pole ============'
    do i = ipole(2)-3,ipole(2)+3
       print '(i5,6f12.5)',i,lonct(i,nj),xlonct(i),lonct(i,nj)+xlonct(i),latct(i,nj),xlatct(i),latct(i,nj)-xlatct(i)
    enddo
    print *

    print *,'============== Ct grid ==============='
    print *,'============== Left edge ============'
    do i = 1,5
       print '(i5,6f12.5)',i,lonct(i,nj),xlonct(i),lonct(i,nj)+xlonct(i),latct(i,nj),xlatct(i),latct(i,nj)-xlatct(i)
    enddo
    print *
    print *,'============== Right edge ==========='
    do i = ni-4,ni
       print '(i5,6f12.5)',i,lonct(i,nj),xlonct(i),lonct(i,nj)+xlonct(i),latct(i,nj),xlatct(i),latct(i,nj)-xlatct(i)
    enddo
    print *


    print *,'============== Cu grid ==============='
    print *,'============== Left pole ============='
    do i = ipole(1)-3,ipole(1)+3
       print '(i5,6f12.5)',i,loncu(i,nj),xloncu(i),loncu(i,nj)+xloncu(i),latcu(i,nj),xlatcu(i),latcu(i,nj)-xlatcu(i)
    enddo
    print *

    print *,'============ Right pole ============'
    do i = ipole(2)-3,ipole(2)+3
       print '(i5,6f12.5)',i,loncu(i,nj),xloncu(i),loncu(i,nj)+xloncu(i),latcu(i,nj),xlatcu(i),latcu(i,nj)-xlatcu(i)
    enddo
    print *

    print *,'============== Cu grid ==============='
    print *,'============== Left edge ============'
    do i = 1,5
       print '(i5,6f12.5)',i,loncu(i,nj),xloncu(i),loncu(i,nj)+xloncu(i),latcu(i,nj),xlatcu(i),latcu(i,nj)-xlatcu(i)
    enddo
    print *
    print *,'============== Right edge ==========='
    do i = ni-4,ni
       print '(i5,6f12.5)',i,loncu(i,nj),xloncu(i),loncu(i,nj)+xloncu(i),latcu(i,nj),xlatcu(i),latcu(i,nj)-xlatcu(i)
    enddo
    print *

  end subroutine checkxlatlon


  subroutine checkpoint

    ! local variables
    integer :: i,j

    ! check
    i = 1; j = nj
    print '(f12.5,a,f12.5)',latbu_vert(i,j,2),'        ',latbu_vert(i,j,1)
    print '(a12,f12.5)','          ',latbu(i,j)
    print '(f12.5,a,f12.5)',latbu_vert(i,j,3),'        ',latbu_vert(i,j,4)
    print *
    print '(f12.5,a,f12.5)',lonbu_vert(i,j,2),'        ',lonbu_vert(i,j,1)
    print '(a12,f12.5)','          ',lonbu(i,j)
    print '(f12.5,a,f12.5)',lonbu_vert(i,j,3),'        ',lonbu_vert(i,j,4)
    print *
    print *
    ! check
    print '(f12.5,a,f12.5)',latcv_vert(i,j,2),'        ',latcv_vert(i,j,1)
    print '(a12,f12.5)','          ',latcv(i,j)
    print '(f12.5,a,f12.5)',latcv_vert(i,j,3),'        ',latcv_vert(i,j,4)
    print *
    print '(f12.5,a,f12.5)',loncv_vert(i,j,2),'        ',loncv_vert(i,j,1)
    print '(a12,f12.5)','          ',loncv(i,j)
    print '(f12.5,a,f12.5)',loncv_vert(i,j,3),'        ',loncv_vert(i,j,4)

    print *
    print *

    i = 1; j = 10
    print '(f12.5,a,f12.5)',latct_vert(i,j,2),'        ',latct_vert(i,j,1)
    print '(a12,f12.5)','          ',latct(i,j)
    print '(f12.5,a,f12.5)',latct_vert(i,j,3),'        ',latct_vert(i,j,4)
    print *
    print '(f12.5,a,f12.5)',lonct_vert(i,j,2),'        ',lonct_vert(i,j,1)
    print '(a12,f12.5)','          ',lonct(i,j)
    print '(f12.5,a,f12.5)',lonct_vert(i,j,3),'        ',lonct_vert(i,j,4)
    print *
    print *
    ! check
    print '(f12.5,a,f12.5)',latcu_vert(i,j,2),'        ',latcu_vert(i,j,1)
    print '(a12,f12.5)','          ',latcu(i,j)
    print '(f12.5,a,f12.5)',latcu_vert(i,j,3),'        ',latcu_vert(i,j,4)
    print *
    print '(f12.5,a,f12.5)',loncu_vert(i,j,2),'        ',loncu_vert(i,j,1)
    print '(a12,f12.5)','          ',loncu(i,j)
    print '(f12.5,a,f12.5)',loncu_vert(i,j,3),'        ',loncu_vert(i,j,4)


    i = ni; j = 10
    print '(f12.5,a,f12.5)',latct_vert(i,j,2),'        ',latct_vert(i,j,1)
    print '(a12,f12.5)','          ',latct(i,j)
    print '(f12.5,a,f12.5)',latct_vert(i,j,3),'        ',latct_vert(i,j,4)
    print *
    print '(f12.5,a,f12.5)',lonct_vert(i,j,2),'        ',lonct_vert(i,j,1)
    print '(a12,f12.5)','          ',lonct(i,j)
    print '(f12.5,a,f12.5)',lonct_vert(i,j,3),'        ',lonct_vert(i,j,4)
    print *
    print *
    ! check
    print '(f12.5,a,f12.5)',latcu_vert(i,j,2),'        ',latcu_vert(i,j,1)
    print '(a12,f12.5)','          ',latcu(i,j)
    print '(f12.5,a,f12.5)',latcu_vert(i,j,3),'        ',latcu_vert(i,j,4)
    print *
    print '(f12.5,a,f12.5)',loncu_vert(i,j,2),'        ',loncu_vert(i,j,1)
    print '(a12,f12.5)','          ',loncu(i,j)
    print '(f12.5,a,f12.5)',loncu_vert(i,j,3),'        ',loncu_vert(i,j,4)

    print *,"latCt minmax ",minval(latct),maxval(latct)
    print *,"latCu minmax ",minval(latcu),maxval(latcu)
    print *,"latCv minmax ",minval(latcv),maxval(latcv)
    print *,"latBu minmax ",minval(latbu),maxval(latbu)

    !   print *,minval(latCt_vert),maxval(latCt_vert)
    !   print *,minval(lonCt_vert),maxval(lonCt_vert)
    !   print *,minval(latBu_vert),maxval(latBu_vert)
    !   print *,minval(lonBu_vert),maxval(lonBu_vert)
  end subroutine checkpoint
end module debugprint