pkg/fizhi/step_physics.F

a456aa407c Andr*0001 #include "FIZHI_OPTIONS.h"
f98d2ec0f4 Andr*0002        subroutine step_physics(uphy,vphy,thphy,sphy,dt,
c2377c66e1 Andr*0003      .  idim1,idim2,jdim1,jdim2,Nrphys,Nsx,Nsy,im1,im2,jm1,jm2,bi,bj,
e337e4ca8c Andr*0004      .                                         duphy,dvphy,dthphy,dsphy)
                0005 c-----------------------------------------------------------------------
                0006 c Routine to step physics state forward using the increments
                0007 c
                0008 c-----------------------------------------------------------------------
                0009        implicit none
                0010        integer im1,im2,jm1,jm2,idim1,idim2,jdim1,jdim2
                0011        integer Nrphys,Nsx,Nsy,bi,bj
c2377c66e1 Andr*0012        _RL uphy(idim1:idim2,jdim1:jdim2,Nrphys,Nsx,Nsy)
                0013        _RL vphy(idim1:idim2,jdim1:jdim2,Nrphys,Nsx,Nsy)
                0014        _RL thphy(idim1:idim2,jdim1:jdim2,Nrphys,Nsx,Nsy)
                0015        _RL sphy(idim1:idim2,jdim1:jdim2,Nrphys,Nsx,Nsy)
                0016        _RL duphy(idim1:idim2,jdim1:jdim2,Nrphys,Nsx,Nsy)
                0017        _RL dvphy(idim1:idim2,jdim1:jdim2,Nrphys,Nsx,Nsy)
                0018        _RL dthphy(idim1:idim2,jdim1:jdim2,Nrphys,Nsx,Nsy)
                0019        _RL dsphy(idim1:idim2,jdim1:jdim2,Nrphys,Nsx,Nsy)
f98d2ec0f4 Andr*0020        _RL dt
e337e4ca8c Andr*0021 c
                0022        integer i,j,L
                0023
                0024        do L = 1,Nrphys
c2377c66e1 Andr*0025         do j = jm1,jm2
                0026         do i = im1,im2
f98d2ec0f4 Andr*0027          uphy(i,j,L,bi,bj) = uphy(i,j,L,bi,bj) + duphy(i,j,L,bi,bj) * dt
                0028          vphy(i,j,L,bi,bj) = vphy(i,j,L,bi,bj) + dvphy(i,j,L,bi,bj) * dt
                0029          thphy(i,j,L,bi,bj)=thphy(i,j,L,bi,bj) + dthphy(i,j,L,bi,bj)* dt
                0030          sphy(i,j,L,bi,bj) = sphy(i,j,L,bi,bj) + dsphy(i,j,L,bi,bj) * dt
                0031         enddo
e337e4ca8c Andr*0032         enddo
                0033        enddo
                0034
                0035        return
                0036        end