PGF90 (Version     12.8)          08/30/2020  00:27:52      page 1

Switches: -noasm -nodclchk -nodebug -nodlines -noline -list
          -idir ../include
          -inform warn -opt 1 -nosave -object -noonetrip
          -depchk on -nostandard     
          -nosymbol -noupcase    

Filename: llwinds.f90

(    1) !&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
(    2)       subroutine llwinds(usd, vsd, ull, vll,lsd)
(    3) !     ******************************************************************
(    4) !     *                                                                *
(    5) !     *   Convert covariant to ll winds on a cubic grid                *
(    6) !     *                                                                *
(    7) !     ******************************************************************
(    8) !-----------------------------------------------------------------------
(    9)        USE mod_graph
(   10)       include 'param_o.h'
(   11) !-----------------------------------------------------------------------
../include/param_o.h
(    1)*      integer,parameter  :: im0=401
(    2)*      integer, parameter :: nm=6
(    3)*      integer, parameter :: lm = 50
(    4)*      integer, parameter :: nsub = 10
(    5)*
(    6)*      integer, parameter :: im=im0, jm=im
(    7)*      integer, parameter :: im1=im-1, jm1=jm-1
(    8)*
(    9)*      integer, parameter :: lm1 = lm-1, lp1 =lm +1  
(   10)*
(   11)*      integer, parameter :: ixm = nsub, jym = ixm, nxy = ixm*jym
(   12)*      integer, parameter :: ildom = (im - 1)/ixm, jldom = (jm - 1)/jym
(   13)*      integer, parameter :: ilm = (im - 1)/ixm +1, jlm = (jm - 1)/jym +1 
(   14)*
(   15)*      logical,parameter::flat=.false.
(   16)*      logical,parameter::hstst=.false.
(   17)*
(   18)*!      integer,parameter::igm=360, jgm = 181
(   19)*      integer,parameter::igm=(im0-1)*4, jgm = (igm/2)+1      
(   20)*      real,parameter::alfa=0.0, beta=0.0, gamm=0.0
(   21)*!      real,parameter::alfa=0.*3.1415926/180.,beta=66.*3.1415926/180. &
(   22)*!	              ,gamm=175.*3.1415926/180.
(   23)*
(   24)*!GSM      integer,parameter::lsm=20
(   25)*      integer,parameter::lsm=31
(   26)*
(   27)*!RESTART
(   28)*
(   29)*      character(len=10):: restartdate='2020080412'    !restart date to create appropriate folder name
(   12)       integer::lsd
(   13) !-----------------------------------------------------------------------
(   14)       real, dimension(0:im+1,0:jm+1,nm,lsd)::  usd, vsd, ull, vll
(   15) !-----------------------------------------------------------------------   
(   16)        dimension  &
(   17)        us(lsd), vs(lsd)
(   18) !-----------------------------------------------------------------------
(   19)       do n=1,nm
(   20)       do j=1,jm1
(   21)       do i=1,im1






PGF90 (Version     12.8)          08/30/2020  00:27:52      page 2

(   22)         us(:) = usd(i,j,n,:)
(   23)         vs(:) = vsd(i,j,n,:)
(   24)         ull(i,j,n,:) = u2us(i,j,n) * us(:) + v2us(i,j,n) * vs(:)
(   25)         vll(i,j,n,:) = u2vs(i,j,n) * us(:) + v2vs(i,j,n) * vs(:)
(   26)       end do
(   27)       end do
(   28)       end do
(   29) !-----------------------------------------------------------------------
(   30)       end subroutine llwinds